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(0
+ ,162556
+ ,1081
+ ,213118
+ ,6282154
+ ,0
+ ,29790
+ ,309
+ ,81767
+ ,4321023
+ ,0
+ ,87550
+ ,458
+ ,153198
+ ,4111912
+ ,1
+ ,84738
+ ,588
+ ,-26007
+ ,223193
+ ,0
+ ,54660
+ ,302
+ ,126942
+ ,1491348
+ ,0
+ ,42634
+ ,156
+ ,157214
+ ,1629616
+ ,1
+ ,40949
+ ,481
+ ,129352
+ ,1398893
+ ,0
+ ,45187
+ ,353
+ ,234817
+ ,1926517
+ ,0
+ ,37704
+ ,452
+ ,60448
+ ,983660
+ ,0
+ ,16275
+ ,109
+ ,47818
+ ,1443586
+ ,1
+ ,25830
+ ,115
+ ,245546
+ ,1073089
+ ,1
+ ,12679
+ ,110
+ ,48020
+ ,984885
+ ,0
+ ,18014
+ ,239
+ ,-1710
+ ,1405225
+ ,1
+ ,43556
+ ,247
+ ,32648
+ ,227132
+ ,0
+ ,24811
+ ,505
+ ,95350
+ ,929118
+ ,1
+ ,6575
+ ,159
+ ,151352
+ ,1071292
+ ,1
+ ,7123
+ ,109
+ ,288170
+ ,638830
+ ,0
+ ,21950
+ ,519
+ ,114337
+ ,856956
+ ,0
+ ,37597
+ ,248
+ ,37884
+ ,992426
+ ,1
+ ,17821
+ ,373
+ ,122844
+ ,444477
+ ,0
+ ,12988
+ ,119
+ ,82340
+ ,857217
+ ,0
+ ,22330
+ ,84
+ ,79801
+ ,711969
+ ,1
+ ,13326
+ ,102
+ ,165548
+ ,702380
+ ,1
+ ,16189
+ ,295
+ ,116384
+ ,358589
+ ,1
+ ,7146
+ ,105
+ ,134028
+ ,297978
+ ,1
+ ,15824
+ ,64
+ ,63838
+ ,585715
+ ,0
+ ,27664
+ ,282
+ ,74996
+ ,657954
+ ,1
+ ,11920
+ ,182
+ ,31080
+ ,209458
+ ,1
+ ,8568
+ ,37
+ ,32168
+ ,786690
+ ,1
+ ,14416
+ ,361
+ ,49857
+ ,439798
+ ,0
+ ,3369
+ ,28
+ ,87161
+ ,688779
+ ,0
+ ,11819
+ ,85
+ ,106113
+ ,574339
+ ,0
+ ,6984
+ ,45
+ ,80570
+ ,741409
+ ,0
+ ,4519
+ ,49
+ ,102129
+ ,597793
+ ,1
+ ,2220
+ ,22
+ ,301670
+ ,644190
+ ,1
+ ,18562
+ ,155
+ ,102313
+ ,377934
+ ,1
+ ,10327
+ ,91
+ ,88577
+ ,640273
+ ,0
+ ,5336
+ ,81
+ ,112477
+ ,697458
+ ,0
+ ,2365
+ ,79
+ ,191778
+ ,550608
+ ,1
+ ,4069
+ ,145
+ ,79804
+ ,207393
+ ,1
+ ,8636
+ ,855
+ ,128294
+ ,301607
+ ,1
+ ,13718
+ ,61
+ ,96448
+ ,345783
+ ,1
+ ,4525
+ ,226
+ ,93811
+ ,501749
+ ,1
+ ,6869
+ ,105
+ ,117520
+ ,379983
+ ,1
+ ,4628
+ ,62
+ ,69159
+ ,387475
+ ,0
+ ,3689
+ ,25
+ ,101792
+ ,377305
+ ,0
+ ,4891
+ ,217
+ ,210568
+ ,370837
+ ,0
+ ,7489
+ ,322
+ ,136996
+ ,430866
+ ,1
+ ,4901
+ ,84
+ ,121920
+ ,469107
+ ,1
+ ,2284
+ ,33
+ ,76403
+ ,194493
+ ,0
+ ,3160
+ ,108
+ ,108094
+ ,530670
+ ,0
+ ,4150
+ ,150
+ ,134759
+ ,518365
+ ,0
+ ,7285
+ ,115
+ ,188873
+ ,491303
+ ,0
+ ,1134
+ ,162
+ ,146216
+ ,527021
+ ,0
+ ,4658
+ ,158
+ ,156608
+ ,233773
+ ,1
+ ,2384
+ ,97
+ ,61348
+ ,405972
+ ,1
+ ,3748
+ ,9
+ ,50350
+ ,652925
+ ,1
+ ,5371
+ ,66
+ ,87720
+ ,446211
+ ,1
+ ,1285
+ ,107
+ ,99489
+ ,341340
+ ,0
+ ,9327
+ ,101
+ ,87419
+ ,387699
+ ,0
+ ,5565
+ ,47
+ ,94355
+ ,493408
+ ,1
+ ,1528
+ ,38
+ ,60326
+ ,146494
+ ,0
+ ,3122
+ ,34
+ ,94670
+ ,414462
+ ,0
+ ,7561
+ ,87
+ ,82425
+ ,364304
+ ,1
+ ,2675
+ ,79
+ ,59017
+ ,355178
+ ,1
+ ,13253
+ ,947
+ ,90829
+ ,357760
+ ,1
+ ,880
+ ,74
+ ,80791
+ ,261216
+ ,0
+ ,2053
+ ,53
+ ,100423
+ ,397144
+ ,1
+ ,1424
+ ,94
+ ,131116
+ ,374943
+ ,0
+ ,4036
+ ,63
+ ,100269
+ ,424898
+ ,0
+ ,3045
+ ,58
+ ,27330
+ ,202055
+ ,1
+ ,5119
+ ,49
+ ,39039
+ ,378525
+ ,1
+ ,1431
+ ,34
+ ,106885
+ ,310768
+ ,1
+ ,554
+ ,11
+ ,79285
+ ,325738
+ ,1
+ ,1975
+ ,35
+ ,118881
+ ,394510
+ ,0
+ ,1765
+ ,20
+ ,77623
+ ,247060
+ ,1
+ ,1012
+ ,47
+ ,114768
+ ,368078
+ ,1
+ ,810
+ ,43
+ ,74015
+ ,236761
+ ,1
+ ,1280
+ ,117
+ ,69465
+ ,312378
+ ,0
+ ,666
+ ,171
+ ,117869
+ ,339836
+ ,1
+ ,1380
+ ,26
+ ,60982
+ ,347385
+ ,0
+ ,4677
+ ,75
+ ,90131
+ ,426280
+ ,1
+ ,876
+ ,59
+ ,138971
+ ,352850
+ ,1
+ ,814
+ ,18
+ ,39625
+ ,301881
+ ,1
+ ,514
+ ,15
+ ,102725
+ ,377516
+ ,0
+ ,5692
+ ,72
+ ,64239
+ ,357312
+ ,1
+ ,3642
+ ,86
+ ,90262
+ ,458343
+ ,1
+ ,540
+ ,14
+ ,103960
+ ,354228
+ ,1
+ ,2099
+ ,64
+ ,106611
+ ,308636
+ ,1
+ ,567
+ ,11
+ ,103345
+ ,386212
+ ,1
+ ,2001
+ ,52
+ ,95551
+ ,393343
+ ,0
+ ,2949
+ ,41
+ ,82903
+ ,378509
+ ,1
+ ,2253
+ ,99
+ ,63593
+ ,452469
+ ,0
+ ,6533
+ ,75
+ ,126910
+ ,364839
+ ,1
+ ,1889
+ ,45
+ ,37527
+ ,358649
+ ,0
+ ,3055
+ ,43
+ ,60247
+ ,376641
+ ,1
+ ,272
+ ,8
+ ,112995
+ ,429112
+ ,0
+ ,1414
+ ,198
+ ,70184
+ ,330546
+ ,1
+ ,2564
+ ,22
+ ,130140
+ ,403560
+ ,0
+ ,1383
+ ,11
+ ,73221
+ ,317892
+ ,0
+ ,1261
+ ,33
+ ,76114
+ ,307528
+ ,1
+ ,975
+ ,23
+ ,90534
+ ,235133
+ ,1
+ ,3366
+ ,80
+ ,108479
+ ,299243
+ ,1
+ ,576
+ ,18
+ ,113761
+ ,314073
+ ,1
+ ,1686
+ ,40
+ ,68696
+ ,368186
+ ,1
+ ,746
+ ,23
+ ,71561
+ ,269661
+ ,0
+ ,3192
+ ,60
+ ,59831
+ ,125390
+ ,0
+ ,2045
+ ,20
+ ,97890
+ ,510834
+ ,1
+ ,5702
+ ,61
+ ,101481
+ ,321896
+ ,0
+ ,1932
+ ,36
+ ,72954
+ ,249898
+ ,1
+ ,936
+ ,30
+ ,67939
+ ,408881
+ ,0
+ ,3437
+ ,47
+ ,48022
+ ,158492
+ ,1
+ ,5131
+ ,71
+ ,86111
+ ,292154
+ ,0
+ ,2397
+ ,14
+ ,74020
+ ,289513
+ ,0
+ ,1389
+ ,9
+ ,57530
+ ,378049
+ ,1
+ ,1503
+ ,39
+ ,56364
+ ,343466
+ ,1
+ ,402
+ ,26
+ ,84990
+ ,332743
+ ,1
+ ,2239
+ ,21
+ ,88590
+ ,442882
+ ,0
+ ,2234
+ ,16
+ ,77200
+ ,214215
+ ,1
+ ,837
+ ,69
+ ,61262
+ ,315688
+ ,1
+ ,10579
+ ,92
+ ,110309
+ ,375195
+ ,1
+ ,875
+ ,14
+ ,67000
+ ,334280
+ ,1
+ ,1585
+ ,107
+ ,93099
+ ,355864
+ ,0
+ ,1659
+ ,29
+ ,107577
+ ,480382
+ ,0
+ ,2647
+ ,37
+ ,62920
+ ,353058
+ ,0
+ ,3294
+ ,23
+ ,75832
+ ,217193
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,94
+ ,7
+ ,60793
+ ,314533
+ ,1
+ ,422
+ ,28
+ ,57935
+ ,318056
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,34
+ ,8
+ ,60630
+ ,314353
+ ,1
+ ,1558
+ ,63
+ ,55637
+ ,369448
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,43
+ ,3
+ ,60887
+ ,312846
+ ,1
+ ,645
+ ,5
+ ,60720
+ ,312075
+ ,1
+ ,316
+ ,9
+ ,60505
+ ,315009
+ ,1
+ ,115
+ ,13
+ ,60945
+ ,318903
+ ,1
+ ,5
+ ,2
+ ,60720
+ ,314887
+ ,1
+ ,897
+ ,5
+ ,60720
+ ,314913
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,389
+ ,14
+ ,58990
+ ,325506
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,1002
+ ,15
+ ,56750
+ ,298568
+ ,1
+ ,36
+ ,3
+ ,60894
+ ,315834
+ ,1
+ ,460
+ ,15
+ ,63346
+ ,329784
+ ,1
+ ,309
+ ,11
+ ,56535
+ ,312878
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,9
+ ,6
+ ,60835
+ ,314987
+ ,1
+ ,271
+ ,2
+ ,60720
+ ,325249
+ ,1
+ ,14
+ ,1
+ ,61016
+ ,315877
+ ,1
+ ,520
+ ,10
+ ,58650
+ ,291650
+ ,1
+ ,1766
+ ,73
+ ,60438
+ ,305959
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,458
+ ,11
+ ,58625
+ ,297765
+ ,1
+ ,20
+ ,3
+ ,60938
+ ,315245
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,98
+ ,2
+ ,61490
+ ,315236
+ ,1
+ ,405
+ ,7
+ ,60845
+ ,336425
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,483
+ ,27
+ ,60830
+ ,306268
+ ,1
+ ,454
+ ,51
+ ,63261
+ ,302187
+ ,1
+ ,47
+ ,3
+ ,60720
+ ,314882
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,757
+ ,19
+ ,45689
+ ,382712
+ ,1
+ ,4655
+ ,393
+ ,60720
+ ,341570
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,36
+ ,4
+ ,61564
+ ,312412
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,203
+ ,9
+ ,61938
+ ,309596
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,126
+ ,10
+ ,60951
+ ,315547
+ ,1
+ ,400
+ ,152
+ ,60720
+ ,313267
+ ,0
+ ,71
+ ,1
+ ,60745
+ ,316176
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,972
+ ,34
+ ,71642
+ ,359335
+ ,0
+ ,531
+ ,10
+ ,71641
+ ,330068
+ ,1
+ ,2461
+ ,57
+ ,55792
+ ,314289
+ ,0
+ ,378
+ ,52
+ ,71873
+ ,297413
+ ,0
+ ,23
+ ,5
+ ,62555
+ ,314806
+ ,0
+ ,638
+ ,14
+ ,60370
+ ,333210
+ ,0
+ ,2300
+ ,29
+ ,64873
+ ,352108
+ ,1
+ ,149
+ ,5
+ ,62041
+ ,313332
+ ,1
+ ,226
+ ,5
+ ,65745
+ ,291787
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,275
+ ,4
+ ,59500
+ ,318745
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,141
+ ,6
+ ,61630
+ ,315366
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,28
+ ,2
+ ,60890
+ ,315688
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,4980
+ ,91
+ ,113521
+ ,409642
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,472
+ ,20
+ ,80045
+ ,269587
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,203
+ ,27
+ ,50804
+ ,300962
+ ,0
+ ,496
+ ,17
+ ,87390
+ ,325479
+ ,0
+ ,10
+ ,2
+ ,61656
+ ,316155
+ ,0
+ ,63
+ ,4
+ ,65688
+ ,318574
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1136
+ ,32
+ ,48522
+ ,343613
+ ,0
+ ,265
+ ,31
+ ,60720
+ ,306948
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,267
+ ,32
+ ,57640
+ ,330059
+ ,1
+ ,474
+ ,20
+ ,61977
+ ,288985
+ ,1
+ ,534
+ ,7
+ ,62620
+ ,304485
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,15
+ ,8
+ ,60831
+ ,315688
+ ,1
+ ,397
+ ,28
+ ,60646
+ ,317736
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,1866
+ ,29
+ ,56225
+ ,322331
+ ,1
+ ,288
+ ,4
+ ,60510
+ ,296656
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,3
+ ,2
+ ,60698
+ ,315354
+ ,1
+ ,468
+ ,21
+ ,60720
+ ,312161
+ ,1
+ ,20
+ ,2
+ ,60805
+ ,315576
+ ,1
+ ,278
+ ,26
+ ,61404
+ ,314922
+ ,1
+ ,61
+ ,14
+ ,60720
+ ,314551
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,192
+ ,4
+ ,65276
+ ,312339
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,317
+ ,9
+ ,63915
+ ,298700
+ ,1
+ ,738
+ ,10
+ ,60720
+ ,321376
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,368
+ ,17
+ ,61686
+ ,303230
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,2
+ ,1
+ ,60743
+ ,315487
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,53
+ ,6
+ ,60349
+ ,315793
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,94
+ ,3
+ ,61360
+ ,312887
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,24
+ ,8
+ ,59818
+ ,315637
+ ,1
+ ,2332
+ ,4
+ ,72680
+ ,324385
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,131
+ ,11
+ ,61808
+ ,308989
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,206
+ ,9
+ ,53110
+ ,296702
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,167
+ ,2
+ ,64245
+ ,307322
+ ,1
+ ,622
+ ,73
+ ,73007
+ ,304376
+ ,1
+ ,2328
+ ,94
+ ,82732
+ ,253588
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,365
+ ,8
+ ,54820
+ ,309560
+ ,1
+ ,364
+ ,35
+ ,47705
+ ,298466
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,226
+ ,12
+ ,72835
+ ,343929
+ ,1
+ ,307
+ ,15
+ ,58856
+ ,331955
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,188
+ ,11
+ ,77655
+ ,381180
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,138
+ ,6
+ ,69817
+ ,331420
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,125
+ ,12
+ ,60798
+ ,310201
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,282
+ ,30
+ ,62452
+ ,320016
+ ,1
+ ,335
+ ,33
+ ,64175
+ ,320398
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1324
+ ,117
+ ,67440
+ ,291841
+ ,1
+ ,176
+ ,28
+ ,68136
+ ,310670
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,249
+ ,72
+ ,56726
+ ,313491
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,333
+ ,13
+ ,70811
+ ,331323
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,601
+ ,6
+ ,60720
+ ,319210
+ ,1
+ ,30
+ ,4
+ ,62045
+ ,318098
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,249
+ ,62
+ ,54323
+ ,292754
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,165
+ ,24
+ ,62841
+ ,325176
+ ,1
+ ,453
+ ,21
+ ,81125
+ ,365959
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,53
+ ,14
+ ,59506
+ ,302409
+ ,0
+ ,382
+ ,21
+ ,59365
+ ,340968
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,30
+ ,4
+ ,60798
+ ,313164
+ ,1
+ ,290
+ ,2
+ ,58790
+ ,301164
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,366
+ ,53
+ ,61808
+ ,344425
+ ,1
+ ,2
+ ,9
+ ,60735
+ ,315394
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,209
+ ,13
+ ,64016
+ ,316647
+ ,1
+ ,384
+ ,22
+ ,54683
+ ,309836
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,365
+ ,83
+ ,87192
+ ,346611
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,49
+ ,8
+ ,64107
+ ,322031
+ ,1
+ ,3
+ ,4
+ ,60761
+ ,315656
+ ,1
+ ,133
+ ,14
+ ,65990
+ ,339445
+ ,1
+ ,32
+ ,1
+ ,59988
+ ,314964
+ ,1
+ ,368
+ ,17
+ ,61167
+ ,297141
+ ,1
+ ,1
+ ,6
+ ,60719
+ ,315372
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,22
+ ,2
+ ,60722
+ ,312502
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,96
+ ,5
+ ,60379
+ ,313729
+ ,1
+ ,1
+ ,2
+ ,60727
+ ,315388
+ ,1
+ ,314
+ ,5
+ ,60720
+ ,315371
+ ,1
+ ,844
+ ,78
+ ,60925
+ ,296139
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,26
+ ,1
+ ,60896
+ ,313880
+ ,1
+ ,125
+ ,13
+ ,59734
+ ,317698
+ ,1
+ ,304
+ ,15
+ ,62969
+ ,295580
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,621
+ ,48
+ ,60720
+ ,308256
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,119
+ ,6
+ ,59118
+ ,303677
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,1595
+ ,17
+ ,60720
+ ,319369
+ ,1
+ ,312
+ ,14
+ ,58598
+ ,318690
+ ,1
+ ,60
+ ,10
+ ,61124
+ ,314049
+ ,1
+ ,587
+ ,12
+ ,59595
+ ,325699
+ ,1
+ ,135
+ ,2
+ ,62065
+ ,314210
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,514
+ ,52
+ ,78780
+ ,322378
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,1
+ ,4
+ ,60722
+ ,315398
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1763
+ ,24
+ ,61600
+ ,308336
+ ,1
+ ,180
+ ,11
+ ,59635
+ ,316386
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,218
+ ,21
+ ,60720
+ ,315553
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,448
+ ,40
+ ,59781
+ ,323361
+ ,1
+ ,227
+ ,9
+ ,76644
+ ,336639
+ ,1
+ ,174
+ ,1
+ ,64820
+ ,307424
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,121
+ ,24
+ ,56178
+ ,295370
+ ,1
+ ,607
+ ,11
+ ,60436
+ ,322340
+ ,1
+ ,2212
+ ,14
+ ,60720
+ ,319864
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,530
+ ,60
+ ,73433
+ ,317291
+ ,1
+ ,571
+ ,80
+ ,41477
+ ,280398
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,78
+ ,16
+ ,62700
+ ,317330
+ ,1
+ ,2489
+ ,40
+ ,67804
+ ,238125
+ ,1
+ ,131
+ ,6
+ ,59661
+ ,327071
+ ,1
+ ,923
+ ,8
+ ,58620
+ ,309038
+ ,1
+ ,72
+ ,3
+ ,60398
+ ,314210
+ ,1
+ ,572
+ ,16
+ ,58580
+ ,307930
+ ,1
+ ,397
+ ,10
+ ,62710
+ ,322327
+ ,1
+ ,450
+ ,8
+ ,59325
+ ,292136
+ ,1
+ ,622
+ ,7
+ ,60950
+ ,263276
+ ,1
+ ,694
+ ,8
+ ,68060
+ ,367655
+ ,0
+ ,3425
+ ,12
+ ,83620
+ ,283910
+ ,1
+ ,562
+ ,13
+ ,58456
+ ,283587
+ ,1
+ ,4917
+ ,42
+ ,52811
+ ,243650
+ ,0
+ ,1442
+ ,118
+ ,121173
+ ,438493
+ ,1
+ ,529
+ ,9
+ ,63870
+ ,296261
+ ,0
+ ,2126
+ ,138
+ ,21001
+ ,230621
+ ,1
+ ,1061
+ ,5
+ ,70415
+ ,304252
+ ,1
+ ,776
+ ,9
+ ,64230
+ ,333505
+ ,1
+ ,611
+ ,8
+ ,59190
+ ,296919
+ ,0
+ ,1526
+ ,25
+ ,69351
+ ,278990
+ ,1
+ ,592
+ ,7
+ ,64270
+ ,276898
+ ,1
+ ,1182
+ ,13
+ ,70694
+ ,327007
+ ,1
+ ,621
+ ,16
+ ,68005
+ ,317046
+ ,1
+ ,989
+ ,11
+ ,58930
+ ,304555
+ ,1
+ ,438
+ ,11
+ ,58320
+ ,298096
+ ,1
+ ,726
+ ,3
+ ,69980
+ ,231861
+ ,1
+ ,1303
+ ,61
+ ,69863
+ ,309422
+ ,0
+ ,7419
+ ,29
+ ,63255
+ ,286963
+ ,0
+ ,1164
+ ,17
+ ,57320
+ ,269753
+ ,0
+ ,3310
+ ,33
+ ,75230
+ ,448243
+ ,1
+ ,1920
+ ,15
+ ,79420
+ ,165404
+ ,1
+ ,965
+ ,3
+ ,73490
+ ,204325
+ ,1
+ ,3256
+ ,66
+ ,35250
+ ,407159
+ ,0
+ ,1135
+ ,17
+ ,62285
+ ,290476
+ ,1
+ ,1270
+ ,26
+ ,69206
+ ,275311
+ ,1
+ ,661
+ ,3
+ ,65920
+ ,246541
+ ,1
+ ,1013
+ ,2
+ ,69770
+ ,253468
+ ,1
+ ,2844
+ ,67
+ ,72683
+ ,240897
+ ,0
+ ,11528
+ ,70
+ ,-14545
+ ,-83265
+ ,1
+ ,6526
+ ,26
+ ,55830
+ ,-42143
+ ,1
+ ,2264
+ ,24
+ ,55174
+ ,272713
+ ,0
+ ,5109
+ ,97
+ ,67038
+ ,215362
+ ,1
+ ,3999
+ ,30
+ ,51252
+ ,42754
+ ,1
+ ,35624
+ ,223
+ ,157278
+ ,306275
+ ,1
+ ,9252
+ ,48
+ ,79510
+ ,253537
+ ,1
+ ,15236
+ ,90
+ ,77440
+ ,372631
+ ,1
+ ,18073
+ ,180
+ ,27284
+ ,-7170)
+ ,dim=c(5
+ ,431)
+ ,dimnames=list(c('Group'
+ ,'Costs'
+ ,'Trades'
+ ,'Dividends'
+ ,'Wealth')
+ ,1:431))
> y <- array(NA,dim=c(5,431),dimnames=list(c('Group','Costs','Trades','Dividends','Wealth'),1:431))
> for (i in 1:dim(x)[1])
+ {
+ for (j in 1:dim(x)[2])
+ {
+ y[i,j] <- as.numeric(x[i,j])
+ }
+ }
> par3 = 'No Linear Trend'
> par2 = 'Do not include Seasonal Dummies'
> par1 = '5'
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Prof. Dr. P. Wessa
> #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #Technical description: Write here your technical program description (don't use hard returns!)
> library(lattice)
> library(lmtest)
Loading required package: zoo
Attaching package: 'zoo'
The following object(s) are masked from package:base :
as.Date.numeric
> n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test
> par1 <- as.numeric(par1)
> x <- t(y)
> k <- length(x[1,])
> n <- length(x[,1])
> x1 <- cbind(x[,par1], x[,1:k!=par1])
> mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1])
> colnames(x1) <- mycolnames #colnames(x)[par1]
> x <- x1
> if (par3 == 'First Differences'){
+ x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep='')))
+ for (i in 1:n-1) {
+ for (j in 1:k) {
+ x2[i,j] <- x[i+1,j] - x[i,j]
+ }
+ }
+ x <- x2
+ }
> if (par2 == 'Include Monthly Dummies'){
+ x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep ='')))
+ for (i in 1:11){
+ x2[seq(i,n,12),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> if (par2 == 'Include Quarterly Dummies'){
+ x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep ='')))
+ for (i in 1:3){
+ x2[seq(i,n,4),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> k <- length(x[1,])
> if (par3 == 'Linear Trend'){
+ x <- cbind(x, c(1:n))
+ colnames(x)[k+1] <- 't'
+ }
> x
Wealth Group Costs Trades Dividends
1 6282154 0 162556 1081 213118
2 4321023 0 29790 309 81767
3 4111912 0 87550 458 153198
4 223193 1 84738 588 -26007
5 1491348 0 54660 302 126942
6 1629616 0 42634 156 157214
7 1398893 1 40949 481 129352
8 1926517 0 45187 353 234817
9 983660 0 37704 452 60448
10 1443586 0 16275 109 47818
11 1073089 1 25830 115 245546
12 984885 1 12679 110 48020
13 1405225 0 18014 239 -1710
14 227132 1 43556 247 32648
15 929118 0 24811 505 95350
16 1071292 1 6575 159 151352
17 638830 1 7123 109 288170
18 856956 0 21950 519 114337
19 992426 0 37597 248 37884
20 444477 1 17821 373 122844
21 857217 0 12988 119 82340
22 711969 0 22330 84 79801
23 702380 1 13326 102 165548
24 358589 1 16189 295 116384
25 297978 1 7146 105 134028
26 585715 1 15824 64 63838
27 657954 0 27664 282 74996
28 209458 1 11920 182 31080
29 786690 1 8568 37 32168
30 439798 1 14416 361 49857
31 688779 0 3369 28 87161
32 574339 0 11819 85 106113
33 741409 0 6984 45 80570
34 597793 0 4519 49 102129
35 644190 1 2220 22 301670
36 377934 1 18562 155 102313
37 640273 1 10327 91 88577
38 697458 0 5336 81 112477
39 550608 0 2365 79 191778
40 207393 1 4069 145 79804
41 301607 1 8636 855 128294
42 345783 1 13718 61 96448
43 501749 1 4525 226 93811
44 379983 1 6869 105 117520
45 387475 1 4628 62 69159
46 377305 0 3689 25 101792
47 370837 0 4891 217 210568
48 430866 0 7489 322 136996
49 469107 1 4901 84 121920
50 194493 1 2284 33 76403
51 530670 0 3160 108 108094
52 518365 0 4150 150 134759
53 491303 0 7285 115 188873
54 527021 0 1134 162 146216
55 233773 0 4658 158 156608
56 405972 1 2384 97 61348
57 652925 1 3748 9 50350
58 446211 1 5371 66 87720
59 341340 1 1285 107 99489
60 387699 0 9327 101 87419
61 493408 0 5565 47 94355
62 146494 1 1528 38 60326
63 414462 0 3122 34 94670
64 364304 0 7561 87 82425
65 355178 1 2675 79 59017
66 357760 1 13253 947 90829
67 261216 1 880 74 80791
68 397144 0 2053 53 100423
69 374943 1 1424 94 131116
70 424898 0 4036 63 100269
71 202055 0 3045 58 27330
72 378525 1 5119 49 39039
73 310768 1 1431 34 106885
74 325738 1 554 11 79285
75 394510 1 1975 35 118881
76 247060 0 1765 20 77623
77 368078 1 1012 47 114768
78 236761 1 810 43 74015
79 312378 1 1280 117 69465
80 339836 0 666 171 117869
81 347385 1 1380 26 60982
82 426280 0 4677 75 90131
83 352850 1 876 59 138971
84 301881 1 814 18 39625
85 377516 1 514 15 102725
86 357312 0 5692 72 64239
87 458343 1 3642 86 90262
88 354228 1 540 14 103960
89 308636 1 2099 64 106611
90 386212 1 567 11 103345
91 393343 1 2001 52 95551
92 378509 0 2949 41 82903
93 452469 1 2253 99 63593
94 364839 0 6533 75 126910
95 358649 1 1889 45 37527
96 376641 0 3055 43 60247
97 429112 1 272 8 112995
98 330546 0 1414 198 70184
99 403560 1 2564 22 130140
100 317892 0 1383 11 73221
101 307528 0 1261 33 76114
102 235133 1 975 23 90534
103 299243 1 3366 80 108479
104 314073 1 576 18 113761
105 368186 1 1686 40 68696
106 269661 1 746 23 71561
107 125390 0 3192 60 59831
108 510834 0 2045 20 97890
109 321896 1 5702 61 101481
110 249898 0 1932 36 72954
111 408881 1 936 30 67939
112 158492 0 3437 47 48022
113 292154 1 5131 71 86111
114 289513 0 2397 14 74020
115 378049 0 1389 9 57530
116 343466 1 1503 39 56364
117 332743 1 402 26 84990
118 442882 1 2239 21 88590
119 214215 0 2234 16 77200
120 315688 1 837 69 61262
121 375195 1 10579 92 110309
122 334280 1 875 14 67000
123 355864 1 1585 107 93099
124 480382 0 1659 29 107577
125 353058 0 2647 37 62920
126 217193 0 3294 23 75832
127 315380 1 0 0 60720
128 314533 1 94 7 60793
129 318056 1 422 28 57935
130 315380 1 0 0 60720
131 314353 1 34 8 60630
132 369448 1 1558 63 55637
133 315380 1 0 0 60720
134 312846 1 43 3 60887
135 312075 1 645 5 60720
136 315009 1 316 9 60505
137 318903 1 115 13 60945
138 314887 1 5 2 60720
139 314913 1 897 5 60720
140 315380 1 0 0 60720
141 325506 1 389 14 58990
142 315380 1 0 0 60720
143 298568 1 1002 15 56750
144 315834 1 36 3 60894
145 329784 1 460 15 63346
146 312878 1 309 11 56535
147 315380 1 0 0 60720
148 314987 1 9 6 60835
149 325249 1 271 2 60720
150 315877 1 14 1 61016
151 291650 1 520 10 58650
152 305959 1 1766 73 60438
153 315380 1 0 0 60720
154 297765 1 458 11 58625
155 315245 1 20 3 60938
156 315380 1 0 0 60720
157 315380 1 0 0 60720
158 315236 1 98 2 61490
159 336425 1 405 7 60845
160 315380 1 0 0 60720
161 315380 1 0 0 60720
162 315380 1 0 0 60720
163 315380 1 0 0 60720
164 306268 1 483 27 60830
165 302187 1 454 51 63261
166 314882 1 47 3 60720
167 315380 1 0 0 60720
168 382712 1 757 19 45689
169 341570 1 4655 393 60720
170 315380 1 0 0 60720
171 315380 1 0 0 60720
172 312412 1 36 4 61564
173 315380 1 0 0 60720
174 309596 1 203 9 61938
175 315380 1 0 0 60720
176 315547 1 126 10 60951
177 313267 1 400 152 60720
178 316176 0 71 1 60745
179 315380 1 0 0 60720
180 315380 1 0 0 60720
181 359335 1 972 34 71642
182 330068 0 531 10 71641
183 314289 1 2461 57 55792
184 297413 0 378 52 71873
185 314806 0 23 5 62555
186 333210 0 638 14 60370
187 352108 0 2300 29 64873
188 313332 1 149 5 62041
189 291787 1 226 5 65745
190 315380 1 0 0 60720
191 318745 1 275 4 59500
192 315380 1 0 0 60720
193 315366 1 141 6 61630
194 315380 1 0 0 60720
195 315688 1 28 2 60890
196 315380 0 0 0 60720
197 409642 0 4980 91 113521
198 315380 0 0 0 60720
199 315380 0 0 0 60720
200 269587 0 472 20 80045
201 315380 0 0 0 60720
202 315380 0 0 0 60720
203 315380 0 0 0 60720
204 300962 0 203 27 50804
205 325479 0 496 17 87390
206 316155 0 10 2 61656
207 318574 0 63 4 65688
208 315380 0 0 0 60720
209 343613 0 1136 32 48522
210 306948 0 265 31 60720
211 315380 1 0 0 60720
212 315380 1 0 0 60720
213 330059 1 267 32 57640
214 288985 1 474 20 61977
215 304485 1 534 7 62620
216 315380 1 0 0 60720
217 315688 1 15 8 60831
218 317736 1 397 28 60646
219 315380 1 0 0 60720
220 322331 1 1866 29 56225
221 296656 1 288 4 60510
222 315380 1 0 0 60720
223 315354 1 3 2 60698
224 312161 1 468 21 60720
225 315576 1 20 2 60805
226 314922 1 278 26 61404
227 314551 1 61 14 60720
228 315380 1 0 0 60720
229 312339 1 192 4 65276
230 315380 1 0 0 60720
231 298700 1 317 9 63915
232 321376 1 738 10 60720
233 315380 1 0 0 60720
234 303230 1 368 17 61686
235 315380 1 0 0 60720
236 315487 1 2 1 60743
237 315380 1 0 0 60720
238 315793 1 53 6 60349
239 315380 1 0 0 60720
240 315380 1 0 0 60720
241 315380 1 0 0 60720
242 312887 1 94 3 61360
243 315380 1 0 0 60720
244 315637 1 24 8 59818
245 324385 1 2332 4 72680
246 315380 0 0 0 60720
247 315380 0 0 0 60720
248 308989 1 131 11 61808
249 315380 0 0 0 60720
250 315380 0 0 0 60720
251 296702 1 206 9 53110
252 315380 0 0 0 60720
253 307322 1 167 2 64245
254 304376 1 622 73 73007
255 253588 1 2328 94 82732
256 315380 1 0 0 60720
257 309560 1 365 8 54820
258 298466 1 364 35 47705
259 315380 0 0 0 60720
260 315380 1 0 0 60720
261 315380 0 0 0 60720
262 315380 0 0 0 60720
263 343929 1 226 12 72835
264 331955 1 307 15 58856
265 315380 0 0 0 60720
266 315380 1 0 0 60720
267 315380 0 0 0 60720
268 381180 1 188 11 77655
269 315380 0 0 0 60720
270 331420 1 138 6 69817
271 315380 0 0 0 60720
272 315380 0 0 0 60720
273 315380 0 0 0 60720
274 310201 1 125 12 60798
275 315380 1 0 0 60720
276 320016 1 282 30 62452
277 320398 1 335 33 64175
278 315380 1 0 0 60720
279 291841 0 1324 117 67440
280 310670 1 176 28 68136
281 315380 0 0 0 60720
282 315380 1 0 0 60720
283 313491 1 249 72 56726
284 315380 1 0 0 60720
285 331323 1 333 13 70811
286 315380 1 0 0 60720
287 319210 0 601 6 60720
288 318098 1 30 4 62045
289 315380 0 0 0 60720
290 292754 1 249 62 54323
291 315380 0 0 0 60720
292 325176 1 165 24 62841
293 365959 1 453 21 81125
294 315380 1 0 0 60720
295 302409 1 53 14 59506
296 340968 0 382 21 59365
297 315380 1 0 0 60720
298 315380 0 0 0 60720
299 315380 1 0 0 60720
300 315380 1 0 0 60720
301 313164 0 30 4 60798
302 301164 1 290 2 58790
303 315380 1 0 0 60720
304 315380 0 0 0 60720
305 344425 1 366 53 61808
306 315394 1 2 9 60735
307 315380 0 0 0 60720
308 316647 0 209 13 64016
309 309836 1 384 22 54683
310 315380 0 0 0 60720
311 315380 0 0 0 60720
312 346611 1 365 83 87192
313 315380 0 0 0 60720
314 322031 0 49 8 64107
315 315656 1 3 4 60761
316 339445 1 133 14 65990
317 314964 1 32 1 59988
318 297141 1 368 17 61167
319 315372 1 1 6 60719
320 315380 1 0 0 60720
321 315380 1 0 0 60720
322 315380 1 0 0 60720
323 315380 1 0 0 60720
324 315380 0 0 0 60720
325 315380 1 0 0 60720
326 312502 1 22 2 60722
327 315380 0 0 0 60720
328 315380 0 0 0 60720
329 315380 1 0 0 60720
330 315380 1 0 0 60720
331 315380 1 0 0 60720
332 315380 1 0 0 60720
333 315380 1 0 0 60720
334 313729 1 96 5 60379
335 315388 1 1 2 60727
336 315371 1 314 5 60720
337 296139 1 844 78 60925
338 315380 1 0 0 60720
339 313880 1 26 1 60896
340 317698 1 125 13 59734
341 295580 1 304 15 62969
342 315380 1 0 0 60720
343 315380 1 0 0 60720
344 315380 1 0 0 60720
345 308256 1 621 48 60720
346 315380 1 0 0 60720
347 303677 1 119 6 59118
348 315380 1 0 0 60720
349 315380 1 0 0 60720
350 319369 1 1595 17 60720
351 318690 1 312 14 58598
352 314049 1 60 10 61124
353 325699 1 587 12 59595
354 314210 1 135 2 62065
355 315380 1 0 0 60720
356 315380 1 0 0 60720
357 322378 1 514 52 78780
358 315380 1 0 0 60720
359 315380 1 0 0 60720
360 315380 1 0 0 60720
361 315398 1 1 4 60722
362 315380 1 0 0 60720
363 315380 1 0 0 60720
364 308336 0 1763 24 61600
365 316386 1 180 11 59635
366 315380 1 0 0 60720
367 315380 1 0 0 60720
368 315380 1 0 0 60720
369 315380 1 0 0 60720
370 315553 1 218 21 60720
371 315380 1 0 0 60720
372 323361 1 448 40 59781
373 336639 1 227 9 76644
374 307424 1 174 1 64820
375 315380 1 0 0 60720
376 315380 1 0 0 60720
377 295370 1 121 24 56178
378 322340 1 607 11 60436
379 319864 1 2212 14 60720
380 315380 1 0 0 60720
381 315380 1 0 0 60720
382 317291 1 530 60 73433
383 280398 1 571 80 41477
384 315380 1 0 0 60720
385 317330 1 78 16 62700
386 238125 1 2489 40 67804
387 327071 1 131 6 59661
388 309038 1 923 8 58620
389 314210 1 72 3 60398
390 307930 1 572 16 58580
391 322327 1 397 10 62710
392 292136 1 450 8 59325
393 263276 1 622 7 60950
394 367655 1 694 8 68060
395 283910 0 3425 12 83620
396 283587 1 562 13 58456
397 243650 1 4917 42 52811
398 438493 0 1442 118 121173
399 296261 1 529 9 63870
400 230621 0 2126 138 21001
401 304252 1 1061 5 70415
402 333505 1 776 9 64230
403 296919 1 611 8 59190
404 278990 0 1526 25 69351
405 276898 1 592 7 64270
406 327007 1 1182 13 70694
407 317046 1 621 16 68005
408 304555 1 989 11 58930
409 298096 1 438 11 58320
410 231861 1 726 3 69980
411 309422 1 1303 61 69863
412 286963 0 7419 29 63255
413 269753 0 1164 17 57320
414 448243 0 3310 33 75230
415 165404 1 1920 15 79420
416 204325 1 965 3 73490
417 407159 1 3256 66 35250
418 290476 0 1135 17 62285
419 275311 1 1270 26 69206
420 246541 1 661 3 65920
421 253468 1 1013 2 69770
422 240897 1 2844 67 72683
423 -83265 0 11528 70 -14545
424 -42143 1 6526 26 55830
425 272713 1 2264 24 55174
426 215362 0 5109 97 67038
427 42754 1 3999 30 51252
428 306275 1 35624 223 157278
429 253537 1 9252 48 79510
430 372631 1 15236 90 77440
431 -7170 1 18073 180 27284
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) Group Costs Trades Dividends
205978.971 -55485.819 29.811 -407.150 1.926
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-2163938 -14344 35201 47931 3195284
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 2.060e+05 3.854e+04 5.345 1.48e-07 ***
Group -5.549e+04 2.784e+04 -1.993 0.0469 *
Costs 2.981e+01 1.578e+00 18.889 < 2e-16 ***
Trades -4.072e+02 1.663e+02 -2.448 0.0148 *
Dividends 1.926e+00 4.105e-01 4.692 3.65e-06 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 253400 on 426 degrees of freedom
Multiple R-squared: 0.6569, Adjusted R-squared: 0.6537
F-statistic: 203.9 on 4 and 426 DF, p-value: < 2.2e-16
> if (n > n25) {
+ kp3 <- k + 3
+ nmkm3 <- n - k - 3
+ gqarr <- array(NA, dim=c(nmkm3-kp3+1,3))
+ numgqtests <- 0
+ numsignificant1 <- 0
+ numsignificant5 <- 0
+ numsignificant10 <- 0
+ for (mypoint in kp3:nmkm3) {
+ j <- 0
+ numgqtests <- numgqtests + 1
+ for (myalt in c('greater', 'two.sided', 'less')) {
+ j <- j + 1
+ gqarr[mypoint-kp3+1,j] <- gqtest(mylm, point=mypoint, alternative=myalt)$p.value
+ }
+ if (gqarr[mypoint-kp3+1,2] < 0.01) numsignificant1 <- numsignificant1 + 1
+ if (gqarr[mypoint-kp3+1,2] < 0.05) numsignificant5 <- numsignificant5 + 1
+ if (gqarr[mypoint-kp3+1,2] < 0.10) numsignificant10 <- numsignificant10 + 1
+ }
+ gqarr
+ }
[,1] [,2] [,3]
[1,] 1.0000000 6.580461e-56 3.290230e-56
[2,] 1.0000000 9.697449e-93 4.848724e-93
[3,] 1.0000000 3.784523e-102 1.892262e-102
[4,] 1.0000000 2.398605e-104 1.199303e-104
[5,] 1.0000000 2.524261e-116 1.262131e-116
[6,] 1.0000000 2.796176e-132 1.398088e-132
[7,] 1.0000000 1.296946e-134 6.484732e-135
[8,] 1.0000000 1.043990e-142 5.219952e-143
[9,] 1.0000000 5.532265e-153 2.766133e-153
[10,] 1.0000000 1.188565e-152 5.942826e-153
[11,] 1.0000000 1.597306e-157 7.986528e-158
[12,] 1.0000000 1.622078e-164 8.110390e-165
[13,] 1.0000000 2.174270e-163 1.087135e-163
[14,] 1.0000000 4.837755e-169 2.418878e-169
[15,] 1.0000000 1.876296e-172 9.381478e-173
[16,] 1.0000000 5.172766e-173 2.586383e-173
[17,] 1.0000000 5.308774e-172 2.654387e-172
[18,] 1.0000000 1.870383e-171 9.351913e-172
[19,] 1.0000000 4.297350e-173 2.148675e-173
[20,] 1.0000000 1.219015e-177 6.095076e-178
[21,] 1.0000000 5.159968e-177 2.579984e-177
[22,] 1.0000000 6.552292e-190 3.276146e-190
[23,] 1.0000000 1.649405e-190 8.247026e-191
[24,] 1.0000000 1.876207e-194 9.381033e-195
[25,] 1.0000000 1.330203e-196 6.651013e-197
[26,] 1.0000000 3.554703e-204 1.777351e-204
[27,] 1.0000000 2.959186e-206 1.479593e-206
[28,] 1.0000000 1.866239e-205 9.331196e-206
[29,] 1.0000000 1.143810e-204 5.719049e-205
[30,] 1.0000000 1.635797e-210 8.178985e-211
[31,] 1.0000000 2.790585e-216 1.395292e-216
[32,] 1.0000000 2.423065e-216 1.211532e-216
[33,] 1.0000000 3.344570e-216 1.672285e-216
[34,] 1.0000000 6.743693e-216 3.371847e-216
[35,] 1.0000000 5.618652e-215 2.809326e-215
[36,] 1.0000000 6.314230e-216 3.157115e-216
[37,] 1.0000000 7.510434e-215 3.755217e-215
[38,] 1.0000000 1.254756e-214 6.273779e-215
[39,] 1.0000000 4.183044e-214 2.091522e-214
[40,] 1.0000000 6.006534e-217 3.003267e-217
[41,] 1.0000000 1.506022e-216 7.530108e-217
[42,] 1.0000000 4.509424e-216 2.254712e-216
[43,] 1.0000000 3.030670e-216 1.515335e-216
[44,] 1.0000000 6.898805e-217 3.449402e-217
[45,] 1.0000000 8.546720e-217 4.273360e-217
[46,] 1.0000000 1.213618e-216 6.068091e-217
[47,] 1.0000000 4.108265e-216 2.054132e-216
[48,] 1.0000000 1.821504e-219 9.107519e-220
[49,] 1.0000000 1.134672e-219 5.673360e-220
[50,] 1.0000000 2.818042e-232 1.409021e-232
[51,] 1.0000000 1.096219e-232 5.481096e-233
[52,] 1.0000000 8.301183e-232 4.150591e-232
[53,] 1.0000000 8.007338e-232 4.003669e-232
[54,] 1.0000000 2.876362e-233 1.438181e-233
[55,] 1.0000000 1.474903e-234 7.374514e-235
[56,] 1.0000000 5.924023e-234 2.962012e-234
[57,] 1.0000000 1.467914e-233 7.339571e-234
[58,] 1.0000000 4.940991e-233 2.470495e-233
[59,] 1.0000000 5.649985e-232 2.824993e-232
[60,] 1.0000000 1.276709e-231 6.383544e-232
[61,] 1.0000000 1.167663e-230 5.838316e-231
[62,] 1.0000000 7.212219e-230 3.606110e-230
[63,] 1.0000000 2.128623e-229 1.064312e-229
[64,] 1.0000000 1.713134e-228 8.565670e-229
[65,] 1.0000000 1.227108e-229 6.135539e-230
[66,] 1.0000000 4.201118e-229 2.100559e-229
[67,] 1.0000000 3.799274e-228 1.899637e-228
[68,] 1.0000000 4.250804e-227 2.125402e-227
[69,] 1.0000000 7.622571e-227 3.811286e-227
[70,] 1.0000000 6.129166e-226 3.064583e-226
[71,] 1.0000000 7.428973e-226 3.714486e-226
[72,] 1.0000000 6.530709e-225 3.265355e-225
[73,] 1.0000000 9.533083e-225 4.766542e-225
[74,] 1.0000000 4.964354e-224 2.482177e-224
[75,] 1.0000000 6.469603e-224 3.234801e-224
[76,] 1.0000000 4.289373e-224 2.144687e-224
[77,] 1.0000000 2.451482e-223 1.225741e-223
[78,] 1.0000000 2.370248e-222 1.185124e-222
[79,] 1.0000000 3.987638e-222 1.993819e-222
[80,] 1.0000000 1.416318e-222 7.081592e-223
[81,] 1.0000000 1.039626e-221 5.198128e-222
[82,] 1.0000000 2.049985e-221 1.024992e-221
[83,] 1.0000000 2.062041e-220 1.031020e-220
[84,] 1.0000000 1.847115e-219 9.235577e-220
[85,] 1.0000000 1.133793e-218 5.668964e-219
[86,] 1.0000000 2.454126e-220 1.227063e-220
[87,] 1.0000000 1.126815e-219 5.634075e-220
[88,] 1.0000000 3.228381e-220 1.614190e-220
[89,] 1.0000000 3.360509e-220 1.680254e-220
[90,] 1.0000000 3.273129e-219 1.636565e-219
[91,] 1.0000000 3.889547e-218 1.944773e-218
[92,] 1.0000000 3.315415e-217 1.657708e-217
[93,] 1.0000000 3.596537e-216 1.798269e-216
[94,] 1.0000000 3.249622e-215 1.624811e-215
[95,] 1.0000000 2.372608e-216 1.186304e-216
[96,] 1.0000000 1.928803e-216 9.644017e-217
[97,] 1.0000000 2.834180e-217 1.417090e-217
[98,] 1.0000000 1.308202e-216 6.541011e-217
[99,] 1.0000000 4.425533e-216 2.212766e-216
[100,] 1.0000000 1.770637e-218 8.853185e-219
[101,] 1.0000000 1.113593e-219 5.567967e-220
[102,] 1.0000000 8.177657e-219 4.088829e-219
[103,] 1.0000000 1.414773e-218 7.073863e-219
[104,] 1.0000000 1.055940e-218 5.279699e-219
[105,] 1.0000000 5.001559e-219 2.500779e-219
[106,] 1.0000000 3.706097e-218 1.853048e-218
[107,] 1.0000000 2.808432e-217 1.404216e-217
[108,] 1.0000000 2.568869e-217 1.284434e-217
[109,] 1.0000000 9.624234e-217 4.812117e-217
[110,] 1.0000000 6.720738e-216 3.360369e-216
[111,] 1.0000000 6.046519e-216 3.023259e-216
[112,] 1.0000000 4.683900e-217 2.341950e-217
[113,] 1.0000000 4.786938e-216 2.393469e-216
[114,] 1.0000000 2.798069e-215 1.399034e-215
[115,] 1.0000000 2.695036e-214 1.347518e-214
[116,] 1.0000000 2.225628e-213 1.112814e-213
[117,] 1.0000000 4.454487e-213 2.227244e-213
[118,] 1.0000000 1.515581e-212 7.577903e-213
[119,] 1.0000000 3.594020e-213 1.797010e-213
[120,] 1.0000000 3.632179e-212 1.816090e-212
[121,] 1.0000000 3.700289e-211 1.850145e-211
[122,] 1.0000000 3.446286e-210 1.723143e-210
[123,] 1.0000000 3.499927e-209 1.749963e-209
[124,] 1.0000000 3.566552e-208 1.783276e-208
[125,] 1.0000000 3.410846e-208 1.705423e-208
[126,] 1.0000000 3.497851e-207 1.748926e-207
[127,] 1.0000000 3.633494e-206 1.816747e-206
[128,] 1.0000000 3.805841e-205 1.902921e-205
[129,] 1.0000000 3.893032e-204 1.946516e-204
[130,] 1.0000000 3.898673e-203 1.949337e-203
[131,] 1.0000000 3.988873e-202 1.994436e-202
[132,] 1.0000000 4.012666e-201 2.006333e-201
[133,] 1.0000000 4.081896e-200 2.040948e-200
[134,] 1.0000000 3.488969e-199 1.744485e-199
[135,] 1.0000000 3.538448e-198 1.769224e-198
[136,] 1.0000000 3.719160e-197 1.859580e-197
[137,] 1.0000000 3.751240e-196 1.875620e-196
[138,] 1.0000000 3.462981e-195 1.731491e-195
[139,] 1.0000000 3.246167e-194 1.623083e-194
[140,] 1.0000000 3.253262e-193 1.626631e-193
[141,] 1.0000000 3.259631e-192 1.629815e-192
[142,] 1.0000000 2.928456e-191 1.464228e-191
[143,] 1.0000000 2.914539e-190 1.457269e-190
[144,] 1.0000000 2.941399e-189 1.470699e-189
[145,] 1.0000000 2.997652e-188 1.498826e-188
[146,] 1.0000000 2.951277e-187 1.475639e-187
[147,] 1.0000000 2.990784e-186 1.495392e-186
[148,] 1.0000000 2.933777e-185 1.466889e-185
[149,] 1.0000000 2.860569e-184 1.430285e-184
[150,] 1.0000000 2.780129e-183 1.390065e-183
[151,] 1.0000000 2.716643e-182 1.358322e-182
[152,] 1.0000000 1.942374e-181 9.711868e-182
[153,] 1.0000000 1.870041e-180 9.350207e-181
[154,] 1.0000000 1.794144e-179 8.970720e-180
[155,] 1.0000000 1.715259e-178 8.576296e-179
[156,] 1.0000000 1.633987e-177 8.169935e-178
[157,] 1.0000000 1.589681e-176 7.948406e-177
[158,] 1.0000000 1.446267e-175 7.231333e-176
[159,] 1.0000000 1.365745e-174 6.828725e-175
[160,] 1.0000000 1.282240e-173 6.411199e-174
[161,] 1.0000000 1.362101e-174 6.810505e-175
[162,] 1.0000000 4.685075e-174 2.342537e-174
[163,] 1.0000000 4.451372e-173 2.225686e-173
[164,] 1.0000000 4.213031e-172 2.106516e-172
[165,] 1.0000000 4.033156e-171 2.016578e-171
[166,] 1.0000000 3.787567e-170 1.893784e-170
[167,] 1.0000000 3.617722e-169 1.808861e-169
[168,] 1.0000000 3.370592e-168 1.685296e-168
[169,] 1.0000000 3.117500e-167 1.558750e-167
[170,] 1.0000000 2.696623e-166 1.348311e-166
[171,] 1.0000000 2.570902e-165 1.285451e-165
[172,] 1.0000000 2.357167e-164 1.178583e-164
[173,] 1.0000000 2.152478e-163 1.076239e-163
[174,] 1.0000000 1.404651e-162 7.023257e-163
[175,] 1.0000000 1.304927e-161 6.524635e-162
[176,] 1.0000000 7.202553e-161 3.601276e-161
[177,] 1.0000000 4.132589e-160 2.066295e-160
[178,] 1.0000000 3.843771e-159 1.921886e-159
[179,] 1.0000000 2.845878e-158 1.422939e-158
[180,] 1.0000000 1.165707e-157 5.828537e-158
[181,] 1.0000000 1.052981e-156 5.264906e-157
[182,] 1.0000000 7.604224e-156 3.802112e-156
[183,] 1.0000000 6.710470e-155 3.355235e-155
[184,] 1.0000000 5.549776e-154 2.774888e-154
[185,] 1.0000000 4.857509e-153 2.428754e-153
[186,] 1.0000000 4.259409e-152 2.129705e-152
[187,] 1.0000000 3.697245e-151 1.848623e-151
[188,] 1.0000000 3.193753e-150 1.596877e-150
[189,] 1.0000000 2.838005e-149 1.419002e-149
[190,] 1.0000000 2.157850e-148 1.078925e-148
[191,] 1.0000000 1.901158e-147 9.505790e-148
[192,] 1.0000000 1.667024e-146 8.335121e-147
[193,] 1.0000000 1.406492e-146 7.032462e-147
[194,] 1.0000000 1.238510e-145 6.192551e-146
[195,] 1.0000000 1.085448e-144 5.427242e-145
[196,] 1.0000000 9.468132e-144 4.734066e-144
[197,] 1.0000000 7.679675e-143 3.839838e-143
[198,] 1.0000000 2.923387e-142 1.461694e-142
[199,] 1.0000000 2.528231e-141 1.264115e-141
[200,] 1.0000000 2.134462e-140 1.067231e-140
[201,] 1.0000000 1.828339e-139 9.141693e-140
[202,] 1.0000000 3.151406e-139 1.575703e-139
[203,] 1.0000000 2.688550e-138 1.344275e-138
[204,] 1.0000000 2.208901e-137 1.104450e-137
[205,] 1.0000000 1.807269e-136 9.036346e-137
[206,] 1.0000000 1.060026e-135 5.300132e-136
[207,] 1.0000000 8.074529e-135 4.037265e-135
[208,] 1.0000000 6.686373e-134 3.343186e-134
[209,] 1.0000000 5.390174e-133 2.695087e-133
[210,] 1.0000000 4.311870e-132 2.155935e-132
[211,] 1.0000000 3.297804e-131 1.648902e-131
[212,] 1.0000000 2.625235e-130 1.312617e-130
[213,] 1.0000000 1.227528e-129 6.137642e-130
[214,] 1.0000000 9.775420e-129 4.887710e-129
[215,] 1.0000000 7.697290e-128 3.848645e-128
[216,] 1.0000000 6.028882e-127 3.014441e-127
[217,] 1.0000000 4.688149e-126 2.344075e-126
[218,] 1.0000000 3.639727e-125 1.819864e-125
[219,] 1.0000000 2.800355e-124 1.400178e-124
[220,] 1.0000000 2.151705e-123 1.075853e-123
[221,] 1.0000000 1.650697e-122 8.253484e-123
[222,] 1.0000000 1.278360e-121 6.391799e-122
[223,] 1.0000000 9.723369e-121 4.861685e-121
[224,] 1.0000000 7.140998e-120 3.570499e-120
[225,] 1.0000000 4.900579e-119 2.450290e-119
[226,] 1.0000000 3.680856e-118 1.840428e-118
[227,] 1.0000000 2.809766e-117 1.404883e-117
[228,] 1.0000000 2.091925e-116 1.045962e-116
[229,] 1.0000000 1.549537e-115 7.747684e-116
[230,] 1.0000000 1.143310e-114 5.716552e-115
[231,] 1.0000000 8.307959e-114 4.153980e-114
[232,] 1.0000000 6.074670e-113 3.037335e-113
[233,] 1.0000000 4.421512e-112 2.210756e-112
[234,] 1.0000000 3.203556e-111 1.601778e-111
[235,] 1.0000000 2.341111e-110 1.170555e-110
[236,] 1.0000000 1.680760e-109 8.403800e-110
[237,] 1.0000000 1.180110e-108 5.900549e-109
[238,] 1.0000000 8.518591e-108 4.259296e-108
[239,] 1.0000000 6.190722e-107 3.095361e-107
[240,] 1.0000000 4.477218e-106 2.238609e-106
[241,] 1.0000000 3.204745e-105 1.602372e-105
[242,] 1.0000000 2.295426e-104 1.147713e-104
[243,] 1.0000000 1.636100e-103 8.180498e-104
[244,] 1.0000000 1.113834e-102 5.569168e-103
[245,] 1.0000000 7.863750e-102 3.931875e-102
[246,] 1.0000000 5.443887e-101 2.721943e-101
[247,] 1.0000000 3.127684e-100 1.563842e-100
[248,] 1.0000000 2.535970e-100 1.267985e-100
[249,] 1.0000000 1.735667e-99 8.678335e-100
[250,] 1.0000000 1.069487e-98 5.347437e-99
[251,] 1.0000000 5.986435e-98 2.993218e-98
[252,] 1.0000000 4.144425e-97 2.072213e-97
[253,] 1.0000000 2.785412e-96 1.392706e-96
[254,] 1.0000000 1.909385e-95 9.546924e-96
[255,] 1.0000000 1.302166e-94 6.510832e-95
[256,] 1.0000000 8.718921e-94 4.359461e-94
[257,] 1.0000000 4.501405e-93 2.250703e-93
[258,] 1.0000000 3.029487e-92 1.514744e-92
[259,] 1.0000000 1.979688e-91 9.898441e-92
[260,] 1.0000000 1.318994e-90 6.594972e-91
[261,] 1.0000000 6.980669e-90 3.490335e-90
[262,] 1.0000000 4.609886e-89 2.304943e-89
[263,] 1.0000000 3.022311e-88 1.511155e-88
[264,] 1.0000000 1.975367e-87 9.876836e-88
[265,] 1.0000000 1.284132e-86 6.420662e-87
[266,] 1.0000000 8.302334e-86 4.151167e-86
[267,] 1.0000000 5.314259e-85 2.657129e-85
[268,] 1.0000000 3.325474e-84 1.662737e-84
[269,] 1.0000000 2.061703e-83 1.030852e-83
[270,] 1.0000000 1.292498e-82 6.462492e-83
[271,] 1.0000000 7.968747e-82 3.984373e-82
[272,] 1.0000000 3.903588e-81 1.951794e-81
[273,] 1.0000000 2.318534e-80 1.159267e-80
[274,] 1.0000000 1.444618e-79 7.223088e-80
[275,] 1.0000000 8.735332e-79 4.367666e-79
[276,] 1.0000000 5.099288e-78 2.549644e-78
[277,] 1.0000000 3.053613e-77 1.526806e-77
[278,] 1.0000000 1.866870e-76 9.334352e-77
[279,] 1.0000000 1.106537e-75 5.532685e-76
[280,] 1.0000000 6.429784e-75 3.214892e-75
[281,] 1.0000000 3.788418e-74 1.894209e-74
[282,] 1.0000000 2.263444e-73 1.131722e-73
[283,] 1.0000000 1.342261e-72 6.711306e-73
[284,] 1.0000000 7.930846e-72 3.965423e-72
[285,] 1.0000000 4.477237e-71 2.238618e-71
[286,] 1.0000000 2.586387e-70 1.293193e-70
[287,] 1.0000000 1.470859e-69 7.354296e-70
[288,] 1.0000000 8.539692e-69 4.269846e-69
[289,] 1.0000000 3.765639e-68 1.882819e-68
[290,] 1.0000000 2.111891e-67 1.055945e-67
[291,] 1.0000000 1.198781e-66 5.993906e-67
[292,] 1.0000000 6.650270e-66 3.325135e-66
[293,] 1.0000000 3.668631e-65 1.834315e-65
[294,] 1.0000000 2.060884e-64 1.030442e-64
[295,] 1.0000000 1.151110e-63 5.755551e-64
[296,] 1.0000000 6.245272e-63 3.122636e-63
[297,] 1.0000000 3.423197e-62 1.711598e-62
[298,] 1.0000000 1.477790e-61 7.388951e-62
[299,] 1.0000000 7.918983e-61 3.959492e-61
[300,] 1.0000000 4.276271e-60 2.138136e-60
[301,] 1.0000000 2.326807e-59 1.163404e-59
[302,] 1.0000000 1.125418e-58 5.627090e-59
[303,] 1.0000000 5.968006e-58 2.984003e-58
[304,] 1.0000000 3.139574e-57 1.569787e-57
[305,] 1.0000000 1.332911e-56 6.664553e-57
[306,] 1.0000000 6.903447e-56 3.451724e-56
[307,] 1.0000000 3.591003e-55 1.795501e-55
[308,] 1.0000000 1.829546e-54 9.147728e-55
[309,] 1.0000000 8.782550e-54 4.391275e-54
[310,] 1.0000000 4.378783e-53 2.189392e-53
[311,] 1.0000000 2.248769e-52 1.124385e-52
[312,] 1.0000000 1.121044e-51 5.605221e-52
[313,] 1.0000000 5.535517e-51 2.767759e-51
[314,] 1.0000000 2.716167e-50 1.358084e-50
[315,] 1.0000000 1.324332e-49 6.621658e-50
[316,] 1.0000000 6.415903e-49 3.207952e-49
[317,] 1.0000000 3.093586e-48 1.546793e-48
[318,] 1.0000000 1.479997e-47 7.399987e-48
[319,] 1.0000000 7.119930e-47 3.559965e-47
[320,] 1.0000000 3.349358e-46 1.674679e-46
[321,] 1.0000000 1.552927e-45 7.764635e-46
[322,] 1.0000000 7.241226e-45 3.620613e-45
[323,] 1.0000000 3.353885e-44 1.676943e-44
[324,] 1.0000000 1.542886e-43 7.714428e-44
[325,] 1.0000000 7.049247e-43 3.524623e-43
[326,] 1.0000000 3.198514e-42 1.599257e-42
[327,] 1.0000000 1.446390e-41 7.231949e-42
[328,] 1.0000000 6.483087e-41 3.241544e-41
[329,] 1.0000000 2.854850e-40 1.427425e-40
[330,] 1.0000000 1.283172e-39 6.415861e-40
[331,] 1.0000000 5.624803e-39 2.812401e-39
[332,] 1.0000000 2.468498e-38 1.234249e-38
[333,] 1.0000000 1.048704e-37 5.243518e-38
[334,] 1.0000000 4.633853e-37 2.316926e-37
[335,] 1.0000000 1.975253e-36 9.876264e-37
[336,] 1.0000000 8.356319e-36 4.178159e-36
[337,] 1.0000000 3.508191e-35 1.754095e-35
[338,] 1.0000000 1.516141e-34 7.580706e-35
[339,] 1.0000000 6.271530e-34 3.135765e-34
[340,] 1.0000000 2.645058e-33 1.322529e-33
[341,] 1.0000000 1.077222e-32 5.386108e-33
[342,] 1.0000000 4.351832e-32 2.175916e-32
[343,] 1.0000000 1.580431e-31 7.902153e-32
[344,] 1.0000000 5.981313e-31 2.990656e-31
[345,] 1.0000000 2.401539e-30 1.200770e-30
[346,] 1.0000000 8.459500e-30 4.229750e-30
[347,] 1.0000000 3.340945e-29 1.670472e-29
[348,] 1.0000000 1.285048e-28 6.425239e-29
[349,] 1.0000000 4.899308e-28 2.449654e-28
[350,] 1.0000000 1.763136e-27 8.815681e-28
[351,] 1.0000000 6.628002e-27 3.314001e-27
[352,] 1.0000000 2.468904e-26 1.234452e-26
[353,] 1.0000000 9.111468e-26 4.555734e-26
[354,] 1.0000000 3.349750e-25 1.674875e-25
[355,] 1.0000000 1.213088e-24 6.065441e-25
[356,] 1.0000000 4.350442e-24 2.175221e-24
[357,] 1.0000000 1.524975e-23 7.624877e-24
[358,] 1.0000000 5.313013e-23 2.656507e-23
[359,] 1.0000000 1.853148e-22 9.265742e-23
[360,] 1.0000000 6.396865e-22 3.198433e-22
[361,] 1.0000000 2.184853e-21 1.092427e-21
[362,] 1.0000000 7.382103e-21 3.691052e-21
[363,] 1.0000000 2.527319e-20 1.263659e-20
[364,] 1.0000000 8.359335e-20 4.179668e-20
[365,] 1.0000000 2.723883e-19 1.361941e-19
[366,] 1.0000000 9.501666e-19 4.750833e-19
[367,] 1.0000000 3.249154e-18 1.624577e-18
[368,] 1.0000000 1.030860e-17 5.154301e-18
[369,] 1.0000000 3.230398e-17 1.615199e-17
[370,] 1.0000000 1.059753e-16 5.298765e-17
[371,] 1.0000000 3.065219e-16 1.532609e-16
[372,] 1.0000000 7.980859e-16 3.990429e-16
[373,] 1.0000000 2.377767e-15 1.188883e-15
[374,] 1.0000000 6.983956e-15 3.491978e-15
[375,] 1.0000000 2.202736e-14 1.101368e-14
[376,] 1.0000000 6.534443e-14 3.267222e-14
[377,] 1.0000000 1.856724e-13 9.283621e-14
[378,] 1.0000000 5.428907e-13 2.714453e-13
[379,] 1.0000000 1.429329e-12 7.146643e-13
[380,] 1.0000000 3.629379e-12 1.814690e-12
[381,] 1.0000000 9.523941e-12 4.761970e-12
[382,] 1.0000000 2.547438e-11 1.273719e-11
[383,] 1.0000000 6.756486e-11 3.378243e-11
[384,] 1.0000000 1.730454e-10 8.652272e-11
[385,] 1.0000000 4.751507e-10 2.375753e-10
[386,] 1.0000000 1.358775e-09 6.793876e-10
[387,] 1.0000000 2.576872e-09 1.288436e-09
[388,] 1.0000000 6.635957e-09 3.317979e-09
[389,] 1.0000000 1.761673e-08 8.808366e-09
[390,] 1.0000000 4.432430e-08 2.216215e-08
[391,] 1.0000000 8.306709e-08 4.153354e-08
[392,] 0.9999999 2.131984e-07 1.065992e-07
[393,] 0.9999997 5.623451e-07 2.811726e-07
[394,] 0.9999993 1.408515e-06 7.042574e-07
[395,] 0.9999986 2.844751e-06 1.422376e-06
[396,] 0.9999969 6.263561e-06 3.131780e-06
[397,] 0.9999928 1.444125e-05 7.220627e-06
[398,] 0.9999828 3.446735e-05 1.723367e-05
[399,] 0.9999633 7.343965e-05 3.671983e-05
[400,] 0.9999209 1.581122e-04 7.905612e-05
[401,] 0.9998541 2.918114e-04 1.459057e-04
[402,] 0.9997261 5.477223e-04 2.738612e-04
[403,] 0.9993926 1.214886e-03 6.074429e-04
[404,] 0.9986910 2.618086e-03 1.309043e-03
[405,] 0.9975113 4.977359e-03 2.488680e-03
[406,] 0.9948779 1.024411e-02 5.122056e-03
[407,] 0.9943533 1.129332e-02 5.646660e-03
[408,] 0.9911964 1.760730e-02 8.803649e-03
[409,] 0.9835631 3.287384e-02 1.643692e-02
[410,] 0.9948437 1.031264e-02 5.156320e-03
[411,] 0.9884542 2.309153e-02 1.154577e-02
[412,] 0.9747599 5.048022e-02 2.524011e-02
[413,] 0.9463041 1.073918e-01 5.369590e-02
[414,] 0.8926270 2.147460e-01 1.073730e-01
[415,] 0.7962125 4.075750e-01 2.037875e-01
[416,] 0.6628701 6.742599e-01 3.371299e-01
> postscript(file="/var/www/html/freestat/rcomp/tmp/1qki21291414557.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/2qki21291414557.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/3jbi51291414557.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/4jbi51291414557.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/5jbi51291414557.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 6
1259827.443 3195283.878 1187360.002 -2163938.266 -465658.107 -86633.872
7 8 9 10 11 12
-25648.977 64892.532 -278718.282 704705.338 -273562.970 408709.063
13 14 15 16 17 18
762830.868 -1184133.917 5446.994 498000.209 -234688.381 -12297.139
19 20 21 22 23 24
-306359.379 -322029.866 153902.340 -279201.028 -122718.005 -378580.851
25 26 27 28 29 30
-280953.361 -133413.771 -402357.714 -282147.420 333879.073 -89503.531
31 32 33 34 35 36
225880.987 -153760.016 190360.572 80331.658 -144589.537 -459875.779
37 38 39 40 41 42
48357.958 148737.922 -63103.461 -159079.702 -5335.757 -374596.864
43 44 45 46 47 48
127681.941 -158893.721 -8951.678 -124535.591 -298183.103 -131141.638
49 50 51 52 53 54
-28126.641 -157816.795 66254.339 -9824.282 -248826.358 71559.952
55 56 57 58 59 60
-348387.682 105736.955 297382.198 -6488.100 4473.442 -223588.738
61 62 63 64 65 66
-41076.166 -150276.044 -53093.255 -190418.162 43429.078 22800.421
67 68 69 70 71 72
-40997.771 -41888.465 -32278.740 -68881.815 -123725.963 20183.871
73 74 75 76 77 78
-74418.783 10492.987 -29593.037 -152906.495 -14508.840 -62936.025
79 80 81 82 83 84
37562.856 -43408.390 48877.645 -62195.598 -67415.436 58126.423
85 86 87 88 89 90
19942.976 -112771.362 60434.071 -6906.064 -83722.349 24236.170
91 92 93 94 95 96
20323.606 -58373.820 152629.298 -249808.092 97881.580 -18948.548
97 98 99 100 101 102
56121.615 27844.604 -65081.319 -65872.062 -69214.165 -109444.008
103 104 105 106 107 108
-127969.710 -65383.903 51398.186 -31544.333 -266560.840 63483.018
109 110 111 112 113 114
-169212.146 -139539.338 111838.134 -223309.564 -148255.495 -124797.099
115 116 117 118 119 120
23515.058 55479.808 17147.822 63553.861 -200546.755 50336.225
121 122 123 124 125 126
-265684.483 34349.869 22362.228 29543.810 -37960.116 -223683.527
127 128 129 130 131 132
47930.758 46990.954 54791.004 47930.758 49320.735 90994.114
133 134 135 136 137 138
47930.758 45014.662 27433.327 42217.914 52885.045 48103.003
139 140 141 142 143 144
22758.922 47930.758 55492.576 47930.758 15002.087 48197.857
145 146 147 148 149 150
49670.808 48756.725 47930.758 49490.851 50535.241 47847.411
151 152 153 154 155 156
16757.609 16128.433 47930.758 25176.204 48001.084 47930.758
157 158 159 160 161 162
47930.758 44196.429 59511.531 47930.758 47930.758 47930.758
163 164 165 166 167 168
47930.758 35201.158 37073.804 47253.085 47930.758 129383.608
169 170 171 172 173 174
95359.943 47930.758 47930.758 43892.484 47930.758 37413.394
175 176 177 178 179 180
47930.758 47968.115 95780.130 -8516.655 47930.758 47930.758
181 182 183 184 185 186
55714.984 -25660.802 6174.204 -37101.259 -10313.459 -2370.308
187 188 189 190 191 192
-35584.633 40932.200 9957.267 47930.758 47076.205 47930.758
193 194 195 196 197 198
44403.489 47930.758 47890.901 -7555.060 -126404.718 -7555.060
199 200 201 202 203 204
-7555.060 -96498.844 -7555.060 -7555.060 -7555.060 2068.080
205 206 207 208 209 210
-56691.365 -8066.752 -14179.696 -7555.060 23336.527 -11265.356
211 212 213 214 215 216
47930.758 47930.758 73611.546 13127.108 20306.970 47930.758
217 218 219 220 221 222
50834.989 49994.478 47930.758 19719.600 22654.245 47930.758
223 224 225 226 227 228
48672.001 39310.301 48181.113 48453.678 50983.381 47930.758
229 230 231 232 233 234
32019.062 47930.758 19310.917 35997.643 47930.758 29871.149
235 236 237 238 239 240
47930.758 48340.985 47930.758 49921.272 47930.758 47930.758
241 242 243 244 245 246
47930.758 42624.223 47930.758 52466.884 -33992.012 -7555.060
247 248 249 250 251 252
-7555.060 40017.496 -7555.060 -7555.060 41434.050 -7555.060
253 254 255 256 257 258
28918.905 24439.537 -87387.956 47930.758 45851.206 59484.659
259 260 261 262 263 264
-7555.060 47930.758 -7555.060 -7555.060 51292.884 65051.344
265 266 267 268 269 270
-7555.060 47930.758 -7555.060 79985.494 -7555.060 44777.497
271 272 273 274 275 276
-7555.060 -7555.060 -7555.060 43760.928 47930.758 53038.424
277 278 279 280 281 282
49743.121 47930.758 -35871.190 35089.843 -7555.060 47930.758
283 284 285 286 287 288
75626.658 47930.758 39802.780 47930.758 -19198.651 48830.870
289 290 291 292 293 294
-7555.060 55446.705 -7555.060 58494.152 54252.289 47930.758
295 296 297 298 299 300
41418.222 17805.179 47930.758 -7555.060 47930.758 47930.758
301 302 303 304 305 306
-9187.034 29601.307 47930.758 -7555.060 85548.186 51520.595
307 308 309 310 311 312
-7555.060 -13574.240 51524.780 -7555.060 -7555.060 51085.001
313 314 315 316 317 318
-7555.060 -5631.489 49666.953 63580.147 48377.897 24781.823
319 320 321 322 323 324
50337.774 47930.758 47930.758 47930.758 47930.758 -7555.060
325 326 327 328 329 330
47930.758 45207.361 -7555.060 -7555.060 47930.758 47930.758
331 332 333 334 335 336
47930.758 47930.758 47930.758 46110.459 48709.764 40596.813
337 338 339 340 341 342
34892.012 47930.758 45723.816 53714.507 20843.505 47930.758
343 344 345 346 347 348
47930.758 47930.758 41837.252 47930.758 38208.833 47930.758
349 350 351 352 353 354
47930.758 11292.552 51727.086 48104.426 47803.348 40959.878
355 356 357 358 359 360
47930.758 47930.758 25991.297 47930.758 47930.758 47930.758
361 362 363 364 365 366
49543.695 47930.758 47930.758 -59079.502 50139.283 47930.758
367 368 369 370 371 372
47930.758 47930.758 47930.758 50155.084 47930.758 60651.036
373 374 375 376 377 378
35414.901 27297.538 47930.758 47930.758 42833.808 41821.079
379 380 381 382 383 384
-7827.369 47930.758 47930.758 33983.667 65563.601 47930.758
385 386 387 388 389 390
50256.107 -100883.026 60199.198 21375.207 46456.029 34065.162
391 392 393 394 395 396
43281.192 17215.934 -20308.732 68636.060 -180351.327 9037.666
397 398 399 400 401 402
-138046.328 4172.197 10638.634 -23000.886 -11465.171 39825.868
403 404 405 406 407 408
17459.373 -95882.735 -12187.230 10402.487 23566.412 15549.014
409 410 411 412 413 414
26690.904 -73845.864 10354.181 -250216.313 -74411.743 12120.541
415 416 417 418 419 420
-189194.453 -115267.526 118575.766 -62387.577 -35757.824 -49407.953
421 422 423 424 425 426
-60797.317 -107098.631 -576390.302 -484134.906 -41774.594 -232554.023
427 428 429 430 431
-313458.634 -1118357.214 -306374.089 -344582.470 -675705.948
> postscript(file="/var/www/html/freestat/rcomp/tmp/6c2hq1291414557.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 1259827.443 NA
1 3195283.878 1259827.443
2 1187360.002 3195283.878
3 -2163938.266 1187360.002
4 -465658.107 -2163938.266
5 -86633.872 -465658.107
6 -25648.977 -86633.872
7 64892.532 -25648.977
8 -278718.282 64892.532
9 704705.338 -278718.282
10 -273562.970 704705.338
11 408709.063 -273562.970
12 762830.868 408709.063
13 -1184133.917 762830.868
14 5446.994 -1184133.917
15 498000.209 5446.994
16 -234688.381 498000.209
17 -12297.139 -234688.381
18 -306359.379 -12297.139
19 -322029.866 -306359.379
20 153902.340 -322029.866
21 -279201.028 153902.340
22 -122718.005 -279201.028
23 -378580.851 -122718.005
24 -280953.361 -378580.851
25 -133413.771 -280953.361
26 -402357.714 -133413.771
27 -282147.420 -402357.714
28 333879.073 -282147.420
29 -89503.531 333879.073
30 225880.987 -89503.531
31 -153760.016 225880.987
32 190360.572 -153760.016
33 80331.658 190360.572
34 -144589.537 80331.658
35 -459875.779 -144589.537
36 48357.958 -459875.779
37 148737.922 48357.958
38 -63103.461 148737.922
39 -159079.702 -63103.461
40 -5335.757 -159079.702
41 -374596.864 -5335.757
42 127681.941 -374596.864
43 -158893.721 127681.941
44 -8951.678 -158893.721
45 -124535.591 -8951.678
46 -298183.103 -124535.591
47 -131141.638 -298183.103
48 -28126.641 -131141.638
49 -157816.795 -28126.641
50 66254.339 -157816.795
51 -9824.282 66254.339
52 -248826.358 -9824.282
53 71559.952 -248826.358
54 -348387.682 71559.952
55 105736.955 -348387.682
56 297382.198 105736.955
57 -6488.100 297382.198
58 4473.442 -6488.100
59 -223588.738 4473.442
60 -41076.166 -223588.738
61 -150276.044 -41076.166
62 -53093.255 -150276.044
63 -190418.162 -53093.255
64 43429.078 -190418.162
65 22800.421 43429.078
66 -40997.771 22800.421
67 -41888.465 -40997.771
68 -32278.740 -41888.465
69 -68881.815 -32278.740
70 -123725.963 -68881.815
71 20183.871 -123725.963
72 -74418.783 20183.871
73 10492.987 -74418.783
74 -29593.037 10492.987
75 -152906.495 -29593.037
76 -14508.840 -152906.495
77 -62936.025 -14508.840
78 37562.856 -62936.025
79 -43408.390 37562.856
80 48877.645 -43408.390
81 -62195.598 48877.645
82 -67415.436 -62195.598
83 58126.423 -67415.436
84 19942.976 58126.423
85 -112771.362 19942.976
86 60434.071 -112771.362
87 -6906.064 60434.071
88 -83722.349 -6906.064
89 24236.170 -83722.349
90 20323.606 24236.170
91 -58373.820 20323.606
92 152629.298 -58373.820
93 -249808.092 152629.298
94 97881.580 -249808.092
95 -18948.548 97881.580
96 56121.615 -18948.548
97 27844.604 56121.615
98 -65081.319 27844.604
99 -65872.062 -65081.319
100 -69214.165 -65872.062
101 -109444.008 -69214.165
102 -127969.710 -109444.008
103 -65383.903 -127969.710
104 51398.186 -65383.903
105 -31544.333 51398.186
106 -266560.840 -31544.333
107 63483.018 -266560.840
108 -169212.146 63483.018
109 -139539.338 -169212.146
110 111838.134 -139539.338
111 -223309.564 111838.134
112 -148255.495 -223309.564
113 -124797.099 -148255.495
114 23515.058 -124797.099
115 55479.808 23515.058
116 17147.822 55479.808
117 63553.861 17147.822
118 -200546.755 63553.861
119 50336.225 -200546.755
120 -265684.483 50336.225
121 34349.869 -265684.483
122 22362.228 34349.869
123 29543.810 22362.228
124 -37960.116 29543.810
125 -223683.527 -37960.116
126 47930.758 -223683.527
127 46990.954 47930.758
128 54791.004 46990.954
129 47930.758 54791.004
130 49320.735 47930.758
131 90994.114 49320.735
132 47930.758 90994.114
133 45014.662 47930.758
134 27433.327 45014.662
135 42217.914 27433.327
136 52885.045 42217.914
137 48103.003 52885.045
138 22758.922 48103.003
139 47930.758 22758.922
140 55492.576 47930.758
141 47930.758 55492.576
142 15002.087 47930.758
143 48197.857 15002.087
144 49670.808 48197.857
145 48756.725 49670.808
146 47930.758 48756.725
147 49490.851 47930.758
148 50535.241 49490.851
149 47847.411 50535.241
150 16757.609 47847.411
151 16128.433 16757.609
152 47930.758 16128.433
153 25176.204 47930.758
154 48001.084 25176.204
155 47930.758 48001.084
156 47930.758 47930.758
157 44196.429 47930.758
158 59511.531 44196.429
159 47930.758 59511.531
160 47930.758 47930.758
161 47930.758 47930.758
162 47930.758 47930.758
163 35201.158 47930.758
164 37073.804 35201.158
165 47253.085 37073.804
166 47930.758 47253.085
167 129383.608 47930.758
168 95359.943 129383.608
169 47930.758 95359.943
170 47930.758 47930.758
171 43892.484 47930.758
172 47930.758 43892.484
173 37413.394 47930.758
174 47930.758 37413.394
175 47968.115 47930.758
176 95780.130 47968.115
177 -8516.655 95780.130
178 47930.758 -8516.655
179 47930.758 47930.758
180 55714.984 47930.758
181 -25660.802 55714.984
182 6174.204 -25660.802
183 -37101.259 6174.204
184 -10313.459 -37101.259
185 -2370.308 -10313.459
186 -35584.633 -2370.308
187 40932.200 -35584.633
188 9957.267 40932.200
189 47930.758 9957.267
190 47076.205 47930.758
191 47930.758 47076.205
192 44403.489 47930.758
193 47930.758 44403.489
194 47890.901 47930.758
195 -7555.060 47890.901
196 -126404.718 -7555.060
197 -7555.060 -126404.718
198 -7555.060 -7555.060
199 -96498.844 -7555.060
200 -7555.060 -96498.844
201 -7555.060 -7555.060
202 -7555.060 -7555.060
203 2068.080 -7555.060
204 -56691.365 2068.080
205 -8066.752 -56691.365
206 -14179.696 -8066.752
207 -7555.060 -14179.696
208 23336.527 -7555.060
209 -11265.356 23336.527
210 47930.758 -11265.356
211 47930.758 47930.758
212 73611.546 47930.758
213 13127.108 73611.546
214 20306.970 13127.108
215 47930.758 20306.970
216 50834.989 47930.758
217 49994.478 50834.989
218 47930.758 49994.478
219 19719.600 47930.758
220 22654.245 19719.600
221 47930.758 22654.245
222 48672.001 47930.758
223 39310.301 48672.001
224 48181.113 39310.301
225 48453.678 48181.113
226 50983.381 48453.678
227 47930.758 50983.381
228 32019.062 47930.758
229 47930.758 32019.062
230 19310.917 47930.758
231 35997.643 19310.917
232 47930.758 35997.643
233 29871.149 47930.758
234 47930.758 29871.149
235 48340.985 47930.758
236 47930.758 48340.985
237 49921.272 47930.758
238 47930.758 49921.272
239 47930.758 47930.758
240 47930.758 47930.758
241 42624.223 47930.758
242 47930.758 42624.223
243 52466.884 47930.758
244 -33992.012 52466.884
245 -7555.060 -33992.012
246 -7555.060 -7555.060
247 40017.496 -7555.060
248 -7555.060 40017.496
249 -7555.060 -7555.060
250 41434.050 -7555.060
251 -7555.060 41434.050
252 28918.905 -7555.060
253 24439.537 28918.905
254 -87387.956 24439.537
255 47930.758 -87387.956
256 45851.206 47930.758
257 59484.659 45851.206
258 -7555.060 59484.659
259 47930.758 -7555.060
260 -7555.060 47930.758
261 -7555.060 -7555.060
262 51292.884 -7555.060
263 65051.344 51292.884
264 -7555.060 65051.344
265 47930.758 -7555.060
266 -7555.060 47930.758
267 79985.494 -7555.060
268 -7555.060 79985.494
269 44777.497 -7555.060
270 -7555.060 44777.497
271 -7555.060 -7555.060
272 -7555.060 -7555.060
273 43760.928 -7555.060
274 47930.758 43760.928
275 53038.424 47930.758
276 49743.121 53038.424
277 47930.758 49743.121
278 -35871.190 47930.758
279 35089.843 -35871.190
280 -7555.060 35089.843
281 47930.758 -7555.060
282 75626.658 47930.758
283 47930.758 75626.658
284 39802.780 47930.758
285 47930.758 39802.780
286 -19198.651 47930.758
287 48830.870 -19198.651
288 -7555.060 48830.870
289 55446.705 -7555.060
290 -7555.060 55446.705
291 58494.152 -7555.060
292 54252.289 58494.152
293 47930.758 54252.289
294 41418.222 47930.758
295 17805.179 41418.222
296 47930.758 17805.179
297 -7555.060 47930.758
298 47930.758 -7555.060
299 47930.758 47930.758
300 -9187.034 47930.758
301 29601.307 -9187.034
302 47930.758 29601.307
303 -7555.060 47930.758
304 85548.186 -7555.060
305 51520.595 85548.186
306 -7555.060 51520.595
307 -13574.240 -7555.060
308 51524.780 -13574.240
309 -7555.060 51524.780
310 -7555.060 -7555.060
311 51085.001 -7555.060
312 -7555.060 51085.001
313 -5631.489 -7555.060
314 49666.953 -5631.489
315 63580.147 49666.953
316 48377.897 63580.147
317 24781.823 48377.897
318 50337.774 24781.823
319 47930.758 50337.774
320 47930.758 47930.758
321 47930.758 47930.758
322 47930.758 47930.758
323 -7555.060 47930.758
324 47930.758 -7555.060
325 45207.361 47930.758
326 -7555.060 45207.361
327 -7555.060 -7555.060
328 47930.758 -7555.060
329 47930.758 47930.758
330 47930.758 47930.758
331 47930.758 47930.758
332 47930.758 47930.758
333 46110.459 47930.758
334 48709.764 46110.459
335 40596.813 48709.764
336 34892.012 40596.813
337 47930.758 34892.012
338 45723.816 47930.758
339 53714.507 45723.816
340 20843.505 53714.507
341 47930.758 20843.505
342 47930.758 47930.758
343 47930.758 47930.758
344 41837.252 47930.758
345 47930.758 41837.252
346 38208.833 47930.758
347 47930.758 38208.833
348 47930.758 47930.758
349 11292.552 47930.758
350 51727.086 11292.552
351 48104.426 51727.086
352 47803.348 48104.426
353 40959.878 47803.348
354 47930.758 40959.878
355 47930.758 47930.758
356 25991.297 47930.758
357 47930.758 25991.297
358 47930.758 47930.758
359 47930.758 47930.758
360 49543.695 47930.758
361 47930.758 49543.695
362 47930.758 47930.758
363 -59079.502 47930.758
364 50139.283 -59079.502
365 47930.758 50139.283
366 47930.758 47930.758
367 47930.758 47930.758
368 47930.758 47930.758
369 50155.084 47930.758
370 47930.758 50155.084
371 60651.036 47930.758
372 35414.901 60651.036
373 27297.538 35414.901
374 47930.758 27297.538
375 47930.758 47930.758
376 42833.808 47930.758
377 41821.079 42833.808
378 -7827.369 41821.079
379 47930.758 -7827.369
380 47930.758 47930.758
381 33983.667 47930.758
382 65563.601 33983.667
383 47930.758 65563.601
384 50256.107 47930.758
385 -100883.026 50256.107
386 60199.198 -100883.026
387 21375.207 60199.198
388 46456.029 21375.207
389 34065.162 46456.029
390 43281.192 34065.162
391 17215.934 43281.192
392 -20308.732 17215.934
393 68636.060 -20308.732
394 -180351.327 68636.060
395 9037.666 -180351.327
396 -138046.328 9037.666
397 4172.197 -138046.328
398 10638.634 4172.197
399 -23000.886 10638.634
400 -11465.171 -23000.886
401 39825.868 -11465.171
402 17459.373 39825.868
403 -95882.735 17459.373
404 -12187.230 -95882.735
405 10402.487 -12187.230
406 23566.412 10402.487
407 15549.014 23566.412
408 26690.904 15549.014
409 -73845.864 26690.904
410 10354.181 -73845.864
411 -250216.313 10354.181
412 -74411.743 -250216.313
413 12120.541 -74411.743
414 -189194.453 12120.541
415 -115267.526 -189194.453
416 118575.766 -115267.526
417 -62387.577 118575.766
418 -35757.824 -62387.577
419 -49407.953 -35757.824
420 -60797.317 -49407.953
421 -107098.631 -60797.317
422 -576390.302 -107098.631
423 -484134.906 -576390.302
424 -41774.594 -484134.906
425 -232554.023 -41774.594
426 -313458.634 -232554.023
427 -1118357.214 -313458.634
428 -306374.089 -1118357.214
429 -344582.470 -306374.089
430 -675705.948 -344582.470
431 NA -675705.948
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] 3195283.878 1259827.443
[2,] 1187360.002 3195283.878
[3,] -2163938.266 1187360.002
[4,] -465658.107 -2163938.266
[5,] -86633.872 -465658.107
[6,] -25648.977 -86633.872
[7,] 64892.532 -25648.977
[8,] -278718.282 64892.532
[9,] 704705.338 -278718.282
[10,] -273562.970 704705.338
[11,] 408709.063 -273562.970
[12,] 762830.868 408709.063
[13,] -1184133.917 762830.868
[14,] 5446.994 -1184133.917
[15,] 498000.209 5446.994
[16,] -234688.381 498000.209
[17,] -12297.139 -234688.381
[18,] -306359.379 -12297.139
[19,] -322029.866 -306359.379
[20,] 153902.340 -322029.866
[21,] -279201.028 153902.340
[22,] -122718.005 -279201.028
[23,] -378580.851 -122718.005
[24,] -280953.361 -378580.851
[25,] -133413.771 -280953.361
[26,] -402357.714 -133413.771
[27,] -282147.420 -402357.714
[28,] 333879.073 -282147.420
[29,] -89503.531 333879.073
[30,] 225880.987 -89503.531
[31,] -153760.016 225880.987
[32,] 190360.572 -153760.016
[33,] 80331.658 190360.572
[34,] -144589.537 80331.658
[35,] -459875.779 -144589.537
[36,] 48357.958 -459875.779
[37,] 148737.922 48357.958
[38,] -63103.461 148737.922
[39,] -159079.702 -63103.461
[40,] -5335.757 -159079.702
[41,] -374596.864 -5335.757
[42,] 127681.941 -374596.864
[43,] -158893.721 127681.941
[44,] -8951.678 -158893.721
[45,] -124535.591 -8951.678
[46,] -298183.103 -124535.591
[47,] -131141.638 -298183.103
[48,] -28126.641 -131141.638
[49,] -157816.795 -28126.641
[50,] 66254.339 -157816.795
[51,] -9824.282 66254.339
[52,] -248826.358 -9824.282
[53,] 71559.952 -248826.358
[54,] -348387.682 71559.952
[55,] 105736.955 -348387.682
[56,] 297382.198 105736.955
[57,] -6488.100 297382.198
[58,] 4473.442 -6488.100
[59,] -223588.738 4473.442
[60,] -41076.166 -223588.738
[61,] -150276.044 -41076.166
[62,] -53093.255 -150276.044
[63,] -190418.162 -53093.255
[64,] 43429.078 -190418.162
[65,] 22800.421 43429.078
[66,] -40997.771 22800.421
[67,] -41888.465 -40997.771
[68,] -32278.740 -41888.465
[69,] -68881.815 -32278.740
[70,] -123725.963 -68881.815
[71,] 20183.871 -123725.963
[72,] -74418.783 20183.871
[73,] 10492.987 -74418.783
[74,] -29593.037 10492.987
[75,] -152906.495 -29593.037
[76,] -14508.840 -152906.495
[77,] -62936.025 -14508.840
[78,] 37562.856 -62936.025
[79,] -43408.390 37562.856
[80,] 48877.645 -43408.390
[81,] -62195.598 48877.645
[82,] -67415.436 -62195.598
[83,] 58126.423 -67415.436
[84,] 19942.976 58126.423
[85,] -112771.362 19942.976
[86,] 60434.071 -112771.362
[87,] -6906.064 60434.071
[88,] -83722.349 -6906.064
[89,] 24236.170 -83722.349
[90,] 20323.606 24236.170
[91,] -58373.820 20323.606
[92,] 152629.298 -58373.820
[93,] -249808.092 152629.298
[94,] 97881.580 -249808.092
[95,] -18948.548 97881.580
[96,] 56121.615 -18948.548
[97,] 27844.604 56121.615
[98,] -65081.319 27844.604
[99,] -65872.062 -65081.319
[100,] -69214.165 -65872.062
[101,] -109444.008 -69214.165
[102,] -127969.710 -109444.008
[103,] -65383.903 -127969.710
[104,] 51398.186 -65383.903
[105,] -31544.333 51398.186
[106,] -266560.840 -31544.333
[107,] 63483.018 -266560.840
[108,] -169212.146 63483.018
[109,] -139539.338 -169212.146
[110,] 111838.134 -139539.338
[111,] -223309.564 111838.134
[112,] -148255.495 -223309.564
[113,] -124797.099 -148255.495
[114,] 23515.058 -124797.099
[115,] 55479.808 23515.058
[116,] 17147.822 55479.808
[117,] 63553.861 17147.822
[118,] -200546.755 63553.861
[119,] 50336.225 -200546.755
[120,] -265684.483 50336.225
[121,] 34349.869 -265684.483
[122,] 22362.228 34349.869
[123,] 29543.810 22362.228
[124,] -37960.116 29543.810
[125,] -223683.527 -37960.116
[126,] 47930.758 -223683.527
[127,] 46990.954 47930.758
[128,] 54791.004 46990.954
[129,] 47930.758 54791.004
[130,] 49320.735 47930.758
[131,] 90994.114 49320.735
[132,] 47930.758 90994.114
[133,] 45014.662 47930.758
[134,] 27433.327 45014.662
[135,] 42217.914 27433.327
[136,] 52885.045 42217.914
[137,] 48103.003 52885.045
[138,] 22758.922 48103.003
[139,] 47930.758 22758.922
[140,] 55492.576 47930.758
[141,] 47930.758 55492.576
[142,] 15002.087 47930.758
[143,] 48197.857 15002.087
[144,] 49670.808 48197.857
[145,] 48756.725 49670.808
[146,] 47930.758 48756.725
[147,] 49490.851 47930.758
[148,] 50535.241 49490.851
[149,] 47847.411 50535.241
[150,] 16757.609 47847.411
[151,] 16128.433 16757.609
[152,] 47930.758 16128.433
[153,] 25176.204 47930.758
[154,] 48001.084 25176.204
[155,] 47930.758 48001.084
[156,] 47930.758 47930.758
[157,] 44196.429 47930.758
[158,] 59511.531 44196.429
[159,] 47930.758 59511.531
[160,] 47930.758 47930.758
[161,] 47930.758 47930.758
[162,] 47930.758 47930.758
[163,] 35201.158 47930.758
[164,] 37073.804 35201.158
[165,] 47253.085 37073.804
[166,] 47930.758 47253.085
[167,] 129383.608 47930.758
[168,] 95359.943 129383.608
[169,] 47930.758 95359.943
[170,] 47930.758 47930.758
[171,] 43892.484 47930.758
[172,] 47930.758 43892.484
[173,] 37413.394 47930.758
[174,] 47930.758 37413.394
[175,] 47968.115 47930.758
[176,] 95780.130 47968.115
[177,] -8516.655 95780.130
[178,] 47930.758 -8516.655
[179,] 47930.758 47930.758
[180,] 55714.984 47930.758
[181,] -25660.802 55714.984
[182,] 6174.204 -25660.802
[183,] -37101.259 6174.204
[184,] -10313.459 -37101.259
[185,] -2370.308 -10313.459
[186,] -35584.633 -2370.308
[187,] 40932.200 -35584.633
[188,] 9957.267 40932.200
[189,] 47930.758 9957.267
[190,] 47076.205 47930.758
[191,] 47930.758 47076.205
[192,] 44403.489 47930.758
[193,] 47930.758 44403.489
[194,] 47890.901 47930.758
[195,] -7555.060 47890.901
[196,] -126404.718 -7555.060
[197,] -7555.060 -126404.718
[198,] -7555.060 -7555.060
[199,] -96498.844 -7555.060
[200,] -7555.060 -96498.844
[201,] -7555.060 -7555.060
[202,] -7555.060 -7555.060
[203,] 2068.080 -7555.060
[204,] -56691.365 2068.080
[205,] -8066.752 -56691.365
[206,] -14179.696 -8066.752
[207,] -7555.060 -14179.696
[208,] 23336.527 -7555.060
[209,] -11265.356 23336.527
[210,] 47930.758 -11265.356
[211,] 47930.758 47930.758
[212,] 73611.546 47930.758
[213,] 13127.108 73611.546
[214,] 20306.970 13127.108
[215,] 47930.758 20306.970
[216,] 50834.989 47930.758
[217,] 49994.478 50834.989
[218,] 47930.758 49994.478
[219,] 19719.600 47930.758
[220,] 22654.245 19719.600
[221,] 47930.758 22654.245
[222,] 48672.001 47930.758
[223,] 39310.301 48672.001
[224,] 48181.113 39310.301
[225,] 48453.678 48181.113
[226,] 50983.381 48453.678
[227,] 47930.758 50983.381
[228,] 32019.062 47930.758
[229,] 47930.758 32019.062
[230,] 19310.917 47930.758
[231,] 35997.643 19310.917
[232,] 47930.758 35997.643
[233,] 29871.149 47930.758
[234,] 47930.758 29871.149
[235,] 48340.985 47930.758
[236,] 47930.758 48340.985
[237,] 49921.272 47930.758
[238,] 47930.758 49921.272
[239,] 47930.758 47930.758
[240,] 47930.758 47930.758
[241,] 42624.223 47930.758
[242,] 47930.758 42624.223
[243,] 52466.884 47930.758
[244,] -33992.012 52466.884
[245,] -7555.060 -33992.012
[246,] -7555.060 -7555.060
[247,] 40017.496 -7555.060
[248,] -7555.060 40017.496
[249,] -7555.060 -7555.060
[250,] 41434.050 -7555.060
[251,] -7555.060 41434.050
[252,] 28918.905 -7555.060
[253,] 24439.537 28918.905
[254,] -87387.956 24439.537
[255,] 47930.758 -87387.956
[256,] 45851.206 47930.758
[257,] 59484.659 45851.206
[258,] -7555.060 59484.659
[259,] 47930.758 -7555.060
[260,] -7555.060 47930.758
[261,] -7555.060 -7555.060
[262,] 51292.884 -7555.060
[263,] 65051.344 51292.884
[264,] -7555.060 65051.344
[265,] 47930.758 -7555.060
[266,] -7555.060 47930.758
[267,] 79985.494 -7555.060
[268,] -7555.060 79985.494
[269,] 44777.497 -7555.060
[270,] -7555.060 44777.497
[271,] -7555.060 -7555.060
[272,] -7555.060 -7555.060
[273,] 43760.928 -7555.060
[274,] 47930.758 43760.928
[275,] 53038.424 47930.758
[276,] 49743.121 53038.424
[277,] 47930.758 49743.121
[278,] -35871.190 47930.758
[279,] 35089.843 -35871.190
[280,] -7555.060 35089.843
[281,] 47930.758 -7555.060
[282,] 75626.658 47930.758
[283,] 47930.758 75626.658
[284,] 39802.780 47930.758
[285,] 47930.758 39802.780
[286,] -19198.651 47930.758
[287,] 48830.870 -19198.651
[288,] -7555.060 48830.870
[289,] 55446.705 -7555.060
[290,] -7555.060 55446.705
[291,] 58494.152 -7555.060
[292,] 54252.289 58494.152
[293,] 47930.758 54252.289
[294,] 41418.222 47930.758
[295,] 17805.179 41418.222
[296,] 47930.758 17805.179
[297,] -7555.060 47930.758
[298,] 47930.758 -7555.060
[299,] 47930.758 47930.758
[300,] -9187.034 47930.758
[301,] 29601.307 -9187.034
[302,] 47930.758 29601.307
[303,] -7555.060 47930.758
[304,] 85548.186 -7555.060
[305,] 51520.595 85548.186
[306,] -7555.060 51520.595
[307,] -13574.240 -7555.060
[308,] 51524.780 -13574.240
[309,] -7555.060 51524.780
[310,] -7555.060 -7555.060
[311,] 51085.001 -7555.060
[312,] -7555.060 51085.001
[313,] -5631.489 -7555.060
[314,] 49666.953 -5631.489
[315,] 63580.147 49666.953
[316,] 48377.897 63580.147
[317,] 24781.823 48377.897
[318,] 50337.774 24781.823
[319,] 47930.758 50337.774
[320,] 47930.758 47930.758
[321,] 47930.758 47930.758
[322,] 47930.758 47930.758
[323,] -7555.060 47930.758
[324,] 47930.758 -7555.060
[325,] 45207.361 47930.758
[326,] -7555.060 45207.361
[327,] -7555.060 -7555.060
[328,] 47930.758 -7555.060
[329,] 47930.758 47930.758
[330,] 47930.758 47930.758
[331,] 47930.758 47930.758
[332,] 47930.758 47930.758
[333,] 46110.459 47930.758
[334,] 48709.764 46110.459
[335,] 40596.813 48709.764
[336,] 34892.012 40596.813
[337,] 47930.758 34892.012
[338,] 45723.816 47930.758
[339,] 53714.507 45723.816
[340,] 20843.505 53714.507
[341,] 47930.758 20843.505
[342,] 47930.758 47930.758
[343,] 47930.758 47930.758
[344,] 41837.252 47930.758
[345,] 47930.758 41837.252
[346,] 38208.833 47930.758
[347,] 47930.758 38208.833
[348,] 47930.758 47930.758
[349,] 11292.552 47930.758
[350,] 51727.086 11292.552
[351,] 48104.426 51727.086
[352,] 47803.348 48104.426
[353,] 40959.878 47803.348
[354,] 47930.758 40959.878
[355,] 47930.758 47930.758
[356,] 25991.297 47930.758
[357,] 47930.758 25991.297
[358,] 47930.758 47930.758
[359,] 47930.758 47930.758
[360,] 49543.695 47930.758
[361,] 47930.758 49543.695
[362,] 47930.758 47930.758
[363,] -59079.502 47930.758
[364,] 50139.283 -59079.502
[365,] 47930.758 50139.283
[366,] 47930.758 47930.758
[367,] 47930.758 47930.758
[368,] 47930.758 47930.758
[369,] 50155.084 47930.758
[370,] 47930.758 50155.084
[371,] 60651.036 47930.758
[372,] 35414.901 60651.036
[373,] 27297.538 35414.901
[374,] 47930.758 27297.538
[375,] 47930.758 47930.758
[376,] 42833.808 47930.758
[377,] 41821.079 42833.808
[378,] -7827.369 41821.079
[379,] 47930.758 -7827.369
[380,] 47930.758 47930.758
[381,] 33983.667 47930.758
[382,] 65563.601 33983.667
[383,] 47930.758 65563.601
[384,] 50256.107 47930.758
[385,] -100883.026 50256.107
[386,] 60199.198 -100883.026
[387,] 21375.207 60199.198
[388,] 46456.029 21375.207
[389,] 34065.162 46456.029
[390,] 43281.192 34065.162
[391,] 17215.934 43281.192
[392,] -20308.732 17215.934
[393,] 68636.060 -20308.732
[394,] -180351.327 68636.060
[395,] 9037.666 -180351.327
[396,] -138046.328 9037.666
[397,] 4172.197 -138046.328
[398,] 10638.634 4172.197
[399,] -23000.886 10638.634
[400,] -11465.171 -23000.886
[401,] 39825.868 -11465.171
[402,] 17459.373 39825.868
[403,] -95882.735 17459.373
[404,] -12187.230 -95882.735
[405,] 10402.487 -12187.230
[406,] 23566.412 10402.487
[407,] 15549.014 23566.412
[408,] 26690.904 15549.014
[409,] -73845.864 26690.904
[410,] 10354.181 -73845.864
[411,] -250216.313 10354.181
[412,] -74411.743 -250216.313
[413,] 12120.541 -74411.743
[414,] -189194.453 12120.541
[415,] -115267.526 -189194.453
[416,] 118575.766 -115267.526
[417,] -62387.577 118575.766
[418,] -35757.824 -62387.577
[419,] -49407.953 -35757.824
[420,] -60797.317 -49407.953
[421,] -107098.631 -60797.317
[422,] -576390.302 -107098.631
[423,] -484134.906 -576390.302
[424,] -41774.594 -484134.906
[425,] -232554.023 -41774.594
[426,] -313458.634 -232554.023
[427,] -1118357.214 -313458.634
[428,] -306374.089 -1118357.214
[429,] -344582.470 -306374.089
[430,] -675705.948 -344582.470
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 3195283.878 1259827.443
2 1187360.002 3195283.878
3 -2163938.266 1187360.002
4 -465658.107 -2163938.266
5 -86633.872 -465658.107
6 -25648.977 -86633.872
7 64892.532 -25648.977
8 -278718.282 64892.532
9 704705.338 -278718.282
10 -273562.970 704705.338
11 408709.063 -273562.970
12 762830.868 408709.063
13 -1184133.917 762830.868
14 5446.994 -1184133.917
15 498000.209 5446.994
16 -234688.381 498000.209
17 -12297.139 -234688.381
18 -306359.379 -12297.139
19 -322029.866 -306359.379
20 153902.340 -322029.866
21 -279201.028 153902.340
22 -122718.005 -279201.028
23 -378580.851 -122718.005
24 -280953.361 -378580.851
25 -133413.771 -280953.361
26 -402357.714 -133413.771
27 -282147.420 -402357.714
28 333879.073 -282147.420
29 -89503.531 333879.073
30 225880.987 -89503.531
31 -153760.016 225880.987
32 190360.572 -153760.016
33 80331.658 190360.572
34 -144589.537 80331.658
35 -459875.779 -144589.537
36 48357.958 -459875.779
37 148737.922 48357.958
38 -63103.461 148737.922
39 -159079.702 -63103.461
40 -5335.757 -159079.702
41 -374596.864 -5335.757
42 127681.941 -374596.864
43 -158893.721 127681.941
44 -8951.678 -158893.721
45 -124535.591 -8951.678
46 -298183.103 -124535.591
47 -131141.638 -298183.103
48 -28126.641 -131141.638
49 -157816.795 -28126.641
50 66254.339 -157816.795
51 -9824.282 66254.339
52 -248826.358 -9824.282
53 71559.952 -248826.358
54 -348387.682 71559.952
55 105736.955 -348387.682
56 297382.198 105736.955
57 -6488.100 297382.198
58 4473.442 -6488.100
59 -223588.738 4473.442
60 -41076.166 -223588.738
61 -150276.044 -41076.166
62 -53093.255 -150276.044
63 -190418.162 -53093.255
64 43429.078 -190418.162
65 22800.421 43429.078
66 -40997.771 22800.421
67 -41888.465 -40997.771
68 -32278.740 -41888.465
69 -68881.815 -32278.740
70 -123725.963 -68881.815
71 20183.871 -123725.963
72 -74418.783 20183.871
73 10492.987 -74418.783
74 -29593.037 10492.987
75 -152906.495 -29593.037
76 -14508.840 -152906.495
77 -62936.025 -14508.840
78 37562.856 -62936.025
79 -43408.390 37562.856
80 48877.645 -43408.390
81 -62195.598 48877.645
82 -67415.436 -62195.598
83 58126.423 -67415.436
84 19942.976 58126.423
85 -112771.362 19942.976
86 60434.071 -112771.362
87 -6906.064 60434.071
88 -83722.349 -6906.064
89 24236.170 -83722.349
90 20323.606 24236.170
91 -58373.820 20323.606
92 152629.298 -58373.820
93 -249808.092 152629.298
94 97881.580 -249808.092
95 -18948.548 97881.580
96 56121.615 -18948.548
97 27844.604 56121.615
98 -65081.319 27844.604
99 -65872.062 -65081.319
100 -69214.165 -65872.062
101 -109444.008 -69214.165
102 -127969.710 -109444.008
103 -65383.903 -127969.710
104 51398.186 -65383.903
105 -31544.333 51398.186
106 -266560.840 -31544.333
107 63483.018 -266560.840
108 -169212.146 63483.018
109 -139539.338 -169212.146
110 111838.134 -139539.338
111 -223309.564 111838.134
112 -148255.495 -223309.564
113 -124797.099 -148255.495
114 23515.058 -124797.099
115 55479.808 23515.058
116 17147.822 55479.808
117 63553.861 17147.822
118 -200546.755 63553.861
119 50336.225 -200546.755
120 -265684.483 50336.225
121 34349.869 -265684.483
122 22362.228 34349.869
123 29543.810 22362.228
124 -37960.116 29543.810
125 -223683.527 -37960.116
126 47930.758 -223683.527
127 46990.954 47930.758
128 54791.004 46990.954
129 47930.758 54791.004
130 49320.735 47930.758
131 90994.114 49320.735
132 47930.758 90994.114
133 45014.662 47930.758
134 27433.327 45014.662
135 42217.914 27433.327
136 52885.045 42217.914
137 48103.003 52885.045
138 22758.922 48103.003
139 47930.758 22758.922
140 55492.576 47930.758
141 47930.758 55492.576
142 15002.087 47930.758
143 48197.857 15002.087
144 49670.808 48197.857
145 48756.725 49670.808
146 47930.758 48756.725
147 49490.851 47930.758
148 50535.241 49490.851
149 47847.411 50535.241
150 16757.609 47847.411
151 16128.433 16757.609
152 47930.758 16128.433
153 25176.204 47930.758
154 48001.084 25176.204
155 47930.758 48001.084
156 47930.758 47930.758
157 44196.429 47930.758
158 59511.531 44196.429
159 47930.758 59511.531
160 47930.758 47930.758
161 47930.758 47930.758
162 47930.758 47930.758
163 35201.158 47930.758
164 37073.804 35201.158
165 47253.085 37073.804
166 47930.758 47253.085
167 129383.608 47930.758
168 95359.943 129383.608
169 47930.758 95359.943
170 47930.758 47930.758
171 43892.484 47930.758
172 47930.758 43892.484
173 37413.394 47930.758
174 47930.758 37413.394
175 47968.115 47930.758
176 95780.130 47968.115
177 -8516.655 95780.130
178 47930.758 -8516.655
179 47930.758 47930.758
180 55714.984 47930.758
181 -25660.802 55714.984
182 6174.204 -25660.802
183 -37101.259 6174.204
184 -10313.459 -37101.259
185 -2370.308 -10313.459
186 -35584.633 -2370.308
187 40932.200 -35584.633
188 9957.267 40932.200
189 47930.758 9957.267
190 47076.205 47930.758
191 47930.758 47076.205
192 44403.489 47930.758
193 47930.758 44403.489
194 47890.901 47930.758
195 -7555.060 47890.901
196 -126404.718 -7555.060
197 -7555.060 -126404.718
198 -7555.060 -7555.060
199 -96498.844 -7555.060
200 -7555.060 -96498.844
201 -7555.060 -7555.060
202 -7555.060 -7555.060
203 2068.080 -7555.060
204 -56691.365 2068.080
205 -8066.752 -56691.365
206 -14179.696 -8066.752
207 -7555.060 -14179.696
208 23336.527 -7555.060
209 -11265.356 23336.527
210 47930.758 -11265.356
211 47930.758 47930.758
212 73611.546 47930.758
213 13127.108 73611.546
214 20306.970 13127.108
215 47930.758 20306.970
216 50834.989 47930.758
217 49994.478 50834.989
218 47930.758 49994.478
219 19719.600 47930.758
220 22654.245 19719.600
221 47930.758 22654.245
222 48672.001 47930.758
223 39310.301 48672.001
224 48181.113 39310.301
225 48453.678 48181.113
226 50983.381 48453.678
227 47930.758 50983.381
228 32019.062 47930.758
229 47930.758 32019.062
230 19310.917 47930.758
231 35997.643 19310.917
232 47930.758 35997.643
233 29871.149 47930.758
234 47930.758 29871.149
235 48340.985 47930.758
236 47930.758 48340.985
237 49921.272 47930.758
238 47930.758 49921.272
239 47930.758 47930.758
240 47930.758 47930.758
241 42624.223 47930.758
242 47930.758 42624.223
243 52466.884 47930.758
244 -33992.012 52466.884
245 -7555.060 -33992.012
246 -7555.060 -7555.060
247 40017.496 -7555.060
248 -7555.060 40017.496
249 -7555.060 -7555.060
250 41434.050 -7555.060
251 -7555.060 41434.050
252 28918.905 -7555.060
253 24439.537 28918.905
254 -87387.956 24439.537
255 47930.758 -87387.956
256 45851.206 47930.758
257 59484.659 45851.206
258 -7555.060 59484.659
259 47930.758 -7555.060
260 -7555.060 47930.758
261 -7555.060 -7555.060
262 51292.884 -7555.060
263 65051.344 51292.884
264 -7555.060 65051.344
265 47930.758 -7555.060
266 -7555.060 47930.758
267 79985.494 -7555.060
268 -7555.060 79985.494
269 44777.497 -7555.060
270 -7555.060 44777.497
271 -7555.060 -7555.060
272 -7555.060 -7555.060
273 43760.928 -7555.060
274 47930.758 43760.928
275 53038.424 47930.758
276 49743.121 53038.424
277 47930.758 49743.121
278 -35871.190 47930.758
279 35089.843 -35871.190
280 -7555.060 35089.843
281 47930.758 -7555.060
282 75626.658 47930.758
283 47930.758 75626.658
284 39802.780 47930.758
285 47930.758 39802.780
286 -19198.651 47930.758
287 48830.870 -19198.651
288 -7555.060 48830.870
289 55446.705 -7555.060
290 -7555.060 55446.705
291 58494.152 -7555.060
292 54252.289 58494.152
293 47930.758 54252.289
294 41418.222 47930.758
295 17805.179 41418.222
296 47930.758 17805.179
297 -7555.060 47930.758
298 47930.758 -7555.060
299 47930.758 47930.758
300 -9187.034 47930.758
301 29601.307 -9187.034
302 47930.758 29601.307
303 -7555.060 47930.758
304 85548.186 -7555.060
305 51520.595 85548.186
306 -7555.060 51520.595
307 -13574.240 -7555.060
308 51524.780 -13574.240
309 -7555.060 51524.780
310 -7555.060 -7555.060
311 51085.001 -7555.060
312 -7555.060 51085.001
313 -5631.489 -7555.060
314 49666.953 -5631.489
315 63580.147 49666.953
316 48377.897 63580.147
317 24781.823 48377.897
318 50337.774 24781.823
319 47930.758 50337.774
320 47930.758 47930.758
321 47930.758 47930.758
322 47930.758 47930.758
323 -7555.060 47930.758
324 47930.758 -7555.060
325 45207.361 47930.758
326 -7555.060 45207.361
327 -7555.060 -7555.060
328 47930.758 -7555.060
329 47930.758 47930.758
330 47930.758 47930.758
331 47930.758 47930.758
332 47930.758 47930.758
333 46110.459 47930.758
334 48709.764 46110.459
335 40596.813 48709.764
336 34892.012 40596.813
337 47930.758 34892.012
338 45723.816 47930.758
339 53714.507 45723.816
340 20843.505 53714.507
341 47930.758 20843.505
342 47930.758 47930.758
343 47930.758 47930.758
344 41837.252 47930.758
345 47930.758 41837.252
346 38208.833 47930.758
347 47930.758 38208.833
348 47930.758 47930.758
349 11292.552 47930.758
350 51727.086 11292.552
351 48104.426 51727.086
352 47803.348 48104.426
353 40959.878 47803.348
354 47930.758 40959.878
355 47930.758 47930.758
356 25991.297 47930.758
357 47930.758 25991.297
358 47930.758 47930.758
359 47930.758 47930.758
360 49543.695 47930.758
361 47930.758 49543.695
362 47930.758 47930.758
363 -59079.502 47930.758
364 50139.283 -59079.502
365 47930.758 50139.283
366 47930.758 47930.758
367 47930.758 47930.758
368 47930.758 47930.758
369 50155.084 47930.758
370 47930.758 50155.084
371 60651.036 47930.758
372 35414.901 60651.036
373 27297.538 35414.901
374 47930.758 27297.538
375 47930.758 47930.758
376 42833.808 47930.758
377 41821.079 42833.808
378 -7827.369 41821.079
379 47930.758 -7827.369
380 47930.758 47930.758
381 33983.667 47930.758
382 65563.601 33983.667
383 47930.758 65563.601
384 50256.107 47930.758
385 -100883.026 50256.107
386 60199.198 -100883.026
387 21375.207 60199.198
388 46456.029 21375.207
389 34065.162 46456.029
390 43281.192 34065.162
391 17215.934 43281.192
392 -20308.732 17215.934
393 68636.060 -20308.732
394 -180351.327 68636.060
395 9037.666 -180351.327
396 -138046.328 9037.666
397 4172.197 -138046.328
398 10638.634 4172.197
399 -23000.886 10638.634
400 -11465.171 -23000.886
401 39825.868 -11465.171
402 17459.373 39825.868
403 -95882.735 17459.373
404 -12187.230 -95882.735
405 10402.487 -12187.230
406 23566.412 10402.487
407 15549.014 23566.412
408 26690.904 15549.014
409 -73845.864 26690.904
410 10354.181 -73845.864
411 -250216.313 10354.181
412 -74411.743 -250216.313
413 12120.541 -74411.743
414 -189194.453 12120.541
415 -115267.526 -189194.453
416 118575.766 -115267.526
417 -62387.577 118575.766
418 -35757.824 -62387.577
419 -49407.953 -35757.824
420 -60797.317 -49407.953
421 -107098.631 -60797.317
422 -576390.302 -107098.631
423 -484134.906 -576390.302
424 -41774.594 -484134.906
425 -232554.023 -41774.594
426 -313458.634 -232554.023
427 -1118357.214 -313458.634
428 -306374.089 -1118357.214
429 -344582.470 -306374.089
430 -675705.948 -344582.470
> 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/74cys1291414557.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/84cys1291414557.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/94cys1291414557.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/10x3fd1291414557.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/11i3e11291414557.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/1234cp1291414557.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/13zdag1291414557.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/14le9m1291414557.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/15oxpa1291414557.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/16ax6y1291414557.tab")
+ }
>
> try(system("convert tmp/1qki21291414557.ps tmp/1qki21291414557.png",intern=TRUE))
character(0)
> try(system("convert tmp/2qki21291414557.ps tmp/2qki21291414557.png",intern=TRUE))
character(0)
> try(system("convert tmp/3jbi51291414557.ps tmp/3jbi51291414557.png",intern=TRUE))
character(0)
> try(system("convert tmp/4jbi51291414557.ps tmp/4jbi51291414557.png",intern=TRUE))
character(0)
> try(system("convert tmp/5jbi51291414557.ps tmp/5jbi51291414557.png",intern=TRUE))
character(0)
> try(system("convert tmp/6c2hq1291414557.ps tmp/6c2hq1291414557.png",intern=TRUE))
character(0)
> try(system("convert tmp/74cys1291414557.ps tmp/74cys1291414557.png",intern=TRUE))
character(0)
> try(system("convert tmp/84cys1291414557.ps tmp/84cys1291414557.png",intern=TRUE))
character(0)
> try(system("convert tmp/94cys1291414557.ps tmp/94cys1291414557.png",intern=TRUE))
character(0)
> try(system("convert tmp/10x3fd1291414557.ps tmp/10x3fd1291414557.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
13.963 3.164 14.377