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(6282154
+ ,807
+ ,1
+ ,213118
+ ,162556
+ ,4321023
+ ,444
+ ,1
+ ,81767
+ ,29790
+ ,4111912
+ ,412
+ ,1
+ ,153198
+ ,87550
+ ,223193
+ ,428
+ ,0
+ ,-26007
+ ,84738
+ ,1491348
+ ,315
+ ,1
+ ,126942
+ ,54660
+ ,1629616
+ ,168
+ ,1
+ ,157214
+ ,42634
+ ,1398893
+ ,263
+ ,0
+ ,129352
+ ,40949
+ ,1926517
+ ,267
+ ,1
+ ,234817
+ ,45187
+ ,983660
+ ,228
+ ,1
+ ,60448
+ ,37704
+ ,1443586
+ ,129
+ ,1
+ ,47818
+ ,16275
+ ,1073089
+ ,104
+ ,0
+ ,245546
+ ,25830
+ ,984885
+ ,122
+ ,0
+ ,48020
+ ,12679
+ ,1405225
+ ,393
+ ,1
+ ,-1710
+ ,18014
+ ,227132
+ ,190
+ ,0
+ ,32648
+ ,43556
+ ,929118
+ ,280
+ ,1
+ ,95350
+ ,24811
+ ,1071292
+ ,63
+ ,0
+ ,151352
+ ,6575
+ ,638830
+ ,102
+ ,0
+ ,288170
+ ,7123
+ ,856956
+ ,265
+ ,1
+ ,114337
+ ,21950
+ ,992426
+ ,234
+ ,1
+ ,37884
+ ,37597
+ ,444477
+ ,277
+ ,0
+ ,122844
+ ,17821
+ ,857217
+ ,73
+ ,1
+ ,82340
+ ,12988
+ ,711969
+ ,67
+ ,1
+ ,79801
+ ,22330
+ ,702380
+ ,103
+ ,0
+ ,165548
+ ,13326
+ ,358589
+ ,290
+ ,0
+ ,116384
+ ,16189
+ ,297978
+ ,83
+ ,0
+ ,134028
+ ,7146
+ ,585715
+ ,56
+ ,0
+ ,63838
+ ,15824
+ ,657954
+ ,236
+ ,1
+ ,74996
+ ,27664
+ ,209458
+ ,73
+ ,0
+ ,31080
+ ,11920
+ ,786690
+ ,34
+ ,0
+ ,32168
+ ,8568
+ ,439798
+ ,139
+ ,0
+ ,49857
+ ,14416
+ ,688779
+ ,26
+ ,1
+ ,87161
+ ,3369
+ ,574339
+ ,70
+ ,1
+ ,106113
+ ,11819
+ ,741409
+ ,40
+ ,1
+ ,80570
+ ,6984
+ ,597793
+ ,42
+ ,1
+ ,102129
+ ,4519
+ ,644190
+ ,12
+ ,0
+ ,301670
+ ,2220
+ ,377934
+ ,211
+ ,0
+ ,102313
+ ,18562
+ ,640273
+ ,74
+ ,0
+ ,88577
+ ,10327
+ ,697458
+ ,80
+ ,1
+ ,112477
+ ,5336
+ ,550608
+ ,83
+ ,1
+ ,191778
+ ,2365
+ ,207393
+ ,131
+ ,0
+ ,79804
+ ,4069
+ ,301607
+ ,203
+ ,0
+ ,128294
+ ,8636
+ ,345783
+ ,56
+ ,0
+ ,96448
+ ,13718
+ ,501749
+ ,89
+ ,0
+ ,93811
+ ,4525
+ ,379983
+ ,88
+ ,0
+ ,117520
+ ,6869
+ ,387475
+ ,39
+ ,0
+ ,69159
+ ,4628
+ ,377305
+ ,25
+ ,1
+ ,101792
+ ,3689
+ ,370837
+ ,49
+ ,1
+ ,210568
+ ,4891
+ ,430866
+ ,149
+ ,1
+ ,136996
+ ,7489
+ ,469107
+ ,58
+ ,0
+ ,121920
+ ,4901
+ ,194493
+ ,41
+ ,0
+ ,76403
+ ,2284
+ ,530670
+ ,90
+ ,1
+ ,108094
+ ,3160
+ ,518365
+ ,136
+ ,1
+ ,134759
+ ,4150
+ ,491303
+ ,97
+ ,1
+ ,188873
+ ,7285
+ ,527021
+ ,63
+ ,1
+ ,146216
+ ,1134
+ ,233773
+ ,114
+ ,1
+ ,156608
+ ,4658
+ ,405972
+ ,77
+ ,0
+ ,61348
+ ,2384
+ ,652925
+ ,6
+ ,0
+ ,50350
+ ,3748
+ ,446211
+ ,47
+ ,0
+ ,87720
+ ,5371
+ ,341340
+ ,51
+ ,0
+ ,99489
+ ,1285
+ ,387699
+ ,85
+ ,1
+ ,87419
+ ,9327
+ ,493408
+ ,43
+ ,1
+ ,94355
+ ,5565
+ ,146494
+ ,32
+ ,0
+ ,60326
+ ,1528
+ ,414462
+ ,25
+ ,1
+ ,94670
+ ,3122
+ ,364304
+ ,77
+ ,1
+ ,82425
+ ,7561
+ ,355178
+ ,54
+ ,0
+ ,59017
+ ,2675
+ ,357760
+ ,251
+ ,0
+ ,90829
+ ,13253
+ ,261216
+ ,15
+ ,0
+ ,80791
+ ,880
+ ,397144
+ ,44
+ ,1
+ ,100423
+ ,2053
+ ,374943
+ ,73
+ ,0
+ ,131116
+ ,1424
+ ,424898
+ ,85
+ ,1
+ ,100269
+ ,4036
+ ,202055
+ ,49
+ ,1
+ ,27330
+ ,3045
+ ,378525
+ ,38
+ ,0
+ ,39039
+ ,5119
+ ,310768
+ ,35
+ ,0
+ ,106885
+ ,1431
+ ,325738
+ ,9
+ ,0
+ ,79285
+ ,554
+ ,394510
+ ,34
+ ,0
+ ,118881
+ ,1975
+ ,247060
+ ,20
+ ,1
+ ,77623
+ ,1765
+ ,368078
+ ,29
+ ,0
+ ,114768
+ ,1012
+ ,236761
+ ,11
+ ,0
+ ,74015
+ ,810
+ ,312378
+ ,52
+ ,0
+ ,69465
+ ,1280
+ ,339836
+ ,13
+ ,1
+ ,117869
+ ,666
+ ,347385
+ ,29
+ ,0
+ ,60982
+ ,1380
+ ,426280
+ ,66
+ ,1
+ ,90131
+ ,4677
+ ,352850
+ ,33
+ ,0
+ ,138971
+ ,876
+ ,301881
+ ,15
+ ,0
+ ,39625
+ ,814
+ ,377516
+ ,15
+ ,0
+ ,102725
+ ,514
+ ,357312
+ ,68
+ ,1
+ ,64239
+ ,5692
+ ,458343
+ ,100
+ ,0
+ ,90262
+ ,3642
+ ,354228
+ ,13
+ ,0
+ ,103960
+ ,540
+ ,308636
+ ,45
+ ,0
+ ,106611
+ ,2099
+ ,386212
+ ,14
+ ,0
+ ,103345
+ ,567
+ ,393343
+ ,36
+ ,0
+ ,95551
+ ,2001
+ ,378509
+ ,40
+ ,1
+ ,82903
+ ,2949
+ ,452469
+ ,68
+ ,0
+ ,63593
+ ,2253
+ ,364839
+ ,29
+ ,1
+ ,126910
+ ,6533
+ ,358649
+ ,43
+ ,0
+ ,37527
+ ,1889
+ ,376641
+ ,30
+ ,1
+ ,60247
+ ,3055
+ ,429112
+ ,9
+ ,0
+ ,112995
+ ,272
+ ,330546
+ ,22
+ ,1
+ ,70184
+ ,1414
+ ,403560
+ ,19
+ ,0
+ ,130140
+ ,2564
+ ,317892
+ ,9
+ ,1
+ ,73221
+ ,1383
+ ,307528
+ ,31
+ ,1
+ ,76114
+ ,1261
+ ,235133
+ ,19
+ ,0
+ ,90534
+ ,975
+ ,299243
+ ,55
+ ,0
+ ,108479
+ ,3366
+ ,314073
+ ,8
+ ,0
+ ,113761
+ ,576
+ ,368186
+ ,28
+ ,0
+ ,68696
+ ,1686
+ ,269661
+ ,29
+ ,0
+ ,71561
+ ,746
+ ,125390
+ ,48
+ ,1
+ ,59831
+ ,3192
+ ,510834
+ ,16
+ ,1
+ ,97890
+ ,2045
+ ,321896
+ ,47
+ ,0
+ ,101481
+ ,5702
+ ,249898
+ ,20
+ ,1
+ ,72954
+ ,1932
+ ,408881
+ ,22
+ ,0
+ ,67939
+ ,936
+ ,158492
+ ,33
+ ,1
+ ,48022
+ ,3437
+ ,292154
+ ,44
+ ,0
+ ,86111
+ ,5131
+ ,289513
+ ,13
+ ,1
+ ,74020
+ ,2397
+ ,378049
+ ,6
+ ,1
+ ,57530
+ ,1389
+ ,343466
+ ,35
+ ,0
+ ,56364
+ ,1503
+ ,332743
+ ,8
+ ,0
+ ,84990
+ ,402
+ ,442882
+ ,17
+ ,0
+ ,88590
+ ,2239
+ ,214215
+ ,11
+ ,1
+ ,77200
+ ,2234
+ ,315688
+ ,21
+ ,0
+ ,61262
+ ,837
+ ,375195
+ ,92
+ ,0
+ ,110309
+ ,10579
+ ,334280
+ ,12
+ ,0
+ ,67000
+ ,875
+ ,355864
+ ,112
+ ,0
+ ,93099
+ ,1585
+ ,480382
+ ,25
+ ,1
+ ,107577
+ ,1659
+ ,353058
+ ,17
+ ,1
+ ,62920
+ ,2647
+ ,217193
+ ,23
+ ,1
+ ,75832
+ ,3294
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,314533
+ ,10
+ ,0
+ ,60793
+ ,94
+ ,318056
+ ,23
+ ,0
+ ,57935
+ ,422
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,314353
+ ,7
+ ,0
+ ,60630
+ ,34
+ ,369448
+ ,25
+ ,0
+ ,55637
+ ,1558
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,312846
+ ,20
+ ,0
+ ,60887
+ ,43
+ ,312075
+ ,4
+ ,0
+ ,60720
+ ,645
+ ,315009
+ ,4
+ ,0
+ ,60505
+ ,316
+ ,318903
+ ,10
+ ,0
+ ,60945
+ ,115
+ ,314887
+ ,1
+ ,0
+ ,60720
+ ,5
+ ,314913
+ ,4
+ ,0
+ ,60720
+ ,897
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325506
+ ,8
+ ,0
+ ,58990
+ ,389
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,298568
+ ,11
+ ,0
+ ,56750
+ ,1002
+ ,315834
+ ,4
+ ,0
+ ,60894
+ ,36
+ ,329784
+ ,15
+ ,0
+ ,63346
+ ,460
+ ,312878
+ ,9
+ ,0
+ ,56535
+ ,309
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,314987
+ ,7
+ ,0
+ ,60835
+ ,9
+ ,325249
+ ,2
+ ,0
+ ,60720
+ ,271
+ ,315877
+ ,0
+ ,0
+ ,61016
+ ,14
+ ,291650
+ ,7
+ ,0
+ ,58650
+ ,520
+ ,305959
+ ,46
+ ,0
+ ,60438
+ ,1766
+ ,315380
+ ,5
+ ,0
+ ,60720
+ ,0
+ ,297765
+ ,7
+ ,0
+ ,58625
+ ,458
+ ,315245
+ ,2
+ ,0
+ ,60938
+ ,20
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315236
+ ,2
+ ,0
+ ,61490
+ ,98
+ ,336425
+ ,5
+ ,0
+ ,60845
+ ,405
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,306268
+ ,7
+ ,0
+ ,60830
+ ,483
+ ,302187
+ ,24
+ ,0
+ ,63261
+ ,454
+ ,314882
+ ,1
+ ,0
+ ,60720
+ ,47
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,382712
+ ,18
+ ,0
+ ,45689
+ ,757
+ ,341570
+ ,55
+ ,0
+ ,60720
+ ,4655
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,312412
+ ,3
+ ,0
+ ,61564
+ ,36
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,309596
+ ,9
+ ,0
+ ,61938
+ ,203
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315547
+ ,8
+ ,0
+ ,60951
+ ,126
+ ,313267
+ ,113
+ ,0
+ ,60720
+ ,400
+ ,316176
+ ,0
+ ,1
+ ,60745
+ ,71
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,359335
+ ,19
+ ,0
+ ,71642
+ ,972
+ ,330068
+ ,11
+ ,1
+ ,71641
+ ,531
+ ,314289
+ ,25
+ ,0
+ ,55792
+ ,2461
+ ,297413
+ ,16
+ ,1
+ ,71873
+ ,378
+ ,314806
+ ,5
+ ,1
+ ,62555
+ ,23
+ ,333210
+ ,11
+ ,1
+ ,60370
+ ,638
+ ,352108
+ ,23
+ ,1
+ ,64873
+ ,2300
+ ,313332
+ ,6
+ ,0
+ ,62041
+ ,149
+ ,291787
+ ,5
+ ,0
+ ,65745
+ ,226
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,318745
+ ,7
+ ,0
+ ,59500
+ ,275
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315366
+ ,7
+ ,0
+ ,61630
+ ,141
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315688
+ ,3
+ ,0
+ ,60890
+ ,28
+ ,315380
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,409642
+ ,89
+ ,1
+ ,113521
+ ,4980
+ ,315380
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,269587
+ ,19
+ ,1
+ ,80045
+ ,472
+ ,315380
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,300962
+ ,12
+ ,1
+ ,50804
+ ,203
+ ,325479
+ ,12
+ ,1
+ ,87390
+ ,496
+ ,316155
+ ,5
+ ,1
+ ,61656
+ ,10
+ ,318574
+ ,2
+ ,1
+ ,65688
+ ,63
+ ,315380
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,343613
+ ,26
+ ,1
+ ,48522
+ ,1136
+ ,306948
+ ,3
+ ,1
+ ,60720
+ ,265
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,330059
+ ,11
+ ,0
+ ,57640
+ ,267
+ ,288985
+ ,10
+ ,0
+ ,61977
+ ,474
+ ,304485
+ ,5
+ ,0
+ ,62620
+ ,534
+ ,315380
+ ,2
+ ,0
+ ,60720
+ ,0
+ ,315688
+ ,6
+ ,0
+ ,60831
+ ,15
+ ,317736
+ ,7
+ ,0
+ ,60646
+ ,397
+ ,315380
+ ,2
+ ,0
+ ,60720
+ ,0
+ ,322331
+ ,28
+ ,0
+ ,56225
+ ,1866
+ ,296656
+ ,3
+ ,0
+ ,60510
+ ,288
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315354
+ ,1
+ ,0
+ ,60698
+ ,3
+ ,312161
+ ,20
+ ,0
+ ,60720
+ ,468
+ ,315576
+ ,1
+ ,0
+ ,60805
+ ,20
+ ,314922
+ ,22
+ ,0
+ ,61404
+ ,278
+ ,314551
+ ,9
+ ,0
+ ,60720
+ ,61
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,312339
+ ,2
+ ,0
+ ,65276
+ ,192
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,298700
+ ,7
+ ,0
+ ,63915
+ ,317
+ ,321376
+ ,9
+ ,0
+ ,60720
+ ,738
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,303230
+ ,13
+ ,0
+ ,61686
+ ,368
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315487
+ ,0
+ ,0
+ ,60743
+ ,2
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315793
+ ,6
+ ,0
+ ,60349
+ ,53
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,312887
+ ,3
+ ,0
+ ,61360
+ ,94
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315637
+ ,7
+ ,0
+ ,59818
+ ,24
+ ,324385
+ ,2
+ ,0
+ ,72680
+ ,2332
+ ,315380
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,308989
+ ,15
+ ,0
+ ,61808
+ ,131
+ ,315380
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,296702
+ ,9
+ ,0
+ ,53110
+ ,206
+ ,315380
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,307322
+ ,1
+ ,0
+ ,64245
+ ,167
+ ,304376
+ ,38
+ ,0
+ ,73007
+ ,622
+ ,253588
+ ,57
+ ,0
+ ,82732
+ ,2328
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,309560
+ ,7
+ ,0
+ ,54820
+ ,365
+ ,298466
+ ,26
+ ,0
+ ,47705
+ ,364
+ ,315380
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,343929
+ ,13
+ ,0
+ ,72835
+ ,226
+ ,331955
+ ,10
+ ,0
+ ,58856
+ ,307
+ ,315380
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,381180
+ ,9
+ ,0
+ ,77655
+ ,188
+ ,315380
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,331420
+ ,26
+ ,0
+ ,69817
+ ,138
+ ,315380
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,310201
+ ,19
+ ,0
+ ,60798
+ ,125
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,320016
+ ,12
+ ,0
+ ,62452
+ ,282
+ ,320398
+ ,23
+ ,0
+ ,64175
+ ,335
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,291841
+ ,29
+ ,1
+ ,67440
+ ,1324
+ ,310670
+ ,8
+ ,0
+ ,68136
+ ,176
+ ,315380
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,313491
+ ,26
+ ,0
+ ,56726
+ ,249
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,331323
+ ,9
+ ,0
+ ,70811
+ ,333
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,319210
+ ,5
+ ,1
+ ,60720
+ ,601
+ ,318098
+ ,3
+ ,0
+ ,62045
+ ,30
+ ,315380
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,292754
+ ,13
+ ,0
+ ,54323
+ ,249
+ ,315380
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,325176
+ ,12
+ ,0
+ ,62841
+ ,165
+ ,365959
+ ,19
+ ,0
+ ,81125
+ ,453
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,302409
+ ,10
+ ,0
+ ,59506
+ ,53
+ ,340968
+ ,9
+ ,1
+ ,59365
+ ,382
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,9
+ ,0
+ ,60720
+ ,0
+ ,313164
+ ,4
+ ,1
+ ,60798
+ ,30
+ ,301164
+ ,1
+ ,0
+ ,58790
+ ,290
+ ,315380
+ ,1
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,344425
+ ,14
+ ,0
+ ,61808
+ ,366
+ ,315394
+ ,12
+ ,0
+ ,60735
+ ,2
+ ,315380
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,316647
+ ,19
+ ,1
+ ,64016
+ ,209
+ ,309836
+ ,17
+ ,0
+ ,54683
+ ,384
+ ,315380
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,346611
+ ,32
+ ,0
+ ,87192
+ ,365
+ ,315380
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,322031
+ ,14
+ ,1
+ ,64107
+ ,49
+ ,315656
+ ,8
+ ,0
+ ,60761
+ ,3
+ ,339445
+ ,4
+ ,0
+ ,65990
+ ,133
+ ,314964
+ ,0
+ ,0
+ ,59988
+ ,32
+ ,297141
+ ,20
+ ,0
+ ,61167
+ ,368
+ ,315372
+ ,5
+ ,0
+ ,60719
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,312502
+ ,1
+ ,0
+ ,60722
+ ,22
+ ,315380
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,313729
+ ,4
+ ,0
+ ,60379
+ ,96
+ ,315388
+ ,1
+ ,0
+ ,60727
+ ,1
+ ,315371
+ ,4
+ ,0
+ ,60720
+ ,314
+ ,296139
+ ,20
+ ,0
+ ,60925
+ ,844
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,313880
+ ,1
+ ,0
+ ,60896
+ ,26
+ ,317698
+ ,10
+ ,0
+ ,59734
+ ,125
+ ,295580
+ ,12
+ ,0
+ ,62969
+ ,304
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,308256
+ ,13
+ ,0
+ ,60720
+ ,621
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,303677
+ ,3
+ ,0
+ ,59118
+ ,119
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,319369
+ ,10
+ ,0
+ ,60720
+ ,1595
+ ,318690
+ ,3
+ ,0
+ ,58598
+ ,312
+ ,314049
+ ,7
+ ,0
+ ,61124
+ ,60
+ ,325699
+ ,10
+ ,0
+ ,59595
+ ,587
+ ,314210
+ ,1
+ ,0
+ ,62065
+ ,135
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,322378
+ ,15
+ ,0
+ ,78780
+ ,514
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315398
+ ,4
+ ,0
+ ,60722
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,308336
+ ,28
+ ,1
+ ,61600
+ ,1763
+ ,316386
+ ,9
+ ,0
+ ,59635
+ ,180
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315553
+ ,7
+ ,0
+ ,60720
+ ,218
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,323361
+ ,7
+ ,0
+ ,59781
+ ,448
+ ,336639
+ ,7
+ ,0
+ ,76644
+ ,227
+ ,307424
+ ,3
+ ,0
+ ,64820
+ ,174
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,295370
+ ,11
+ ,0
+ ,56178
+ ,121
+ ,322340
+ ,7
+ ,0
+ ,60436
+ ,607
+ ,319864
+ ,10
+ ,0
+ ,60720
+ ,2212
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,317291
+ ,18
+ ,0
+ ,73433
+ ,530
+ ,280398
+ ,14
+ ,0
+ ,41477
+ ,571
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,317330
+ ,12
+ ,0
+ ,62700
+ ,78
+ ,238125
+ ,29
+ ,0
+ ,67804
+ ,2489
+ ,327071
+ ,3
+ ,0
+ ,59661
+ ,131
+ ,309038
+ ,6
+ ,0
+ ,58620
+ ,923
+ ,314210
+ ,3
+ ,0
+ ,60398
+ ,72
+ ,307930
+ ,8
+ ,0
+ ,58580
+ ,572
+ ,322327
+ ,10
+ ,0
+ ,62710
+ ,397
+ ,292136
+ ,6
+ ,0
+ ,59325
+ ,450
+ ,263276
+ ,8
+ ,0
+ ,60950
+ ,622
+ ,367655
+ ,6
+ ,0
+ ,68060
+ ,694
+ ,283910
+ ,9
+ ,1
+ ,83620
+ ,3425
+ ,283587
+ ,8
+ ,0
+ ,58456
+ ,562
+ ,243650
+ ,26
+ ,0
+ ,52811
+ ,4917
+ ,438493
+ ,239
+ ,1
+ ,121173
+ ,1442
+ ,296261
+ ,7
+ ,0
+ ,63870
+ ,529
+ ,230621
+ ,41
+ ,1
+ ,21001
+ ,2126
+ ,304252
+ ,3
+ ,0
+ ,70415
+ ,1061
+ ,333505
+ ,8
+ ,0
+ ,64230
+ ,776
+ ,296919
+ ,6
+ ,0
+ ,59190
+ ,611
+ ,278990
+ ,21
+ ,1
+ ,69351
+ ,1526
+ ,276898
+ ,7
+ ,0
+ ,64270
+ ,592
+ ,327007
+ ,11
+ ,0
+ ,70694
+ ,1182
+ ,317046
+ ,11
+ ,0
+ ,68005
+ ,621
+ ,304555
+ ,12
+ ,0
+ ,58930
+ ,989
+ ,298096
+ ,9
+ ,0
+ ,58320
+ ,438
+ ,231861
+ ,3
+ ,0
+ ,69980
+ ,726
+ ,309422
+ ,57
+ ,0
+ ,69863
+ ,1303
+ ,286963
+ ,21
+ ,1
+ ,63255
+ ,7419
+ ,269753
+ ,15
+ ,1
+ ,57320
+ ,1164
+ ,448243
+ ,32
+ ,1
+ ,75230
+ ,3310
+ ,165404
+ ,11
+ ,0
+ ,79420
+ ,1920
+ ,204325
+ ,2
+ ,0
+ ,73490
+ ,965
+ ,407159
+ ,23
+ ,0
+ ,35250
+ ,3256
+ ,290476
+ ,20
+ ,1
+ ,62285
+ ,1135
+ ,275311
+ ,24
+ ,0
+ ,69206
+ ,1270
+ ,246541
+ ,1
+ ,0
+ ,65920
+ ,661
+ ,253468
+ ,1
+ ,0
+ ,69770
+ ,1013
+ ,240897
+ ,74
+ ,0
+ ,72683
+ ,2844
+ ,-83265
+ ,68
+ ,1
+ ,-14545
+ ,11528
+ ,-42143
+ ,20
+ ,0
+ ,55830
+ ,6526
+ ,272713
+ ,20
+ ,0
+ ,55174
+ ,2264
+ ,215362
+ ,82
+ ,1
+ ,67038
+ ,5109
+ ,42754
+ ,21
+ ,0
+ ,51252
+ ,3999
+ ,306275
+ ,244
+ ,0
+ ,157278
+ ,35624
+ ,253537
+ ,32
+ ,0
+ ,79510
+ ,9252
+ ,372631
+ ,86
+ ,0
+ ,77440
+ ,15236
+ ,-7170
+ ,69
+ ,0
+ ,27284
+ ,18073)
+ ,dim=c(5
+ ,431)
+ ,dimnames=list(c('wealth'
+ ,'Orders'
+ ,'Group'
+ ,'Dividends'
+ ,'Costs
')
+ ,1:431))
> y <- array(NA,dim=c(5,431),dimnames=list(c('wealth','Orders','Group','Dividends','Costs
'),1:431))
> for (i in 1:dim(x)[1])
+ {
+ for (j in 1:dim(x)[2])
+ {
+ y[i,j] <- as.numeric(x[i,j])
+ }
+ }
> par3 = 'No Linear Trend'
> par2 = 'Do not include Seasonal Dummies'
> par1 = '1'
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Prof. Dr. P. Wessa
> #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #Technical description: Write here your technical program description (don't use hard returns!)
> library(lattice)
> library(lmtest)
Loading required package: zoo
Attaching package: 'zoo'
The following object(s) are masked from package:base :
as.Date.numeric
> n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test
> par1 <- as.numeric(par1)
> x <- t(y)
> k <- length(x[1,])
> n <- length(x[,1])
> x1 <- cbind(x[,par1], x[,1:k!=par1])
> mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1])
> colnames(x1) <- mycolnames #colnames(x)[par1]
> x <- x1
> if (par3 == 'First Differences'){
+ x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep='')))
+ for (i in 1:n-1) {
+ for (j in 1:k) {
+ x2[i,j] <- x[i+1,j] - x[i,j]
+ }
+ }
+ x <- x2
+ }
> if (par2 == 'Include Monthly Dummies'){
+ x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep ='')))
+ for (i in 1:11){
+ x2[seq(i,n,12),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> if (par2 == 'Include Quarterly Dummies'){
+ x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep ='')))
+ for (i in 1:3){
+ x2[seq(i,n,4),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> k <- length(x[1,])
> if (par3 == 'Linear Trend'){
+ x <- cbind(x, c(1:n))
+ colnames(x)[k+1] <- 't'
+ }
> x
wealth Orders Group Dividends Costs\r
1 6282154 807 1 213118 162556
2 4321023 444 1 81767 29790
3 4111912 412 1 153198 87550
4 223193 428 0 -26007 84738
5 1491348 315 1 126942 54660
6 1629616 168 1 157214 42634
7 1398893 263 0 129352 40949
8 1926517 267 1 234817 45187
9 983660 228 1 60448 37704
10 1443586 129 1 47818 16275
11 1073089 104 0 245546 25830
12 984885 122 0 48020 12679
13 1405225 393 1 -1710 18014
14 227132 190 0 32648 43556
15 929118 280 1 95350 24811
16 1071292 63 0 151352 6575
17 638830 102 0 288170 7123
18 856956 265 1 114337 21950
19 992426 234 1 37884 37597
20 444477 277 0 122844 17821
21 857217 73 1 82340 12988
22 711969 67 1 79801 22330
23 702380 103 0 165548 13326
24 358589 290 0 116384 16189
25 297978 83 0 134028 7146
26 585715 56 0 63838 15824
27 657954 236 1 74996 27664
28 209458 73 0 31080 11920
29 786690 34 0 32168 8568
30 439798 139 0 49857 14416
31 688779 26 1 87161 3369
32 574339 70 1 106113 11819
33 741409 40 1 80570 6984
34 597793 42 1 102129 4519
35 644190 12 0 301670 2220
36 377934 211 0 102313 18562
37 640273 74 0 88577 10327
38 697458 80 1 112477 5336
39 550608 83 1 191778 2365
40 207393 131 0 79804 4069
41 301607 203 0 128294 8636
42 345783 56 0 96448 13718
43 501749 89 0 93811 4525
44 379983 88 0 117520 6869
45 387475 39 0 69159 4628
46 377305 25 1 101792 3689
47 370837 49 1 210568 4891
48 430866 149 1 136996 7489
49 469107 58 0 121920 4901
50 194493 41 0 76403 2284
51 530670 90 1 108094 3160
52 518365 136 1 134759 4150
53 491303 97 1 188873 7285
54 527021 63 1 146216 1134
55 233773 114 1 156608 4658
56 405972 77 0 61348 2384
57 652925 6 0 50350 3748
58 446211 47 0 87720 5371
59 341340 51 0 99489 1285
60 387699 85 1 87419 9327
61 493408 43 1 94355 5565
62 146494 32 0 60326 1528
63 414462 25 1 94670 3122
64 364304 77 1 82425 7561
65 355178 54 0 59017 2675
66 357760 251 0 90829 13253
67 261216 15 0 80791 880
68 397144 44 1 100423 2053
69 374943 73 0 131116 1424
70 424898 85 1 100269 4036
71 202055 49 1 27330 3045
72 378525 38 0 39039 5119
73 310768 35 0 106885 1431
74 325738 9 0 79285 554
75 394510 34 0 118881 1975
76 247060 20 1 77623 1765
77 368078 29 0 114768 1012
78 236761 11 0 74015 810
79 312378 52 0 69465 1280
80 339836 13 1 117869 666
81 347385 29 0 60982 1380
82 426280 66 1 90131 4677
83 352850 33 0 138971 876
84 301881 15 0 39625 814
85 377516 15 0 102725 514
86 357312 68 1 64239 5692
87 458343 100 0 90262 3642
88 354228 13 0 103960 540
89 308636 45 0 106611 2099
90 386212 14 0 103345 567
91 393343 36 0 95551 2001
92 378509 40 1 82903 2949
93 452469 68 0 63593 2253
94 364839 29 1 126910 6533
95 358649 43 0 37527 1889
96 376641 30 1 60247 3055
97 429112 9 0 112995 272
98 330546 22 1 70184 1414
99 403560 19 0 130140 2564
100 317892 9 1 73221 1383
101 307528 31 1 76114 1261
102 235133 19 0 90534 975
103 299243 55 0 108479 3366
104 314073 8 0 113761 576
105 368186 28 0 68696 1686
106 269661 29 0 71561 746
107 125390 48 1 59831 3192
108 510834 16 1 97890 2045
109 321896 47 0 101481 5702
110 249898 20 1 72954 1932
111 408881 22 0 67939 936
112 158492 33 1 48022 3437
113 292154 44 0 86111 5131
114 289513 13 1 74020 2397
115 378049 6 1 57530 1389
116 343466 35 0 56364 1503
117 332743 8 0 84990 402
118 442882 17 0 88590 2239
119 214215 11 1 77200 2234
120 315688 21 0 61262 837
121 375195 92 0 110309 10579
122 334280 12 0 67000 875
123 355864 112 0 93099 1585
124 480382 25 1 107577 1659
125 353058 17 1 62920 2647
126 217193 23 1 75832 3294
127 315380 0 0 60720 0
128 314533 10 0 60793 94
129 318056 23 0 57935 422
130 315380 0 0 60720 0
131 314353 7 0 60630 34
132 369448 25 0 55637 1558
133 315380 0 0 60720 0
134 312846 20 0 60887 43
135 312075 4 0 60720 645
136 315009 4 0 60505 316
137 318903 10 0 60945 115
138 314887 1 0 60720 5
139 314913 4 0 60720 897
140 315380 0 0 60720 0
141 325506 8 0 58990 389
142 315380 0 0 60720 0
143 298568 11 0 56750 1002
144 315834 4 0 60894 36
145 329784 15 0 63346 460
146 312878 9 0 56535 309
147 315380 0 0 60720 0
148 314987 7 0 60835 9
149 325249 2 0 60720 271
150 315877 0 0 61016 14
151 291650 7 0 58650 520
152 305959 46 0 60438 1766
153 315380 5 0 60720 0
154 297765 7 0 58625 458
155 315245 2 0 60938 20
156 315380 0 0 60720 0
157 315380 0 0 60720 0
158 315236 2 0 61490 98
159 336425 5 0 60845 405
160 315380 0 0 60720 0
161 315380 0 0 60720 0
162 315380 0 0 60720 0
163 315380 0 0 60720 0
164 306268 7 0 60830 483
165 302187 24 0 63261 454
166 314882 1 0 60720 47
167 315380 0 0 60720 0
168 382712 18 0 45689 757
169 341570 55 0 60720 4655
170 315380 0 0 60720 0
171 315380 0 0 60720 0
172 312412 3 0 61564 36
173 315380 0 0 60720 0
174 309596 9 0 61938 203
175 315380 0 0 60720 0
176 315547 8 0 60951 126
177 313267 113 0 60720 400
178 316176 0 1 60745 71
179 315380 0 0 60720 0
180 315380 0 0 60720 0
181 359335 19 0 71642 972
182 330068 11 1 71641 531
183 314289 25 0 55792 2461
184 297413 16 1 71873 378
185 314806 5 1 62555 23
186 333210 11 1 60370 638
187 352108 23 1 64873 2300
188 313332 6 0 62041 149
189 291787 5 0 65745 226
190 315380 0 0 60720 0
191 318745 7 0 59500 275
192 315380 0 0 60720 0
193 315366 7 0 61630 141
194 315380 0 0 60720 0
195 315688 3 0 60890 28
196 315380 0 1 60720 0
197 409642 89 1 113521 4980
198 315380 0 1 60720 0
199 315380 0 1 60720 0
200 269587 19 1 80045 472
201 315380 0 1 60720 0
202 315380 0 1 60720 0
203 315380 0 1 60720 0
204 300962 12 1 50804 203
205 325479 12 1 87390 496
206 316155 5 1 61656 10
207 318574 2 1 65688 63
208 315380 0 1 60720 0
209 343613 26 1 48522 1136
210 306948 3 1 60720 265
211 315380 0 0 60720 0
212 315380 0 0 60720 0
213 330059 11 0 57640 267
214 288985 10 0 61977 474
215 304485 5 0 62620 534
216 315380 2 0 60720 0
217 315688 6 0 60831 15
218 317736 7 0 60646 397
219 315380 2 0 60720 0
220 322331 28 0 56225 1866
221 296656 3 0 60510 288
222 315380 0 0 60720 0
223 315354 1 0 60698 3
224 312161 20 0 60720 468
225 315576 1 0 60805 20
226 314922 22 0 61404 278
227 314551 9 0 60720 61
228 315380 0 0 60720 0
229 312339 2 0 65276 192
230 315380 0 0 60720 0
231 298700 7 0 63915 317
232 321376 9 0 60720 738
233 315380 0 0 60720 0
234 303230 13 0 61686 368
235 315380 0 0 60720 0
236 315487 0 0 60743 2
237 315380 0 0 60720 0
238 315793 6 0 60349 53
239 315380 0 0 60720 0
240 315380 0 0 60720 0
241 315380 0 0 60720 0
242 312887 3 0 61360 94
243 315380 0 0 60720 0
244 315637 7 0 59818 24
245 324385 2 0 72680 2332
246 315380 0 1 60720 0
247 315380 0 1 60720 0
248 308989 15 0 61808 131
249 315380 0 1 60720 0
250 315380 0 1 60720 0
251 296702 9 0 53110 206
252 315380 0 1 60720 0
253 307322 1 0 64245 167
254 304376 38 0 73007 622
255 253588 57 0 82732 2328
256 315380 0 0 60720 0
257 309560 7 0 54820 365
258 298466 26 0 47705 364
259 315380 0 1 60720 0
260 315380 0 0 60720 0
261 315380 0 1 60720 0
262 315380 0 1 60720 0
263 343929 13 0 72835 226
264 331955 10 0 58856 307
265 315380 0 1 60720 0
266 315380 0 0 60720 0
267 315380 0 1 60720 0
268 381180 9 0 77655 188
269 315380 0 1 60720 0
270 331420 26 0 69817 138
271 315380 0 1 60720 0
272 315380 0 1 60720 0
273 315380 0 1 60720 0
274 310201 19 0 60798 125
275 315380 0 0 60720 0
276 320016 12 0 62452 282
277 320398 23 0 64175 335
278 315380 0 0 60720 0
279 291841 29 1 67440 1324
280 310670 8 0 68136 176
281 315380 0 1 60720 0
282 315380 0 0 60720 0
283 313491 26 0 56726 249
284 315380 0 0 60720 0
285 331323 9 0 70811 333
286 315380 0 0 60720 0
287 319210 5 1 60720 601
288 318098 3 0 62045 30
289 315380 0 1 60720 0
290 292754 13 0 54323 249
291 315380 0 1 60720 0
292 325176 12 0 62841 165
293 365959 19 0 81125 453
294 315380 0 0 60720 0
295 302409 10 0 59506 53
296 340968 9 1 59365 382
297 315380 0 0 60720 0
298 315380 0 1 60720 0
299 315380 0 0 60720 0
300 315380 9 0 60720 0
301 313164 4 1 60798 30
302 301164 1 0 58790 290
303 315380 1 0 60720 0
304 315380 0 1 60720 0
305 344425 14 0 61808 366
306 315394 12 0 60735 2
307 315380 0 1 60720 0
308 316647 19 1 64016 209
309 309836 17 0 54683 384
310 315380 0 1 60720 0
311 315380 0 1 60720 0
312 346611 32 0 87192 365
313 315380 0 1 60720 0
314 322031 14 1 64107 49
315 315656 8 0 60761 3
316 339445 4 0 65990 133
317 314964 0 0 59988 32
318 297141 20 0 61167 368
319 315372 5 0 60719 1
320 315380 0 0 60720 0
321 315380 0 0 60720 0
322 315380 0 0 60720 0
323 315380 0 0 60720 0
324 315380 0 1 60720 0
325 315380 0 0 60720 0
326 312502 1 0 60722 22
327 315380 0 1 60720 0
328 315380 0 1 60720 0
329 315380 0 0 60720 0
330 315380 0 0 60720 0
331 315380 0 0 60720 0
332 315380 0 0 60720 0
333 315380 0 0 60720 0
334 313729 4 0 60379 96
335 315388 1 0 60727 1
336 315371 4 0 60720 314
337 296139 20 0 60925 844
338 315380 0 0 60720 0
339 313880 1 0 60896 26
340 317698 10 0 59734 125
341 295580 12 0 62969 304
342 315380 0 0 60720 0
343 315380 0 0 60720 0
344 315380 0 0 60720 0
345 308256 13 0 60720 621
346 315380 0 0 60720 0
347 303677 3 0 59118 119
348 315380 0 0 60720 0
349 315380 0 0 60720 0
350 319369 10 0 60720 1595
351 318690 3 0 58598 312
352 314049 7 0 61124 60
353 325699 10 0 59595 587
354 314210 1 0 62065 135
355 315380 0 0 60720 0
356 315380 0 0 60720 0
357 322378 15 0 78780 514
358 315380 0 0 60720 0
359 315380 0 0 60720 0
360 315380 0 0 60720 0
361 315398 4 0 60722 1
362 315380 0 0 60720 0
363 315380 0 0 60720 0
364 308336 28 1 61600 1763
365 316386 9 0 59635 180
366 315380 0 0 60720 0
367 315380 0 0 60720 0
368 315380 0 0 60720 0
369 315380 0 0 60720 0
370 315553 7 0 60720 218
371 315380 0 0 60720 0
372 323361 7 0 59781 448
373 336639 7 0 76644 227
374 307424 3 0 64820 174
375 315380 0 0 60720 0
376 315380 0 0 60720 0
377 295370 11 0 56178 121
378 322340 7 0 60436 607
379 319864 10 0 60720 2212
380 315380 0 0 60720 0
381 315380 0 0 60720 0
382 317291 18 0 73433 530
383 280398 14 0 41477 571
384 315380 0 0 60720 0
385 317330 12 0 62700 78
386 238125 29 0 67804 2489
387 327071 3 0 59661 131
388 309038 6 0 58620 923
389 314210 3 0 60398 72
390 307930 8 0 58580 572
391 322327 10 0 62710 397
392 292136 6 0 59325 450
393 263276 8 0 60950 622
394 367655 6 0 68060 694
395 283910 9 1 83620 3425
396 283587 8 0 58456 562
397 243650 26 0 52811 4917
398 438493 239 1 121173 1442
399 296261 7 0 63870 529
400 230621 41 1 21001 2126
401 304252 3 0 70415 1061
402 333505 8 0 64230 776
403 296919 6 0 59190 611
404 278990 21 1 69351 1526
405 276898 7 0 64270 592
406 327007 11 0 70694 1182
407 317046 11 0 68005 621
408 304555 12 0 58930 989
409 298096 9 0 58320 438
410 231861 3 0 69980 726
411 309422 57 0 69863 1303
412 286963 21 1 63255 7419
413 269753 15 1 57320 1164
414 448243 32 1 75230 3310
415 165404 11 0 79420 1920
416 204325 2 0 73490 965
417 407159 23 0 35250 3256
418 290476 20 1 62285 1135
419 275311 24 0 69206 1270
420 246541 1 0 65920 661
421 253468 1 0 69770 1013
422 240897 74 0 72683 2844
423 -83265 68 1 -14545 11528
424 -42143 20 0 55830 6526
425 272713 20 0 55174 2264
426 215362 82 1 67038 5109
427 42754 21 0 51252 3999
428 306275 244 0 157278 35624
429 253537 32 0 79510 9252
430 372631 86 0 77440 15236
431 -7170 69 0 27284 18073
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) Orders Group Dividends `Costs\r`
1.573e+05 9.510e+02 4.756e+04 1.518e+00 2.204e+01
> (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) 1.573e+05 3.090e+04 5.089 5.4e-07 ***
Orders 9.510e+02 3.506e+02 2.713 0.006947 **
Group 4.756e+04 2.793e+04 1.703 0.089345 .
Dividends 1.518e+00 4.116e-01 3.689 0.000254 ***
`Costs\r` 2.204e+01 2.150e+00 10.253 < 2e-16 ***
---
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/1xpby1291221210.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/2xpby1291221210.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/38ga11291221210.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/48ga11291221210.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/58ga11291221210.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/61p941291221210.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/7by871291221210.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/8by871291221210.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/9by871291221210.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/10mq8s1291221210.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/11pqoy1291221210.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/12brn41291221210.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/13hskx1291221210.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/14sjji1291221210.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/15e2zo1291221210.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/16abxx1291221210.tab")
+ }
>
> try(system("convert tmp/1xpby1291221210.ps tmp/1xpby1291221210.png",intern=TRUE))
character(0)
> try(system("convert tmp/2xpby1291221210.ps tmp/2xpby1291221210.png",intern=TRUE))
character(0)
> try(system("convert tmp/38ga11291221210.ps tmp/38ga11291221210.png",intern=TRUE))
character(0)
> try(system("convert tmp/48ga11291221210.ps tmp/48ga11291221210.png",intern=TRUE))
character(0)
> try(system("convert tmp/58ga11291221210.ps tmp/58ga11291221210.png",intern=TRUE))
character(0)
> try(system("convert tmp/61p941291221210.ps tmp/61p941291221210.png",intern=TRUE))
character(0)
> try(system("convert tmp/7by871291221210.ps tmp/7by871291221210.png",intern=TRUE))
character(0)
> try(system("convert tmp/8by871291221210.ps tmp/8by871291221210.png",intern=TRUE))
character(0)
> try(system("convert tmp/9by871291221210.ps tmp/9by871291221210.png",intern=TRUE))
character(0)
> try(system("convert tmp/10mq8s1291221210.ps tmp/10mq8s1291221210.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
12.309 2.246 29.216