R version 2.9.0 (2009-04-17)
Copyright (C) 2009 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.
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
+ ,807
+ ,213118
+ ,6282154
+ ,0
+ ,29790
+ ,444
+ ,81767
+ ,4321023
+ ,0
+ ,87550
+ ,412
+ ,153198
+ ,4111912
+ ,1
+ ,84738
+ ,428
+ ,-26007
+ ,223193
+ ,0
+ ,54660
+ ,315
+ ,126942
+ ,1491348
+ ,0
+ ,42634
+ ,168
+ ,157214
+ ,1629616
+ ,1
+ ,40949
+ ,263
+ ,129352
+ ,1398893
+ ,0
+ ,45187
+ ,267
+ ,234817
+ ,1926517
+ ,0
+ ,37704
+ ,228
+ ,60448
+ ,983660
+ ,0
+ ,16275
+ ,129
+ ,47818
+ ,1443586
+ ,1
+ ,25830
+ ,104
+ ,245546
+ ,1073089
+ ,1
+ ,12679
+ ,122
+ ,48020
+ ,984885
+ ,0
+ ,18014
+ ,393
+ ,-1710
+ ,1405225
+ ,1
+ ,43556
+ ,190
+ ,32648
+ ,227132
+ ,0
+ ,24811
+ ,280
+ ,95350
+ ,929118
+ ,1
+ ,6575
+ ,63
+ ,151352
+ ,1071292
+ ,1
+ ,7123
+ ,102
+ ,288170
+ ,638830
+ ,0
+ ,21950
+ ,265
+ ,114337
+ ,856956
+ ,0
+ ,37597
+ ,234
+ ,37884
+ ,992426
+ ,1
+ ,17821
+ ,277
+ ,122844
+ ,444477
+ ,0
+ ,12988
+ ,73
+ ,82340
+ ,857217
+ ,0
+ ,22330
+ ,67
+ ,79801
+ ,711969
+ ,1
+ ,13326
+ ,103
+ ,165548
+ ,702380
+ ,1
+ ,16189
+ ,290
+ ,116384
+ ,358589
+ ,1
+ ,7146
+ ,83
+ ,134028
+ ,297978
+ ,1
+ ,15824
+ ,56
+ ,63838
+ ,585715
+ ,0
+ ,27664
+ ,236
+ ,74996
+ ,657954
+ ,1
+ ,11920
+ ,73
+ ,31080
+ ,209458
+ ,1
+ ,8568
+ ,34
+ ,32168
+ ,786690
+ ,1
+ ,14416
+ ,139
+ ,49857
+ ,439798
+ ,0
+ ,3369
+ ,26
+ ,87161
+ ,688779
+ ,0
+ ,11819
+ ,70
+ ,106113
+ ,574339
+ ,0
+ ,6984
+ ,40
+ ,80570
+ ,741409
+ ,0
+ ,4519
+ ,42
+ ,102129
+ ,597793
+ ,1
+ ,2220
+ ,12
+ ,301670
+ ,644190
+ ,1
+ ,18562
+ ,211
+ ,102313
+ ,377934
+ ,1
+ ,10327
+ ,74
+ ,88577
+ ,640273
+ ,0
+ ,5336
+ ,80
+ ,112477
+ ,697458
+ ,0
+ ,2365
+ ,83
+ ,191778
+ ,550608
+ ,1
+ ,4069
+ ,131
+ ,79804
+ ,207393
+ ,1
+ ,8636
+ ,203
+ ,128294
+ ,301607
+ ,1
+ ,13718
+ ,56
+ ,96448
+ ,345783
+ ,1
+ ,4525
+ ,89
+ ,93811
+ ,501749
+ ,1
+ ,6869
+ ,88
+ ,117520
+ ,379983
+ ,1
+ ,4628
+ ,39
+ ,69159
+ ,387475
+ ,0
+ ,3689
+ ,25
+ ,101792
+ ,377305
+ ,0
+ ,4891
+ ,49
+ ,210568
+ ,370837
+ ,0
+ ,7489
+ ,149
+ ,136996
+ ,430866
+ ,1
+ ,4901
+ ,58
+ ,121920
+ ,469107
+ ,1
+ ,2284
+ ,41
+ ,76403
+ ,194493
+ ,0
+ ,3160
+ ,90
+ ,108094
+ ,530670
+ ,0
+ ,4150
+ ,136
+ ,134759
+ ,518365
+ ,0
+ ,7285
+ ,97
+ ,188873
+ ,491303
+ ,0
+ ,1134
+ ,63
+ ,146216
+ ,527021
+ ,0
+ ,4658
+ ,114
+ ,156608
+ ,233773
+ ,1
+ ,2384
+ ,77
+ ,61348
+ ,405972
+ ,1
+ ,3748
+ ,6
+ ,50350
+ ,652925
+ ,1
+ ,5371
+ ,47
+ ,87720
+ ,446211
+ ,1
+ ,1285
+ ,51
+ ,99489
+ ,341340
+ ,0
+ ,9327
+ ,85
+ ,87419
+ ,387699
+ ,0
+ ,5565
+ ,43
+ ,94355
+ ,493408
+ ,1
+ ,1528
+ ,32
+ ,60326
+ ,146494
+ ,0
+ ,3122
+ ,25
+ ,94670
+ ,414462
+ ,0
+ ,7561
+ ,77
+ ,82425
+ ,364304
+ ,1
+ ,2675
+ ,54
+ ,59017
+ ,355178
+ ,1
+ ,13253
+ ,251
+ ,90829
+ ,357760
+ ,1
+ ,880
+ ,15
+ ,80791
+ ,261216
+ ,0
+ ,2053
+ ,44
+ ,100423
+ ,397144
+ ,1
+ ,1424
+ ,73
+ ,131116
+ ,374943
+ ,0
+ ,4036
+ ,85
+ ,100269
+ ,424898
+ ,0
+ ,3045
+ ,49
+ ,27330
+ ,202055
+ ,1
+ ,5119
+ ,38
+ ,39039
+ ,378525
+ ,1
+ ,1431
+ ,35
+ ,106885
+ ,310768
+ ,1
+ ,554
+ ,9
+ ,79285
+ ,325738
+ ,1
+ ,1975
+ ,34
+ ,118881
+ ,394510
+ ,0
+ ,1765
+ ,20
+ ,77623
+ ,247060
+ ,1
+ ,1012
+ ,29
+ ,114768
+ ,368078
+ ,1
+ ,810
+ ,11
+ ,74015
+ ,236761
+ ,1
+ ,1280
+ ,52
+ ,69465
+ ,312378
+ ,0
+ ,666
+ ,13
+ ,117869
+ ,339836
+ ,1
+ ,1380
+ ,29
+ ,60982
+ ,347385
+ ,0
+ ,4677
+ ,66
+ ,90131
+ ,426280
+ ,1
+ ,876
+ ,33
+ ,138971
+ ,352850
+ ,1
+ ,814
+ ,15
+ ,39625
+ ,301881
+ ,1
+ ,514
+ ,15
+ ,102725
+ ,377516
+ ,0
+ ,5692
+ ,68
+ ,64239
+ ,357312
+ ,1
+ ,3642
+ ,100
+ ,90262
+ ,458343
+ ,1
+ ,540
+ ,13
+ ,103960
+ ,354228
+ ,1
+ ,2099
+ ,45
+ ,106611
+ ,308636
+ ,1
+ ,567
+ ,14
+ ,103345
+ ,386212
+ ,1
+ ,2001
+ ,36
+ ,95551
+ ,393343
+ ,0
+ ,2949
+ ,40
+ ,82903
+ ,378509
+ ,1
+ ,2253
+ ,68
+ ,63593
+ ,452469
+ ,0
+ ,6533
+ ,29
+ ,126910
+ ,364839
+ ,1
+ ,1889
+ ,43
+ ,37527
+ ,358649
+ ,0
+ ,3055
+ ,30
+ ,60247
+ ,376641
+ ,1
+ ,272
+ ,9
+ ,112995
+ ,429112
+ ,0
+ ,1414
+ ,22
+ ,70184
+ ,330546
+ ,1
+ ,2564
+ ,19
+ ,130140
+ ,403560
+ ,0
+ ,1383
+ ,9
+ ,73221
+ ,317892
+ ,0
+ ,1261
+ ,31
+ ,76114
+ ,307528
+ ,1
+ ,975
+ ,19
+ ,90534
+ ,235133
+ ,1
+ ,3366
+ ,55
+ ,108479
+ ,299243
+ ,1
+ ,576
+ ,8
+ ,113761
+ ,314073
+ ,1
+ ,1686
+ ,28
+ ,68696
+ ,368186
+ ,1
+ ,746
+ ,29
+ ,71561
+ ,269661
+ ,0
+ ,3192
+ ,48
+ ,59831
+ ,125390
+ ,0
+ ,2045
+ ,16
+ ,97890
+ ,510834
+ ,1
+ ,5702
+ ,47
+ ,101481
+ ,321896
+ ,0
+ ,1932
+ ,20
+ ,72954
+ ,249898
+ ,1
+ ,936
+ ,22
+ ,67939
+ ,408881
+ ,0
+ ,3437
+ ,33
+ ,48022
+ ,158492
+ ,1
+ ,5131
+ ,44
+ ,86111
+ ,292154
+ ,0
+ ,2397
+ ,13
+ ,74020
+ ,289513
+ ,0
+ ,1389
+ ,6
+ ,57530
+ ,378049
+ ,1
+ ,1503
+ ,35
+ ,56364
+ ,343466
+ ,1
+ ,402
+ ,8
+ ,84990
+ ,332743
+ ,1
+ ,2239
+ ,17
+ ,88590
+ ,442882
+ ,0
+ ,2234
+ ,11
+ ,77200
+ ,214215
+ ,1
+ ,837
+ ,21
+ ,61262
+ ,315688
+ ,1
+ ,10579
+ ,92
+ ,110309
+ ,375195
+ ,1
+ ,875
+ ,12
+ ,67000
+ ,334280
+ ,1
+ ,1585
+ ,112
+ ,93099
+ ,355864
+ ,0
+ ,1659
+ ,25
+ ,107577
+ ,480382
+ ,0
+ ,2647
+ ,17
+ ,62920
+ ,353058
+ ,0
+ ,3294
+ ,23
+ ,75832
+ ,217193
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,94
+ ,10
+ ,60793
+ ,314533
+ ,1
+ ,422
+ ,23
+ ,57935
+ ,318056
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,34
+ ,7
+ ,60630
+ ,314353
+ ,1
+ ,1558
+ ,25
+ ,55637
+ ,369448
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,43
+ ,20
+ ,60887
+ ,312846
+ ,1
+ ,645
+ ,4
+ ,60720
+ ,312075
+ ,1
+ ,316
+ ,4
+ ,60505
+ ,315009
+ ,1
+ ,115
+ ,10
+ ,60945
+ ,318903
+ ,1
+ ,5
+ ,1
+ ,60720
+ ,314887
+ ,1
+ ,897
+ ,4
+ ,60720
+ ,314913
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,389
+ ,8
+ ,58990
+ ,325506
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,1002
+ ,11
+ ,56750
+ ,298568
+ ,1
+ ,36
+ ,4
+ ,60894
+ ,315834
+ ,1
+ ,460
+ ,15
+ ,63346
+ ,329784
+ ,1
+ ,309
+ ,9
+ ,56535
+ ,312878
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,9
+ ,7
+ ,60835
+ ,314987
+ ,1
+ ,271
+ ,2
+ ,60720
+ ,325249
+ ,1
+ ,14
+ ,0
+ ,61016
+ ,315877
+ ,1
+ ,520
+ ,7
+ ,58650
+ ,291650
+ ,1
+ ,1766
+ ,46
+ ,60438
+ ,305959
+ ,1
+ ,0
+ ,5
+ ,60720
+ ,315380
+ ,1
+ ,458
+ ,7
+ ,58625
+ ,297765
+ ,1
+ ,20
+ ,2
+ ,60938
+ ,315245
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,98
+ ,2
+ ,61490
+ ,315236
+ ,1
+ ,405
+ ,5
+ ,60845
+ ,336425
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,483
+ ,7
+ ,60830
+ ,306268
+ ,1
+ ,454
+ ,24
+ ,63261
+ ,302187
+ ,1
+ ,47
+ ,1
+ ,60720
+ ,314882
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,757
+ ,18
+ ,45689
+ ,382712
+ ,1
+ ,4655
+ ,55
+ ,60720
+ ,341570
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,36
+ ,3
+ ,61564
+ ,312412
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,203
+ ,9
+ ,61938
+ ,309596
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,126
+ ,8
+ ,60951
+ ,315547
+ ,1
+ ,400
+ ,113
+ ,60720
+ ,313267
+ ,0
+ ,71
+ ,0
+ ,60745
+ ,316176
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,972
+ ,19
+ ,71642
+ ,359335
+ ,0
+ ,531
+ ,11
+ ,71641
+ ,330068
+ ,1
+ ,2461
+ ,25
+ ,55792
+ ,314289
+ ,0
+ ,378
+ ,16
+ ,71873
+ ,297413
+ ,0
+ ,23
+ ,5
+ ,62555
+ ,314806
+ ,0
+ ,638
+ ,11
+ ,60370
+ ,333210
+ ,0
+ ,2300
+ ,23
+ ,64873
+ ,352108
+ ,1
+ ,149
+ ,6
+ ,62041
+ ,313332
+ ,1
+ ,226
+ ,5
+ ,65745
+ ,291787
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,275
+ ,7
+ ,59500
+ ,318745
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,141
+ ,7
+ ,61630
+ ,315366
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,28
+ ,3
+ ,60890
+ ,315688
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,4980
+ ,89
+ ,113521
+ ,409642
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,472
+ ,19
+ ,80045
+ ,269587
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,203
+ ,12
+ ,50804
+ ,300962
+ ,0
+ ,496
+ ,12
+ ,87390
+ ,325479
+ ,0
+ ,10
+ ,5
+ ,61656
+ ,316155
+ ,0
+ ,63
+ ,2
+ ,65688
+ ,318574
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1136
+ ,26
+ ,48522
+ ,343613
+ ,0
+ ,265
+ ,3
+ ,60720
+ ,306948
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,267
+ ,11
+ ,57640
+ ,330059
+ ,1
+ ,474
+ ,10
+ ,61977
+ ,288985
+ ,1
+ ,534
+ ,5
+ ,62620
+ ,304485
+ ,1
+ ,0
+ ,2
+ ,60720
+ ,315380
+ ,1
+ ,15
+ ,6
+ ,60831
+ ,315688
+ ,1
+ ,397
+ ,7
+ ,60646
+ ,317736
+ ,1
+ ,0
+ ,2
+ ,60720
+ ,315380
+ ,1
+ ,1866
+ ,28
+ ,56225
+ ,322331
+ ,1
+ ,288
+ ,3
+ ,60510
+ ,296656
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,3
+ ,1
+ ,60698
+ ,315354
+ ,1
+ ,468
+ ,20
+ ,60720
+ ,312161
+ ,1
+ ,20
+ ,1
+ ,60805
+ ,315576
+ ,1
+ ,278
+ ,22
+ ,61404
+ ,314922
+ ,1
+ ,61
+ ,9
+ ,60720
+ ,314551
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,192
+ ,2
+ ,65276
+ ,312339
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,317
+ ,7
+ ,63915
+ ,298700
+ ,1
+ ,738
+ ,9
+ ,60720
+ ,321376
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,368
+ ,13
+ ,61686
+ ,303230
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,2
+ ,0
+ ,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
+ ,7
+ ,59818
+ ,315637
+ ,1
+ ,2332
+ ,2
+ ,72680
+ ,324385
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,131
+ ,15
+ ,61808
+ ,308989
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,206
+ ,9
+ ,53110
+ ,296702
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,167
+ ,1
+ ,64245
+ ,307322
+ ,1
+ ,622
+ ,38
+ ,73007
+ ,304376
+ ,1
+ ,2328
+ ,57
+ ,82732
+ ,253588
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,365
+ ,7
+ ,54820
+ ,309560
+ ,1
+ ,364
+ ,26
+ ,47705
+ ,298466
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,226
+ ,13
+ ,72835
+ ,343929
+ ,1
+ ,307
+ ,10
+ ,58856
+ ,331955
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,188
+ ,9
+ ,77655
+ ,381180
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,138
+ ,26
+ ,69817
+ ,331420
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,125
+ ,19
+ ,60798
+ ,310201
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,282
+ ,12
+ ,62452
+ ,320016
+ ,1
+ ,335
+ ,23
+ ,64175
+ ,320398
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1324
+ ,29
+ ,67440
+ ,291841
+ ,1
+ ,176
+ ,8
+ ,68136
+ ,310670
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,249
+ ,26
+ ,56726
+ ,313491
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,333
+ ,9
+ ,70811
+ ,331323
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,601
+ ,5
+ ,60720
+ ,319210
+ ,1
+ ,30
+ ,3
+ ,62045
+ ,318098
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,249
+ ,13
+ ,54323
+ ,292754
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,165
+ ,12
+ ,62841
+ ,325176
+ ,1
+ ,453
+ ,19
+ ,81125
+ ,365959
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,53
+ ,10
+ ,59506
+ ,302409
+ ,0
+ ,382
+ ,9
+ ,59365
+ ,340968
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,9
+ ,60720
+ ,315380
+ ,0
+ ,30
+ ,4
+ ,60798
+ ,313164
+ ,1
+ ,290
+ ,1
+ ,58790
+ ,301164
+ ,1
+ ,0
+ ,1
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,366
+ ,14
+ ,61808
+ ,344425
+ ,1
+ ,2
+ ,12
+ ,60735
+ ,315394
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,209
+ ,19
+ ,64016
+ ,316647
+ ,1
+ ,384
+ ,17
+ ,54683
+ ,309836
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,365
+ ,32
+ ,87192
+ ,346611
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,49
+ ,14
+ ,64107
+ ,322031
+ ,1
+ ,3
+ ,8
+ ,60761
+ ,315656
+ ,1
+ ,133
+ ,4
+ ,65990
+ ,339445
+ ,1
+ ,32
+ ,0
+ ,59988
+ ,314964
+ ,1
+ ,368
+ ,20
+ ,61167
+ ,297141
+ ,1
+ ,1
+ ,5
+ ,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
+ ,1
+ ,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
+ ,4
+ ,60379
+ ,313729
+ ,1
+ ,1
+ ,1
+ ,60727
+ ,315388
+ ,1
+ ,314
+ ,4
+ ,60720
+ ,315371
+ ,1
+ ,844
+ ,20
+ ,60925
+ ,296139
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,26
+ ,1
+ ,60896
+ ,313880
+ ,1
+ ,125
+ ,10
+ ,59734
+ ,317698
+ ,1
+ ,304
+ ,12
+ ,62969
+ ,295580
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,621
+ ,13
+ ,60720
+ ,308256
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,119
+ ,3
+ ,59118
+ ,303677
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,1595
+ ,10
+ ,60720
+ ,319369
+ ,1
+ ,312
+ ,3
+ ,58598
+ ,318690
+ ,1
+ ,60
+ ,7
+ ,61124
+ ,314049
+ ,1
+ ,587
+ ,10
+ ,59595
+ ,325699
+ ,1
+ ,135
+ ,1
+ ,62065
+ ,314210
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,514
+ ,15
+ ,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
+ ,28
+ ,61600
+ ,308336
+ ,1
+ ,180
+ ,9
+ ,59635
+ ,316386
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,218
+ ,7
+ ,60720
+ ,315553
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,448
+ ,7
+ ,59781
+ ,323361
+ ,1
+ ,227
+ ,7
+ ,76644
+ ,336639
+ ,1
+ ,174
+ ,3
+ ,64820
+ ,307424
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,121
+ ,11
+ ,56178
+ ,295370
+ ,1
+ ,607
+ ,7
+ ,60436
+ ,322340
+ ,1
+ ,2212
+ ,10
+ ,60720
+ ,319864
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,530
+ ,18
+ ,73433
+ ,317291
+ ,1
+ ,571
+ ,14
+ ,41477
+ ,280398
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,78
+ ,12
+ ,62700
+ ,317330
+ ,1
+ ,2489
+ ,29
+ ,67804
+ ,238125
+ ,1
+ ,131
+ ,3
+ ,59661
+ ,327071
+ ,1
+ ,923
+ ,6
+ ,58620
+ ,309038
+ ,1
+ ,72
+ ,3
+ ,60398
+ ,314210
+ ,1
+ ,572
+ ,8
+ ,58580
+ ,307930
+ ,1
+ ,397
+ ,10
+ ,62710
+ ,322327
+ ,1
+ ,450
+ ,6
+ ,59325
+ ,292136
+ ,1
+ ,622
+ ,8
+ ,60950
+ ,263276
+ ,1
+ ,694
+ ,6
+ ,68060
+ ,367655
+ ,0
+ ,3425
+ ,9
+ ,83620
+ ,283910
+ ,1
+ ,562
+ ,8
+ ,58456
+ ,283587
+ ,1
+ ,4917
+ ,26
+ ,52811
+ ,243650
+ ,0
+ ,1442
+ ,239
+ ,121173
+ ,438493
+ ,1
+ ,529
+ ,7
+ ,63870
+ ,296261
+ ,0
+ ,2126
+ ,41
+ ,21001
+ ,230621
+ ,1
+ ,1061
+ ,3
+ ,70415
+ ,304252
+ ,1
+ ,776
+ ,8
+ ,64230
+ ,333505
+ ,1
+ ,611
+ ,6
+ ,59190
+ ,296919
+ ,0
+ ,1526
+ ,21
+ ,69351
+ ,278990
+ ,1
+ ,592
+ ,7
+ ,64270
+ ,276898
+ ,1
+ ,1182
+ ,11
+ ,70694
+ ,327007
+ ,1
+ ,621
+ ,11
+ ,68005
+ ,317046
+ ,1
+ ,989
+ ,12
+ ,58930
+ ,304555
+ ,1
+ ,438
+ ,9
+ ,58320
+ ,298096
+ ,1
+ ,726
+ ,3
+ ,69980
+ ,231861
+ ,1
+ ,1303
+ ,57
+ ,69863
+ ,309422
+ ,0
+ ,7419
+ ,21
+ ,63255
+ ,286963
+ ,0
+ ,1164
+ ,15
+ ,57320
+ ,269753
+ ,0
+ ,3310
+ ,32
+ ,75230
+ ,448243
+ ,1
+ ,1920
+ ,11
+ ,79420
+ ,165404
+ ,1
+ ,965
+ ,2
+ ,73490
+ ,204325
+ ,1
+ ,3256
+ ,23
+ ,35250
+ ,407159
+ ,0
+ ,1135
+ ,20
+ ,62285
+ ,290476
+ ,1
+ ,1270
+ ,24
+ ,69206
+ ,275311
+ ,1
+ ,661
+ ,1
+ ,65920
+ ,246541
+ ,1
+ ,1013
+ ,1
+ ,69770
+ ,253468
+ ,1
+ ,2844
+ ,74
+ ,72683
+ ,240897
+ ,0
+ ,11528
+ ,68
+ ,-14545
+ ,-83265
+ ,1
+ ,6526
+ ,20
+ ,55830
+ ,-42143
+ ,1
+ ,2264
+ ,20
+ ,55174
+ ,272713
+ ,0
+ ,5109
+ ,82
+ ,67038
+ ,215362
+ ,1
+ ,3999
+ ,21
+ ,51252
+ ,42754
+ ,1
+ ,35624
+ ,244
+ ,157278
+ ,306275
+ ,1
+ ,9252
+ ,32
+ ,79510
+ ,253537
+ ,1
+ ,15236
+ ,86
+ ,77440
+ ,372631
+ ,1
+ ,18073
+ ,69
+ ,27284
+ ,-7170)
+ ,dim=c(5
+ ,431)
+ ,dimnames=list(c('group'
+ ,'costs'
+ ,'orders'
+ ,'dividends'
+ ,'wealth
')
+ ,1:431))
> y <- array(NA,dim=c(5,431),dimnames=list(c('group','costs','orders','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\r group costs orders dividends
1 6282154 0 162556 807 213118
2 4321023 0 29790 444 81767
3 4111912 0 87550 412 153198
4 223193 1 84738 428 -26007
5 1491348 0 54660 315 126942
6 1629616 0 42634 168 157214
7 1398893 1 40949 263 129352
8 1926517 0 45187 267 234817
9 983660 0 37704 228 60448
10 1443586 0 16275 129 47818
11 1073089 1 25830 104 245546
12 984885 1 12679 122 48020
13 1405225 0 18014 393 -1710
14 227132 1 43556 190 32648
15 929118 0 24811 280 95350
16 1071292 1 6575 63 151352
17 638830 1 7123 102 288170
18 856956 0 21950 265 114337
19 992426 0 37597 234 37884
20 444477 1 17821 277 122844
21 857217 0 12988 73 82340
22 711969 0 22330 67 79801
23 702380 1 13326 103 165548
24 358589 1 16189 290 116384
25 297978 1 7146 83 134028
26 585715 1 15824 56 63838
27 657954 0 27664 236 74996
28 209458 1 11920 73 31080
29 786690 1 8568 34 32168
30 439798 1 14416 139 49857
31 688779 0 3369 26 87161
32 574339 0 11819 70 106113
33 741409 0 6984 40 80570
34 597793 0 4519 42 102129
35 644190 1 2220 12 301670
36 377934 1 18562 211 102313
37 640273 1 10327 74 88577
38 697458 0 5336 80 112477
39 550608 0 2365 83 191778
40 207393 1 4069 131 79804
41 301607 1 8636 203 128294
42 345783 1 13718 56 96448
43 501749 1 4525 89 93811
44 379983 1 6869 88 117520
45 387475 1 4628 39 69159
46 377305 0 3689 25 101792
47 370837 0 4891 49 210568
48 430866 0 7489 149 136996
49 469107 1 4901 58 121920
50 194493 1 2284 41 76403
51 530670 0 3160 90 108094
52 518365 0 4150 136 134759
53 491303 0 7285 97 188873
54 527021 0 1134 63 146216
55 233773 0 4658 114 156608
56 405972 1 2384 77 61348
57 652925 1 3748 6 50350
58 446211 1 5371 47 87720
59 341340 1 1285 51 99489
60 387699 0 9327 85 87419
61 493408 0 5565 43 94355
62 146494 1 1528 32 60326
63 414462 0 3122 25 94670
64 364304 0 7561 77 82425
65 355178 1 2675 54 59017
66 357760 1 13253 251 90829
67 261216 1 880 15 80791
68 397144 0 2053 44 100423
69 374943 1 1424 73 131116
70 424898 0 4036 85 100269
71 202055 0 3045 49 27330
72 378525 1 5119 38 39039
73 310768 1 1431 35 106885
74 325738 1 554 9 79285
75 394510 1 1975 34 118881
76 247060 0 1765 20 77623
77 368078 1 1012 29 114768
78 236761 1 810 11 74015
79 312378 1 1280 52 69465
80 339836 0 666 13 117869
81 347385 1 1380 29 60982
82 426280 0 4677 66 90131
83 352850 1 876 33 138971
84 301881 1 814 15 39625
85 377516 1 514 15 102725
86 357312 0 5692 68 64239
87 458343 1 3642 100 90262
88 354228 1 540 13 103960
89 308636 1 2099 45 106611
90 386212 1 567 14 103345
91 393343 1 2001 36 95551
92 378509 0 2949 40 82903
93 452469 1 2253 68 63593
94 364839 0 6533 29 126910
95 358649 1 1889 43 37527
96 376641 0 3055 30 60247
97 429112 1 272 9 112995
98 330546 0 1414 22 70184
99 403560 1 2564 19 130140
100 317892 0 1383 9 73221
101 307528 0 1261 31 76114
102 235133 1 975 19 90534
103 299243 1 3366 55 108479
104 314073 1 576 8 113761
105 368186 1 1686 28 68696
106 269661 1 746 29 71561
107 125390 0 3192 48 59831
108 510834 0 2045 16 97890
109 321896 1 5702 47 101481
110 249898 0 1932 20 72954
111 408881 1 936 22 67939
112 158492 0 3437 33 48022
113 292154 1 5131 44 86111
114 289513 0 2397 13 74020
115 378049 0 1389 6 57530
116 343466 1 1503 35 56364
117 332743 1 402 8 84990
118 442882 1 2239 17 88590
119 214215 0 2234 11 77200
120 315688 1 837 21 61262
121 375195 1 10579 92 110309
122 334280 1 875 12 67000
123 355864 1 1585 112 93099
124 480382 0 1659 25 107577
125 353058 0 2647 17 62920
126 217193 0 3294 23 75832
127 315380 1 0 0 60720
128 314533 1 94 10 60793
129 318056 1 422 23 57935
130 315380 1 0 0 60720
131 314353 1 34 7 60630
132 369448 1 1558 25 55637
133 315380 1 0 0 60720
134 312846 1 43 20 60887
135 312075 1 645 4 60720
136 315009 1 316 4 60505
137 318903 1 115 10 60945
138 314887 1 5 1 60720
139 314913 1 897 4 60720
140 315380 1 0 0 60720
141 325506 1 389 8 58990
142 315380 1 0 0 60720
143 298568 1 1002 11 56750
144 315834 1 36 4 60894
145 329784 1 460 15 63346
146 312878 1 309 9 56535
147 315380 1 0 0 60720
148 314987 1 9 7 60835
149 325249 1 271 2 60720
150 315877 1 14 0 61016
151 291650 1 520 7 58650
152 305959 1 1766 46 60438
153 315380 1 0 5 60720
154 297765 1 458 7 58625
155 315245 1 20 2 60938
156 315380 1 0 0 60720
157 315380 1 0 0 60720
158 315236 1 98 2 61490
159 336425 1 405 5 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 7 60830
165 302187 1 454 24 63261
166 314882 1 47 1 60720
167 315380 1 0 0 60720
168 382712 1 757 18 45689
169 341570 1 4655 55 60720
170 315380 1 0 0 60720
171 315380 1 0 0 60720
172 312412 1 36 3 61564
173 315380 1 0 0 60720
174 309596 1 203 9 61938
175 315380 1 0 0 60720
176 315547 1 126 8 60951
177 313267 1 400 113 60720
178 316176 0 71 0 60745
179 315380 1 0 0 60720
180 315380 1 0 0 60720
181 359335 1 972 19 71642
182 330068 0 531 11 71641
183 314289 1 2461 25 55792
184 297413 0 378 16 71873
185 314806 0 23 5 62555
186 333210 0 638 11 60370
187 352108 0 2300 23 64873
188 313332 1 149 6 62041
189 291787 1 226 5 65745
190 315380 1 0 0 60720
191 318745 1 275 7 59500
192 315380 1 0 0 60720
193 315366 1 141 7 61630
194 315380 1 0 0 60720
195 315688 1 28 3 60890
196 315380 0 0 0 60720
197 409642 0 4980 89 113521
198 315380 0 0 0 60720
199 315380 0 0 0 60720
200 269587 0 472 19 80045
201 315380 0 0 0 60720
202 315380 0 0 0 60720
203 315380 0 0 0 60720
204 300962 0 203 12 50804
205 325479 0 496 12 87390
206 316155 0 10 5 61656
207 318574 0 63 2 65688
208 315380 0 0 0 60720
209 343613 0 1136 26 48522
210 306948 0 265 3 60720
211 315380 1 0 0 60720
212 315380 1 0 0 60720
213 330059 1 267 11 57640
214 288985 1 474 10 61977
215 304485 1 534 5 62620
216 315380 1 0 2 60720
217 315688 1 15 6 60831
218 317736 1 397 7 60646
219 315380 1 0 2 60720
220 322331 1 1866 28 56225
221 296656 1 288 3 60510
222 315380 1 0 0 60720
223 315354 1 3 1 60698
224 312161 1 468 20 60720
225 315576 1 20 1 60805
226 314922 1 278 22 61404
227 314551 1 61 9 60720
228 315380 1 0 0 60720
229 312339 1 192 2 65276
230 315380 1 0 0 60720
231 298700 1 317 7 63915
232 321376 1 738 9 60720
233 315380 1 0 0 60720
234 303230 1 368 13 61686
235 315380 1 0 0 60720
236 315487 1 2 0 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 7 59818
245 324385 1 2332 2 72680
246 315380 0 0 0 60720
247 315380 0 0 0 60720
248 308989 1 131 15 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 1 64245
254 304376 1 622 38 73007
255 253588 1 2328 57 82732
256 315380 1 0 0 60720
257 309560 1 365 7 54820
258 298466 1 364 26 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 13 72835
264 331955 1 307 10 58856
265 315380 0 0 0 60720
266 315380 1 0 0 60720
267 315380 0 0 0 60720
268 381180 1 188 9 77655
269 315380 0 0 0 60720
270 331420 1 138 26 69817
271 315380 0 0 0 60720
272 315380 0 0 0 60720
273 315380 0 0 0 60720
274 310201 1 125 19 60798
275 315380 1 0 0 60720
276 320016 1 282 12 62452
277 320398 1 335 23 64175
278 315380 1 0 0 60720
279 291841 0 1324 29 67440
280 310670 1 176 8 68136
281 315380 0 0 0 60720
282 315380 1 0 0 60720
283 313491 1 249 26 56726
284 315380 1 0 0 60720
285 331323 1 333 9 70811
286 315380 1 0 0 60720
287 319210 0 601 5 60720
288 318098 1 30 3 62045
289 315380 0 0 0 60720
290 292754 1 249 13 54323
291 315380 0 0 0 60720
292 325176 1 165 12 62841
293 365959 1 453 19 81125
294 315380 1 0 0 60720
295 302409 1 53 10 59506
296 340968 0 382 9 59365
297 315380 1 0 0 60720
298 315380 0 0 0 60720
299 315380 1 0 0 60720
300 315380 1 0 9 60720
301 313164 0 30 4 60798
302 301164 1 290 1 58790
303 315380 1 0 1 60720
304 315380 0 0 0 60720
305 344425 1 366 14 61808
306 315394 1 2 12 60735
307 315380 0 0 0 60720
308 316647 0 209 19 64016
309 309836 1 384 17 54683
310 315380 0 0 0 60720
311 315380 0 0 0 60720
312 346611 1 365 32 87192
313 315380 0 0 0 60720
314 322031 0 49 14 64107
315 315656 1 3 8 60761
316 339445 1 133 4 65990
317 314964 1 32 0 59988
318 297141 1 368 20 61167
319 315372 1 1 5 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 1 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 4 60379
335 315388 1 1 1 60727
336 315371 1 314 4 60720
337 296139 1 844 20 60925
338 315380 1 0 0 60720
339 313880 1 26 1 60896
340 317698 1 125 10 59734
341 295580 1 304 12 62969
342 315380 1 0 0 60720
343 315380 1 0 0 60720
344 315380 1 0 0 60720
345 308256 1 621 13 60720
346 315380 1 0 0 60720
347 303677 1 119 3 59118
348 315380 1 0 0 60720
349 315380 1 0 0 60720
350 319369 1 1595 10 60720
351 318690 1 312 3 58598
352 314049 1 60 7 61124
353 325699 1 587 10 59595
354 314210 1 135 1 62065
355 315380 1 0 0 60720
356 315380 1 0 0 60720
357 322378 1 514 15 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 28 61600
365 316386 1 180 9 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 7 60720
371 315380 1 0 0 60720
372 323361 1 448 7 59781
373 336639 1 227 7 76644
374 307424 1 174 3 64820
375 315380 1 0 0 60720
376 315380 1 0 0 60720
377 295370 1 121 11 56178
378 322340 1 607 7 60436
379 319864 1 2212 10 60720
380 315380 1 0 0 60720
381 315380 1 0 0 60720
382 317291 1 530 18 73433
383 280398 1 571 14 41477
384 315380 1 0 0 60720
385 317330 1 78 12 62700
386 238125 1 2489 29 67804
387 327071 1 131 3 59661
388 309038 1 923 6 58620
389 314210 1 72 3 60398
390 307930 1 572 8 58580
391 322327 1 397 10 62710
392 292136 1 450 6 59325
393 263276 1 622 8 60950
394 367655 1 694 6 68060
395 283910 0 3425 9 83620
396 283587 1 562 8 58456
397 243650 1 4917 26 52811
398 438493 0 1442 239 121173
399 296261 1 529 7 63870
400 230621 0 2126 41 21001
401 304252 1 1061 3 70415
402 333505 1 776 8 64230
403 296919 1 611 6 59190
404 278990 0 1526 21 69351
405 276898 1 592 7 64270
406 327007 1 1182 11 70694
407 317046 1 621 11 68005
408 304555 1 989 12 58930
409 298096 1 438 9 58320
410 231861 1 726 3 69980
411 309422 1 1303 57 69863
412 286963 0 7419 21 63255
413 269753 0 1164 15 57320
414 448243 0 3310 32 75230
415 165404 1 1920 11 79420
416 204325 1 965 2 73490
417 407159 1 3256 23 35250
418 290476 0 1135 20 62285
419 275311 1 1270 24 69206
420 246541 1 661 1 65920
421 253468 1 1013 1 69770
422 240897 1 2844 74 72683
423 -83265 0 11528 68 -14545
424 -42143 1 6526 20 55830
425 272713 1 2264 20 55174
426 215362 0 5109 82 67038
427 42754 1 3999 21 51252
428 306275 1 35624 244 157278
429 253537 1 9252 32 79510
430 372631 1 15236 86 77440
431 -7170 1 18073 69 27284
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) group costs orders dividends
204815.678 -47556.317 22.044 950.981 1.518
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-2169550 -26507 35623 65929 2913139
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 2.048e+05 3.849e+04 5.322 1.66e-07 ***
group -4.756e+04 2.793e+04 -1.703 0.089345 .
costs 2.204e+01 2.150e+00 10.253 < 2e-16 ***
orders 9.510e+02 3.506e+02 2.713 0.006947 **
dividends 1.518e+00 4.116e-01 3.689 0.000254 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 253000 on 426 degrees of freedom
Multiple R-squared: 0.658, Adjusted R-squared: 0.6547
F-statistic: 204.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 8.929724e-44 4.464862e-44
[2,] 1.0000000 1.710130e-43 8.550649e-44
[3,] 1.0000000 1.472796e-56 7.363978e-57
[4,] 1.0000000 4.388800e-60 2.194400e-60
[5,] 1.0000000 6.105175e-74 3.052588e-74
[6,] 1.0000000 8.510729e-98 4.255365e-98
[7,] 1.0000000 1.060595e-99 5.302973e-100
[8,] 1.0000000 8.189103e-111 4.094551e-111
[9,] 1.0000000 1.455076e-120 7.275382e-121
[10,] 1.0000000 3.847186e-126 1.923593e-126
[11,] 1.0000000 1.133439e-134 5.667195e-135
[12,] 1.0000000 1.932293e-139 9.661465e-140
[13,] 1.0000000 1.562321e-142 7.811607e-143
[14,] 1.0000000 6.183486e-148 3.091743e-148
[15,] 1.0000000 1.613759e-149 8.068793e-150
[16,] 1.0000000 2.276656e-150 1.138328e-150
[17,] 1.0000000 2.100393e-152 1.050196e-152
[18,] 1.0000000 4.716624e-152 2.358312e-152
[19,] 1.0000000 3.782349e-157 1.891175e-157
[20,] 1.0000000 2.669433e-162 1.334716e-162
[21,] 1.0000000 7.852758e-163 3.926379e-163
[22,] 1.0000000 2.042638e-179 1.021319e-179
[23,] 1.0000000 9.704061e-180 4.852030e-180
[24,] 1.0000000 8.316716e-184 4.158358e-184
[25,] 1.0000000 2.684244e-185 1.342122e-185
[26,] 1.0000000 2.821387e-193 1.410693e-193
[27,] 1.0000000 3.137464e-195 1.568732e-195
[28,] 1.0000000 3.129579e-194 1.564790e-194
[29,] 1.0000000 6.222566e-194 3.111283e-194
[30,] 1.0000000 1.561918e-200 7.809590e-201
[31,] 1.0000000 4.094542e-206 2.047271e-206
[32,] 1.0000000 1.694215e-206 8.471075e-207
[33,] 1.0000000 9.272081e-207 4.636041e-207
[34,] 1.0000000 2.923527e-207 1.461763e-207
[35,] 1.0000000 1.286388e-206 6.431941e-207
[36,] 1.0000000 1.569793e-207 7.848964e-208
[37,] 1.0000000 1.744535e-206 8.722676e-207
[38,] 1.0000000 7.524827e-207 3.762413e-207
[39,] 1.0000000 5.955238e-206 2.977619e-206
[40,] 1.0000000 5.505826e-208 2.752913e-208
[41,] 1.0000000 1.962118e-208 9.810589e-209
[42,] 1.0000000 3.513849e-208 1.756924e-208
[43,] 1.0000000 1.164878e-208 5.824391e-209
[44,] 1.0000000 4.067565e-209 2.033783e-209
[45,] 1.0000000 3.537442e-209 1.768721e-209
[46,] 1.0000000 3.931178e-209 1.965589e-209
[47,] 1.0000000 1.181033e-208 5.905163e-209
[48,] 1.0000000 9.862900e-213 4.931450e-213
[49,] 1.0000000 8.382714e-213 4.191357e-213
[50,] 1.0000000 4.190671e-227 2.095336e-227
[51,] 1.0000000 8.613017e-228 4.306509e-228
[52,] 1.0000000 5.429237e-227 2.714619e-227
[53,] 1.0000000 8.915173e-227 4.457587e-227
[54,] 1.0000000 2.834720e-228 1.417360e-228
[55,] 1.0000000 6.806395e-230 3.403198e-230
[56,] 1.0000000 2.652220e-229 1.326110e-229
[57,] 1.0000000 1.001598e-228 5.007990e-229
[58,] 1.0000000 3.221348e-228 1.610674e-228
[59,] 1.0000000 6.432763e-228 3.216382e-228
[60,] 1.0000000 1.089958e-227 5.449789e-228
[61,] 1.0000000 1.117477e-226 5.587384e-227
[62,] 1.0000000 5.704019e-226 2.852010e-226
[63,] 1.0000000 2.524624e-225 1.262312e-225
[64,] 1.0000000 1.877316e-224 9.386578e-225
[65,] 1.0000000 8.390189e-226 4.195095e-226
[66,] 1.0000000 2.237458e-225 1.118729e-225
[67,] 1.0000000 1.134072e-224 5.670358e-225
[68,] 1.0000000 1.010447e-223 5.052234e-224
[69,] 1.0000000 2.371207e-223 1.185603e-223
[70,] 1.0000000 1.581880e-222 7.909402e-223
[71,] 1.0000000 1.595748e-222 7.978741e-223
[72,] 1.0000000 1.272961e-221 6.364805e-222
[73,] 1.0000000 4.949137e-221 2.474569e-221
[74,] 1.0000000 1.947523e-220 9.737613e-221
[75,] 1.0000000 3.071400e-220 1.535700e-220
[76,] 1.0000000 2.706247e-220 1.353123e-220
[77,] 1.0000000 1.000773e-219 5.003864e-220
[78,] 1.0000000 7.101104e-219 3.550552e-219
[79,] 1.0000000 1.642677e-218 8.213386e-219
[80,] 1.0000000 1.183264e-218 5.916319e-219
[81,] 1.0000000 6.844389e-218 3.422195e-218
[82,] 1.0000000 1.390227e-217 6.951136e-218
[83,] 1.0000000 1.030231e-216 5.151155e-217
[84,] 1.0000000 7.673293e-216 3.836647e-216
[85,] 1.0000000 5.126159e-215 2.563079e-215
[86,] 1.0000000 1.372404e-216 6.862019e-217
[87,] 1.0000000 1.071483e-215 5.357414e-216
[88,] 1.0000000 2.983363e-216 1.491681e-216
[89,] 1.0000000 3.073925e-216 1.536962e-216
[90,] 1.0000000 2.120654e-215 1.060327e-215
[91,] 1.0000000 2.348028e-214 1.174014e-214
[92,] 1.0000000 1.824227e-213 9.121133e-214
[93,] 1.0000000 2.064612e-212 1.032306e-212
[94,] 1.0000000 1.946156e-211 9.730779e-212
[95,] 1.0000000 1.420118e-212 7.100592e-213
[96,] 1.0000000 1.331733e-212 6.658663e-213
[97,] 1.0000000 1.913082e-213 9.565411e-214
[98,] 1.0000000 6.987772e-213 3.493886e-213
[99,] 1.0000000 2.149489e-212 1.074745e-212
[100,] 1.0000000 1.128984e-214 5.644918e-215
[101,] 1.0000000 5.620141e-216 2.810071e-216
[102,] 1.0000000 4.368869e-215 2.184434e-215
[103,] 1.0000000 9.236168e-215 4.618084e-215
[104,] 1.0000000 5.210613e-215 2.605307e-215
[105,] 1.0000000 2.884043e-215 1.442021e-215
[106,] 1.0000000 2.258819e-214 1.129409e-214
[107,] 1.0000000 1.881515e-213 9.407576e-214
[108,] 1.0000000 1.568024e-213 7.840122e-214
[109,] 1.0000000 5.293107e-213 2.646554e-213
[110,] 1.0000000 3.201947e-212 1.600974e-212
[111,] 1.0000000 2.139583e-212 1.069792e-212
[112,] 1.0000000 2.138771e-213 1.069386e-213
[113,] 1.0000000 1.872131e-212 9.360656e-213
[114,] 1.0000000 1.081818e-211 5.409088e-212
[115,] 1.0000000 8.438831e-211 4.219416e-211
[116,] 1.0000000 7.663780e-210 3.831890e-210
[117,] 1.0000000 1.241985e-209 6.209925e-210
[118,] 1.0000000 3.889246e-209 1.944623e-209
[119,] 1.0000000 1.130349e-209 5.651744e-210
[120,] 1.0000000 8.681176e-209 4.340588e-209
[121,] 1.0000000 7.245568e-208 3.622784e-208
[122,] 1.0000000 5.985974e-207 2.992987e-207
[123,] 1.0000000 4.768247e-206 2.384124e-206
[124,] 1.0000000 4.015630e-205 2.007815e-205
[125,] 1.0000000 2.980651e-205 1.490326e-205
[126,] 1.0000000 2.466344e-204 1.233172e-204
[127,] 1.0000000 2.287487e-203 1.143744e-203
[128,] 1.0000000 2.017381e-202 1.008690e-202
[129,] 1.0000000 1.742127e-201 8.710637e-202
[130,] 1.0000000 1.520844e-200 7.604218e-201
[131,] 1.0000000 1.309392e-199 6.546959e-200
[132,] 1.0000000 1.140130e-198 5.700652e-199
[133,] 1.0000000 9.838625e-198 4.919313e-198
[134,] 1.0000000 7.395459e-197 3.697730e-197
[135,] 1.0000000 6.439118e-196 3.219559e-196
[136,] 1.0000000 6.165250e-195 3.082625e-195
[137,] 1.0000000 5.479407e-194 2.739703e-194
[138,] 1.0000000 4.652079e-193 2.326040e-193
[139,] 1.0000000 3.960683e-192 1.980342e-192
[140,] 1.0000000 3.485863e-191 1.742932e-191
[141,] 1.0000000 3.160173e-190 1.580086e-190
[142,] 1.0000000 2.531805e-189 1.265903e-189
[143,] 1.0000000 2.241069e-188 1.120534e-188
[144,] 1.0000000 2.085650e-187 1.042825e-187
[145,] 1.0000000 2.103201e-186 1.051600e-186
[146,] 1.0000000 1.884981e-185 9.424907e-186
[147,] 1.0000000 1.768562e-184 8.842812e-185
[148,] 1.0000000 1.574150e-183 7.870749e-184
[149,] 1.0000000 1.387979e-182 6.939896e-183
[150,] 1.0000000 1.223813e-181 6.119066e-182
[151,] 1.0000000 1.095283e-180 5.476414e-181
[152,] 1.0000000 7.241611e-180 3.620805e-180
[153,] 1.0000000 6.380519e-179 3.190260e-179
[154,] 1.0000000 5.617652e-178 2.808826e-178
[155,] 1.0000000 4.941539e-177 2.470769e-177
[156,] 1.0000000 4.342227e-176 2.171114e-176
[157,] 1.0000000 4.015465e-175 2.007733e-175
[158,] 1.0000000 3.670896e-174 1.835448e-174
[159,] 1.0000000 3.222541e-173 1.611271e-173
[160,] 1.0000000 2.810352e-172 1.405176e-172
[161,] 1.0000000 3.212993e-173 1.606497e-173
[162,] 1.0000000 5.545004e-173 2.772502e-173
[163,] 1.0000000 4.904097e-172 2.452049e-172
[164,] 1.0000000 4.329083e-171 2.164542e-171
[165,] 1.0000000 3.912148e-170 1.956074e-170
[166,] 1.0000000 3.438805e-169 1.719403e-169
[167,] 1.0000000 3.151383e-168 1.575692e-168
[168,] 1.0000000 2.756594e-167 1.378297e-167
[169,] 1.0000000 2.443510e-166 1.221755e-166
[170,] 1.0000000 2.316665e-165 1.158332e-165
[171,] 1.0000000 2.180605e-164 1.090303e-164
[172,] 1.0000000 1.880971e-163 9.404853e-164
[173,] 1.0000000 1.618505e-162 8.092524e-163
[174,] 1.0000000 1.035890e-161 5.179450e-162
[175,] 1.0000000 9.663808e-161 4.831904e-161
[176,] 1.0000000 5.433790e-160 2.716895e-160
[177,] 1.0000000 3.059999e-159 1.530000e-159
[178,] 1.0000000 2.826486e-158 1.413243e-158
[179,] 1.0000000 2.071980e-157 1.035990e-157
[180,] 1.0000000 8.547270e-157 4.273635e-157
[181,] 1.0000000 7.404772e-156 3.702386e-156
[182,] 1.0000000 5.203063e-155 2.601532e-155
[183,] 1.0000000 4.353463e-154 2.176731e-154
[184,] 1.0000000 3.454640e-153 1.727320e-153
[185,] 1.0000000 2.873879e-152 1.436940e-152
[186,] 1.0000000 2.429975e-151 1.214988e-151
[187,] 1.0000000 2.009169e-150 1.004585e-150
[188,] 1.0000000 1.663834e-149 8.319169e-150
[189,] 1.0000000 1.463920e-148 7.319598e-149
[190,] 1.0000000 1.055848e-147 5.279239e-148
[191,] 1.0000000 9.200750e-147 4.600375e-147
[192,] 1.0000000 7.980481e-146 3.990240e-146
[193,] 1.0000000 6.884047e-146 3.442024e-146
[194,] 1.0000000 5.994865e-145 2.997433e-145
[195,] 1.0000000 5.196791e-144 2.598396e-144
[196,] 1.0000000 4.484394e-143 2.242197e-143
[197,] 1.0000000 3.623067e-142 1.811533e-142
[198,] 1.0000000 1.398935e-141 6.994675e-142
[199,] 1.0000000 1.199565e-140 5.997823e-141
[200,] 1.0000000 1.006251e-139 5.031253e-140
[201,] 1.0000000 8.529467e-139 4.264734e-139
[202,] 1.0000000 1.461810e-138 7.309050e-139
[203,] 1.0000000 1.238013e-137 6.190064e-138
[204,] 1.0000000 9.757329e-137 4.878665e-137
[205,] 1.0000000 7.665771e-136 3.832885e-136
[206,] 1.0000000 4.442382e-135 2.221191e-135
[207,] 1.0000000 3.322306e-134 1.661153e-134
[208,] 1.0000000 2.682326e-133 1.341163e-133
[209,] 1.0000000 2.086211e-132 1.043106e-132
[210,] 1.0000000 1.622627e-131 8.113137e-132
[211,] 1.0000000 1.219075e-130 6.095377e-131
[212,] 1.0000000 9.381112e-130 4.690556e-130
[213,] 1.0000000 4.307616e-129 2.153808e-129
[214,] 1.0000000 3.352272e-128 1.676136e-128
[215,] 1.0000000 2.551118e-127 1.275559e-127
[216,] 1.0000000 1.935829e-126 9.679143e-127
[217,] 1.0000000 1.488273e-125 7.441367e-126
[218,] 1.0000000 1.120540e-124 5.602700e-125
[219,] 1.0000000 8.556595e-124 4.278297e-124
[220,] 1.0000000 6.450997e-123 3.225498e-123
[221,] 1.0000000 4.798742e-122 2.399371e-122
[222,] 1.0000000 3.627399e-121 1.813699e-121
[223,] 1.0000000 2.678659e-120 1.339330e-120
[224,] 1.0000000 1.937733e-119 9.688663e-120
[225,] 1.0000000 1.301685e-118 6.508424e-119
[226,] 1.0000000 9.507079e-118 4.753540e-118
[227,] 1.0000000 7.176897e-117 3.588449e-117
[228,] 1.0000000 5.200686e-116 2.600343e-116
[229,] 1.0000000 3.752854e-115 1.876427e-115
[230,] 1.0000000 2.698230e-114 1.349115e-114
[231,] 1.0000000 1.922043e-113 9.610217e-114
[232,] 1.0000000 1.370863e-112 6.854316e-113
[233,] 1.0000000 9.738268e-112 4.869134e-112
[234,] 1.0000000 6.889927e-111 3.444964e-111
[235,] 1.0000000 4.938277e-110 2.469138e-110
[236,] 1.0000000 3.465376e-109 1.732688e-109
[237,] 1.0000000 2.393388e-108 1.196694e-108
[238,] 1.0000000 1.723177e-107 8.615886e-108
[239,] 1.0000000 1.247045e-106 6.235225e-107
[240,] 1.0000000 8.981264e-106 4.490632e-106
[241,] 1.0000000 6.373885e-105 3.186943e-105
[242,] 1.0000000 4.546689e-104 2.273345e-104
[243,] 1.0000000 3.227538e-103 1.613769e-103
[244,] 1.0000000 2.168018e-102 1.084009e-102
[245,] 1.0000000 1.524638e-101 7.623188e-102
[246,] 1.0000000 1.038565e-100 5.192823e-101
[247,] 1.0000000 5.954458e-100 2.977229e-100
[248,] 1.0000000 5.159507e-100 2.579754e-100
[249,] 1.0000000 3.470232e-99 1.735116e-99
[250,] 1.0000000 2.118732e-98 1.059366e-98
[251,] 1.0000000 1.164181e-97 5.820905e-98
[252,] 1.0000000 8.038866e-97 4.019433e-97
[253,] 1.0000000 5.315619e-96 2.657809e-96
[254,] 1.0000000 3.635138e-95 1.817569e-95
[255,] 1.0000000 2.473386e-94 1.236693e-94
[256,] 1.0000000 1.636538e-93 8.182691e-94
[257,] 1.0000000 8.323230e-93 4.161615e-93
[258,] 1.0000000 5.590268e-92 2.795134e-92
[259,] 1.0000000 3.599310e-91 1.799655e-91
[260,] 1.0000000 2.393728e-90 1.196864e-90
[261,] 1.0000000 1.244427e-89 6.222137e-90
[262,] 1.0000000 8.203833e-89 4.101917e-89
[263,] 1.0000000 5.354740e-88 2.677370e-88
[264,] 1.0000000 3.494243e-87 1.747121e-87
[265,] 1.0000000 2.268154e-86 1.134077e-86
[266,] 1.0000000 1.464461e-85 7.322305e-86
[267,] 1.0000000 9.324322e-85 4.662161e-85
[268,] 1.0000000 5.758235e-84 2.879117e-84
[269,] 1.0000000 3.526792e-83 1.763396e-83
[270,] 1.0000000 2.197237e-82 1.098618e-82
[271,] 1.0000000 1.337620e-81 6.688100e-82
[272,] 1.0000000 7.175639e-81 3.587819e-81
[273,] 1.0000000 4.282514e-80 2.141257e-80
[274,] 1.0000000 2.664198e-79 1.332099e-79
[275,] 1.0000000 1.593667e-78 7.968334e-79
[276,] 1.0000000 8.914938e-78 4.457469e-78
[277,] 1.0000000 5.281743e-77 2.640872e-77
[278,] 1.0000000 3.196466e-76 1.598233e-76
[279,] 1.0000000 1.875505e-75 9.377524e-76
[280,] 1.0000000 1.096519e-74 5.482597e-75
[281,] 1.0000000 6.395781e-74 3.197890e-74
[282,] 1.0000000 3.819060e-73 1.909530e-73
[283,] 1.0000000 2.241909e-72 1.120954e-72
[284,] 1.0000000 1.324254e-71 6.621270e-72
[285,] 1.0000000 7.362502e-71 3.681251e-71
[286,] 1.0000000 4.208816e-70 2.104408e-70
[287,] 1.0000000 2.371638e-69 1.185819e-69
[288,] 1.0000000 1.367691e-68 6.838456e-69
[289,] 1.0000000 5.971010e-68 2.985505e-68
[290,] 1.0000000 3.319440e-67 1.659720e-67
[291,] 1.0000000 1.885851e-66 9.429255e-67
[292,] 1.0000000 1.037447e-65 5.187236e-66
[293,] 1.0000000 5.698661e-65 2.849331e-65
[294,] 1.0000000 3.199901e-64 1.599951e-64
[295,] 1.0000000 1.773839e-63 8.869194e-64
[296,] 1.0000000 9.556411e-63 4.778206e-63
[297,] 1.0000000 5.247472e-62 2.623736e-62
[298,] 1.0000000 2.133624e-61 1.066812e-61
[299,] 1.0000000 1.136810e-60 5.684050e-61
[300,] 1.0000000 6.149827e-60 3.074914e-60
[301,] 1.0000000 3.339058e-59 1.669529e-59
[302,] 1.0000000 1.596429e-58 7.982144e-59
[303,] 1.0000000 8.488333e-58 4.244166e-58
[304,] 1.0000000 4.479978e-57 2.239989e-57
[305,] 1.0000000 2.100607e-56 1.050303e-56
[306,] 1.0000000 1.094907e-55 5.474537e-56
[307,] 1.0000000 5.705074e-55 2.852537e-55
[308,] 1.0000000 2.890970e-54 1.445485e-54
[309,] 1.0000000 1.365145e-53 6.825723e-54
[310,] 1.0000000 6.768636e-53 3.384318e-53
[311,] 1.0000000 3.465018e-52 1.732509e-52
[312,] 1.0000000 1.715208e-51 8.576040e-52
[313,] 1.0000000 8.423252e-51 4.211626e-51
[314,] 1.0000000 4.111810e-50 2.055905e-50
[315,] 1.0000000 1.995060e-49 9.975302e-50
[316,] 1.0000000 9.621166e-49 4.810583e-49
[317,] 1.0000000 4.681033e-48 2.340517e-48
[318,] 1.0000000 2.230061e-47 1.115031e-47
[319,] 1.0000000 1.066973e-46 5.334867e-47
[320,] 1.0000000 5.075600e-46 2.537800e-46
[321,] 1.0000000 2.385956e-45 1.192978e-45
[322,] 1.0000000 1.108746e-44 5.543729e-45
[323,] 1.0000000 5.119272e-44 2.559636e-44
[324,] 1.0000000 2.348372e-43 1.174186e-43
[325,] 1.0000000 1.070245e-42 5.351224e-43
[326,] 1.0000000 4.845424e-42 2.422712e-42
[327,] 1.0000000 2.183570e-41 1.091785e-41
[328,] 1.0000000 9.759100e-41 4.879550e-41
[329,] 1.0000000 4.293055e-40 2.146527e-40
[330,] 1.0000000 1.960476e-39 9.802379e-40
[331,] 1.0000000 8.581514e-39 4.290757e-39
[332,] 1.0000000 3.757470e-38 1.878735e-38
[333,] 1.0000000 1.584330e-37 7.921651e-38
[334,] 1.0000000 6.971298e-37 3.485649e-37
[335,] 1.0000000 2.969279e-36 1.484640e-36
[336,] 1.0000000 1.255578e-35 6.277890e-36
[337,] 1.0000000 5.270589e-35 2.635294e-35
[338,] 1.0000000 2.225750e-34 1.112875e-34
[339,] 1.0000000 9.205314e-34 4.602657e-34
[340,] 1.0000000 3.864565e-33 1.932283e-33
[341,] 1.0000000 1.574401e-32 7.872007e-33
[342,] 1.0000000 6.364734e-32 3.182367e-32
[343,] 1.0000000 2.313304e-31 1.156652e-31
[344,] 1.0000000 8.664499e-31 4.332249e-31
[345,] 1.0000000 3.458423e-30 1.729211e-30
[346,] 1.0000000 1.218564e-29 6.092819e-30
[347,] 1.0000000 4.808086e-29 2.404043e-29
[348,] 1.0000000 1.853821e-28 9.269104e-29
[349,] 1.0000000 7.087858e-28 3.543929e-28
[350,] 1.0000000 2.703424e-27 1.351712e-27
[351,] 1.0000000 1.018846e-26 5.094228e-27
[352,] 1.0000000 3.806428e-26 1.903214e-26
[353,] 1.0000000 1.409578e-25 7.047889e-26
[354,] 1.0000000 5.181652e-25 2.590826e-25
[355,] 1.0000000 1.884555e-24 9.422777e-25
[356,] 1.0000000 6.791133e-24 3.395566e-24
[357,] 1.0000000 2.417417e-23 1.208709e-23
[358,] 1.0000000 8.387977e-23 4.193988e-23
[359,] 1.0000000 2.942338e-22 1.471169e-22
[360,] 1.0000000 1.022085e-21 5.110426e-22
[361,] 1.0000000 3.515363e-21 1.757682e-21
[362,] 1.0000000 1.196928e-20 5.984640e-21
[363,] 1.0000000 4.013921e-20 2.006960e-20
[364,] 1.0000000 1.338684e-19 6.693420e-20
[365,] 1.0000000 4.103891e-19 2.051945e-19
[366,] 1.0000000 1.417499e-18 7.087495e-19
[367,] 1.0000000 4.830696e-18 2.415348e-18
[368,] 1.0000000 1.545519e-17 7.727594e-18
[369,] 1.0000000 4.888215e-17 2.444108e-17
[370,] 1.0000000 1.575527e-16 7.877633e-17
[371,] 1.0000000 4.574695e-16 2.287347e-16
[372,] 1.0000000 1.225756e-15 6.128782e-16
[373,] 1.0000000 3.696777e-15 1.848389e-15
[374,] 1.0000000 1.100606e-14 5.503029e-15
[375,] 1.0000000 3.499313e-14 1.749657e-14
[376,] 1.0000000 8.946392e-14 4.473196e-14
[377,] 1.0000000 2.565646e-13 1.282823e-13
[378,] 1.0000000 7.391409e-13 3.695705e-13
[379,] 1.0000000 1.996102e-12 9.980509e-13
[380,] 1.0000000 5.092436e-12 2.546218e-12
[381,] 1.0000000 1.350930e-11 6.754652e-12
[382,] 1.0000000 3.638123e-11 1.819062e-11
[383,] 1.0000000 9.534922e-11 4.767461e-11
[384,] 1.0000000 2.439596e-10 1.219798e-10
[385,] 1.0000000 6.685122e-10 3.342561e-10
[386,] 1.0000000 1.890270e-09 9.451348e-10
[387,] 1.0000000 3.651164e-09 1.825582e-09
[388,] 1.0000000 8.801424e-09 4.400712e-09
[389,] 1.0000000 2.319614e-08 1.159807e-08
[390,] 1.0000000 5.791891e-08 2.895945e-08
[391,] 0.9999999 1.200941e-07 6.004706e-08
[392,] 0.9999998 3.061282e-07 1.530641e-07
[393,] 0.9999997 6.854829e-07 3.427415e-07
[394,] 0.9999991 1.714987e-06 8.574935e-07
[395,] 0.9999982 3.530216e-06 1.765108e-06
[396,] 0.9999960 7.901662e-06 3.950831e-06
[397,] 0.9999908 1.848380e-05 9.241898e-06
[398,] 0.9999780 4.393486e-05 2.196743e-05
[399,] 0.9999532 9.353515e-05 4.676758e-05
[400,] 0.9999007 1.985446e-04 9.927231e-05
[401,] 0.9998100 3.799579e-04 1.899790e-04
[402,] 0.9996361 7.278703e-04 3.639352e-04
[403,] 0.9992154 1.569141e-03 7.845706e-04
[404,] 0.9983970 3.206047e-03 1.603024e-03
[405,] 0.9968513 6.297453e-03 3.148727e-03
[406,] 0.9936497 1.270069e-02 6.350347e-03
[407,] 0.9937102 1.257957e-02 6.289785e-03
[408,] 0.9899248 2.015036e-02 1.007518e-02
[409,] 0.9814504 3.709920e-02 1.854960e-02
[410,] 0.9948193 1.036136e-02 5.180679e-03
[411,] 0.9885426 2.291473e-02 1.145737e-02
[412,] 0.9748603 5.027941e-02 2.513971e-02
[413,] 0.9465972 1.068056e-01 5.340278e-02
[414,] 0.8937335 2.125330e-01 1.062665e-01
[415,] 0.7986671 4.026659e-01 2.013329e-01
[416,] 0.6517514 6.964971e-01 3.482486e-01
> postscript(file="/var/www/html/rcomp/tmp/11m531291219010.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/rcomp/tmp/21m531291219010.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/rcomp/tmp/31m531291219010.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/rcomp/tmp/4tv4o1291219010.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/rcomp/tmp/5tv4o1291219010.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
1402961.479 2913138.769 1352753.355 -2169549.620 -410679.272 86519.825
7 8 9 10 11 12
-107543.117 115170.675 -360898.474 684728.075 -125278.630 359203.174
13 14 15 16 17 18
432172.621 -1120524.455 -233672.311 479383.913 -209978.041 -257330.290
19 20 21 22 23 24
-321220.597 -555562.003 171657.086 -169963.511 -97939.098 -608028.961
25 26 27 28 29 30
-299233.585 -70546.516 -494980.590 -327174.386 359384.866 -243129.762
31 32 33 34 35 36
252635.029 -118693.439 222270.011 98356.917 -31446.011 -544502.394
37 38 39 40 41 42
50507.376 128163.400 -76450.561 -285308.135 -433861.615 -313566.191
43 44 45 46 47 48
17670.150 -190813.141 -13896.029 -87156.361 -308099.912 -288734.002
49 50 51 52 53 54
-36458.009 -168107.866 6487.705 -111871.450 -253114.256 15295.099
55 56 57 58 59 60
-419913.622 29789.549 330892.738 -7327.732 -43800.516 -236281.701
61 62 63 64 65 66
-18233.594 -166473.087 -26687.129 -205557.027 -2007.474 -468248.697
67 68 69 70 71 72
-52372.158 -47243.728 -82202.671 -101959.036 -157977.521 13012.828
73 74 75 76 77 78
-73604.619 27328.585 -19117.089 -133538.174 -13321.181 -61192.121
79 80 81 82 83 84
-28017.638 -70984.910 39537.333 -81245.903 -66102.720 52250.348
85 86 87 88 89 90
38693.345 -135178.356 -11343.467 14859.058 -99555.689 46230.652
91 92 93 94 95 96
12662.905 -55225.344 84324.520 -204255.765 61879.176 -15521.425
97 98 99 100 101 102
85736.827 -32922.067 -25880.808 -37141.056 -70129.747 -99146.107
103 104 105 106 107 108
-149224.227 -36215.540 42831.684 -40273.045 -286278.469 97095.955
109 110 111 112 113 114
-159832.658 -127292.519 106914.799 -226382.812 -150798.315 -92889.550
115 116 117 118 119 120
49560.374 34212.597 29973.288 85593.041 -167520.659 26992.853
121 122 123 124 125 126
-270239.165 44593.971 -84197.513 51886.210 -21806.226 -197243.845
127 128 129 130 131 132
65929.058 53389.293 41658.487 65929.058 57632.350 69595.804
133 134 135 136 137 138
65929.058 43174.000 44601.864 55114.720 57065.590 64374.858
139 140 141 142 143 144
41884.819 65929.058 62498.832 65929.058 22596.029 61521.371
145 146 147 148 149 150
51941.111 54410.801 65929.058 58506.192 67922.218 65668.025
151 152 153 154 155 156
27222.295 -25738.300 61174.154 34741.970 63120.229 65929.058
157 158 159 160 161 162
65929.058 60553.705 73101.614 65929.058 65929.058 65929.058
163 164 165 166 167 168
65929.058 39346.008 16046.603 63444.017 65929.058 122277.892
169 170 171 172 173 174
-62798.915 65929.058 65929.058 58033.087 65929.058 45262.036
175 176 177 178 179 180
65929.058 55359.960 -52462.303 17565.671 65929.058 65929.058
181 182 183 184 185 186
53805.874 -5686.749 -5704.113 -40076.194 9750.736 12209.395
187 188 189 190 191 192
-23778.155 52884.959 24969.757 65929.058 58427.473 65929.058
193 194 195 196 197 198
54768.353 65929.058 62508.776 18372.740 -161948.935 18372.740
199 200 201 202 203 204
18372.740 -85234.859 18372.740 18372.740 18372.740 3123.603
205 206 207 208 209 210
-34367.008 12751.263 10733.070 18372.740 15358.754 1246.183
211 212 213 214 215 216
65929.058 65929.058 68937.951 17666.963 35622.966 64027.096
217 218 219 220 221 222
60031.984 52989.146 64027.096 11943.597 38322.337 65929.058
223 224 225 226 227 228
64919.348 33373.929 64604.144 37382.774 55196.557 65929.058
229 230 231 232 233 234
49836.275 65929.058 30753.308 47097.884 65929.058 31837.494
235 236 237 238 239 240
65929.058 65957.049 65929.058 60031.142 65929.058 65929.058
241 242 243 244 245 246
65929.058 57539.279 65929.058 60369.653 3466.936 18372.740
247 248 249 250 251 252
18372.740 40733.687 18372.740 18372.740 45705.516 18372.740
253 254 255 256 257 258
47886.726 -13578.909 -134807.851 65929.058 54364.204 36026.365
259 260 261 262 263 264
18372.740 65929.058 18372.740 18372.740 58739.116 69056.918
265 266 267 268 269 270
18372.740 65929.058 18372.740 93313.467 18372.740 40389.473
271 272 273 274 275 276
18372.740 18372.740 18372.740 39807.516 65929.058 50307.221
277 278 279 280 281 282
36444.067 65929.058 -72133.754 38471.735 18372.740 65929.058
283 284 285 286 287 288
39889.761 65929.058 50651.402 65929.058 4199.495 63121.044
289 290 291 292 293 294
18372.740 35164.002 18372.740 57455.728 57472.520 65929.058
295 296 297 298 299 300
44123.151 29038.476 65929.058 18372.740 65929.058 57370.231
301 302 303 304 305 306
11573.074 47299.698 64978.077 18372.740 71940.367 54464.426
307 308 309 310 311 312
18372.740 -8040.394 44919.571 18372.740 18372.740 18490.060
313 314 315 316 317 318
18372.740 5487.357 58468.829 75256.829 65919.056 19879.629
319 320 321 322 323 324
61145.628 65929.058 65929.058 65929.058 65929.058 18372.740
325 326 327 328 329 330
65929.058 61612.076 18372.740 18372.740 65929.058 65929.058
331 332 333 334 335 336
65929.058 65929.058 65929.058 58875.669 64953.405 55194.372
337 338 339 340 341 342
8752.196 65929.058 62637.715 57478.821 24601.292 65929.058
343 344 345 346 347 348
65929.058 65929.058 32753.089 65929.058 51182.227 65929.058
349 350 351 352 353 354
65929.058 25248.341 62730.287 56005.167 55506.617 58790.037
355 356 357 358 359 360
65929.058 65929.058 19911.198 65929.058 65929.058 65929.058
361 362 363 364 365 366
62118.054 65929.058 65929.058 -55498.104 56055.704 65929.058
367 368 369 370 371 372
65929.058 65929.058 65929.058 54639.637 65929.058 58803.246
373 374 375 376 377 378
51349.727 45059.431 65929.058 65929.058 39687.115 53282.786
379 380 381 382 383 384
12142.297 65929.058 65929.058 19736.941 34263.075 65929.058
385 386 387 388 389 390
51741.622 -104527.163 73487.260 36723.162 60807.854 41511.317
391 392 393 394 395 396
51593.419 29177.487 -7843.262 86055.384 -131925.428 17577.025
397 398 399 400 401 402
-126907.668 -209372.064 23709.340 -91936.031 13839.628 54010.943
403 404 405 406 407 408
30616.402 -84731.246 2350.256 25895.871 32384.186 24608.709
409 410 411 412 413 414
34074.969 -50506.225 -36839.836 -197406.940 -62015.747 25808.647
415 416 417 418 419 420
-165224.220 -87688.976 102732.060 -52946.773 -37843.496 -26327.090
421 422 423 424 425 426
-33004.999 -159782.676 -584784.875 -447047.074 -37244.261 -281840.274
427 428 429 430 431
-300445.488 -1107109.319 -258823.827 -319850.175 -669870.658
> postscript(file="/var/www/html/rcomp/tmp/64m3r1291219010.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 1402961.479 NA
1 2913138.769 1402961.479
2 1352753.355 2913138.769
3 -2169549.620 1352753.355
4 -410679.272 -2169549.620
5 86519.825 -410679.272
6 -107543.117 86519.825
7 115170.675 -107543.117
8 -360898.474 115170.675
9 684728.075 -360898.474
10 -125278.630 684728.075
11 359203.174 -125278.630
12 432172.621 359203.174
13 -1120524.455 432172.621
14 -233672.311 -1120524.455
15 479383.913 -233672.311
16 -209978.041 479383.913
17 -257330.290 -209978.041
18 -321220.597 -257330.290
19 -555562.003 -321220.597
20 171657.086 -555562.003
21 -169963.511 171657.086
22 -97939.098 -169963.511
23 -608028.961 -97939.098
24 -299233.585 -608028.961
25 -70546.516 -299233.585
26 -494980.590 -70546.516
27 -327174.386 -494980.590
28 359384.866 -327174.386
29 -243129.762 359384.866
30 252635.029 -243129.762
31 -118693.439 252635.029
32 222270.011 -118693.439
33 98356.917 222270.011
34 -31446.011 98356.917
35 -544502.394 -31446.011
36 50507.376 -544502.394
37 128163.400 50507.376
38 -76450.561 128163.400
39 -285308.135 -76450.561
40 -433861.615 -285308.135
41 -313566.191 -433861.615
42 17670.150 -313566.191
43 -190813.141 17670.150
44 -13896.029 -190813.141
45 -87156.361 -13896.029
46 -308099.912 -87156.361
47 -288734.002 -308099.912
48 -36458.009 -288734.002
49 -168107.866 -36458.009
50 6487.705 -168107.866
51 -111871.450 6487.705
52 -253114.256 -111871.450
53 15295.099 -253114.256
54 -419913.622 15295.099
55 29789.549 -419913.622
56 330892.738 29789.549
57 -7327.732 330892.738
58 -43800.516 -7327.732
59 -236281.701 -43800.516
60 -18233.594 -236281.701
61 -166473.087 -18233.594
62 -26687.129 -166473.087
63 -205557.027 -26687.129
64 -2007.474 -205557.027
65 -468248.697 -2007.474
66 -52372.158 -468248.697
67 -47243.728 -52372.158
68 -82202.671 -47243.728
69 -101959.036 -82202.671
70 -157977.521 -101959.036
71 13012.828 -157977.521
72 -73604.619 13012.828
73 27328.585 -73604.619
74 -19117.089 27328.585
75 -133538.174 -19117.089
76 -13321.181 -133538.174
77 -61192.121 -13321.181
78 -28017.638 -61192.121
79 -70984.910 -28017.638
80 39537.333 -70984.910
81 -81245.903 39537.333
82 -66102.720 -81245.903
83 52250.348 -66102.720
84 38693.345 52250.348
85 -135178.356 38693.345
86 -11343.467 -135178.356
87 14859.058 -11343.467
88 -99555.689 14859.058
89 46230.652 -99555.689
90 12662.905 46230.652
91 -55225.344 12662.905
92 84324.520 -55225.344
93 -204255.765 84324.520
94 61879.176 -204255.765
95 -15521.425 61879.176
96 85736.827 -15521.425
97 -32922.067 85736.827
98 -25880.808 -32922.067
99 -37141.056 -25880.808
100 -70129.747 -37141.056
101 -99146.107 -70129.747
102 -149224.227 -99146.107
103 -36215.540 -149224.227
104 42831.684 -36215.540
105 -40273.045 42831.684
106 -286278.469 -40273.045
107 97095.955 -286278.469
108 -159832.658 97095.955
109 -127292.519 -159832.658
110 106914.799 -127292.519
111 -226382.812 106914.799
112 -150798.315 -226382.812
113 -92889.550 -150798.315
114 49560.374 -92889.550
115 34212.597 49560.374
116 29973.288 34212.597
117 85593.041 29973.288
118 -167520.659 85593.041
119 26992.853 -167520.659
120 -270239.165 26992.853
121 44593.971 -270239.165
122 -84197.513 44593.971
123 51886.210 -84197.513
124 -21806.226 51886.210
125 -197243.845 -21806.226
126 65929.058 -197243.845
127 53389.293 65929.058
128 41658.487 53389.293
129 65929.058 41658.487
130 57632.350 65929.058
131 69595.804 57632.350
132 65929.058 69595.804
133 43174.000 65929.058
134 44601.864 43174.000
135 55114.720 44601.864
136 57065.590 55114.720
137 64374.858 57065.590
138 41884.819 64374.858
139 65929.058 41884.819
140 62498.832 65929.058
141 65929.058 62498.832
142 22596.029 65929.058
143 61521.371 22596.029
144 51941.111 61521.371
145 54410.801 51941.111
146 65929.058 54410.801
147 58506.192 65929.058
148 67922.218 58506.192
149 65668.025 67922.218
150 27222.295 65668.025
151 -25738.300 27222.295
152 61174.154 -25738.300
153 34741.970 61174.154
154 63120.229 34741.970
155 65929.058 63120.229
156 65929.058 65929.058
157 60553.705 65929.058
158 73101.614 60553.705
159 65929.058 73101.614
160 65929.058 65929.058
161 65929.058 65929.058
162 65929.058 65929.058
163 39346.008 65929.058
164 16046.603 39346.008
165 63444.017 16046.603
166 65929.058 63444.017
167 122277.892 65929.058
168 -62798.915 122277.892
169 65929.058 -62798.915
170 65929.058 65929.058
171 58033.087 65929.058
172 65929.058 58033.087
173 45262.036 65929.058
174 65929.058 45262.036
175 55359.960 65929.058
176 -52462.303 55359.960
177 17565.671 -52462.303
178 65929.058 17565.671
179 65929.058 65929.058
180 53805.874 65929.058
181 -5686.749 53805.874
182 -5704.113 -5686.749
183 -40076.194 -5704.113
184 9750.736 -40076.194
185 12209.395 9750.736
186 -23778.155 12209.395
187 52884.959 -23778.155
188 24969.757 52884.959
189 65929.058 24969.757
190 58427.473 65929.058
191 65929.058 58427.473
192 54768.353 65929.058
193 65929.058 54768.353
194 62508.776 65929.058
195 18372.740 62508.776
196 -161948.935 18372.740
197 18372.740 -161948.935
198 18372.740 18372.740
199 -85234.859 18372.740
200 18372.740 -85234.859
201 18372.740 18372.740
202 18372.740 18372.740
203 3123.603 18372.740
204 -34367.008 3123.603
205 12751.263 -34367.008
206 10733.070 12751.263
207 18372.740 10733.070
208 15358.754 18372.740
209 1246.183 15358.754
210 65929.058 1246.183
211 65929.058 65929.058
212 68937.951 65929.058
213 17666.963 68937.951
214 35622.966 17666.963
215 64027.096 35622.966
216 60031.984 64027.096
217 52989.146 60031.984
218 64027.096 52989.146
219 11943.597 64027.096
220 38322.337 11943.597
221 65929.058 38322.337
222 64919.348 65929.058
223 33373.929 64919.348
224 64604.144 33373.929
225 37382.774 64604.144
226 55196.557 37382.774
227 65929.058 55196.557
228 49836.275 65929.058
229 65929.058 49836.275
230 30753.308 65929.058
231 47097.884 30753.308
232 65929.058 47097.884
233 31837.494 65929.058
234 65929.058 31837.494
235 65957.049 65929.058
236 65929.058 65957.049
237 60031.142 65929.058
238 65929.058 60031.142
239 65929.058 65929.058
240 65929.058 65929.058
241 57539.279 65929.058
242 65929.058 57539.279
243 60369.653 65929.058
244 3466.936 60369.653
245 18372.740 3466.936
246 18372.740 18372.740
247 40733.687 18372.740
248 18372.740 40733.687
249 18372.740 18372.740
250 45705.516 18372.740
251 18372.740 45705.516
252 47886.726 18372.740
253 -13578.909 47886.726
254 -134807.851 -13578.909
255 65929.058 -134807.851
256 54364.204 65929.058
257 36026.365 54364.204
258 18372.740 36026.365
259 65929.058 18372.740
260 18372.740 65929.058
261 18372.740 18372.740
262 58739.116 18372.740
263 69056.918 58739.116
264 18372.740 69056.918
265 65929.058 18372.740
266 18372.740 65929.058
267 93313.467 18372.740
268 18372.740 93313.467
269 40389.473 18372.740
270 18372.740 40389.473
271 18372.740 18372.740
272 18372.740 18372.740
273 39807.516 18372.740
274 65929.058 39807.516
275 50307.221 65929.058
276 36444.067 50307.221
277 65929.058 36444.067
278 -72133.754 65929.058
279 38471.735 -72133.754
280 18372.740 38471.735
281 65929.058 18372.740
282 39889.761 65929.058
283 65929.058 39889.761
284 50651.402 65929.058
285 65929.058 50651.402
286 4199.495 65929.058
287 63121.044 4199.495
288 18372.740 63121.044
289 35164.002 18372.740
290 18372.740 35164.002
291 57455.728 18372.740
292 57472.520 57455.728
293 65929.058 57472.520
294 44123.151 65929.058
295 29038.476 44123.151
296 65929.058 29038.476
297 18372.740 65929.058
298 65929.058 18372.740
299 57370.231 65929.058
300 11573.074 57370.231
301 47299.698 11573.074
302 64978.077 47299.698
303 18372.740 64978.077
304 71940.367 18372.740
305 54464.426 71940.367
306 18372.740 54464.426
307 -8040.394 18372.740
308 44919.571 -8040.394
309 18372.740 44919.571
310 18372.740 18372.740
311 18490.060 18372.740
312 18372.740 18490.060
313 5487.357 18372.740
314 58468.829 5487.357
315 75256.829 58468.829
316 65919.056 75256.829
317 19879.629 65919.056
318 61145.628 19879.629
319 65929.058 61145.628
320 65929.058 65929.058
321 65929.058 65929.058
322 65929.058 65929.058
323 18372.740 65929.058
324 65929.058 18372.740
325 61612.076 65929.058
326 18372.740 61612.076
327 18372.740 18372.740
328 65929.058 18372.740
329 65929.058 65929.058
330 65929.058 65929.058
331 65929.058 65929.058
332 65929.058 65929.058
333 58875.669 65929.058
334 64953.405 58875.669
335 55194.372 64953.405
336 8752.196 55194.372
337 65929.058 8752.196
338 62637.715 65929.058
339 57478.821 62637.715
340 24601.292 57478.821
341 65929.058 24601.292
342 65929.058 65929.058
343 65929.058 65929.058
344 32753.089 65929.058
345 65929.058 32753.089
346 51182.227 65929.058
347 65929.058 51182.227
348 65929.058 65929.058
349 25248.341 65929.058
350 62730.287 25248.341
351 56005.167 62730.287
352 55506.617 56005.167
353 58790.037 55506.617
354 65929.058 58790.037
355 65929.058 65929.058
356 19911.198 65929.058
357 65929.058 19911.198
358 65929.058 65929.058
359 65929.058 65929.058
360 62118.054 65929.058
361 65929.058 62118.054
362 65929.058 65929.058
363 -55498.104 65929.058
364 56055.704 -55498.104
365 65929.058 56055.704
366 65929.058 65929.058
367 65929.058 65929.058
368 65929.058 65929.058
369 54639.637 65929.058
370 65929.058 54639.637
371 58803.246 65929.058
372 51349.727 58803.246
373 45059.431 51349.727
374 65929.058 45059.431
375 65929.058 65929.058
376 39687.115 65929.058
377 53282.786 39687.115
378 12142.297 53282.786
379 65929.058 12142.297
380 65929.058 65929.058
381 19736.941 65929.058
382 34263.075 19736.941
383 65929.058 34263.075
384 51741.622 65929.058
385 -104527.163 51741.622
386 73487.260 -104527.163
387 36723.162 73487.260
388 60807.854 36723.162
389 41511.317 60807.854
390 51593.419 41511.317
391 29177.487 51593.419
392 -7843.262 29177.487
393 86055.384 -7843.262
394 -131925.428 86055.384
395 17577.025 -131925.428
396 -126907.668 17577.025
397 -209372.064 -126907.668
398 23709.340 -209372.064
399 -91936.031 23709.340
400 13839.628 -91936.031
401 54010.943 13839.628
402 30616.402 54010.943
403 -84731.246 30616.402
404 2350.256 -84731.246
405 25895.871 2350.256
406 32384.186 25895.871
407 24608.709 32384.186
408 34074.969 24608.709
409 -50506.225 34074.969
410 -36839.836 -50506.225
411 -197406.940 -36839.836
412 -62015.747 -197406.940
413 25808.647 -62015.747
414 -165224.220 25808.647
415 -87688.976 -165224.220
416 102732.060 -87688.976
417 -52946.773 102732.060
418 -37843.496 -52946.773
419 -26327.090 -37843.496
420 -33004.999 -26327.090
421 -159782.676 -33004.999
422 -584784.875 -159782.676
423 -447047.074 -584784.875
424 -37244.261 -447047.074
425 -281840.274 -37244.261
426 -300445.488 -281840.274
427 -1107109.319 -300445.488
428 -258823.827 -1107109.319
429 -319850.175 -258823.827
430 -669870.658 -319850.175
431 NA -669870.658
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] 2913138.769 1402961.479
[2,] 1352753.355 2913138.769
[3,] -2169549.620 1352753.355
[4,] -410679.272 -2169549.620
[5,] 86519.825 -410679.272
[6,] -107543.117 86519.825
[7,] 115170.675 -107543.117
[8,] -360898.474 115170.675
[9,] 684728.075 -360898.474
[10,] -125278.630 684728.075
[11,] 359203.174 -125278.630
[12,] 432172.621 359203.174
[13,] -1120524.455 432172.621
[14,] -233672.311 -1120524.455
[15,] 479383.913 -233672.311
[16,] -209978.041 479383.913
[17,] -257330.290 -209978.041
[18,] -321220.597 -257330.290
[19,] -555562.003 -321220.597
[20,] 171657.086 -555562.003
[21,] -169963.511 171657.086
[22,] -97939.098 -169963.511
[23,] -608028.961 -97939.098
[24,] -299233.585 -608028.961
[25,] -70546.516 -299233.585
[26,] -494980.590 -70546.516
[27,] -327174.386 -494980.590
[28,] 359384.866 -327174.386
[29,] -243129.762 359384.866
[30,] 252635.029 -243129.762
[31,] -118693.439 252635.029
[32,] 222270.011 -118693.439
[33,] 98356.917 222270.011
[34,] -31446.011 98356.917
[35,] -544502.394 -31446.011
[36,] 50507.376 -544502.394
[37,] 128163.400 50507.376
[38,] -76450.561 128163.400
[39,] -285308.135 -76450.561
[40,] -433861.615 -285308.135
[41,] -313566.191 -433861.615
[42,] 17670.150 -313566.191
[43,] -190813.141 17670.150
[44,] -13896.029 -190813.141
[45,] -87156.361 -13896.029
[46,] -308099.912 -87156.361
[47,] -288734.002 -308099.912
[48,] -36458.009 -288734.002
[49,] -168107.866 -36458.009
[50,] 6487.705 -168107.866
[51,] -111871.450 6487.705
[52,] -253114.256 -111871.450
[53,] 15295.099 -253114.256
[54,] -419913.622 15295.099
[55,] 29789.549 -419913.622
[56,] 330892.738 29789.549
[57,] -7327.732 330892.738
[58,] -43800.516 -7327.732
[59,] -236281.701 -43800.516
[60,] -18233.594 -236281.701
[61,] -166473.087 -18233.594
[62,] -26687.129 -166473.087
[63,] -205557.027 -26687.129
[64,] -2007.474 -205557.027
[65,] -468248.697 -2007.474
[66,] -52372.158 -468248.697
[67,] -47243.728 -52372.158
[68,] -82202.671 -47243.728
[69,] -101959.036 -82202.671
[70,] -157977.521 -101959.036
[71,] 13012.828 -157977.521
[72,] -73604.619 13012.828
[73,] 27328.585 -73604.619
[74,] -19117.089 27328.585
[75,] -133538.174 -19117.089
[76,] -13321.181 -133538.174
[77,] -61192.121 -13321.181
[78,] -28017.638 -61192.121
[79,] -70984.910 -28017.638
[80,] 39537.333 -70984.910
[81,] -81245.903 39537.333
[82,] -66102.720 -81245.903
[83,] 52250.348 -66102.720
[84,] 38693.345 52250.348
[85,] -135178.356 38693.345
[86,] -11343.467 -135178.356
[87,] 14859.058 -11343.467
[88,] -99555.689 14859.058
[89,] 46230.652 -99555.689
[90,] 12662.905 46230.652
[91,] -55225.344 12662.905
[92,] 84324.520 -55225.344
[93,] -204255.765 84324.520
[94,] 61879.176 -204255.765
[95,] -15521.425 61879.176
[96,] 85736.827 -15521.425
[97,] -32922.067 85736.827
[98,] -25880.808 -32922.067
[99,] -37141.056 -25880.808
[100,] -70129.747 -37141.056
[101,] -99146.107 -70129.747
[102,] -149224.227 -99146.107
[103,] -36215.540 -149224.227
[104,] 42831.684 -36215.540
[105,] -40273.045 42831.684
[106,] -286278.469 -40273.045
[107,] 97095.955 -286278.469
[108,] -159832.658 97095.955
[109,] -127292.519 -159832.658
[110,] 106914.799 -127292.519
[111,] -226382.812 106914.799
[112,] -150798.315 -226382.812
[113,] -92889.550 -150798.315
[114,] 49560.374 -92889.550
[115,] 34212.597 49560.374
[116,] 29973.288 34212.597
[117,] 85593.041 29973.288
[118,] -167520.659 85593.041
[119,] 26992.853 -167520.659
[120,] -270239.165 26992.853
[121,] 44593.971 -270239.165
[122,] -84197.513 44593.971
[123,] 51886.210 -84197.513
[124,] -21806.226 51886.210
[125,] -197243.845 -21806.226
[126,] 65929.058 -197243.845
[127,] 53389.293 65929.058
[128,] 41658.487 53389.293
[129,] 65929.058 41658.487
[130,] 57632.350 65929.058
[131,] 69595.804 57632.350
[132,] 65929.058 69595.804
[133,] 43174.000 65929.058
[134,] 44601.864 43174.000
[135,] 55114.720 44601.864
[136,] 57065.590 55114.720
[137,] 64374.858 57065.590
[138,] 41884.819 64374.858
[139,] 65929.058 41884.819
[140,] 62498.832 65929.058
[141,] 65929.058 62498.832
[142,] 22596.029 65929.058
[143,] 61521.371 22596.029
[144,] 51941.111 61521.371
[145,] 54410.801 51941.111
[146,] 65929.058 54410.801
[147,] 58506.192 65929.058
[148,] 67922.218 58506.192
[149,] 65668.025 67922.218
[150,] 27222.295 65668.025
[151,] -25738.300 27222.295
[152,] 61174.154 -25738.300
[153,] 34741.970 61174.154
[154,] 63120.229 34741.970
[155,] 65929.058 63120.229
[156,] 65929.058 65929.058
[157,] 60553.705 65929.058
[158,] 73101.614 60553.705
[159,] 65929.058 73101.614
[160,] 65929.058 65929.058
[161,] 65929.058 65929.058
[162,] 65929.058 65929.058
[163,] 39346.008 65929.058
[164,] 16046.603 39346.008
[165,] 63444.017 16046.603
[166,] 65929.058 63444.017
[167,] 122277.892 65929.058
[168,] -62798.915 122277.892
[169,] 65929.058 -62798.915
[170,] 65929.058 65929.058
[171,] 58033.087 65929.058
[172,] 65929.058 58033.087
[173,] 45262.036 65929.058
[174,] 65929.058 45262.036
[175,] 55359.960 65929.058
[176,] -52462.303 55359.960
[177,] 17565.671 -52462.303
[178,] 65929.058 17565.671
[179,] 65929.058 65929.058
[180,] 53805.874 65929.058
[181,] -5686.749 53805.874
[182,] -5704.113 -5686.749
[183,] -40076.194 -5704.113
[184,] 9750.736 -40076.194
[185,] 12209.395 9750.736
[186,] -23778.155 12209.395
[187,] 52884.959 -23778.155
[188,] 24969.757 52884.959
[189,] 65929.058 24969.757
[190,] 58427.473 65929.058
[191,] 65929.058 58427.473
[192,] 54768.353 65929.058
[193,] 65929.058 54768.353
[194,] 62508.776 65929.058
[195,] 18372.740 62508.776
[196,] -161948.935 18372.740
[197,] 18372.740 -161948.935
[198,] 18372.740 18372.740
[199,] -85234.859 18372.740
[200,] 18372.740 -85234.859
[201,] 18372.740 18372.740
[202,] 18372.740 18372.740
[203,] 3123.603 18372.740
[204,] -34367.008 3123.603
[205,] 12751.263 -34367.008
[206,] 10733.070 12751.263
[207,] 18372.740 10733.070
[208,] 15358.754 18372.740
[209,] 1246.183 15358.754
[210,] 65929.058 1246.183
[211,] 65929.058 65929.058
[212,] 68937.951 65929.058
[213,] 17666.963 68937.951
[214,] 35622.966 17666.963
[215,] 64027.096 35622.966
[216,] 60031.984 64027.096
[217,] 52989.146 60031.984
[218,] 64027.096 52989.146
[219,] 11943.597 64027.096
[220,] 38322.337 11943.597
[221,] 65929.058 38322.337
[222,] 64919.348 65929.058
[223,] 33373.929 64919.348
[224,] 64604.144 33373.929
[225,] 37382.774 64604.144
[226,] 55196.557 37382.774
[227,] 65929.058 55196.557
[228,] 49836.275 65929.058
[229,] 65929.058 49836.275
[230,] 30753.308 65929.058
[231,] 47097.884 30753.308
[232,] 65929.058 47097.884
[233,] 31837.494 65929.058
[234,] 65929.058 31837.494
[235,] 65957.049 65929.058
[236,] 65929.058 65957.049
[237,] 60031.142 65929.058
[238,] 65929.058 60031.142
[239,] 65929.058 65929.058
[240,] 65929.058 65929.058
[241,] 57539.279 65929.058
[242,] 65929.058 57539.279
[243,] 60369.653 65929.058
[244,] 3466.936 60369.653
[245,] 18372.740 3466.936
[246,] 18372.740 18372.740
[247,] 40733.687 18372.740
[248,] 18372.740 40733.687
[249,] 18372.740 18372.740
[250,] 45705.516 18372.740
[251,] 18372.740 45705.516
[252,] 47886.726 18372.740
[253,] -13578.909 47886.726
[254,] -134807.851 -13578.909
[255,] 65929.058 -134807.851
[256,] 54364.204 65929.058
[257,] 36026.365 54364.204
[258,] 18372.740 36026.365
[259,] 65929.058 18372.740
[260,] 18372.740 65929.058
[261,] 18372.740 18372.740
[262,] 58739.116 18372.740
[263,] 69056.918 58739.116
[264,] 18372.740 69056.918
[265,] 65929.058 18372.740
[266,] 18372.740 65929.058
[267,] 93313.467 18372.740
[268,] 18372.740 93313.467
[269,] 40389.473 18372.740
[270,] 18372.740 40389.473
[271,] 18372.740 18372.740
[272,] 18372.740 18372.740
[273,] 39807.516 18372.740
[274,] 65929.058 39807.516
[275,] 50307.221 65929.058
[276,] 36444.067 50307.221
[277,] 65929.058 36444.067
[278,] -72133.754 65929.058
[279,] 38471.735 -72133.754
[280,] 18372.740 38471.735
[281,] 65929.058 18372.740
[282,] 39889.761 65929.058
[283,] 65929.058 39889.761
[284,] 50651.402 65929.058
[285,] 65929.058 50651.402
[286,] 4199.495 65929.058
[287,] 63121.044 4199.495
[288,] 18372.740 63121.044
[289,] 35164.002 18372.740
[290,] 18372.740 35164.002
[291,] 57455.728 18372.740
[292,] 57472.520 57455.728
[293,] 65929.058 57472.520
[294,] 44123.151 65929.058
[295,] 29038.476 44123.151
[296,] 65929.058 29038.476
[297,] 18372.740 65929.058
[298,] 65929.058 18372.740
[299,] 57370.231 65929.058
[300,] 11573.074 57370.231
[301,] 47299.698 11573.074
[302,] 64978.077 47299.698
[303,] 18372.740 64978.077
[304,] 71940.367 18372.740
[305,] 54464.426 71940.367
[306,] 18372.740 54464.426
[307,] -8040.394 18372.740
[308,] 44919.571 -8040.394
[309,] 18372.740 44919.571
[310,] 18372.740 18372.740
[311,] 18490.060 18372.740
[312,] 18372.740 18490.060
[313,] 5487.357 18372.740
[314,] 58468.829 5487.357
[315,] 75256.829 58468.829
[316,] 65919.056 75256.829
[317,] 19879.629 65919.056
[318,] 61145.628 19879.629
[319,] 65929.058 61145.628
[320,] 65929.058 65929.058
[321,] 65929.058 65929.058
[322,] 65929.058 65929.058
[323,] 18372.740 65929.058
[324,] 65929.058 18372.740
[325,] 61612.076 65929.058
[326,] 18372.740 61612.076
[327,] 18372.740 18372.740
[328,] 65929.058 18372.740
[329,] 65929.058 65929.058
[330,] 65929.058 65929.058
[331,] 65929.058 65929.058
[332,] 65929.058 65929.058
[333,] 58875.669 65929.058
[334,] 64953.405 58875.669
[335,] 55194.372 64953.405
[336,] 8752.196 55194.372
[337,] 65929.058 8752.196
[338,] 62637.715 65929.058
[339,] 57478.821 62637.715
[340,] 24601.292 57478.821
[341,] 65929.058 24601.292
[342,] 65929.058 65929.058
[343,] 65929.058 65929.058
[344,] 32753.089 65929.058
[345,] 65929.058 32753.089
[346,] 51182.227 65929.058
[347,] 65929.058 51182.227
[348,] 65929.058 65929.058
[349,] 25248.341 65929.058
[350,] 62730.287 25248.341
[351,] 56005.167 62730.287
[352,] 55506.617 56005.167
[353,] 58790.037 55506.617
[354,] 65929.058 58790.037
[355,] 65929.058 65929.058
[356,] 19911.198 65929.058
[357,] 65929.058 19911.198
[358,] 65929.058 65929.058
[359,] 65929.058 65929.058
[360,] 62118.054 65929.058
[361,] 65929.058 62118.054
[362,] 65929.058 65929.058
[363,] -55498.104 65929.058
[364,] 56055.704 -55498.104
[365,] 65929.058 56055.704
[366,] 65929.058 65929.058
[367,] 65929.058 65929.058
[368,] 65929.058 65929.058
[369,] 54639.637 65929.058
[370,] 65929.058 54639.637
[371,] 58803.246 65929.058
[372,] 51349.727 58803.246
[373,] 45059.431 51349.727
[374,] 65929.058 45059.431
[375,] 65929.058 65929.058
[376,] 39687.115 65929.058
[377,] 53282.786 39687.115
[378,] 12142.297 53282.786
[379,] 65929.058 12142.297
[380,] 65929.058 65929.058
[381,] 19736.941 65929.058
[382,] 34263.075 19736.941
[383,] 65929.058 34263.075
[384,] 51741.622 65929.058
[385,] -104527.163 51741.622
[386,] 73487.260 -104527.163
[387,] 36723.162 73487.260
[388,] 60807.854 36723.162
[389,] 41511.317 60807.854
[390,] 51593.419 41511.317
[391,] 29177.487 51593.419
[392,] -7843.262 29177.487
[393,] 86055.384 -7843.262
[394,] -131925.428 86055.384
[395,] 17577.025 -131925.428
[396,] -126907.668 17577.025
[397,] -209372.064 -126907.668
[398,] 23709.340 -209372.064
[399,] -91936.031 23709.340
[400,] 13839.628 -91936.031
[401,] 54010.943 13839.628
[402,] 30616.402 54010.943
[403,] -84731.246 30616.402
[404,] 2350.256 -84731.246
[405,] 25895.871 2350.256
[406,] 32384.186 25895.871
[407,] 24608.709 32384.186
[408,] 34074.969 24608.709
[409,] -50506.225 34074.969
[410,] -36839.836 -50506.225
[411,] -197406.940 -36839.836
[412,] -62015.747 -197406.940
[413,] 25808.647 -62015.747
[414,] -165224.220 25808.647
[415,] -87688.976 -165224.220
[416,] 102732.060 -87688.976
[417,] -52946.773 102732.060
[418,] -37843.496 -52946.773
[419,] -26327.090 -37843.496
[420,] -33004.999 -26327.090
[421,] -159782.676 -33004.999
[422,] -584784.875 -159782.676
[423,] -447047.074 -584784.875
[424,] -37244.261 -447047.074
[425,] -281840.274 -37244.261
[426,] -300445.488 -281840.274
[427,] -1107109.319 -300445.488
[428,] -258823.827 -1107109.319
[429,] -319850.175 -258823.827
[430,] -669870.658 -319850.175
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 2913138.769 1402961.479
2 1352753.355 2913138.769
3 -2169549.620 1352753.355
4 -410679.272 -2169549.620
5 86519.825 -410679.272
6 -107543.117 86519.825
7 115170.675 -107543.117
8 -360898.474 115170.675
9 684728.075 -360898.474
10 -125278.630 684728.075
11 359203.174 -125278.630
12 432172.621 359203.174
13 -1120524.455 432172.621
14 -233672.311 -1120524.455
15 479383.913 -233672.311
16 -209978.041 479383.913
17 -257330.290 -209978.041
18 -321220.597 -257330.290
19 -555562.003 -321220.597
20 171657.086 -555562.003
21 -169963.511 171657.086
22 -97939.098 -169963.511
23 -608028.961 -97939.098
24 -299233.585 -608028.961
25 -70546.516 -299233.585
26 -494980.590 -70546.516
27 -327174.386 -494980.590
28 359384.866 -327174.386
29 -243129.762 359384.866
30 252635.029 -243129.762
31 -118693.439 252635.029
32 222270.011 -118693.439
33 98356.917 222270.011
34 -31446.011 98356.917
35 -544502.394 -31446.011
36 50507.376 -544502.394
37 128163.400 50507.376
38 -76450.561 128163.400
39 -285308.135 -76450.561
40 -433861.615 -285308.135
41 -313566.191 -433861.615
42 17670.150 -313566.191
43 -190813.141 17670.150
44 -13896.029 -190813.141
45 -87156.361 -13896.029
46 -308099.912 -87156.361
47 -288734.002 -308099.912
48 -36458.009 -288734.002
49 -168107.866 -36458.009
50 6487.705 -168107.866
51 -111871.450 6487.705
52 -253114.256 -111871.450
53 15295.099 -253114.256
54 -419913.622 15295.099
55 29789.549 -419913.622
56 330892.738 29789.549
57 -7327.732 330892.738
58 -43800.516 -7327.732
59 -236281.701 -43800.516
60 -18233.594 -236281.701
61 -166473.087 -18233.594
62 -26687.129 -166473.087
63 -205557.027 -26687.129
64 -2007.474 -205557.027
65 -468248.697 -2007.474
66 -52372.158 -468248.697
67 -47243.728 -52372.158
68 -82202.671 -47243.728
69 -101959.036 -82202.671
70 -157977.521 -101959.036
71 13012.828 -157977.521
72 -73604.619 13012.828
73 27328.585 -73604.619
74 -19117.089 27328.585
75 -133538.174 -19117.089
76 -13321.181 -133538.174
77 -61192.121 -13321.181
78 -28017.638 -61192.121
79 -70984.910 -28017.638
80 39537.333 -70984.910
81 -81245.903 39537.333
82 -66102.720 -81245.903
83 52250.348 -66102.720
84 38693.345 52250.348
85 -135178.356 38693.345
86 -11343.467 -135178.356
87 14859.058 -11343.467
88 -99555.689 14859.058
89 46230.652 -99555.689
90 12662.905 46230.652
91 -55225.344 12662.905
92 84324.520 -55225.344
93 -204255.765 84324.520
94 61879.176 -204255.765
95 -15521.425 61879.176
96 85736.827 -15521.425
97 -32922.067 85736.827
98 -25880.808 -32922.067
99 -37141.056 -25880.808
100 -70129.747 -37141.056
101 -99146.107 -70129.747
102 -149224.227 -99146.107
103 -36215.540 -149224.227
104 42831.684 -36215.540
105 -40273.045 42831.684
106 -286278.469 -40273.045
107 97095.955 -286278.469
108 -159832.658 97095.955
109 -127292.519 -159832.658
110 106914.799 -127292.519
111 -226382.812 106914.799
112 -150798.315 -226382.812
113 -92889.550 -150798.315
114 49560.374 -92889.550
115 34212.597 49560.374
116 29973.288 34212.597
117 85593.041 29973.288
118 -167520.659 85593.041
119 26992.853 -167520.659
120 -270239.165 26992.853
121 44593.971 -270239.165
122 -84197.513 44593.971
123 51886.210 -84197.513
124 -21806.226 51886.210
125 -197243.845 -21806.226
126 65929.058 -197243.845
127 53389.293 65929.058
128 41658.487 53389.293
129 65929.058 41658.487
130 57632.350 65929.058
131 69595.804 57632.350
132 65929.058 69595.804
133 43174.000 65929.058
134 44601.864 43174.000
135 55114.720 44601.864
136 57065.590 55114.720
137 64374.858 57065.590
138 41884.819 64374.858
139 65929.058 41884.819
140 62498.832 65929.058
141 65929.058 62498.832
142 22596.029 65929.058
143 61521.371 22596.029
144 51941.111 61521.371
145 54410.801 51941.111
146 65929.058 54410.801
147 58506.192 65929.058
148 67922.218 58506.192
149 65668.025 67922.218
150 27222.295 65668.025
151 -25738.300 27222.295
152 61174.154 -25738.300
153 34741.970 61174.154
154 63120.229 34741.970
155 65929.058 63120.229
156 65929.058 65929.058
157 60553.705 65929.058
158 73101.614 60553.705
159 65929.058 73101.614
160 65929.058 65929.058
161 65929.058 65929.058
162 65929.058 65929.058
163 39346.008 65929.058
164 16046.603 39346.008
165 63444.017 16046.603
166 65929.058 63444.017
167 122277.892 65929.058
168 -62798.915 122277.892
169 65929.058 -62798.915
170 65929.058 65929.058
171 58033.087 65929.058
172 65929.058 58033.087
173 45262.036 65929.058
174 65929.058 45262.036
175 55359.960 65929.058
176 -52462.303 55359.960
177 17565.671 -52462.303
178 65929.058 17565.671
179 65929.058 65929.058
180 53805.874 65929.058
181 -5686.749 53805.874
182 -5704.113 -5686.749
183 -40076.194 -5704.113
184 9750.736 -40076.194
185 12209.395 9750.736
186 -23778.155 12209.395
187 52884.959 -23778.155
188 24969.757 52884.959
189 65929.058 24969.757
190 58427.473 65929.058
191 65929.058 58427.473
192 54768.353 65929.058
193 65929.058 54768.353
194 62508.776 65929.058
195 18372.740 62508.776
196 -161948.935 18372.740
197 18372.740 -161948.935
198 18372.740 18372.740
199 -85234.859 18372.740
200 18372.740 -85234.859
201 18372.740 18372.740
202 18372.740 18372.740
203 3123.603 18372.740
204 -34367.008 3123.603
205 12751.263 -34367.008
206 10733.070 12751.263
207 18372.740 10733.070
208 15358.754 18372.740
209 1246.183 15358.754
210 65929.058 1246.183
211 65929.058 65929.058
212 68937.951 65929.058
213 17666.963 68937.951
214 35622.966 17666.963
215 64027.096 35622.966
216 60031.984 64027.096
217 52989.146 60031.984
218 64027.096 52989.146
219 11943.597 64027.096
220 38322.337 11943.597
221 65929.058 38322.337
222 64919.348 65929.058
223 33373.929 64919.348
224 64604.144 33373.929
225 37382.774 64604.144
226 55196.557 37382.774
227 65929.058 55196.557
228 49836.275 65929.058
229 65929.058 49836.275
230 30753.308 65929.058
231 47097.884 30753.308
232 65929.058 47097.884
233 31837.494 65929.058
234 65929.058 31837.494
235 65957.049 65929.058
236 65929.058 65957.049
237 60031.142 65929.058
238 65929.058 60031.142
239 65929.058 65929.058
240 65929.058 65929.058
241 57539.279 65929.058
242 65929.058 57539.279
243 60369.653 65929.058
244 3466.936 60369.653
245 18372.740 3466.936
246 18372.740 18372.740
247 40733.687 18372.740
248 18372.740 40733.687
249 18372.740 18372.740
250 45705.516 18372.740
251 18372.740 45705.516
252 47886.726 18372.740
253 -13578.909 47886.726
254 -134807.851 -13578.909
255 65929.058 -134807.851
256 54364.204 65929.058
257 36026.365 54364.204
258 18372.740 36026.365
259 65929.058 18372.740
260 18372.740 65929.058
261 18372.740 18372.740
262 58739.116 18372.740
263 69056.918 58739.116
264 18372.740 69056.918
265 65929.058 18372.740
266 18372.740 65929.058
267 93313.467 18372.740
268 18372.740 93313.467
269 40389.473 18372.740
270 18372.740 40389.473
271 18372.740 18372.740
272 18372.740 18372.740
273 39807.516 18372.740
274 65929.058 39807.516
275 50307.221 65929.058
276 36444.067 50307.221
277 65929.058 36444.067
278 -72133.754 65929.058
279 38471.735 -72133.754
280 18372.740 38471.735
281 65929.058 18372.740
282 39889.761 65929.058
283 65929.058 39889.761
284 50651.402 65929.058
285 65929.058 50651.402
286 4199.495 65929.058
287 63121.044 4199.495
288 18372.740 63121.044
289 35164.002 18372.740
290 18372.740 35164.002
291 57455.728 18372.740
292 57472.520 57455.728
293 65929.058 57472.520
294 44123.151 65929.058
295 29038.476 44123.151
296 65929.058 29038.476
297 18372.740 65929.058
298 65929.058 18372.740
299 57370.231 65929.058
300 11573.074 57370.231
301 47299.698 11573.074
302 64978.077 47299.698
303 18372.740 64978.077
304 71940.367 18372.740
305 54464.426 71940.367
306 18372.740 54464.426
307 -8040.394 18372.740
308 44919.571 -8040.394
309 18372.740 44919.571
310 18372.740 18372.740
311 18490.060 18372.740
312 18372.740 18490.060
313 5487.357 18372.740
314 58468.829 5487.357
315 75256.829 58468.829
316 65919.056 75256.829
317 19879.629 65919.056
318 61145.628 19879.629
319 65929.058 61145.628
320 65929.058 65929.058
321 65929.058 65929.058
322 65929.058 65929.058
323 18372.740 65929.058
324 65929.058 18372.740
325 61612.076 65929.058
326 18372.740 61612.076
327 18372.740 18372.740
328 65929.058 18372.740
329 65929.058 65929.058
330 65929.058 65929.058
331 65929.058 65929.058
332 65929.058 65929.058
333 58875.669 65929.058
334 64953.405 58875.669
335 55194.372 64953.405
336 8752.196 55194.372
337 65929.058 8752.196
338 62637.715 65929.058
339 57478.821 62637.715
340 24601.292 57478.821
341 65929.058 24601.292
342 65929.058 65929.058
343 65929.058 65929.058
344 32753.089 65929.058
345 65929.058 32753.089
346 51182.227 65929.058
347 65929.058 51182.227
348 65929.058 65929.058
349 25248.341 65929.058
350 62730.287 25248.341
351 56005.167 62730.287
352 55506.617 56005.167
353 58790.037 55506.617
354 65929.058 58790.037
355 65929.058 65929.058
356 19911.198 65929.058
357 65929.058 19911.198
358 65929.058 65929.058
359 65929.058 65929.058
360 62118.054 65929.058
361 65929.058 62118.054
362 65929.058 65929.058
363 -55498.104 65929.058
364 56055.704 -55498.104
365 65929.058 56055.704
366 65929.058 65929.058
367 65929.058 65929.058
368 65929.058 65929.058
369 54639.637 65929.058
370 65929.058 54639.637
371 58803.246 65929.058
372 51349.727 58803.246
373 45059.431 51349.727
374 65929.058 45059.431
375 65929.058 65929.058
376 39687.115 65929.058
377 53282.786 39687.115
378 12142.297 53282.786
379 65929.058 12142.297
380 65929.058 65929.058
381 19736.941 65929.058
382 34263.075 19736.941
383 65929.058 34263.075
384 51741.622 65929.058
385 -104527.163 51741.622
386 73487.260 -104527.163
387 36723.162 73487.260
388 60807.854 36723.162
389 41511.317 60807.854
390 51593.419 41511.317
391 29177.487 51593.419
392 -7843.262 29177.487
393 86055.384 -7843.262
394 -131925.428 86055.384
395 17577.025 -131925.428
396 -126907.668 17577.025
397 -209372.064 -126907.668
398 23709.340 -209372.064
399 -91936.031 23709.340
400 13839.628 -91936.031
401 54010.943 13839.628
402 30616.402 54010.943
403 -84731.246 30616.402
404 2350.256 -84731.246
405 25895.871 2350.256
406 32384.186 25895.871
407 24608.709 32384.186
408 34074.969 24608.709
409 -50506.225 34074.969
410 -36839.836 -50506.225
411 -197406.940 -36839.836
412 -62015.747 -197406.940
413 25808.647 -62015.747
414 -165224.220 25808.647
415 -87688.976 -165224.220
416 102732.060 -87688.976
417 -52946.773 102732.060
418 -37843.496 -52946.773
419 -26327.090 -37843.496
420 -33004.999 -26327.090
421 -159782.676 -33004.999
422 -584784.875 -159782.676
423 -447047.074 -584784.875
424 -37244.261 -447047.074
425 -281840.274 -37244.261
426 -300445.488 -281840.274
427 -1107109.319 -300445.488
428 -258823.827 -1107109.319
429 -319850.175 -258823.827
430 -669870.658 -319850.175
> 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/rcomp/tmp/74m3r1291219010.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/rcomp/tmp/8fw2c1291219010.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/rcomp/tmp/9fw2c1291219010.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/rcomp/tmp/1085kx1291219010.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/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/www/html/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/rcomp/tmp/11t50l1291219010.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/rcomp/tmp/12w6h81291219010.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/rcomp/tmp/13byeh1291219010.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/rcomp/tmp/14egd51291219010.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/rcomp/tmp/15hhbb1291219010.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/rcomp/tmp/163hsh1291219010.tab")
+ }
>
> try(system("convert tmp/11m531291219010.ps tmp/11m531291219010.png",intern=TRUE))
character(0)
> try(system("convert tmp/21m531291219010.ps tmp/21m531291219010.png",intern=TRUE))
character(0)
> try(system("convert tmp/31m531291219010.ps tmp/31m531291219010.png",intern=TRUE))
character(0)
> try(system("convert tmp/4tv4o1291219010.ps tmp/4tv4o1291219010.png",intern=TRUE))
character(0)
> try(system("convert tmp/5tv4o1291219010.ps tmp/5tv4o1291219010.png",intern=TRUE))
character(0)
> try(system("convert tmp/64m3r1291219010.ps tmp/64m3r1291219010.png",intern=TRUE))
character(0)
> try(system("convert tmp/74m3r1291219010.ps tmp/74m3r1291219010.png",intern=TRUE))
character(0)
> try(system("convert tmp/8fw2c1291219010.ps tmp/8fw2c1291219010.png",intern=TRUE))
character(0)
> try(system("convert tmp/9fw2c1291219010.ps tmp/9fw2c1291219010.png",intern=TRUE))
character(0)
> try(system("convert tmp/1085kx1291219010.ps tmp/1085kx1291219010.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
12.354 2.383 27.160