R version 2.8.0 (2008-10-20)
Copyright (C) 2008 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
Natural language support but running in an English locale
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> x <- array(list(1
+ ,162556
+ ,1081
+ ,213118
+ ,6282154
+ ,1
+ ,29790
+ ,309
+ ,81767
+ ,4321023
+ ,1
+ ,87550
+ ,458
+ ,153198
+ ,4111912
+ ,0
+ ,84738
+ ,588
+ ,-26007
+ ,223193
+ ,1
+ ,54660
+ ,302
+ ,126942
+ ,1491348
+ ,1
+ ,42634
+ ,156
+ ,157214
+ ,1629616
+ ,0
+ ,40949
+ ,481
+ ,129352
+ ,1398893
+ ,1
+ ,45187
+ ,353
+ ,234817
+ ,1926517
+ ,1
+ ,37704
+ ,452
+ ,60448
+ ,983660
+ ,1
+ ,16275
+ ,109
+ ,47818
+ ,1443586
+ ,0
+ ,25830
+ ,115
+ ,245546
+ ,1073089
+ ,0
+ ,12679
+ ,110
+ ,48020
+ ,984885
+ ,1
+ ,18014
+ ,239
+ ,-1710
+ ,1405225
+ ,0
+ ,43556
+ ,247
+ ,32648
+ ,227132
+ ,1
+ ,24811
+ ,505
+ ,95350
+ ,929118
+ ,0
+ ,6575
+ ,159
+ ,151352
+ ,1071292
+ ,0
+ ,7123
+ ,109
+ ,288170
+ ,638830
+ ,1
+ ,21950
+ ,519
+ ,114337
+ ,856956
+ ,1
+ ,37597
+ ,248
+ ,37884
+ ,992426
+ ,0
+ ,17821
+ ,373
+ ,122844
+ ,444477
+ ,1
+ ,12988
+ ,119
+ ,82340
+ ,857217
+ ,1
+ ,22330
+ ,84
+ ,79801
+ ,711969
+ ,0
+ ,13326
+ ,102
+ ,165548
+ ,702380
+ ,0
+ ,16189
+ ,295
+ ,116384
+ ,358589
+ ,0
+ ,7146
+ ,105
+ ,134028
+ ,297978
+ ,0
+ ,15824
+ ,64
+ ,63838
+ ,585715
+ ,1
+ ,27664
+ ,282
+ ,74996
+ ,657954
+ ,0
+ ,11920
+ ,182
+ ,31080
+ ,209458
+ ,0
+ ,8568
+ ,37
+ ,32168
+ ,786690
+ ,0
+ ,14416
+ ,361
+ ,49857
+ ,439798
+ ,1
+ ,3369
+ ,28
+ ,87161
+ ,688779
+ ,1
+ ,11819
+ ,85
+ ,106113
+ ,574339
+ ,1
+ ,6984
+ ,45
+ ,80570
+ ,741409
+ ,1
+ ,4519
+ ,49
+ ,102129
+ ,597793
+ ,0
+ ,2220
+ ,22
+ ,301670
+ ,644190
+ ,0
+ ,18562
+ ,155
+ ,102313
+ ,377934
+ ,0
+ ,10327
+ ,91
+ ,88577
+ ,640273
+ ,1
+ ,5336
+ ,81
+ ,112477
+ ,697458
+ ,1
+ ,2365
+ ,79
+ ,191778
+ ,550608
+ ,0
+ ,4069
+ ,145
+ ,79804
+ ,207393
+ ,0
+ ,8636
+ ,855
+ ,128294
+ ,301607
+ ,0
+ ,13718
+ ,61
+ ,96448
+ ,345783
+ ,0
+ ,4525
+ ,226
+ ,93811
+ ,501749
+ ,0
+ ,6869
+ ,105
+ ,117520
+ ,379983
+ ,0
+ ,4628
+ ,62
+ ,69159
+ ,387475
+ ,1
+ ,3689
+ ,25
+ ,101792
+ ,377305
+ ,1
+ ,4891
+ ,217
+ ,210568
+ ,370837
+ ,1
+ ,7489
+ ,322
+ ,136996
+ ,430866
+ ,0
+ ,4901
+ ,84
+ ,121920
+ ,469107
+ ,0
+ ,2284
+ ,33
+ ,76403
+ ,194493
+ ,1
+ ,3160
+ ,108
+ ,108094
+ ,530670
+ ,1
+ ,4150
+ ,150
+ ,134759
+ ,518365
+ ,1
+ ,7285
+ ,115
+ ,188873
+ ,491303
+ ,1
+ ,1134
+ ,162
+ ,146216
+ ,527021
+ ,1
+ ,4658
+ ,158
+ ,156608
+ ,233773
+ ,0
+ ,2384
+ ,97
+ ,61348
+ ,405972
+ ,0
+ ,3748
+ ,9
+ ,50350
+ ,652925
+ ,0
+ ,5371
+ ,66
+ ,87720
+ ,446211
+ ,0
+ ,1285
+ ,107
+ ,99489
+ ,341340
+ ,1
+ ,9327
+ ,101
+ ,87419
+ ,387699
+ ,1
+ ,5565
+ ,47
+ ,94355
+ ,493408
+ ,0
+ ,1528
+ ,38
+ ,60326
+ ,146494
+ ,1
+ ,3122
+ ,34
+ ,94670
+ ,414462
+ ,1
+ ,7561
+ ,87
+ ,82425
+ ,364304
+ ,0
+ ,2675
+ ,79
+ ,59017
+ ,355178
+ ,0
+ ,13253
+ ,947
+ ,90829
+ ,357760
+ ,0
+ ,880
+ ,74
+ ,80791
+ ,261216
+ ,1
+ ,2053
+ ,53
+ ,100423
+ ,397144
+ ,0
+ ,1424
+ ,94
+ ,131116
+ ,374943
+ ,1
+ ,4036
+ ,63
+ ,100269
+ ,424898
+ ,1
+ ,3045
+ ,58
+ ,27330
+ ,202055
+ ,0
+ ,5119
+ ,49
+ ,39039
+ ,378525
+ ,0
+ ,1431
+ ,34
+ ,106885
+ ,310768
+ ,0
+ ,554
+ ,11
+ ,79285
+ ,325738
+ ,0
+ ,1975
+ ,35
+ ,118881
+ ,394510
+ ,1
+ ,1765
+ ,20
+ ,77623
+ ,247060
+ ,0
+ ,1012
+ ,47
+ ,114768
+ ,368078
+ ,0
+ ,810
+ ,43
+ ,74015
+ ,236761
+ ,0
+ ,1280
+ ,117
+ ,69465
+ ,312378
+ ,1
+ ,666
+ ,171
+ ,117869
+ ,339836
+ ,0
+ ,1380
+ ,26
+ ,60982
+ ,347385
+ ,1
+ ,4677
+ ,75
+ ,90131
+ ,426280
+ ,0
+ ,876
+ ,59
+ ,138971
+ ,352850
+ ,0
+ ,814
+ ,18
+ ,39625
+ ,301881
+ ,0
+ ,514
+ ,15
+ ,102725
+ ,377516
+ ,1
+ ,5692
+ ,72
+ ,64239
+ ,357312
+ ,0
+ ,3642
+ ,86
+ ,90262
+ ,458343
+ ,0
+ ,540
+ ,14
+ ,103960
+ ,354228
+ ,0
+ ,2099
+ ,64
+ ,106611
+ ,308636
+ ,0
+ ,567
+ ,11
+ ,103345
+ ,386212
+ ,0
+ ,2001
+ ,52
+ ,95551
+ ,393343
+ ,1
+ ,2949
+ ,41
+ ,82903
+ ,378509
+ ,0
+ ,2253
+ ,99
+ ,63593
+ ,452469
+ ,1
+ ,6533
+ ,75
+ ,126910
+ ,364839
+ ,0
+ ,1889
+ ,45
+ ,37527
+ ,358649
+ ,1
+ ,3055
+ ,43
+ ,60247
+ ,376641
+ ,0
+ ,272
+ ,8
+ ,112995
+ ,429112
+ ,1
+ ,1414
+ ,198
+ ,70184
+ ,330546
+ ,0
+ ,2564
+ ,22
+ ,130140
+ ,403560
+ ,1
+ ,1383
+ ,11
+ ,73221
+ ,317892
+ ,1
+ ,1261
+ ,33
+ ,76114
+ ,307528
+ ,0
+ ,975
+ ,23
+ ,90534
+ ,235133
+ ,0
+ ,3366
+ ,80
+ ,108479
+ ,299243
+ ,0
+ ,576
+ ,18
+ ,113761
+ ,314073
+ ,0
+ ,1686
+ ,40
+ ,68696
+ ,368186
+ ,0
+ ,746
+ ,23
+ ,71561
+ ,269661
+ ,1
+ ,3192
+ ,60
+ ,59831
+ ,125390
+ ,1
+ ,2045
+ ,20
+ ,97890
+ ,510834
+ ,0
+ ,5702
+ ,61
+ ,101481
+ ,321896
+ ,1
+ ,1932
+ ,36
+ ,72954
+ ,249898
+ ,0
+ ,936
+ ,30
+ ,67939
+ ,408881
+ ,1
+ ,3437
+ ,47
+ ,48022
+ ,158492
+ ,0
+ ,5131
+ ,71
+ ,86111
+ ,292154
+ ,1
+ ,2397
+ ,14
+ ,74020
+ ,289513
+ ,1
+ ,1389
+ ,9
+ ,57530
+ ,378049
+ ,0
+ ,1503
+ ,39
+ ,56364
+ ,343466
+ ,0
+ ,402
+ ,26
+ ,84990
+ ,332743
+ ,0
+ ,2239
+ ,21
+ ,88590
+ ,442882
+ ,1
+ ,2234
+ ,16
+ ,77200
+ ,214215
+ ,0
+ ,837
+ ,69
+ ,61262
+ ,315688
+ ,0
+ ,10579
+ ,92
+ ,110309
+ ,375195
+ ,0
+ ,875
+ ,14
+ ,67000
+ ,334280
+ ,0
+ ,1585
+ ,107
+ ,93099
+ ,355864
+ ,1
+ ,1659
+ ,29
+ ,107577
+ ,480382
+ ,1
+ ,2647
+ ,37
+ ,62920
+ ,353058
+ ,1
+ ,3294
+ ,23
+ ,75832
+ ,217193
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,94
+ ,7
+ ,60793
+ ,314533
+ ,0
+ ,422
+ ,28
+ ,57935
+ ,318056
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,34
+ ,8
+ ,60630
+ ,314353
+ ,0
+ ,1558
+ ,63
+ ,55637
+ ,369448
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,43
+ ,3
+ ,60887
+ ,312846
+ ,0
+ ,645
+ ,5
+ ,60720
+ ,312075
+ ,0
+ ,316
+ ,9
+ ,60505
+ ,315009
+ ,0
+ ,115
+ ,13
+ ,60945
+ ,318903
+ ,0
+ ,5
+ ,2
+ ,60720
+ ,314887
+ ,0
+ ,897
+ ,5
+ ,60720
+ ,314913
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,389
+ ,14
+ ,58990
+ ,325506
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1002
+ ,15
+ ,56750
+ ,298568
+ ,0
+ ,36
+ ,3
+ ,60894
+ ,315834
+ ,0
+ ,460
+ ,15
+ ,63346
+ ,329784
+ ,0
+ ,309
+ ,11
+ ,56535
+ ,312878
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,9
+ ,6
+ ,60835
+ ,314987
+ ,0
+ ,271
+ ,2
+ ,60720
+ ,325249
+ ,0
+ ,14
+ ,1
+ ,61016
+ ,315877
+ ,0
+ ,520
+ ,10
+ ,58650
+ ,291650
+ ,0
+ ,1766
+ ,73
+ ,60438
+ ,305959
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,458
+ ,11
+ ,58625
+ ,297765
+ ,0
+ ,20
+ ,3
+ ,60938
+ ,315245
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,98
+ ,2
+ ,61490
+ ,315236
+ ,0
+ ,405
+ ,7
+ ,60845
+ ,336425
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,483
+ ,27
+ ,60830
+ ,306268
+ ,0
+ ,454
+ ,51
+ ,63261
+ ,302187
+ ,0
+ ,47
+ ,3
+ ,60720
+ ,314882
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,757
+ ,19
+ ,45689
+ ,382712
+ ,0
+ ,4655
+ ,393
+ ,60720
+ ,341570
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,36
+ ,4
+ ,61564
+ ,312412
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,203
+ ,9
+ ,61938
+ ,309596
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,126
+ ,10
+ ,60951
+ ,315547
+ ,0
+ ,400
+ ,152
+ ,60720
+ ,313267
+ ,1
+ ,71
+ ,1
+ ,60745
+ ,316176
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,972
+ ,34
+ ,71642
+ ,359335
+ ,1
+ ,531
+ ,10
+ ,71641
+ ,330068
+ ,0
+ ,2461
+ ,57
+ ,55792
+ ,314289
+ ,1
+ ,378
+ ,52
+ ,71873
+ ,297413
+ ,1
+ ,23
+ ,5
+ ,62555
+ ,314806
+ ,1
+ ,638
+ ,14
+ ,60370
+ ,333210
+ ,1
+ ,2300
+ ,29
+ ,64873
+ ,352108
+ ,0
+ ,149
+ ,5
+ ,62041
+ ,313332
+ ,0
+ ,226
+ ,5
+ ,65745
+ ,291787
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,275
+ ,4
+ ,59500
+ ,318745
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,141
+ ,6
+ ,61630
+ ,315366
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,28
+ ,2
+ ,60890
+ ,315688
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,4980
+ ,91
+ ,113521
+ ,409642
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,472
+ ,20
+ ,80045
+ ,269587
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,203
+ ,27
+ ,50804
+ ,300962
+ ,1
+ ,496
+ ,17
+ ,87390
+ ,325479
+ ,1
+ ,10
+ ,2
+ ,61656
+ ,316155
+ ,1
+ ,63
+ ,4
+ ,65688
+ ,318574
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,1136
+ ,32
+ ,48522
+ ,343613
+ ,1
+ ,265
+ ,31
+ ,60720
+ ,306948
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,267
+ ,32
+ ,57640
+ ,330059
+ ,0
+ ,474
+ ,20
+ ,61977
+ ,288985
+ ,0
+ ,534
+ ,7
+ ,62620
+ ,304485
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,15
+ ,8
+ ,60831
+ ,315688
+ ,0
+ ,397
+ ,28
+ ,60646
+ ,317736
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1866
+ ,29
+ ,56225
+ ,322331
+ ,0
+ ,288
+ ,4
+ ,60510
+ ,296656
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,3
+ ,2
+ ,60698
+ ,315354
+ ,0
+ ,468
+ ,21
+ ,60720
+ ,312161
+ ,0
+ ,20
+ ,2
+ ,60805
+ ,315576
+ ,0
+ ,278
+ ,26
+ ,61404
+ ,314922
+ ,0
+ ,61
+ ,14
+ ,60720
+ ,314551
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,192
+ ,4
+ ,65276
+ ,312339
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,317
+ ,9
+ ,63915
+ ,298700
+ ,0
+ ,738
+ ,10
+ ,60720
+ ,321376
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,368
+ ,17
+ ,61686
+ ,303230
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,2
+ ,1
+ ,60743
+ ,315487
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,53
+ ,6
+ ,60349
+ ,315793
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,94
+ ,3
+ ,61360
+ ,312887
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,24
+ ,8
+ ,59818
+ ,315637
+ ,0
+ ,2332
+ ,4
+ ,72680
+ ,324385
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,131
+ ,11
+ ,61808
+ ,308989
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,206
+ ,9
+ ,53110
+ ,296702
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,167
+ ,2
+ ,64245
+ ,307322
+ ,0
+ ,622
+ ,73
+ ,73007
+ ,304376
+ ,0
+ ,2328
+ ,94
+ ,82732
+ ,253588
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,365
+ ,8
+ ,54820
+ ,309560
+ ,0
+ ,364
+ ,35
+ ,47705
+ ,298466
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,226
+ ,12
+ ,72835
+ ,343929
+ ,0
+ ,307
+ ,15
+ ,58856
+ ,331955
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,188
+ ,11
+ ,77655
+ ,381180
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,138
+ ,6
+ ,69817
+ ,331420
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,125
+ ,12
+ ,60798
+ ,310201
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,282
+ ,30
+ ,62452
+ ,320016
+ ,0
+ ,335
+ ,33
+ ,64175
+ ,320398
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,1324
+ ,117
+ ,67440
+ ,291841
+ ,0
+ ,176
+ ,28
+ ,68136
+ ,310670
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,249
+ ,72
+ ,56726
+ ,313491
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,333
+ ,13
+ ,70811
+ ,331323
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,601
+ ,6
+ ,60720
+ ,319210
+ ,0
+ ,30
+ ,4
+ ,62045
+ ,318098
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,249
+ ,62
+ ,54323
+ ,292754
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,165
+ ,24
+ ,62841
+ ,325176
+ ,0
+ ,453
+ ,21
+ ,81125
+ ,365959
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,53
+ ,14
+ ,59506
+ ,302409
+ ,1
+ ,382
+ ,21
+ ,59365
+ ,340968
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,30
+ ,4
+ ,60798
+ ,313164
+ ,0
+ ,290
+ ,2
+ ,58790
+ ,301164
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,366
+ ,53
+ ,61808
+ ,344425
+ ,0
+ ,2
+ ,9
+ ,60735
+ ,315394
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,209
+ ,13
+ ,64016
+ ,316647
+ ,0
+ ,384
+ ,22
+ ,54683
+ ,309836
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,365
+ ,83
+ ,87192
+ ,346611
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,49
+ ,8
+ ,64107
+ ,322031
+ ,0
+ ,3
+ ,4
+ ,60761
+ ,315656
+ ,0
+ ,133
+ ,14
+ ,65990
+ ,339445
+ ,0
+ ,32
+ ,1
+ ,59988
+ ,314964
+ ,0
+ ,368
+ ,17
+ ,61167
+ ,297141
+ ,0
+ ,1
+ ,6
+ ,60719
+ ,315372
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,22
+ ,2
+ ,60722
+ ,312502
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,96
+ ,5
+ ,60379
+ ,313729
+ ,0
+ ,1
+ ,2
+ ,60727
+ ,315388
+ ,0
+ ,314
+ ,5
+ ,60720
+ ,315371
+ ,0
+ ,844
+ ,78
+ ,60925
+ ,296139
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,26
+ ,1
+ ,60896
+ ,313880
+ ,0
+ ,125
+ ,13
+ ,59734
+ ,317698
+ ,0
+ ,304
+ ,15
+ ,62969
+ ,295580
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,621
+ ,48
+ ,60720
+ ,308256
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,119
+ ,6
+ ,59118
+ ,303677
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1595
+ ,17
+ ,60720
+ ,319369
+ ,0
+ ,312
+ ,14
+ ,58598
+ ,318690
+ ,0
+ ,60
+ ,10
+ ,61124
+ ,314049
+ ,0
+ ,587
+ ,12
+ ,59595
+ ,325699
+ ,0
+ ,135
+ ,2
+ ,62065
+ ,314210
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,514
+ ,52
+ ,78780
+ ,322378
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,4
+ ,60722
+ ,315398
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,1763
+ ,24
+ ,61600
+ ,308336
+ ,0
+ ,180
+ ,11
+ ,59635
+ ,316386
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,218
+ ,21
+ ,60720
+ ,315553
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,448
+ ,40
+ ,59781
+ ,323361
+ ,0
+ ,227
+ ,9
+ ,76644
+ ,336639
+ ,0
+ ,174
+ ,1
+ ,64820
+ ,307424
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,121
+ ,24
+ ,56178
+ ,295370
+ ,0
+ ,607
+ ,11
+ ,60436
+ ,322340
+ ,0
+ ,2212
+ ,14
+ ,60720
+ ,319864
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,530
+ ,60
+ ,73433
+ ,317291
+ ,0
+ ,571
+ ,80
+ ,41477
+ ,280398
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,78
+ ,16
+ ,62700
+ ,317330
+ ,0
+ ,2489
+ ,40
+ ,67804
+ ,238125
+ ,0
+ ,131
+ ,6
+ ,59661
+ ,327071
+ ,0
+ ,923
+ ,8
+ ,58620
+ ,309038
+ ,0
+ ,72
+ ,3
+ ,60398
+ ,314210
+ ,0
+ ,572
+ ,16
+ ,58580
+ ,307930
+ ,0
+ ,397
+ ,10
+ ,62710
+ ,322327
+ ,0
+ ,450
+ ,8
+ ,59325
+ ,292136
+ ,0
+ ,622
+ ,7
+ ,60950
+ ,263276
+ ,0
+ ,694
+ ,8
+ ,68060
+ ,367655
+ ,1
+ ,3425
+ ,12
+ ,83620
+ ,283910
+ ,0
+ ,562
+ ,13
+ ,58456
+ ,283587
+ ,0
+ ,4917
+ ,42
+ ,52811
+ ,243650
+ ,1
+ ,1442
+ ,118
+ ,121173
+ ,438493
+ ,0
+ ,529
+ ,9
+ ,63870
+ ,296261
+ ,1
+ ,2126
+ ,138
+ ,21001
+ ,230621
+ ,0
+ ,1061
+ ,5
+ ,70415
+ ,304252
+ ,0
+ ,776
+ ,9
+ ,64230
+ ,333505
+ ,0
+ ,611
+ ,8
+ ,59190
+ ,296919
+ ,1
+ ,1526
+ ,25
+ ,69351
+ ,278990
+ ,0
+ ,592
+ ,7
+ ,64270
+ ,276898
+ ,0
+ ,1182
+ ,13
+ ,70694
+ ,327007
+ ,0
+ ,621
+ ,16
+ ,68005
+ ,317046
+ ,0
+ ,989
+ ,11
+ ,58930
+ ,304555
+ ,0
+ ,438
+ ,11
+ ,58320
+ ,298096
+ ,0
+ ,726
+ ,3
+ ,69980
+ ,231861
+ ,0
+ ,1303
+ ,61
+ ,69863
+ ,309422
+ ,1
+ ,7419
+ ,29
+ ,63255
+ ,286963
+ ,1
+ ,1164
+ ,17
+ ,57320
+ ,269753
+ ,1
+ ,3310
+ ,33
+ ,75230
+ ,448243
+ ,0
+ ,1920
+ ,15
+ ,79420
+ ,165404
+ ,0
+ ,965
+ ,3
+ ,73490
+ ,204325
+ ,0
+ ,3256
+ ,66
+ ,35250
+ ,407159
+ ,1
+ ,1135
+ ,17
+ ,62285
+ ,290476
+ ,0
+ ,1270
+ ,26
+ ,69206
+ ,275311
+ ,0
+ ,661
+ ,3
+ ,65920
+ ,246541
+ ,0
+ ,1013
+ ,2
+ ,69770
+ ,253468
+ ,0
+ ,2844
+ ,67
+ ,72683
+ ,240897
+ ,1
+ ,11528
+ ,70
+ ,-14545
+ ,-83265
+ ,0
+ ,6526
+ ,26
+ ,55830
+ ,-42143
+ ,0
+ ,2264
+ ,24
+ ,55174
+ ,272713
+ ,1
+ ,5109
+ ,97
+ ,67038
+ ,215362
+ ,0
+ ,3999
+ ,30
+ ,51252
+ ,42754
+ ,0
+ ,35624
+ ,223
+ ,157278
+ ,306275
+ ,0
+ ,9252
+ ,48
+ ,79510
+ ,253537
+ ,0
+ ,15236
+ ,90
+ ,77440
+ ,372631
+ ,0
+ ,18073
+ ,180
+ ,27284
+ ,-7170)
+ ,dim=c(5
+ ,431)
+ ,dimnames=list(c('Group'
+ ,'Costs'
+ ,'Trades'
+ ,'Dividends'
+ ,'Wealth
')
+ ,1:431))
> y <- array(NA,dim=c(5,431),dimnames=list(c('Group','Costs','Trades','Dividends','Wealth
'),1:431))
> for (i in 1:dim(x)[1])
+ {
+ for (j in 1:dim(x)[2])
+ {
+ y[i,j] <- as.numeric(x[i,j])
+ }
+ }
> par3 = 'Linear Trend'
> par2 = 'Do not include Seasonal Dummies'
> par1 = '5'
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Prof. Dr. P. Wessa
> #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #Technical description: Write here your technical program description (don't use hard returns!)
> library(lattice)
> library(lmtest)
Loading required package: zoo
Attaching package: 'zoo'
The following object(s) are masked from package:base :
as.Date.numeric
> n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test
> par1 <- as.numeric(par1)
> x <- t(y)
> k <- length(x[1,])
> n <- length(x[,1])
> x1 <- cbind(x[,par1], x[,1:k!=par1])
> mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1])
> colnames(x1) <- mycolnames #colnames(x)[par1]
> x <- x1
> if (par3 == 'First Differences'){
+ x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep='')))
+ for (i in 1:n-1) {
+ for (j in 1:k) {
+ x2[i,j] <- x[i+1,j] - x[i,j]
+ }
+ }
+ x <- x2
+ }
> if (par2 == 'Include Monthly Dummies'){
+ x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep ='')))
+ for (i in 1:11){
+ x2[seq(i,n,12),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> if (par2 == 'Include Quarterly Dummies'){
+ x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep ='')))
+ for (i in 1:3){
+ x2[seq(i,n,4),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> k <- length(x[1,])
> if (par3 == 'Linear Trend'){
+ x <- cbind(x, c(1:n))
+ colnames(x)[k+1] <- 't'
+ }
> x
Wealth\r Group Costs Trades Dividends t
1 6282154 1 162556 1081 213118 1
2 4321023 1 29790 309 81767 2
3 4111912 1 87550 458 153198 3
4 223193 0 84738 588 -26007 4
5 1491348 1 54660 302 126942 5
6 1629616 1 42634 156 157214 6
7 1398893 0 40949 481 129352 7
8 1926517 1 45187 353 234817 8
9 983660 1 37704 452 60448 9
10 1443586 1 16275 109 47818 10
11 1073089 0 25830 115 245546 11
12 984885 0 12679 110 48020 12
13 1405225 1 18014 239 -1710 13
14 227132 0 43556 247 32648 14
15 929118 1 24811 505 95350 15
16 1071292 0 6575 159 151352 16
17 638830 0 7123 109 288170 17
18 856956 1 21950 519 114337 18
19 992426 1 37597 248 37884 19
20 444477 0 17821 373 122844 20
21 857217 1 12988 119 82340 21
22 711969 1 22330 84 79801 22
23 702380 0 13326 102 165548 23
24 358589 0 16189 295 116384 24
25 297978 0 7146 105 134028 25
26 585715 0 15824 64 63838 26
27 657954 1 27664 282 74996 27
28 209458 0 11920 182 31080 28
29 786690 0 8568 37 32168 29
30 439798 0 14416 361 49857 30
31 688779 1 3369 28 87161 31
32 574339 1 11819 85 106113 32
33 741409 1 6984 45 80570 33
34 597793 1 4519 49 102129 34
35 644190 0 2220 22 301670 35
36 377934 0 18562 155 102313 36
37 640273 0 10327 91 88577 37
38 697458 1 5336 81 112477 38
39 550608 1 2365 79 191778 39
40 207393 0 4069 145 79804 40
41 301607 0 8636 855 128294 41
42 345783 0 13718 61 96448 42
43 501749 0 4525 226 93811 43
44 379983 0 6869 105 117520 44
45 387475 0 4628 62 69159 45
46 377305 1 3689 25 101792 46
47 370837 1 4891 217 210568 47
48 430866 1 7489 322 136996 48
49 469107 0 4901 84 121920 49
50 194493 0 2284 33 76403 50
51 530670 1 3160 108 108094 51
52 518365 1 4150 150 134759 52
53 491303 1 7285 115 188873 53
54 527021 1 1134 162 146216 54
55 233773 1 4658 158 156608 55
56 405972 0 2384 97 61348 56
57 652925 0 3748 9 50350 57
58 446211 0 5371 66 87720 58
59 341340 0 1285 107 99489 59
60 387699 1 9327 101 87419 60
61 493408 1 5565 47 94355 61
62 146494 0 1528 38 60326 62
63 414462 1 3122 34 94670 63
64 364304 1 7561 87 82425 64
65 355178 0 2675 79 59017 65
66 357760 0 13253 947 90829 66
67 261216 0 880 74 80791 67
68 397144 1 2053 53 100423 68
69 374943 0 1424 94 131116 69
70 424898 1 4036 63 100269 70
71 202055 1 3045 58 27330 71
72 378525 0 5119 49 39039 72
73 310768 0 1431 34 106885 73
74 325738 0 554 11 79285 74
75 394510 0 1975 35 118881 75
76 247060 1 1765 20 77623 76
77 368078 0 1012 47 114768 77
78 236761 0 810 43 74015 78
79 312378 0 1280 117 69465 79
80 339836 1 666 171 117869 80
81 347385 0 1380 26 60982 81
82 426280 1 4677 75 90131 82
83 352850 0 876 59 138971 83
84 301881 0 814 18 39625 84
85 377516 0 514 15 102725 85
86 357312 1 5692 72 64239 86
87 458343 0 3642 86 90262 87
88 354228 0 540 14 103960 88
89 308636 0 2099 64 106611 89
90 386212 0 567 11 103345 90
91 393343 0 2001 52 95551 91
92 378509 1 2949 41 82903 92
93 452469 0 2253 99 63593 93
94 364839 1 6533 75 126910 94
95 358649 0 1889 45 37527 95
96 376641 1 3055 43 60247 96
97 429112 0 272 8 112995 97
98 330546 1 1414 198 70184 98
99 403560 0 2564 22 130140 99
100 317892 1 1383 11 73221 100
101 307528 1 1261 33 76114 101
102 235133 0 975 23 90534 102
103 299243 0 3366 80 108479 103
104 314073 0 576 18 113761 104
105 368186 0 1686 40 68696 105
106 269661 0 746 23 71561 106
107 125390 1 3192 60 59831 107
108 510834 1 2045 20 97890 108
109 321896 0 5702 61 101481 109
110 249898 1 1932 36 72954 110
111 408881 0 936 30 67939 111
112 158492 1 3437 47 48022 112
113 292154 0 5131 71 86111 113
114 289513 1 2397 14 74020 114
115 378049 1 1389 9 57530 115
116 343466 0 1503 39 56364 116
117 332743 0 402 26 84990 117
118 442882 0 2239 21 88590 118
119 214215 1 2234 16 77200 119
120 315688 0 837 69 61262 120
121 375195 0 10579 92 110309 121
122 334280 0 875 14 67000 122
123 355864 0 1585 107 93099 123
124 480382 1 1659 29 107577 124
125 353058 1 2647 37 62920 125
126 217193 1 3294 23 75832 126
127 315380 0 0 0 60720 127
128 314533 0 94 7 60793 128
129 318056 0 422 28 57935 129
130 315380 0 0 0 60720 130
131 314353 0 34 8 60630 131
132 369448 0 1558 63 55637 132
133 315380 0 0 0 60720 133
134 312846 0 43 3 60887 134
135 312075 0 645 5 60720 135
136 315009 0 316 9 60505 136
137 318903 0 115 13 60945 137
138 314887 0 5 2 60720 138
139 314913 0 897 5 60720 139
140 315380 0 0 0 60720 140
141 325506 0 389 14 58990 141
142 315380 0 0 0 60720 142
143 298568 0 1002 15 56750 143
144 315834 0 36 3 60894 144
145 329784 0 460 15 63346 145
146 312878 0 309 11 56535 146
147 315380 0 0 0 60720 147
148 314987 0 9 6 60835 148
149 325249 0 271 2 60720 149
150 315877 0 14 1 61016 150
151 291650 0 520 10 58650 151
152 305959 0 1766 73 60438 152
153 315380 0 0 0 60720 153
154 297765 0 458 11 58625 154
155 315245 0 20 3 60938 155
156 315380 0 0 0 60720 156
157 315380 0 0 0 60720 157
158 315236 0 98 2 61490 158
159 336425 0 405 7 60845 159
160 315380 0 0 0 60720 160
161 315380 0 0 0 60720 161
162 315380 0 0 0 60720 162
163 315380 0 0 0 60720 163
164 306268 0 483 27 60830 164
165 302187 0 454 51 63261 165
166 314882 0 47 3 60720 166
167 315380 0 0 0 60720 167
168 382712 0 757 19 45689 168
169 341570 0 4655 393 60720 169
170 315380 0 0 0 60720 170
171 315380 0 0 0 60720 171
172 312412 0 36 4 61564 172
173 315380 0 0 0 60720 173
174 309596 0 203 9 61938 174
175 315380 0 0 0 60720 175
176 315547 0 126 10 60951 176
177 313267 0 400 152 60720 177
178 316176 1 71 1 60745 178
179 315380 0 0 0 60720 179
180 315380 0 0 0 60720 180
181 359335 0 972 34 71642 181
182 330068 1 531 10 71641 182
183 314289 0 2461 57 55792 183
184 297413 1 378 52 71873 184
185 314806 1 23 5 62555 185
186 333210 1 638 14 60370 186
187 352108 1 2300 29 64873 187
188 313332 0 149 5 62041 188
189 291787 0 226 5 65745 189
190 315380 0 0 0 60720 190
191 318745 0 275 4 59500 191
192 315380 0 0 0 60720 192
193 315366 0 141 6 61630 193
194 315380 0 0 0 60720 194
195 315688 0 28 2 60890 195
196 315380 1 0 0 60720 196
197 409642 1 4980 91 113521 197
198 315380 1 0 0 60720 198
199 315380 1 0 0 60720 199
200 269587 1 472 20 80045 200
201 315380 1 0 0 60720 201
202 315380 1 0 0 60720 202
203 315380 1 0 0 60720 203
204 300962 1 203 27 50804 204
205 325479 1 496 17 87390 205
206 316155 1 10 2 61656 206
207 318574 1 63 4 65688 207
208 315380 1 0 0 60720 208
209 343613 1 1136 32 48522 209
210 306948 1 265 31 60720 210
211 315380 0 0 0 60720 211
212 315380 0 0 0 60720 212
213 330059 0 267 32 57640 213
214 288985 0 474 20 61977 214
215 304485 0 534 7 62620 215
216 315380 0 0 0 60720 216
217 315688 0 15 8 60831 217
218 317736 0 397 28 60646 218
219 315380 0 0 0 60720 219
220 322331 0 1866 29 56225 220
221 296656 0 288 4 60510 221
222 315380 0 0 0 60720 222
223 315354 0 3 2 60698 223
224 312161 0 468 21 60720 224
225 315576 0 20 2 60805 225
226 314922 0 278 26 61404 226
227 314551 0 61 14 60720 227
228 315380 0 0 0 60720 228
229 312339 0 192 4 65276 229
230 315380 0 0 0 60720 230
231 298700 0 317 9 63915 231
232 321376 0 738 10 60720 232
233 315380 0 0 0 60720 233
234 303230 0 368 17 61686 234
235 315380 0 0 0 60720 235
236 315487 0 2 1 60743 236
237 315380 0 0 0 60720 237
238 315793 0 53 6 60349 238
239 315380 0 0 0 60720 239
240 315380 0 0 0 60720 240
241 315380 0 0 0 60720 241
242 312887 0 94 3 61360 242
243 315380 0 0 0 60720 243
244 315637 0 24 8 59818 244
245 324385 0 2332 4 72680 245
246 315380 1 0 0 60720 246
247 315380 1 0 0 60720 247
248 308989 0 131 11 61808 248
249 315380 1 0 0 60720 249
250 315380 1 0 0 60720 250
251 296702 0 206 9 53110 251
252 315380 1 0 0 60720 252
253 307322 0 167 2 64245 253
254 304376 0 622 73 73007 254
255 253588 0 2328 94 82732 255
256 315380 0 0 0 60720 256
257 309560 0 365 8 54820 257
258 298466 0 364 35 47705 258
259 315380 1 0 0 60720 259
260 315380 0 0 0 60720 260
261 315380 1 0 0 60720 261
262 315380 1 0 0 60720 262
263 343929 0 226 12 72835 263
264 331955 0 307 15 58856 264
265 315380 1 0 0 60720 265
266 315380 0 0 0 60720 266
267 315380 1 0 0 60720 267
268 381180 0 188 11 77655 268
269 315380 1 0 0 60720 269
270 331420 0 138 6 69817 270
271 315380 1 0 0 60720 271
272 315380 1 0 0 60720 272
273 315380 1 0 0 60720 273
274 310201 0 125 12 60798 274
275 315380 0 0 0 60720 275
276 320016 0 282 30 62452 276
277 320398 0 335 33 64175 277
278 315380 0 0 0 60720 278
279 291841 1 1324 117 67440 279
280 310670 0 176 28 68136 280
281 315380 1 0 0 60720 281
282 315380 0 0 0 60720 282
283 313491 0 249 72 56726 283
284 315380 0 0 0 60720 284
285 331323 0 333 13 70811 285
286 315380 0 0 0 60720 286
287 319210 1 601 6 60720 287
288 318098 0 30 4 62045 288
289 315380 1 0 0 60720 289
290 292754 0 249 62 54323 290
291 315380 1 0 0 60720 291
292 325176 0 165 24 62841 292
293 365959 0 453 21 81125 293
294 315380 0 0 0 60720 294
295 302409 0 53 14 59506 295
296 340968 1 382 21 59365 296
297 315380 0 0 0 60720 297
298 315380 1 0 0 60720 298
299 315380 0 0 0 60720 299
300 315380 0 0 0 60720 300
301 313164 1 30 4 60798 301
302 301164 0 290 2 58790 302
303 315380 0 0 0 60720 303
304 315380 1 0 0 60720 304
305 344425 0 366 53 61808 305
306 315394 0 2 9 60735 306
307 315380 1 0 0 60720 307
308 316647 1 209 13 64016 308
309 309836 0 384 22 54683 309
310 315380 1 0 0 60720 310
311 315380 1 0 0 60720 311
312 346611 0 365 83 87192 312
313 315380 1 0 0 60720 313
314 322031 1 49 8 64107 314
315 315656 0 3 4 60761 315
316 339445 0 133 14 65990 316
317 314964 0 32 1 59988 317
318 297141 0 368 17 61167 318
319 315372 0 1 6 60719 319
320 315380 0 0 0 60720 320
321 315380 0 0 0 60720 321
322 315380 0 0 0 60720 322
323 315380 0 0 0 60720 323
324 315380 1 0 0 60720 324
325 315380 0 0 0 60720 325
326 312502 0 22 2 60722 326
327 315380 1 0 0 60720 327
328 315380 1 0 0 60720 328
329 315380 0 0 0 60720 329
330 315380 0 0 0 60720 330
331 315380 0 0 0 60720 331
332 315380 0 0 0 60720 332
333 315380 0 0 0 60720 333
334 313729 0 96 5 60379 334
335 315388 0 1 2 60727 335
336 315371 0 314 5 60720 336
337 296139 0 844 78 60925 337
338 315380 0 0 0 60720 338
339 313880 0 26 1 60896 339
340 317698 0 125 13 59734 340
341 295580 0 304 15 62969 341
342 315380 0 0 0 60720 342
343 315380 0 0 0 60720 343
344 315380 0 0 0 60720 344
345 308256 0 621 48 60720 345
346 315380 0 0 0 60720 346
347 303677 0 119 6 59118 347
348 315380 0 0 0 60720 348
349 315380 0 0 0 60720 349
350 319369 0 1595 17 60720 350
351 318690 0 312 14 58598 351
352 314049 0 60 10 61124 352
353 325699 0 587 12 59595 353
354 314210 0 135 2 62065 354
355 315380 0 0 0 60720 355
356 315380 0 0 0 60720 356
357 322378 0 514 52 78780 357
358 315380 0 0 0 60720 358
359 315380 0 0 0 60720 359
360 315380 0 0 0 60720 360
361 315398 0 1 4 60722 361
362 315380 0 0 0 60720 362
363 315380 0 0 0 60720 363
364 308336 1 1763 24 61600 364
365 316386 0 180 11 59635 365
366 315380 0 0 0 60720 366
367 315380 0 0 0 60720 367
368 315380 0 0 0 60720 368
369 315380 0 0 0 60720 369
370 315553 0 218 21 60720 370
371 315380 0 0 0 60720 371
372 323361 0 448 40 59781 372
373 336639 0 227 9 76644 373
374 307424 0 174 1 64820 374
375 315380 0 0 0 60720 375
376 315380 0 0 0 60720 376
377 295370 0 121 24 56178 377
378 322340 0 607 11 60436 378
379 319864 0 2212 14 60720 379
380 315380 0 0 0 60720 380
381 315380 0 0 0 60720 381
382 317291 0 530 60 73433 382
383 280398 0 571 80 41477 383
384 315380 0 0 0 60720 384
385 317330 0 78 16 62700 385
386 238125 0 2489 40 67804 386
387 327071 0 131 6 59661 387
388 309038 0 923 8 58620 388
389 314210 0 72 3 60398 389
390 307930 0 572 16 58580 390
391 322327 0 397 10 62710 391
392 292136 0 450 8 59325 392
393 263276 0 622 7 60950 393
394 367655 0 694 8 68060 394
395 283910 1 3425 12 83620 395
396 283587 0 562 13 58456 396
397 243650 0 4917 42 52811 397
398 438493 1 1442 118 121173 398
399 296261 0 529 9 63870 399
400 230621 1 2126 138 21001 400
401 304252 0 1061 5 70415 401
402 333505 0 776 9 64230 402
403 296919 0 611 8 59190 403
404 278990 1 1526 25 69351 404
405 276898 0 592 7 64270 405
406 327007 0 1182 13 70694 406
407 317046 0 621 16 68005 407
408 304555 0 989 11 58930 408
409 298096 0 438 11 58320 409
410 231861 0 726 3 69980 410
411 309422 0 1303 61 69863 411
412 286963 1 7419 29 63255 412
413 269753 1 1164 17 57320 413
414 448243 1 3310 33 75230 414
415 165404 0 1920 15 79420 415
416 204325 0 965 3 73490 416
417 407159 0 3256 66 35250 417
418 290476 1 1135 17 62285 418
419 275311 0 1270 26 69206 419
420 246541 0 661 3 65920 420
421 253468 0 1013 2 69770 421
422 240897 0 2844 67 72683 422
423 -83265 1 11528 70 -14545 423
424 -42143 0 6526 26 55830 424
425 272713 0 2264 24 55174 425
426 215362 1 5109 97 67038 426
427 42754 0 3999 30 51252 427
428 306275 0 35624 223 157278 428
429 253537 0 9252 48 79510 429
430 372631 0 15236 90 77440 430
431 -7170 0 18073 180 27284 431
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) Group Costs Trades Dividends t
175770.913 53021.362 29.801 -431.749 1.838 -78.375
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-2175875 -15297 35999 51627 3187769
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.758e+05 4.820e+04 3.647 0.000299 ***
Group 5.302e+04 2.809e+04 1.888 0.059763 .
Costs 2.980e+01 1.579e+00 18.870 < 2e-16 ***
Trades -4.317e+02 1.703e+02 -2.536 0.011575 *
Dividends 1.838e+00 4.307e-01 4.267 2.45e-05 ***
t -7.837e+01 1.145e+02 -0.684 0.494218
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 253500 on 425 degrees of freedom
Multiple R-squared: 0.6573, Adjusted R-squared: 0.6532
F-statistic: 163 on 5 and 425 DF, p-value: < 2.2e-16
> if (n > n25) {
+ kp3 <- k + 3
+ nmkm3 <- n - k - 3
+ gqarr <- array(NA, dim=c(nmkm3-kp3+1,3))
+ numgqtests <- 0
+ numsignificant1 <- 0
+ numsignificant5 <- 0
+ numsignificant10 <- 0
+ for (mypoint in kp3:nmkm3) {
+ j <- 0
+ numgqtests <- numgqtests + 1
+ for (myalt in c('greater', 'two.sided', 'less')) {
+ j <- j + 1
+ gqarr[mypoint-kp3+1,j] <- gqtest(mylm, point=mypoint, alternative=myalt)$p.value
+ }
+ if (gqarr[mypoint-kp3+1,2] < 0.01) numsignificant1 <- numsignificant1 + 1
+ if (gqarr[mypoint-kp3+1,2] < 0.05) numsignificant5 <- numsignificant5 + 1
+ if (gqarr[mypoint-kp3+1,2] < 0.10) numsignificant10 <- numsignificant10 + 1
+ }
+ gqarr
+ }
[,1] [,2] [,3]
[1,] 1.0000000 1.515478e-18 7.577388e-19
[2,] 1.0000000 4.419952e-64 2.209976e-64
[3,] 1.0000000 4.935590e-79 2.467795e-79
[4,] 1.0000000 9.104659e-109 4.552329e-109
[5,] 1.0000000 1.680338e-129 8.401691e-130
[6,] 1.0000000 3.121883e-133 1.560941e-133
[7,] 1.0000000 3.217971e-136 1.608985e-136
[8,] 1.0000000 3.820751e-151 1.910375e-151
[9,] 1.0000000 4.497188e-150 2.248594e-150
[10,] 1.0000000 9.318449e-153 4.659224e-153
[11,] 1.0000000 1.161379e-158 5.806893e-159
[12,] 1.0000000 6.858018e-159 3.429009e-159
[13,] 1.0000000 1.208447e-164 6.042236e-165
[14,] 1.0000000 6.215520e-166 3.107760e-166
[15,] 1.0000000 1.817677e-169 9.088387e-170
[16,] 1.0000000 1.297091e-169 6.485457e-170
[17,] 1.0000000 1.694399e-170 8.471993e-171
[18,] 1.0000000 1.406974e-174 7.034868e-175
[19,] 1.0000000 1.561645e-176 7.808225e-177
[20,] 1.0000000 5.798162e-178 2.899081e-178
[21,] 1.0000000 1.993953e-192 9.969767e-193
[22,] 1.0000000 1.315769e-193 6.578846e-194
[23,] 1.0000000 4.356192e-197 2.178096e-197
[24,] 1.0000000 6.875675e-198 3.437838e-198
[25,] 1.0000000 5.143499e-205 2.571749e-205
[26,] 1.0000000 1.611490e-206 8.057452e-207
[27,] 1.0000000 2.024035e-205 1.012018e-205
[28,] 1.0000000 1.343738e-204 6.718692e-205
[29,] 1.0000000 4.749207e-211 2.374604e-211
[30,] 1.0000000 1.690935e-216 8.454676e-217
[31,] 1.0000000 8.165241e-216 4.082621e-216
[32,] 1.0000000 1.032994e-216 5.164970e-217
[33,] 1.0000000 1.818056e-216 9.090279e-217
[34,] 1.0000000 1.311054e-215 6.555271e-216
[35,] 1.0000000 3.547422e-217 1.773711e-217
[36,] 1.0000000 2.262280e-216 1.131140e-216
[37,] 1.0000000 1.882959e-216 9.414796e-217
[38,] 1.0000000 1.865443e-215 9.327213e-216
[39,] 1.0000000 5.660137e-217 2.830068e-217
[40,] 1.0000000 5.208572e-216 2.604286e-216
[41,] 1.0000000 5.039992e-216 2.519996e-216
[42,] 1.0000000 4.853354e-217 2.426677e-217
[43,] 1.0000000 1.602300e-217 8.011499e-218
[44,] 1.0000000 3.142225e-217 1.571112e-217
[45,] 1.0000000 1.847783e-216 9.238914e-217
[46,] 1.0000000 5.500221e-216 2.750110e-216
[47,] 1.0000000 1.261439e-218 6.307196e-219
[48,] 1.0000000 1.472648e-219 7.363239e-220
[49,] 1.0000000 3.768798e-232 1.884399e-232
[50,] 1.0000000 1.082690e-232 5.413451e-233
[51,] 1.0000000 3.711894e-232 1.855947e-232
[52,] 1.0000000 1.143142e-231 5.715709e-232
[53,] 1.0000000 8.744860e-233 4.372430e-233
[54,] 1.0000000 8.764647e-235 4.382323e-235
[55,] 1.0000000 5.489724e-234 2.744862e-234
[56,] 1.0000000 3.315814e-233 1.657907e-233
[57,] 1.0000000 6.575947e-233 3.287974e-233
[58,] 1.0000000 4.504210e-232 2.252105e-232
[59,] 1.0000000 4.223847e-232 2.111924e-232
[60,] 1.0000000 4.861336e-231 2.430668e-231
[61,] 1.0000000 1.962121e-230 9.810604e-231
[62,] 1.0000000 8.842435e-230 4.421218e-230
[63,] 1.0000000 5.884206e-229 2.942103e-229
[64,] 1.0000000 7.544807e-230 3.772404e-230
[65,] 1.0000000 1.554422e-229 7.772108e-230
[66,] 1.0000000 8.834832e-229 4.417416e-229
[67,] 1.0000000 7.675549e-228 3.837774e-228
[68,] 1.0000000 1.076490e-227 5.382448e-228
[69,] 1.0000000 6.193443e-227 3.096721e-227
[70,] 1.0000000 3.306115e-227 1.653058e-227
[71,] 1.0000000 2.062552e-226 1.031276e-226
[72,] 1.0000000 3.482571e-226 1.741286e-226
[73,] 1.0000000 1.898070e-225 9.490351e-226
[74,] 1.0000000 4.012425e-225 2.006213e-225
[75,] 1.0000000 2.208483e-225 1.104242e-225
[76,] 1.0000000 1.344902e-224 6.724510e-225
[77,] 1.0000000 1.083559e-223 5.417793e-224
[78,] 1.0000000 3.431825e-223 1.715913e-223
[79,] 1.0000000 1.620583e-223 8.102917e-224
[80,] 1.0000000 9.558939e-223 4.779470e-223
[81,] 1.0000000 1.350018e-222 6.750091e-223
[82,] 1.0000000 1.218524e-221 6.092618e-222
[83,] 1.0000000 1.130777e-220 5.653884e-221
[84,] 1.0000000 8.945680e-220 4.472840e-220
[85,] 1.0000000 3.543613e-221 1.771806e-221
[86,] 1.0000000 1.740925e-220 8.704624e-221
[87,] 1.0000000 9.237789e-221 4.618895e-221
[88,] 1.0000000 1.509665e-220 7.548325e-221
[89,] 1.0000000 1.426899e-219 7.134495e-220
[90,] 1.0000000 1.709761e-218 8.548805e-219
[91,] 1.0000000 1.347585e-217 6.737925e-218
[92,] 1.0000000 1.457848e-216 7.289241e-217
[93,] 1.0000000 1.225221e-215 6.126103e-216
[94,] 1.0000000 4.867287e-217 2.433643e-217
[95,] 1.0000000 2.612817e-217 1.306409e-217
[96,] 1.0000000 2.196936e-218 1.098468e-218
[97,] 1.0000000 1.293048e-217 6.465242e-218
[98,] 1.0000000 2.754812e-217 1.377406e-217
[99,] 1.0000000 2.822821e-220 1.411410e-220
[100,] 1.0000000 3.184031e-221 1.592015e-221
[101,] 1.0000000 1.683356e-220 8.416779e-221
[102,] 1.0000000 1.590035e-220 7.950175e-221
[103,] 1.0000000 2.031209e-220 1.015604e-220
[104,] 1.0000000 3.355586e-221 1.677793e-221
[105,] 1.0000000 1.581953e-220 7.909765e-221
[106,] 1.0000000 8.652811e-220 4.326406e-220
[107,] 1.0000000 1.638948e-219 8.194741e-220
[108,] 1.0000000 9.051606e-219 4.525803e-219
[109,] 1.0000000 4.368337e-218 2.184169e-218
[110,] 1.0000000 7.049943e-218 3.524971e-218
[111,] 1.0000000 1.475061e-219 7.375304e-220
[112,] 1.0000000 1.445986e-218 7.229931e-219
[113,] 1.0000000 1.334807e-217 6.674035e-218
[114,] 1.0000000 1.379293e-216 6.896466e-217
[115,] 1.0000000 8.159952e-216 4.079976e-216
[116,] 1.0000000 2.805980e-215 1.402990e-215
[117,] 1.0000000 1.621156e-214 8.105779e-215
[118,] 1.0000000 1.033328e-215 5.166640e-216
[119,] 1.0000000 1.005872e-214 5.029362e-215
[120,] 1.0000000 9.861832e-214 4.930916e-214
[121,] 1.0000000 9.986274e-213 4.993137e-213
[122,] 1.0000000 1.000083e-211 5.000417e-212
[123,] 1.0000000 9.970742e-211 4.985371e-211
[124,] 1.0000000 1.941811e-210 9.709054e-211
[125,] 1.0000000 2.007670e-209 1.003835e-209
[126,] 1.0000000 2.051664e-208 1.025832e-208
[127,] 1.0000000 2.227086e-207 1.113543e-207
[128,] 1.0000000 2.372285e-206 1.186142e-206
[129,] 1.0000000 2.501003e-205 1.250502e-205
[130,] 1.0000000 2.606333e-204 1.303166e-204
[131,] 1.0000000 2.893069e-203 1.446534e-203
[132,] 1.0000000 3.029840e-202 1.514920e-202
[133,] 1.0000000 3.152845e-201 1.576422e-201
[134,] 1.0000000 3.308396e-200 1.654198e-200
[135,] 1.0000000 3.584480e-199 1.792240e-199
[136,] 1.0000000 3.753738e-198 1.876869e-198
[137,] 1.0000000 4.017008e-197 2.008504e-197
[138,] 1.0000000 4.294421e-196 2.147211e-196
[139,] 1.0000000 4.483868e-195 2.241934e-195
[140,] 1.0000000 4.636728e-194 2.318364e-194
[141,] 1.0000000 4.902323e-193 2.451162e-193
[142,] 1.0000000 5.089050e-192 2.544525e-192
[143,] 1.0000000 4.536544e-191 2.268272e-191
[144,] 1.0000000 4.810067e-190 2.405033e-190
[145,] 1.0000000 4.943833e-189 2.471917e-189
[146,] 1.0000000 4.693623e-188 2.346812e-188
[147,] 1.0000000 4.768557e-187 2.384279e-187
[148,] 1.0000000 4.850010e-186 2.425005e-186
[149,] 1.0000000 4.916451e-185 2.458226e-185
[150,] 1.0000000 4.938567e-184 2.469284e-184
[151,] 1.0000000 4.691499e-183 2.345749e-183
[152,] 1.0000000 4.718413e-182 2.359206e-182
[153,] 1.0000000 4.727624e-181 2.363812e-181
[154,] 1.0000000 4.718569e-180 2.359285e-180
[155,] 1.0000000 4.690923e-179 2.345461e-179
[156,] 1.0000000 4.393029e-178 2.196514e-178
[157,] 1.0000000 3.414431e-177 1.707215e-177
[158,] 1.0000000 3.341140e-176 1.670570e-176
[159,] 1.0000000 3.260951e-175 1.630475e-175
[160,] 1.0000000 1.036608e-175 5.183042e-176
[161,] 1.0000000 5.933276e-175 2.966638e-175
[162,] 1.0000000 5.881166e-174 2.940583e-174
[163,] 1.0000000 5.803678e-173 2.901839e-173
[164,] 1.0000000 5.532135e-172 2.766068e-172
[165,] 1.0000000 5.408149e-171 2.704074e-171
[166,] 1.0000000 5.001847e-170 2.500923e-170
[167,] 1.0000000 4.842125e-169 2.421062e-169
[168,] 1.0000000 4.683089e-168 2.341544e-168
[169,] 1.0000000 4.458892e-167 2.229446e-167
[170,] 1.0000000 4.033280e-166 2.016640e-166
[171,] 1.0000000 3.835111e-165 1.917555e-165
[172,] 1.0000000 3.629539e-164 1.814769e-164
[173,] 1.0000000 3.132239e-163 1.566119e-163
[174,] 1.0000000 2.565970e-162 1.282985e-162
[175,] 1.0000000 2.019233e-161 1.009617e-161
[176,] 1.0000000 7.887108e-161 3.943554e-161
[177,] 1.0000000 6.739485e-160 3.369742e-160
[178,] 1.0000000 5.611098e-159 2.805549e-159
[179,] 1.0000000 3.006365e-158 1.503183e-158
[180,] 1.0000000 2.722168e-157 1.361084e-157
[181,] 1.0000000 1.485235e-156 7.426177e-157
[182,] 1.0000000 1.358335e-155 6.791673e-156
[183,] 1.0000000 1.256463e-154 6.282317e-155
[184,] 1.0000000 1.139527e-153 5.697634e-154
[185,] 1.0000000 1.024080e-152 5.120398e-153
[186,] 1.0000000 9.201480e-152 4.600740e-152
[187,] 1.0000000 8.234211e-151 4.117106e-151
[188,] 1.0000000 7.049946e-150 3.524973e-150
[189,] 1.0000000 4.862527e-149 2.431264e-149
[190,] 1.0000000 4.177865e-148 2.088933e-148
[191,] 1.0000000 3.580136e-147 1.790068e-147
[192,] 1.0000000 1.421982e-147 7.109909e-148
[193,] 1.0000000 1.216127e-146 6.080634e-147
[194,] 1.0000000 1.036096e-145 5.180481e-146
[195,] 1.0000000 8.791750e-145 4.395875e-145
[196,] 1.0000000 7.669867e-144 3.834933e-144
[197,] 1.0000000 1.804365e-143 9.021823e-144
[198,] 1.0000000 1.489763e-142 7.448815e-143
[199,] 1.0000000 1.124895e-141 5.624473e-142
[200,] 1.0000000 9.285318e-141 4.642659e-141
[201,] 1.0000000 2.830105e-140 1.415053e-140
[202,] 1.0000000 2.233113e-139 1.116556e-139
[203,] 1.0000000 1.897353e-138 9.486763e-139
[204,] 1.0000000 1.605142e-137 8.025712e-138
[205,] 1.0000000 1.236876e-136 6.184382e-137
[206,] 1.0000000 7.594920e-136 3.797460e-136
[207,] 1.0000000 5.802934e-135 2.901467e-135
[208,] 1.0000000 4.828425e-134 2.414213e-134
[209,] 1.0000000 4.001215e-133 2.000607e-133
[210,] 1.0000000 3.358999e-132 1.679500e-132
[211,] 1.0000000 2.759562e-131 1.379781e-131
[212,] 1.0000000 1.859136e-130 9.295678e-131
[213,] 1.0000000 1.320228e-129 6.601142e-130
[214,] 1.0000000 1.073671e-128 5.368357e-129
[215,] 1.0000000 8.694309e-128 4.347155e-128
[216,] 1.0000000 7.040859e-127 3.520430e-127
[217,] 1.0000000 5.653002e-126 2.826501e-126
[218,] 1.0000000 4.522004e-125 2.261002e-125
[219,] 1.0000000 3.591051e-124 1.795525e-124
[220,] 1.0000000 2.842281e-123 1.421141e-123
[221,] 1.0000000 2.009232e-122 1.004616e-122
[222,] 1.0000000 1.574298e-121 7.871490e-122
[223,] 1.0000000 9.592876e-121 4.796438e-121
[224,] 1.0000000 7.566092e-120 3.783046e-120
[225,] 1.0000000 5.845772e-119 2.922886e-119
[226,] 1.0000000 4.101848e-118 2.050924e-118
[227,] 1.0000000 3.136694e-117 1.568347e-117
[228,] 1.0000000 2.388046e-116 1.194023e-116
[229,] 1.0000000 1.808837e-115 9.044183e-116
[230,] 1.0000000 1.374516e-114 6.872582e-115
[231,] 1.0000000 1.031393e-113 5.156967e-114
[232,] 1.0000000 7.701906e-113 3.850953e-113
[233,] 1.0000000 5.723502e-112 2.861751e-112
[234,] 1.0000000 4.151705e-111 2.075853e-111
[235,] 1.0000000 3.054121e-110 1.527061e-110
[236,] 1.0000000 2.262317e-109 1.131158e-109
[237,] 1.0000000 1.577420e-108 7.887100e-109
[238,] 1.0000000 1.117470e-107 5.587352e-108
[239,] 1.0000000 7.885128e-107 3.942564e-107
[240,] 1.0000000 5.358671e-106 2.679335e-106
[241,] 1.0000000 3.747670e-105 1.873835e-105
[242,] 1.0000000 2.609784e-104 1.304892e-104
[243,] 1.0000000 1.870328e-103 9.351641e-104
[244,] 1.0000000 1.291012e-102 6.455060e-103
[245,] 1.0000000 7.857566e-102 3.928783e-102
[246,] 1.0000000 3.186289e-101 1.593144e-101
[247,] 1.0000000 6.151534e-102 3.075767e-102
[248,] 1.0000000 4.270225e-101 2.135112e-101
[249,] 1.0000000 3.086918e-100 1.543459e-100
[250,] 1.0000000 2.205365e-99 1.102682e-99
[251,] 1.0000000 1.507519e-98 7.537593e-99
[252,] 1.0000000 1.031568e-97 5.157840e-98
[253,] 1.0000000 6.987961e-97 3.493980e-97
[254,] 1.0000000 4.711436e-96 2.355718e-96
[255,] 1.0000000 3.116187e-95 1.558094e-95
[256,] 1.0000000 2.116939e-94 1.058470e-94
[257,] 1.0000000 1.410379e-93 7.051893e-94
[258,] 1.0000000 9.423939e-93 4.711970e-93
[259,] 1.0000000 6.219654e-92 3.109827e-92
[260,] 1.0000000 4.156511e-91 2.078255e-91
[261,] 1.0000000 2.721886e-90 1.360943e-90
[262,] 1.0000000 1.693425e-89 8.467125e-90
[263,] 1.0000000 1.097962e-88 5.489808e-89
[264,] 1.0000000 7.080714e-88 3.540357e-88
[265,] 1.0000000 4.540766e-87 2.270383e-87
[266,] 1.0000000 2.829782e-86 1.414891e-86
[267,] 1.0000000 1.817427e-85 9.087134e-86
[268,] 1.0000000 1.161385e-84 5.806923e-85
[269,] 1.0000000 7.189676e-84 3.594838e-84
[270,] 1.0000000 4.546009e-83 2.273005e-83
[271,] 1.0000000 1.398810e-82 6.994051e-83
[272,] 1.0000000 6.275686e-82 3.137843e-82
[273,] 1.0000000 3.863207e-81 1.931604e-81
[274,] 1.0000000 2.398864e-80 1.199432e-80
[275,] 1.0000000 1.468775e-79 7.343877e-80
[276,] 1.0000000 9.028335e-79 4.514168e-79
[277,] 1.0000000 5.052306e-78 2.526153e-78
[278,] 1.0000000 3.068809e-77 1.534405e-77
[279,] 1.0000000 1.890460e-76 9.452298e-77
[280,] 1.0000000 1.129239e-75 5.646194e-76
[281,] 1.0000000 6.704396e-75 3.352198e-75
[282,] 1.0000000 3.492248e-74 1.746124e-74
[283,] 1.0000000 2.052099e-73 1.026050e-73
[284,] 1.0000000 1.216532e-72 6.082660e-73
[285,] 1.0000000 6.933466e-72 3.466733e-72
[286,] 1.0000000 4.045770e-71 2.022885e-71
[287,] 1.0000000 2.113795e-70 1.056898e-70
[288,] 1.0000000 1.209689e-69 6.048443e-70
[289,] 1.0000000 6.968100e-69 3.484050e-69
[290,] 1.0000000 3.988206e-68 1.994103e-68
[291,] 1.0000000 2.275480e-67 1.137740e-67
[292,] 1.0000000 1.291287e-66 6.456433e-67
[293,] 1.0000000 7.145543e-66 3.572771e-66
[294,] 1.0000000 3.816662e-65 1.908331e-65
[295,] 1.0000000 2.130522e-64 1.065261e-64
[296,] 1.0000000 1.179190e-63 5.895950e-64
[297,] 1.0000000 6.724510e-63 3.362255e-63
[298,] 1.0000000 3.674526e-62 1.837263e-62
[299,] 1.0000000 2.003543e-61 1.001771e-61
[300,] 1.0000000 1.020572e-60 5.102861e-61
[301,] 1.0000000 5.729089e-60 2.864545e-60
[302,] 1.0000000 3.052548e-59 1.526274e-59
[303,] 1.0000000 1.610343e-58 8.051713e-59
[304,] 1.0000000 3.626974e-58 1.813487e-58
[305,] 1.0000000 1.893155e-57 9.465775e-58
[306,] 1.0000000 9.417253e-57 4.708626e-57
[307,] 1.0000000 4.936750e-56 2.468375e-56
[308,] 1.0000000 2.663227e-55 1.331614e-55
[309,] 1.0000000 1.401157e-54 7.005784e-55
[310,] 1.0000000 6.002863e-54 3.001431e-54
[311,] 1.0000000 3.062509e-53 1.531254e-53
[312,] 1.0000000 1.568720e-52 7.843602e-53
[313,] 1.0000000 7.992001e-52 3.996001e-52
[314,] 1.0000000 4.049453e-51 2.024727e-51
[315,] 1.0000000 2.040594e-50 1.020297e-50
[316,] 1.0000000 9.858699e-50 4.929349e-50
[317,] 1.0000000 4.911602e-49 2.455801e-49
[318,] 1.0000000 2.388471e-48 1.194235e-48
[319,] 1.0000000 1.118447e-47 5.592235e-48
[320,] 1.0000000 5.091928e-47 2.545964e-47
[321,] 1.0000000 2.471670e-46 1.235835e-46
[322,] 1.0000000 1.192980e-45 5.964902e-46
[323,] 1.0000000 5.725285e-45 2.862643e-45
[324,] 1.0000000 2.731922e-44 1.365961e-44
[325,] 1.0000000 1.296081e-43 6.480404e-44
[326,] 1.0000000 6.078515e-43 3.039257e-43
[327,] 1.0000000 2.843432e-42 1.421716e-42
[328,] 1.0000000 1.333776e-41 6.668881e-42
[329,] 1.0000000 4.466083e-41 2.233042e-41
[330,] 1.0000000 2.055124e-40 1.027562e-40
[331,] 1.0000000 9.298580e-40 4.649290e-40
[332,] 1.0000000 4.273578e-39 2.136789e-39
[333,] 1.0000000 1.500689e-38 7.503443e-39
[334,] 1.0000000 6.723012e-38 3.361506e-38
[335,] 1.0000000 2.993973e-37 1.496986e-37
[336,] 1.0000000 1.325334e-36 6.626668e-37
[337,] 1.0000000 5.176679e-36 2.588339e-36
[338,] 1.0000000 2.260424e-35 1.130212e-35
[339,] 1.0000000 9.323828e-35 4.661914e-35
[340,] 1.0000000 4.016530e-34 2.008265e-34
[341,] 1.0000000 1.719686e-33 8.598431e-34
[342,] 1.0000000 7.574332e-33 3.787166e-33
[343,] 1.0000000 3.295568e-32 1.647784e-32
[344,] 1.0000000 1.356348e-31 6.781740e-32
[345,] 1.0000000 5.881285e-31 2.940642e-31
[346,] 1.0000000 2.407009e-30 1.203504e-30
[347,] 1.0000000 9.976592e-30 4.988296e-30
[348,] 1.0000000 4.108175e-29 2.054088e-29
[349,] 1.0000000 9.444924e-29 4.722462e-29
[350,] 1.0000000 3.834291e-28 1.917146e-28
[351,] 1.0000000 1.546504e-27 7.732519e-28
[352,] 1.0000000 6.196694e-27 3.098347e-27
[353,] 1.0000000 2.453025e-26 1.226513e-26
[354,] 1.0000000 9.696166e-26 4.848083e-26
[355,] 1.0000000 3.806424e-25 1.903212e-25
[356,] 1.0000000 1.336261e-24 6.681303e-25
[357,] 1.0000000 5.190094e-24 2.595047e-24
[358,] 1.0000000 1.991091e-23 9.955454e-24
[359,] 1.0000000 7.583449e-23 3.791725e-23
[360,] 1.0000000 2.867061e-22 1.433531e-22
[361,] 1.0000000 1.075795e-21 5.378977e-22
[362,] 1.0000000 3.946850e-21 1.973425e-21
[363,] 1.0000000 1.457539e-20 7.287697e-21
[364,] 1.0000000 5.355144e-20 2.677572e-20
[365,] 1.0000000 1.762319e-19 8.811597e-20
[366,] 1.0000000 5.903226e-19 2.951613e-19
[367,] 1.0000000 2.103424e-18 1.051712e-18
[368,] 1.0000000 7.430136e-18 3.715068e-18
[369,] 1.0000000 2.492973e-17 1.246486e-17
[370,] 1.0000000 8.711215e-17 4.355607e-17
[371,] 1.0000000 2.980005e-16 1.490002e-16
[372,] 1.0000000 1.015336e-15 5.076679e-16
[373,] 1.0000000 3.421759e-15 1.710880e-15
[374,] 1.0000000 9.315506e-15 4.657753e-15
[375,] 1.0000000 3.054829e-14 1.527415e-14
[376,] 1.0000000 9.975289e-14 4.987645e-14
[377,] 1.0000000 3.197171e-13 1.598586e-13
[378,] 1.0000000 5.175394e-13 2.587697e-13
[379,] 1.0000000 1.648750e-12 8.243748e-13
[380,] 1.0000000 5.197907e-12 2.598953e-12
[381,] 1.0000000 1.616802e-11 8.084009e-12
[382,] 1.0000000 4.954103e-11 2.477052e-11
[383,] 1.0000000 1.492874e-10 7.464372e-11
[384,] 1.0000000 4.398662e-10 2.199331e-10
[385,] 1.0000000 1.144726e-09 5.723628e-10
[386,] 1.0000000 2.983455e-09 1.491728e-09
[387,] 1.0000000 5.922027e-09 2.961013e-09
[388,] 1.0000000 1.656015e-08 8.280077e-09
[389,] 1.0000000 4.396224e-08 2.198112e-08
[390,] 1.0000000 7.238622e-08 3.619311e-08
[391,] 0.9999999 1.964606e-07 9.823030e-08
[392,] 0.9999998 4.450165e-07 2.225082e-07
[393,] 0.9999994 1.161090e-06 5.805450e-07
[394,] 0.9999985 2.984062e-06 1.492031e-06
[395,] 0.9999962 7.664921e-06 3.832461e-06
[396,] 0.9999928 1.430934e-05 7.154669e-06
[397,] 0.9999827 3.460435e-05 1.730218e-05
[398,] 0.9999577 8.460541e-05 4.230270e-05
[399,] 0.9998991 2.018981e-04 1.009490e-04
[400,] 0.9997735 4.530900e-04 2.265450e-04
[401,] 0.9995066 9.868163e-04 4.934082e-04
[402,] 0.9990026 1.994881e-03 9.974404e-04
[403,] 0.9979303 4.139357e-03 2.069679e-03
[404,] 0.9958082 8.383675e-03 4.191837e-03
[405,] 0.9915924 1.681520e-02 8.407600e-03
[406,] 0.9890337 2.193266e-02 1.096633e-02
[407,] 0.9856072 2.878567e-02 1.439284e-02
[408,] 0.9775658 4.486845e-02 2.243423e-02
[409,] 0.9881766 2.364676e-02 1.182338e-02
[410,] 0.9751101 4.977978e-02 2.488989e-02
[411,] 0.9514808 9.703838e-02 4.851919e-02
[412,] 0.9068904 1.862192e-01 9.310962e-02
[413,] 0.8341798 3.316403e-01 1.658202e-01
[414,] 0.7693017 4.613966e-01 2.306983e-01
> postscript(file="/var/www/html/freestat/rcomp/tmp/1fdkq1291223383.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> plot(x[,1], type='l', main='Actuals and Interpolation', ylab='value of Actuals and Interpolation (dots)', xlab='time or index')
> points(x[,1]-mysum$resid)
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/freestat/rcomp/tmp/2fdkq1291223383.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> plot(mysum$resid, type='b', pch=19, main='Residuals', ylab='value of Residuals', xlab='time or index')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/freestat/rcomp/tmp/3fdkq1291223383.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> hist(mysum$resid, main='Residual Histogram', xlab='values of Residuals')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/freestat/rcomp/tmp/484kt1291223383.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> densityplot(~mysum$resid,col='black',main='Residual Density Plot', xlab='values of Residuals')
> dev.off()
null device
1
> postscript(file="/var/www/html/freestat/rcomp/tmp/584kt1291223383.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
1284205.5414 3187769.1732 1190500.6719 -2175874.6365 -468857.7504
6 7 8 9 10
-90790.2427 -26679.9252 72632.9742 -283972.6976 689755.2837
11 12 13 14 15
-273154.5235 391457.0462 746948.1834 -1198904.2326 4923.7717
16 17 18 19 20
491349.8015 -230376.0355 -10589.6859 -317846.2350 -325510.8117
21 22 23 24 25
143082.0979 -290933.1323 -128896.5763 -384255.3887 -289754.1051
26 27 28 29 30
-149268.0248 -409198.1778 -297881.3181 314718.7937 -98990.2136
31 32 33 34 35
213934.4307 -162462.4642 178442.5490 70473.2483 -139861.8286
36 37 38 39 40
-469274.7105 36163.3018 140904.3362 -63919.8915 -170551.2589
41 42 43 44 45
5074.2454 -386406.3673 119682.3108 -167669.1703 -23009.4475
46 47 48 49 50
-134082.1738 -293289.1576 -130071.2149 -36657.4800 -171578.8841
51 52 53 54 55
59693.6577 -12903.4165 -247866.8332 69915.9152 -349096.0264
56 57 58 59 60
92688.0486 281287.5672 -17778.3211 -4730.0110 -231383.0298
61 62 63 64 65
-49544.9854 -164405.2161 -61722.1836 -198703.3558 30439.4419
66 67 68 69 70
34164.2920 -52044.9352 -49159.9099 -38217.7441 -75743.8977
71 72 73 74 75
-137098.0592 5261.7630 -83664.6369 -1691.8845 -37590.3118
76 77 78 79 80
-162383.4934 -22428.0277 -74484.0321 27515.6407 -45306.8495
81 82 83 84 85
35998.9557 -68712.6880 -72428.5580 43390.1431 10792.7087
86 87 88 89 90
-121329.8845 52116.0046 -15736.2398 -90993.6726 15434.7945
91 92 93 94 95
11934.0118 -65600.8623 142726.8062 -252110.6991 84496.8895
96 97 98 99 100
-27816.9201 48646.1016 23808.4614 -70514.9164 -74082.9420
101 102 103 104 105
-76550.7077 -118139.3378 -133571.9433 -71993.9074 41430.7063
106 107 108 109 110
-41607.6573 -274184.3298 58310.5508 -175406.9829 -146369.5250
111 112 113 114 115
102020.2593 -231903.5840 -155256.9593 -131755.8649 15041.9739
116 117 118 119 120
45256.5727 9205.4005 55904.1942 -206784.6764 41591.1305
121 122 123 124 125
-269345.1492 24916.8122 17612.2972 26700.1079 -44470.8119
126 127 128 129 130
-229310.8781 37980.4618 37298.6434 45444.0513 38215.5862
131 132 133 134 135
39873.1091 82551.4642 38450.7106 35702.0045 18239.6111
136 137 138 139 140
33178.5791 44059.3673 39064.0778 13881.2785 38999.3343
141 142 143 144 145
46834.7729 39156.0840 6333.6478 39669.4956 41737.3438
146 147 148 149 150
40198.9096 39547.9580 41344.2876 42361.1560 39750.6735
151 152 153 154 155
8756.4069 9926.2905 40018.2069 17431.8785 40338.5765
156 157 158 159 160
40253.3314 40331.7062 36794.0965 52256.6186 40566.8306
161 162 163 164 165
40645.2054 40723.5803 40801.9551 28829.5624 31585.8041
166 167 168 169 170
40433.6828 41115.4543 121791.5073 98416.2146 41350.5788
171 172 173 174 175
41428.9536 37642.5124 41585.7032 31477.9707 41742.4528
176 177 178 179 180
42125.9025 93491.6577 -11977.8421 42055.9521 42134.3269
181 182 183 184 185
51809.8333 -27618.0822 1604.0199 -37849.6760 -12967.9289
186 187 188 189 190
-4912.1108 -37263.5792 36004.1943 5436.2354 42918.0750
191 192 193 194 195
42135.1276 43074.8247 39855.4995 43231.5743 43334.6202
196 197 198 199 200
-9633.0377 -121441.9611 -9476.2880 -9397.9132 -96056.2200
201 202 203 204 205
-9241.1636 -9162.7888 -9084.4140 405.7495 -55280.3682
206 207 208 209 210
-9228.8433 -14856.8323 -8692.5399 21996.6312 -11480.8201
211 212 213 214 215
44563.9461 44642.3210 70918.7810 10603.4828 17599.4573
216 217 218 219 220
44955.8202 48145.1926 47862.5583 45190.9447 17392.7647
221 222 223 224 225
20153.9316 45426.0691 46292.9672 37483.7138 45968.4724
226 227 228 229 230
46965.4284 49215.5708 45896.3180 30566.5692 46053.0676
231 232 233 234 235
18018.9923 34530.2267 46288.1921 28814.3868 46444.9417
236 237 238 239 240
46960.1975 46601.6913 48785.8793 46758.4409 46836.8158
241 242 243 244 245
46915.1906 41818.4282 47071.9402 51803.6457 -33513.3776
246 247 248 249 250
-5714.2970 -5635.9222 39918.7656 -5479.1725 -5400.7977
251 252 253 254 255
40752.2201 -5244.0481 29206.7092 27332.3264 -83022.1079
256 257 258 259 260
48090.8128 45767.9891 59514.2771 -4695.4244 48404.3120
261 262 263 264 265
-4538.6748 -4460.2999 53371.2611 66045.5447 -4225.1755
266 267 268 269 270
48874.5609 -4068.4259 82857.3389 -3911.6762 46988.9064
271 272 273 274 275
-3754.9266 -3676.5518 -3598.1770 45635.0936 49579.9342
276 277 278 279 280
55660.1016 52669.9998 49815.0587 -27957.7533 38477.7869
281 282 283 284 285
-2971.1785 50128.5579 79323.0019 50285.3076 43451.9650
286 287 288 289 290
50442.0572 -13990.7895 51714.8857 -2344.1800 59233.0149
291 292 293 294 295
-2187.4303 62255.4666 59639.3015 51069.0557 44872.3750
296 297 298 299 300
23965.2363 51304.1802 -1638.8066 51460.9298 51539.3046
301 302 303 304 305
-2930.0513 33247.9534 51774.4290 -1168.5578 90952.3650
306 307 308 309 310
55822.1264 -933.4333 -6260.6167 55849.5040 -698.3089
311 312 313 314 315
-619.9341 60022.2974 -463.1844 2035.8088 54553.1757
316 317 318 319 320
69254.8372 53278.9571 30262.6125 55582.9558 53106.8009
321 322 323 324 325
53185.1757 53263.5505 53341.9253 398.9385 53498.6750
326 327 328 329 330
50903.2519 634.0630 712.4378 53812.1742 53890.5490
331 332 333 334 335
53968.9238 54047.2987 54125.6735 52477.5435 55111.2563
336 337 338 339 340
47153.0532 43345.8869 54517.5475 52429.4205 60691.8419
341 342 343 344 345
28236.5456 54831.0468 54909.4216 54987.7964 50159.7441
346 347 348 349 350
55144.5460 45508.0229 55301.2957 55379.6705 19254.3082
351 352 353 354 355
59492.5143 56070.8167 55767.3742 48970.2762 55849.9194
356 357 358 359 360
55928.2942 36949.9413 56085.0438 56163.4186 56241.7934
361 362 363 364 365
58031.6874 56398.5431 56476.9179 -47304.2524 59018.5915
366 367 368 369 370
56712.0423 56790.4171 56868.7920 56947.1668 59768.6666
371 372 373 374 375
57103.9164 70807.9878 46377.8269 37095.0566 57417.4157
376 377 378 379 380
57495.7905 52666.8413 51794.5121 2339.7637 57809.2897
381 382 383 384 385
57887.6645 46625.4127 75947.9325 58122.7890 61096.1226
386 387 388 389 390
-88897.9241 70681.5594 31901.0989 57085.9673 44937.4683
391 392 393 394 395
54448.0084 28112.8887 -9212.4280 80465.3240 -164475.8443
396 397 398 399 400
20295.3487 -126452.0213 26193.5371 22511.8604 -9189.2937
401 402 403 404 405
1051.0932 51968.6029 29208.1466 -80264.3333 143.0929
406 407 408 409 410
23533.3418 36605.7253 27744.3089 38904.9583 -60715.3444
411 412 413 414 415
24985.3375 -234351.8343 -59353.2269 29258.0085 -174529.2171
416 417 418 419 420
-101353.6737 130756.7812 -46498.0682 -21419.1056 -35853.6695
421 422 423 424 425
-46844.9274 -91192.4427 -565498.6666 -470534.4341 -28246.5365
426 427 428 429 430
-213608.3783 -299955.6190 -1090321.6984 -289716.9918 -326935.9312
431
-660177.0992
> postscript(file="/var/www/html/freestat/rcomp/tmp/684kt1291223383.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 1284205.5414 NA
1 3187769.1732 1284205.5414
2 1190500.6719 3187769.1732
3 -2175874.6365 1190500.6719
4 -468857.7504 -2175874.6365
5 -90790.2427 -468857.7504
6 -26679.9252 -90790.2427
7 72632.9742 -26679.9252
8 -283972.6976 72632.9742
9 689755.2837 -283972.6976
10 -273154.5235 689755.2837
11 391457.0462 -273154.5235
12 746948.1834 391457.0462
13 -1198904.2326 746948.1834
14 4923.7717 -1198904.2326
15 491349.8015 4923.7717
16 -230376.0355 491349.8015
17 -10589.6859 -230376.0355
18 -317846.2350 -10589.6859
19 -325510.8117 -317846.2350
20 143082.0979 -325510.8117
21 -290933.1323 143082.0979
22 -128896.5763 -290933.1323
23 -384255.3887 -128896.5763
24 -289754.1051 -384255.3887
25 -149268.0248 -289754.1051
26 -409198.1778 -149268.0248
27 -297881.3181 -409198.1778
28 314718.7937 -297881.3181
29 -98990.2136 314718.7937
30 213934.4307 -98990.2136
31 -162462.4642 213934.4307
32 178442.5490 -162462.4642
33 70473.2483 178442.5490
34 -139861.8286 70473.2483
35 -469274.7105 -139861.8286
36 36163.3018 -469274.7105
37 140904.3362 36163.3018
38 -63919.8915 140904.3362
39 -170551.2589 -63919.8915
40 5074.2454 -170551.2589
41 -386406.3673 5074.2454
42 119682.3108 -386406.3673
43 -167669.1703 119682.3108
44 -23009.4475 -167669.1703
45 -134082.1738 -23009.4475
46 -293289.1576 -134082.1738
47 -130071.2149 -293289.1576
48 -36657.4800 -130071.2149
49 -171578.8841 -36657.4800
50 59693.6577 -171578.8841
51 -12903.4165 59693.6577
52 -247866.8332 -12903.4165
53 69915.9152 -247866.8332
54 -349096.0264 69915.9152
55 92688.0486 -349096.0264
56 281287.5672 92688.0486
57 -17778.3211 281287.5672
58 -4730.0110 -17778.3211
59 -231383.0298 -4730.0110
60 -49544.9854 -231383.0298
61 -164405.2161 -49544.9854
62 -61722.1836 -164405.2161
63 -198703.3558 -61722.1836
64 30439.4419 -198703.3558
65 34164.2920 30439.4419
66 -52044.9352 34164.2920
67 -49159.9099 -52044.9352
68 -38217.7441 -49159.9099
69 -75743.8977 -38217.7441
70 -137098.0592 -75743.8977
71 5261.7630 -137098.0592
72 -83664.6369 5261.7630
73 -1691.8845 -83664.6369
74 -37590.3118 -1691.8845
75 -162383.4934 -37590.3118
76 -22428.0277 -162383.4934
77 -74484.0321 -22428.0277
78 27515.6407 -74484.0321
79 -45306.8495 27515.6407
80 35998.9557 -45306.8495
81 -68712.6880 35998.9557
82 -72428.5580 -68712.6880
83 43390.1431 -72428.5580
84 10792.7087 43390.1431
85 -121329.8845 10792.7087
86 52116.0046 -121329.8845
87 -15736.2398 52116.0046
88 -90993.6726 -15736.2398
89 15434.7945 -90993.6726
90 11934.0118 15434.7945
91 -65600.8623 11934.0118
92 142726.8062 -65600.8623
93 -252110.6991 142726.8062
94 84496.8895 -252110.6991
95 -27816.9201 84496.8895
96 48646.1016 -27816.9201
97 23808.4614 48646.1016
98 -70514.9164 23808.4614
99 -74082.9420 -70514.9164
100 -76550.7077 -74082.9420
101 -118139.3378 -76550.7077
102 -133571.9433 -118139.3378
103 -71993.9074 -133571.9433
104 41430.7063 -71993.9074
105 -41607.6573 41430.7063
106 -274184.3298 -41607.6573
107 58310.5508 -274184.3298
108 -175406.9829 58310.5508
109 -146369.5250 -175406.9829
110 102020.2593 -146369.5250
111 -231903.5840 102020.2593
112 -155256.9593 -231903.5840
113 -131755.8649 -155256.9593
114 15041.9739 -131755.8649
115 45256.5727 15041.9739
116 9205.4005 45256.5727
117 55904.1942 9205.4005
118 -206784.6764 55904.1942
119 41591.1305 -206784.6764
120 -269345.1492 41591.1305
121 24916.8122 -269345.1492
122 17612.2972 24916.8122
123 26700.1079 17612.2972
124 -44470.8119 26700.1079
125 -229310.8781 -44470.8119
126 37980.4618 -229310.8781
127 37298.6434 37980.4618
128 45444.0513 37298.6434
129 38215.5862 45444.0513
130 39873.1091 38215.5862
131 82551.4642 39873.1091
132 38450.7106 82551.4642
133 35702.0045 38450.7106
134 18239.6111 35702.0045
135 33178.5791 18239.6111
136 44059.3673 33178.5791
137 39064.0778 44059.3673
138 13881.2785 39064.0778
139 38999.3343 13881.2785
140 46834.7729 38999.3343
141 39156.0840 46834.7729
142 6333.6478 39156.0840
143 39669.4956 6333.6478
144 41737.3438 39669.4956
145 40198.9096 41737.3438
146 39547.9580 40198.9096
147 41344.2876 39547.9580
148 42361.1560 41344.2876
149 39750.6735 42361.1560
150 8756.4069 39750.6735
151 9926.2905 8756.4069
152 40018.2069 9926.2905
153 17431.8785 40018.2069
154 40338.5765 17431.8785
155 40253.3314 40338.5765
156 40331.7062 40253.3314
157 36794.0965 40331.7062
158 52256.6186 36794.0965
159 40566.8306 52256.6186
160 40645.2054 40566.8306
161 40723.5803 40645.2054
162 40801.9551 40723.5803
163 28829.5624 40801.9551
164 31585.8041 28829.5624
165 40433.6828 31585.8041
166 41115.4543 40433.6828
167 121791.5073 41115.4543
168 98416.2146 121791.5073
169 41350.5788 98416.2146
170 41428.9536 41350.5788
171 37642.5124 41428.9536
172 41585.7032 37642.5124
173 31477.9707 41585.7032
174 41742.4528 31477.9707
175 42125.9025 41742.4528
176 93491.6577 42125.9025
177 -11977.8421 93491.6577
178 42055.9521 -11977.8421
179 42134.3269 42055.9521
180 51809.8333 42134.3269
181 -27618.0822 51809.8333
182 1604.0199 -27618.0822
183 -37849.6760 1604.0199
184 -12967.9289 -37849.6760
185 -4912.1108 -12967.9289
186 -37263.5792 -4912.1108
187 36004.1943 -37263.5792
188 5436.2354 36004.1943
189 42918.0750 5436.2354
190 42135.1276 42918.0750
191 43074.8247 42135.1276
192 39855.4995 43074.8247
193 43231.5743 39855.4995
194 43334.6202 43231.5743
195 -9633.0377 43334.6202
196 -121441.9611 -9633.0377
197 -9476.2880 -121441.9611
198 -9397.9132 -9476.2880
199 -96056.2200 -9397.9132
200 -9241.1636 -96056.2200
201 -9162.7888 -9241.1636
202 -9084.4140 -9162.7888
203 405.7495 -9084.4140
204 -55280.3682 405.7495
205 -9228.8433 -55280.3682
206 -14856.8323 -9228.8433
207 -8692.5399 -14856.8323
208 21996.6312 -8692.5399
209 -11480.8201 21996.6312
210 44563.9461 -11480.8201
211 44642.3210 44563.9461
212 70918.7810 44642.3210
213 10603.4828 70918.7810
214 17599.4573 10603.4828
215 44955.8202 17599.4573
216 48145.1926 44955.8202
217 47862.5583 48145.1926
218 45190.9447 47862.5583
219 17392.7647 45190.9447
220 20153.9316 17392.7647
221 45426.0691 20153.9316
222 46292.9672 45426.0691
223 37483.7138 46292.9672
224 45968.4724 37483.7138
225 46965.4284 45968.4724
226 49215.5708 46965.4284
227 45896.3180 49215.5708
228 30566.5692 45896.3180
229 46053.0676 30566.5692
230 18018.9923 46053.0676
231 34530.2267 18018.9923
232 46288.1921 34530.2267
233 28814.3868 46288.1921
234 46444.9417 28814.3868
235 46960.1975 46444.9417
236 46601.6913 46960.1975
237 48785.8793 46601.6913
238 46758.4409 48785.8793
239 46836.8158 46758.4409
240 46915.1906 46836.8158
241 41818.4282 46915.1906
242 47071.9402 41818.4282
243 51803.6457 47071.9402
244 -33513.3776 51803.6457
245 -5714.2970 -33513.3776
246 -5635.9222 -5714.2970
247 39918.7656 -5635.9222
248 -5479.1725 39918.7656
249 -5400.7977 -5479.1725
250 40752.2201 -5400.7977
251 -5244.0481 40752.2201
252 29206.7092 -5244.0481
253 27332.3264 29206.7092
254 -83022.1079 27332.3264
255 48090.8128 -83022.1079
256 45767.9891 48090.8128
257 59514.2771 45767.9891
258 -4695.4244 59514.2771
259 48404.3120 -4695.4244
260 -4538.6748 48404.3120
261 -4460.2999 -4538.6748
262 53371.2611 -4460.2999
263 66045.5447 53371.2611
264 -4225.1755 66045.5447
265 48874.5609 -4225.1755
266 -4068.4259 48874.5609
267 82857.3389 -4068.4259
268 -3911.6762 82857.3389
269 46988.9064 -3911.6762
270 -3754.9266 46988.9064
271 -3676.5518 -3754.9266
272 -3598.1770 -3676.5518
273 45635.0936 -3598.1770
274 49579.9342 45635.0936
275 55660.1016 49579.9342
276 52669.9998 55660.1016
277 49815.0587 52669.9998
278 -27957.7533 49815.0587
279 38477.7869 -27957.7533
280 -2971.1785 38477.7869
281 50128.5579 -2971.1785
282 79323.0019 50128.5579
283 50285.3076 79323.0019
284 43451.9650 50285.3076
285 50442.0572 43451.9650
286 -13990.7895 50442.0572
287 51714.8857 -13990.7895
288 -2344.1800 51714.8857
289 59233.0149 -2344.1800
290 -2187.4303 59233.0149
291 62255.4666 -2187.4303
292 59639.3015 62255.4666
293 51069.0557 59639.3015
294 44872.3750 51069.0557
295 23965.2363 44872.3750
296 51304.1802 23965.2363
297 -1638.8066 51304.1802
298 51460.9298 -1638.8066
299 51539.3046 51460.9298
300 -2930.0513 51539.3046
301 33247.9534 -2930.0513
302 51774.4290 33247.9534
303 -1168.5578 51774.4290
304 90952.3650 -1168.5578
305 55822.1264 90952.3650
306 -933.4333 55822.1264
307 -6260.6167 -933.4333
308 55849.5040 -6260.6167
309 -698.3089 55849.5040
310 -619.9341 -698.3089
311 60022.2974 -619.9341
312 -463.1844 60022.2974
313 2035.8088 -463.1844
314 54553.1757 2035.8088
315 69254.8372 54553.1757
316 53278.9571 69254.8372
317 30262.6125 53278.9571
318 55582.9558 30262.6125
319 53106.8009 55582.9558
320 53185.1757 53106.8009
321 53263.5505 53185.1757
322 53341.9253 53263.5505
323 398.9385 53341.9253
324 53498.6750 398.9385
325 50903.2519 53498.6750
326 634.0630 50903.2519
327 712.4378 634.0630
328 53812.1742 712.4378
329 53890.5490 53812.1742
330 53968.9238 53890.5490
331 54047.2987 53968.9238
332 54125.6735 54047.2987
333 52477.5435 54125.6735
334 55111.2563 52477.5435
335 47153.0532 55111.2563
336 43345.8869 47153.0532
337 54517.5475 43345.8869
338 52429.4205 54517.5475
339 60691.8419 52429.4205
340 28236.5456 60691.8419
341 54831.0468 28236.5456
342 54909.4216 54831.0468
343 54987.7964 54909.4216
344 50159.7441 54987.7964
345 55144.5460 50159.7441
346 45508.0229 55144.5460
347 55301.2957 45508.0229
348 55379.6705 55301.2957
349 19254.3082 55379.6705
350 59492.5143 19254.3082
351 56070.8167 59492.5143
352 55767.3742 56070.8167
353 48970.2762 55767.3742
354 55849.9194 48970.2762
355 55928.2942 55849.9194
356 36949.9413 55928.2942
357 56085.0438 36949.9413
358 56163.4186 56085.0438
359 56241.7934 56163.4186
360 58031.6874 56241.7934
361 56398.5431 58031.6874
362 56476.9179 56398.5431
363 -47304.2524 56476.9179
364 59018.5915 -47304.2524
365 56712.0423 59018.5915
366 56790.4171 56712.0423
367 56868.7920 56790.4171
368 56947.1668 56868.7920
369 59768.6666 56947.1668
370 57103.9164 59768.6666
371 70807.9878 57103.9164
372 46377.8269 70807.9878
373 37095.0566 46377.8269
374 57417.4157 37095.0566
375 57495.7905 57417.4157
376 52666.8413 57495.7905
377 51794.5121 52666.8413
378 2339.7637 51794.5121
379 57809.2897 2339.7637
380 57887.6645 57809.2897
381 46625.4127 57887.6645
382 75947.9325 46625.4127
383 58122.7890 75947.9325
384 61096.1226 58122.7890
385 -88897.9241 61096.1226
386 70681.5594 -88897.9241
387 31901.0989 70681.5594
388 57085.9673 31901.0989
389 44937.4683 57085.9673
390 54448.0084 44937.4683
391 28112.8887 54448.0084
392 -9212.4280 28112.8887
393 80465.3240 -9212.4280
394 -164475.8443 80465.3240
395 20295.3487 -164475.8443
396 -126452.0213 20295.3487
397 26193.5371 -126452.0213
398 22511.8604 26193.5371
399 -9189.2937 22511.8604
400 1051.0932 -9189.2937
401 51968.6029 1051.0932
402 29208.1466 51968.6029
403 -80264.3333 29208.1466
404 143.0929 -80264.3333
405 23533.3418 143.0929
406 36605.7253 23533.3418
407 27744.3089 36605.7253
408 38904.9583 27744.3089
409 -60715.3444 38904.9583
410 24985.3375 -60715.3444
411 -234351.8343 24985.3375
412 -59353.2269 -234351.8343
413 29258.0085 -59353.2269
414 -174529.2171 29258.0085
415 -101353.6737 -174529.2171
416 130756.7812 -101353.6737
417 -46498.0682 130756.7812
418 -21419.1056 -46498.0682
419 -35853.6695 -21419.1056
420 -46844.9274 -35853.6695
421 -91192.4427 -46844.9274
422 -565498.6666 -91192.4427
423 -470534.4341 -565498.6666
424 -28246.5365 -470534.4341
425 -213608.3783 -28246.5365
426 -299955.6190 -213608.3783
427 -1090321.6984 -299955.6190
428 -289716.9918 -1090321.6984
429 -326935.9312 -289716.9918
430 -660177.0992 -326935.9312
431 NA -660177.0992
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] 3187769.1732 1284205.5414
[2,] 1190500.6719 3187769.1732
[3,] -2175874.6365 1190500.6719
[4,] -468857.7504 -2175874.6365
[5,] -90790.2427 -468857.7504
[6,] -26679.9252 -90790.2427
[7,] 72632.9742 -26679.9252
[8,] -283972.6976 72632.9742
[9,] 689755.2837 -283972.6976
[10,] -273154.5235 689755.2837
[11,] 391457.0462 -273154.5235
[12,] 746948.1834 391457.0462
[13,] -1198904.2326 746948.1834
[14,] 4923.7717 -1198904.2326
[15,] 491349.8015 4923.7717
[16,] -230376.0355 491349.8015
[17,] -10589.6859 -230376.0355
[18,] -317846.2350 -10589.6859
[19,] -325510.8117 -317846.2350
[20,] 143082.0979 -325510.8117
[21,] -290933.1323 143082.0979
[22,] -128896.5763 -290933.1323
[23,] -384255.3887 -128896.5763
[24,] -289754.1051 -384255.3887
[25,] -149268.0248 -289754.1051
[26,] -409198.1778 -149268.0248
[27,] -297881.3181 -409198.1778
[28,] 314718.7937 -297881.3181
[29,] -98990.2136 314718.7937
[30,] 213934.4307 -98990.2136
[31,] -162462.4642 213934.4307
[32,] 178442.5490 -162462.4642
[33,] 70473.2483 178442.5490
[34,] -139861.8286 70473.2483
[35,] -469274.7105 -139861.8286
[36,] 36163.3018 -469274.7105
[37,] 140904.3362 36163.3018
[38,] -63919.8915 140904.3362
[39,] -170551.2589 -63919.8915
[40,] 5074.2454 -170551.2589
[41,] -386406.3673 5074.2454
[42,] 119682.3108 -386406.3673
[43,] -167669.1703 119682.3108
[44,] -23009.4475 -167669.1703
[45,] -134082.1738 -23009.4475
[46,] -293289.1576 -134082.1738
[47,] -130071.2149 -293289.1576
[48,] -36657.4800 -130071.2149
[49,] -171578.8841 -36657.4800
[50,] 59693.6577 -171578.8841
[51,] -12903.4165 59693.6577
[52,] -247866.8332 -12903.4165
[53,] 69915.9152 -247866.8332
[54,] -349096.0264 69915.9152
[55,] 92688.0486 -349096.0264
[56,] 281287.5672 92688.0486
[57,] -17778.3211 281287.5672
[58,] -4730.0110 -17778.3211
[59,] -231383.0298 -4730.0110
[60,] -49544.9854 -231383.0298
[61,] -164405.2161 -49544.9854
[62,] -61722.1836 -164405.2161
[63,] -198703.3558 -61722.1836
[64,] 30439.4419 -198703.3558
[65,] 34164.2920 30439.4419
[66,] -52044.9352 34164.2920
[67,] -49159.9099 -52044.9352
[68,] -38217.7441 -49159.9099
[69,] -75743.8977 -38217.7441
[70,] -137098.0592 -75743.8977
[71,] 5261.7630 -137098.0592
[72,] -83664.6369 5261.7630
[73,] -1691.8845 -83664.6369
[74,] -37590.3118 -1691.8845
[75,] -162383.4934 -37590.3118
[76,] -22428.0277 -162383.4934
[77,] -74484.0321 -22428.0277
[78,] 27515.6407 -74484.0321
[79,] -45306.8495 27515.6407
[80,] 35998.9557 -45306.8495
[81,] -68712.6880 35998.9557
[82,] -72428.5580 -68712.6880
[83,] 43390.1431 -72428.5580
[84,] 10792.7087 43390.1431
[85,] -121329.8845 10792.7087
[86,] 52116.0046 -121329.8845
[87,] -15736.2398 52116.0046
[88,] -90993.6726 -15736.2398
[89,] 15434.7945 -90993.6726
[90,] 11934.0118 15434.7945
[91,] -65600.8623 11934.0118
[92,] 142726.8062 -65600.8623
[93,] -252110.6991 142726.8062
[94,] 84496.8895 -252110.6991
[95,] -27816.9201 84496.8895
[96,] 48646.1016 -27816.9201
[97,] 23808.4614 48646.1016
[98,] -70514.9164 23808.4614
[99,] -74082.9420 -70514.9164
[100,] -76550.7077 -74082.9420
[101,] -118139.3378 -76550.7077
[102,] -133571.9433 -118139.3378
[103,] -71993.9074 -133571.9433
[104,] 41430.7063 -71993.9074
[105,] -41607.6573 41430.7063
[106,] -274184.3298 -41607.6573
[107,] 58310.5508 -274184.3298
[108,] -175406.9829 58310.5508
[109,] -146369.5250 -175406.9829
[110,] 102020.2593 -146369.5250
[111,] -231903.5840 102020.2593
[112,] -155256.9593 -231903.5840
[113,] -131755.8649 -155256.9593
[114,] 15041.9739 -131755.8649
[115,] 45256.5727 15041.9739
[116,] 9205.4005 45256.5727
[117,] 55904.1942 9205.4005
[118,] -206784.6764 55904.1942
[119,] 41591.1305 -206784.6764
[120,] -269345.1492 41591.1305
[121,] 24916.8122 -269345.1492
[122,] 17612.2972 24916.8122
[123,] 26700.1079 17612.2972
[124,] -44470.8119 26700.1079
[125,] -229310.8781 -44470.8119
[126,] 37980.4618 -229310.8781
[127,] 37298.6434 37980.4618
[128,] 45444.0513 37298.6434
[129,] 38215.5862 45444.0513
[130,] 39873.1091 38215.5862
[131,] 82551.4642 39873.1091
[132,] 38450.7106 82551.4642
[133,] 35702.0045 38450.7106
[134,] 18239.6111 35702.0045
[135,] 33178.5791 18239.6111
[136,] 44059.3673 33178.5791
[137,] 39064.0778 44059.3673
[138,] 13881.2785 39064.0778
[139,] 38999.3343 13881.2785
[140,] 46834.7729 38999.3343
[141,] 39156.0840 46834.7729
[142,] 6333.6478 39156.0840
[143,] 39669.4956 6333.6478
[144,] 41737.3438 39669.4956
[145,] 40198.9096 41737.3438
[146,] 39547.9580 40198.9096
[147,] 41344.2876 39547.9580
[148,] 42361.1560 41344.2876
[149,] 39750.6735 42361.1560
[150,] 8756.4069 39750.6735
[151,] 9926.2905 8756.4069
[152,] 40018.2069 9926.2905
[153,] 17431.8785 40018.2069
[154,] 40338.5765 17431.8785
[155,] 40253.3314 40338.5765
[156,] 40331.7062 40253.3314
[157,] 36794.0965 40331.7062
[158,] 52256.6186 36794.0965
[159,] 40566.8306 52256.6186
[160,] 40645.2054 40566.8306
[161,] 40723.5803 40645.2054
[162,] 40801.9551 40723.5803
[163,] 28829.5624 40801.9551
[164,] 31585.8041 28829.5624
[165,] 40433.6828 31585.8041
[166,] 41115.4543 40433.6828
[167,] 121791.5073 41115.4543
[168,] 98416.2146 121791.5073
[169,] 41350.5788 98416.2146
[170,] 41428.9536 41350.5788
[171,] 37642.5124 41428.9536
[172,] 41585.7032 37642.5124
[173,] 31477.9707 41585.7032
[174,] 41742.4528 31477.9707
[175,] 42125.9025 41742.4528
[176,] 93491.6577 42125.9025
[177,] -11977.8421 93491.6577
[178,] 42055.9521 -11977.8421
[179,] 42134.3269 42055.9521
[180,] 51809.8333 42134.3269
[181,] -27618.0822 51809.8333
[182,] 1604.0199 -27618.0822
[183,] -37849.6760 1604.0199
[184,] -12967.9289 -37849.6760
[185,] -4912.1108 -12967.9289
[186,] -37263.5792 -4912.1108
[187,] 36004.1943 -37263.5792
[188,] 5436.2354 36004.1943
[189,] 42918.0750 5436.2354
[190,] 42135.1276 42918.0750
[191,] 43074.8247 42135.1276
[192,] 39855.4995 43074.8247
[193,] 43231.5743 39855.4995
[194,] 43334.6202 43231.5743
[195,] -9633.0377 43334.6202
[196,] -121441.9611 -9633.0377
[197,] -9476.2880 -121441.9611
[198,] -9397.9132 -9476.2880
[199,] -96056.2200 -9397.9132
[200,] -9241.1636 -96056.2200
[201,] -9162.7888 -9241.1636
[202,] -9084.4140 -9162.7888
[203,] 405.7495 -9084.4140
[204,] -55280.3682 405.7495
[205,] -9228.8433 -55280.3682
[206,] -14856.8323 -9228.8433
[207,] -8692.5399 -14856.8323
[208,] 21996.6312 -8692.5399
[209,] -11480.8201 21996.6312
[210,] 44563.9461 -11480.8201
[211,] 44642.3210 44563.9461
[212,] 70918.7810 44642.3210
[213,] 10603.4828 70918.7810
[214,] 17599.4573 10603.4828
[215,] 44955.8202 17599.4573
[216,] 48145.1926 44955.8202
[217,] 47862.5583 48145.1926
[218,] 45190.9447 47862.5583
[219,] 17392.7647 45190.9447
[220,] 20153.9316 17392.7647
[221,] 45426.0691 20153.9316
[222,] 46292.9672 45426.0691
[223,] 37483.7138 46292.9672
[224,] 45968.4724 37483.7138
[225,] 46965.4284 45968.4724
[226,] 49215.5708 46965.4284
[227,] 45896.3180 49215.5708
[228,] 30566.5692 45896.3180
[229,] 46053.0676 30566.5692
[230,] 18018.9923 46053.0676
[231,] 34530.2267 18018.9923
[232,] 46288.1921 34530.2267
[233,] 28814.3868 46288.1921
[234,] 46444.9417 28814.3868
[235,] 46960.1975 46444.9417
[236,] 46601.6913 46960.1975
[237,] 48785.8793 46601.6913
[238,] 46758.4409 48785.8793
[239,] 46836.8158 46758.4409
[240,] 46915.1906 46836.8158
[241,] 41818.4282 46915.1906
[242,] 47071.9402 41818.4282
[243,] 51803.6457 47071.9402
[244,] -33513.3776 51803.6457
[245,] -5714.2970 -33513.3776
[246,] -5635.9222 -5714.2970
[247,] 39918.7656 -5635.9222
[248,] -5479.1725 39918.7656
[249,] -5400.7977 -5479.1725
[250,] 40752.2201 -5400.7977
[251,] -5244.0481 40752.2201
[252,] 29206.7092 -5244.0481
[253,] 27332.3264 29206.7092
[254,] -83022.1079 27332.3264
[255,] 48090.8128 -83022.1079
[256,] 45767.9891 48090.8128
[257,] 59514.2771 45767.9891
[258,] -4695.4244 59514.2771
[259,] 48404.3120 -4695.4244
[260,] -4538.6748 48404.3120
[261,] -4460.2999 -4538.6748
[262,] 53371.2611 -4460.2999
[263,] 66045.5447 53371.2611
[264,] -4225.1755 66045.5447
[265,] 48874.5609 -4225.1755
[266,] -4068.4259 48874.5609
[267,] 82857.3389 -4068.4259
[268,] -3911.6762 82857.3389
[269,] 46988.9064 -3911.6762
[270,] -3754.9266 46988.9064
[271,] -3676.5518 -3754.9266
[272,] -3598.1770 -3676.5518
[273,] 45635.0936 -3598.1770
[274,] 49579.9342 45635.0936
[275,] 55660.1016 49579.9342
[276,] 52669.9998 55660.1016
[277,] 49815.0587 52669.9998
[278,] -27957.7533 49815.0587
[279,] 38477.7869 -27957.7533
[280,] -2971.1785 38477.7869
[281,] 50128.5579 -2971.1785
[282,] 79323.0019 50128.5579
[283,] 50285.3076 79323.0019
[284,] 43451.9650 50285.3076
[285,] 50442.0572 43451.9650
[286,] -13990.7895 50442.0572
[287,] 51714.8857 -13990.7895
[288,] -2344.1800 51714.8857
[289,] 59233.0149 -2344.1800
[290,] -2187.4303 59233.0149
[291,] 62255.4666 -2187.4303
[292,] 59639.3015 62255.4666
[293,] 51069.0557 59639.3015
[294,] 44872.3750 51069.0557
[295,] 23965.2363 44872.3750
[296,] 51304.1802 23965.2363
[297,] -1638.8066 51304.1802
[298,] 51460.9298 -1638.8066
[299,] 51539.3046 51460.9298
[300,] -2930.0513 51539.3046
[301,] 33247.9534 -2930.0513
[302,] 51774.4290 33247.9534
[303,] -1168.5578 51774.4290
[304,] 90952.3650 -1168.5578
[305,] 55822.1264 90952.3650
[306,] -933.4333 55822.1264
[307,] -6260.6167 -933.4333
[308,] 55849.5040 -6260.6167
[309,] -698.3089 55849.5040
[310,] -619.9341 -698.3089
[311,] 60022.2974 -619.9341
[312,] -463.1844 60022.2974
[313,] 2035.8088 -463.1844
[314,] 54553.1757 2035.8088
[315,] 69254.8372 54553.1757
[316,] 53278.9571 69254.8372
[317,] 30262.6125 53278.9571
[318,] 55582.9558 30262.6125
[319,] 53106.8009 55582.9558
[320,] 53185.1757 53106.8009
[321,] 53263.5505 53185.1757
[322,] 53341.9253 53263.5505
[323,] 398.9385 53341.9253
[324,] 53498.6750 398.9385
[325,] 50903.2519 53498.6750
[326,] 634.0630 50903.2519
[327,] 712.4378 634.0630
[328,] 53812.1742 712.4378
[329,] 53890.5490 53812.1742
[330,] 53968.9238 53890.5490
[331,] 54047.2987 53968.9238
[332,] 54125.6735 54047.2987
[333,] 52477.5435 54125.6735
[334,] 55111.2563 52477.5435
[335,] 47153.0532 55111.2563
[336,] 43345.8869 47153.0532
[337,] 54517.5475 43345.8869
[338,] 52429.4205 54517.5475
[339,] 60691.8419 52429.4205
[340,] 28236.5456 60691.8419
[341,] 54831.0468 28236.5456
[342,] 54909.4216 54831.0468
[343,] 54987.7964 54909.4216
[344,] 50159.7441 54987.7964
[345,] 55144.5460 50159.7441
[346,] 45508.0229 55144.5460
[347,] 55301.2957 45508.0229
[348,] 55379.6705 55301.2957
[349,] 19254.3082 55379.6705
[350,] 59492.5143 19254.3082
[351,] 56070.8167 59492.5143
[352,] 55767.3742 56070.8167
[353,] 48970.2762 55767.3742
[354,] 55849.9194 48970.2762
[355,] 55928.2942 55849.9194
[356,] 36949.9413 55928.2942
[357,] 56085.0438 36949.9413
[358,] 56163.4186 56085.0438
[359,] 56241.7934 56163.4186
[360,] 58031.6874 56241.7934
[361,] 56398.5431 58031.6874
[362,] 56476.9179 56398.5431
[363,] -47304.2524 56476.9179
[364,] 59018.5915 -47304.2524
[365,] 56712.0423 59018.5915
[366,] 56790.4171 56712.0423
[367,] 56868.7920 56790.4171
[368,] 56947.1668 56868.7920
[369,] 59768.6666 56947.1668
[370,] 57103.9164 59768.6666
[371,] 70807.9878 57103.9164
[372,] 46377.8269 70807.9878
[373,] 37095.0566 46377.8269
[374,] 57417.4157 37095.0566
[375,] 57495.7905 57417.4157
[376,] 52666.8413 57495.7905
[377,] 51794.5121 52666.8413
[378,] 2339.7637 51794.5121
[379,] 57809.2897 2339.7637
[380,] 57887.6645 57809.2897
[381,] 46625.4127 57887.6645
[382,] 75947.9325 46625.4127
[383,] 58122.7890 75947.9325
[384,] 61096.1226 58122.7890
[385,] -88897.9241 61096.1226
[386,] 70681.5594 -88897.9241
[387,] 31901.0989 70681.5594
[388,] 57085.9673 31901.0989
[389,] 44937.4683 57085.9673
[390,] 54448.0084 44937.4683
[391,] 28112.8887 54448.0084
[392,] -9212.4280 28112.8887
[393,] 80465.3240 -9212.4280
[394,] -164475.8443 80465.3240
[395,] 20295.3487 -164475.8443
[396,] -126452.0213 20295.3487
[397,] 26193.5371 -126452.0213
[398,] 22511.8604 26193.5371
[399,] -9189.2937 22511.8604
[400,] 1051.0932 -9189.2937
[401,] 51968.6029 1051.0932
[402,] 29208.1466 51968.6029
[403,] -80264.3333 29208.1466
[404,] 143.0929 -80264.3333
[405,] 23533.3418 143.0929
[406,] 36605.7253 23533.3418
[407,] 27744.3089 36605.7253
[408,] 38904.9583 27744.3089
[409,] -60715.3444 38904.9583
[410,] 24985.3375 -60715.3444
[411,] -234351.8343 24985.3375
[412,] -59353.2269 -234351.8343
[413,] 29258.0085 -59353.2269
[414,] -174529.2171 29258.0085
[415,] -101353.6737 -174529.2171
[416,] 130756.7812 -101353.6737
[417,] -46498.0682 130756.7812
[418,] -21419.1056 -46498.0682
[419,] -35853.6695 -21419.1056
[420,] -46844.9274 -35853.6695
[421,] -91192.4427 -46844.9274
[422,] -565498.6666 -91192.4427
[423,] -470534.4341 -565498.6666
[424,] -28246.5365 -470534.4341
[425,] -213608.3783 -28246.5365
[426,] -299955.6190 -213608.3783
[427,] -1090321.6984 -299955.6190
[428,] -289716.9918 -1090321.6984
[429,] -326935.9312 -289716.9918
[430,] -660177.0992 -326935.9312
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 3187769.1732 1284205.5414
2 1190500.6719 3187769.1732
3 -2175874.6365 1190500.6719
4 -468857.7504 -2175874.6365
5 -90790.2427 -468857.7504
6 -26679.9252 -90790.2427
7 72632.9742 -26679.9252
8 -283972.6976 72632.9742
9 689755.2837 -283972.6976
10 -273154.5235 689755.2837
11 391457.0462 -273154.5235
12 746948.1834 391457.0462
13 -1198904.2326 746948.1834
14 4923.7717 -1198904.2326
15 491349.8015 4923.7717
16 -230376.0355 491349.8015
17 -10589.6859 -230376.0355
18 -317846.2350 -10589.6859
19 -325510.8117 -317846.2350
20 143082.0979 -325510.8117
21 -290933.1323 143082.0979
22 -128896.5763 -290933.1323
23 -384255.3887 -128896.5763
24 -289754.1051 -384255.3887
25 -149268.0248 -289754.1051
26 -409198.1778 -149268.0248
27 -297881.3181 -409198.1778
28 314718.7937 -297881.3181
29 -98990.2136 314718.7937
30 213934.4307 -98990.2136
31 -162462.4642 213934.4307
32 178442.5490 -162462.4642
33 70473.2483 178442.5490
34 -139861.8286 70473.2483
35 -469274.7105 -139861.8286
36 36163.3018 -469274.7105
37 140904.3362 36163.3018
38 -63919.8915 140904.3362
39 -170551.2589 -63919.8915
40 5074.2454 -170551.2589
41 -386406.3673 5074.2454
42 119682.3108 -386406.3673
43 -167669.1703 119682.3108
44 -23009.4475 -167669.1703
45 -134082.1738 -23009.4475
46 -293289.1576 -134082.1738
47 -130071.2149 -293289.1576
48 -36657.4800 -130071.2149
49 -171578.8841 -36657.4800
50 59693.6577 -171578.8841
51 -12903.4165 59693.6577
52 -247866.8332 -12903.4165
53 69915.9152 -247866.8332
54 -349096.0264 69915.9152
55 92688.0486 -349096.0264
56 281287.5672 92688.0486
57 -17778.3211 281287.5672
58 -4730.0110 -17778.3211
59 -231383.0298 -4730.0110
60 -49544.9854 -231383.0298
61 -164405.2161 -49544.9854
62 -61722.1836 -164405.2161
63 -198703.3558 -61722.1836
64 30439.4419 -198703.3558
65 34164.2920 30439.4419
66 -52044.9352 34164.2920
67 -49159.9099 -52044.9352
68 -38217.7441 -49159.9099
69 -75743.8977 -38217.7441
70 -137098.0592 -75743.8977
71 5261.7630 -137098.0592
72 -83664.6369 5261.7630
73 -1691.8845 -83664.6369
74 -37590.3118 -1691.8845
75 -162383.4934 -37590.3118
76 -22428.0277 -162383.4934
77 -74484.0321 -22428.0277
78 27515.6407 -74484.0321
79 -45306.8495 27515.6407
80 35998.9557 -45306.8495
81 -68712.6880 35998.9557
82 -72428.5580 -68712.6880
83 43390.1431 -72428.5580
84 10792.7087 43390.1431
85 -121329.8845 10792.7087
86 52116.0046 -121329.8845
87 -15736.2398 52116.0046
88 -90993.6726 -15736.2398
89 15434.7945 -90993.6726
90 11934.0118 15434.7945
91 -65600.8623 11934.0118
92 142726.8062 -65600.8623
93 -252110.6991 142726.8062
94 84496.8895 -252110.6991
95 -27816.9201 84496.8895
96 48646.1016 -27816.9201
97 23808.4614 48646.1016
98 -70514.9164 23808.4614
99 -74082.9420 -70514.9164
100 -76550.7077 -74082.9420
101 -118139.3378 -76550.7077
102 -133571.9433 -118139.3378
103 -71993.9074 -133571.9433
104 41430.7063 -71993.9074
105 -41607.6573 41430.7063
106 -274184.3298 -41607.6573
107 58310.5508 -274184.3298
108 -175406.9829 58310.5508
109 -146369.5250 -175406.9829
110 102020.2593 -146369.5250
111 -231903.5840 102020.2593
112 -155256.9593 -231903.5840
113 -131755.8649 -155256.9593
114 15041.9739 -131755.8649
115 45256.5727 15041.9739
116 9205.4005 45256.5727
117 55904.1942 9205.4005
118 -206784.6764 55904.1942
119 41591.1305 -206784.6764
120 -269345.1492 41591.1305
121 24916.8122 -269345.1492
122 17612.2972 24916.8122
123 26700.1079 17612.2972
124 -44470.8119 26700.1079
125 -229310.8781 -44470.8119
126 37980.4618 -229310.8781
127 37298.6434 37980.4618
128 45444.0513 37298.6434
129 38215.5862 45444.0513
130 39873.1091 38215.5862
131 82551.4642 39873.1091
132 38450.7106 82551.4642
133 35702.0045 38450.7106
134 18239.6111 35702.0045
135 33178.5791 18239.6111
136 44059.3673 33178.5791
137 39064.0778 44059.3673
138 13881.2785 39064.0778
139 38999.3343 13881.2785
140 46834.7729 38999.3343
141 39156.0840 46834.7729
142 6333.6478 39156.0840
143 39669.4956 6333.6478
144 41737.3438 39669.4956
145 40198.9096 41737.3438
146 39547.9580 40198.9096
147 41344.2876 39547.9580
148 42361.1560 41344.2876
149 39750.6735 42361.1560
150 8756.4069 39750.6735
151 9926.2905 8756.4069
152 40018.2069 9926.2905
153 17431.8785 40018.2069
154 40338.5765 17431.8785
155 40253.3314 40338.5765
156 40331.7062 40253.3314
157 36794.0965 40331.7062
158 52256.6186 36794.0965
159 40566.8306 52256.6186
160 40645.2054 40566.8306
161 40723.5803 40645.2054
162 40801.9551 40723.5803
163 28829.5624 40801.9551
164 31585.8041 28829.5624
165 40433.6828 31585.8041
166 41115.4543 40433.6828
167 121791.5073 41115.4543
168 98416.2146 121791.5073
169 41350.5788 98416.2146
170 41428.9536 41350.5788
171 37642.5124 41428.9536
172 41585.7032 37642.5124
173 31477.9707 41585.7032
174 41742.4528 31477.9707
175 42125.9025 41742.4528
176 93491.6577 42125.9025
177 -11977.8421 93491.6577
178 42055.9521 -11977.8421
179 42134.3269 42055.9521
180 51809.8333 42134.3269
181 -27618.0822 51809.8333
182 1604.0199 -27618.0822
183 -37849.6760 1604.0199
184 -12967.9289 -37849.6760
185 -4912.1108 -12967.9289
186 -37263.5792 -4912.1108
187 36004.1943 -37263.5792
188 5436.2354 36004.1943
189 42918.0750 5436.2354
190 42135.1276 42918.0750
191 43074.8247 42135.1276
192 39855.4995 43074.8247
193 43231.5743 39855.4995
194 43334.6202 43231.5743
195 -9633.0377 43334.6202
196 -121441.9611 -9633.0377
197 -9476.2880 -121441.9611
198 -9397.9132 -9476.2880
199 -96056.2200 -9397.9132
200 -9241.1636 -96056.2200
201 -9162.7888 -9241.1636
202 -9084.4140 -9162.7888
203 405.7495 -9084.4140
204 -55280.3682 405.7495
205 -9228.8433 -55280.3682
206 -14856.8323 -9228.8433
207 -8692.5399 -14856.8323
208 21996.6312 -8692.5399
209 -11480.8201 21996.6312
210 44563.9461 -11480.8201
211 44642.3210 44563.9461
212 70918.7810 44642.3210
213 10603.4828 70918.7810
214 17599.4573 10603.4828
215 44955.8202 17599.4573
216 48145.1926 44955.8202
217 47862.5583 48145.1926
218 45190.9447 47862.5583
219 17392.7647 45190.9447
220 20153.9316 17392.7647
221 45426.0691 20153.9316
222 46292.9672 45426.0691
223 37483.7138 46292.9672
224 45968.4724 37483.7138
225 46965.4284 45968.4724
226 49215.5708 46965.4284
227 45896.3180 49215.5708
228 30566.5692 45896.3180
229 46053.0676 30566.5692
230 18018.9923 46053.0676
231 34530.2267 18018.9923
232 46288.1921 34530.2267
233 28814.3868 46288.1921
234 46444.9417 28814.3868
235 46960.1975 46444.9417
236 46601.6913 46960.1975
237 48785.8793 46601.6913
238 46758.4409 48785.8793
239 46836.8158 46758.4409
240 46915.1906 46836.8158
241 41818.4282 46915.1906
242 47071.9402 41818.4282
243 51803.6457 47071.9402
244 -33513.3776 51803.6457
245 -5714.2970 -33513.3776
246 -5635.9222 -5714.2970
247 39918.7656 -5635.9222
248 -5479.1725 39918.7656
249 -5400.7977 -5479.1725
250 40752.2201 -5400.7977
251 -5244.0481 40752.2201
252 29206.7092 -5244.0481
253 27332.3264 29206.7092
254 -83022.1079 27332.3264
255 48090.8128 -83022.1079
256 45767.9891 48090.8128
257 59514.2771 45767.9891
258 -4695.4244 59514.2771
259 48404.3120 -4695.4244
260 -4538.6748 48404.3120
261 -4460.2999 -4538.6748
262 53371.2611 -4460.2999
263 66045.5447 53371.2611
264 -4225.1755 66045.5447
265 48874.5609 -4225.1755
266 -4068.4259 48874.5609
267 82857.3389 -4068.4259
268 -3911.6762 82857.3389
269 46988.9064 -3911.6762
270 -3754.9266 46988.9064
271 -3676.5518 -3754.9266
272 -3598.1770 -3676.5518
273 45635.0936 -3598.1770
274 49579.9342 45635.0936
275 55660.1016 49579.9342
276 52669.9998 55660.1016
277 49815.0587 52669.9998
278 -27957.7533 49815.0587
279 38477.7869 -27957.7533
280 -2971.1785 38477.7869
281 50128.5579 -2971.1785
282 79323.0019 50128.5579
283 50285.3076 79323.0019
284 43451.9650 50285.3076
285 50442.0572 43451.9650
286 -13990.7895 50442.0572
287 51714.8857 -13990.7895
288 -2344.1800 51714.8857
289 59233.0149 -2344.1800
290 -2187.4303 59233.0149
291 62255.4666 -2187.4303
292 59639.3015 62255.4666
293 51069.0557 59639.3015
294 44872.3750 51069.0557
295 23965.2363 44872.3750
296 51304.1802 23965.2363
297 -1638.8066 51304.1802
298 51460.9298 -1638.8066
299 51539.3046 51460.9298
300 -2930.0513 51539.3046
301 33247.9534 -2930.0513
302 51774.4290 33247.9534
303 -1168.5578 51774.4290
304 90952.3650 -1168.5578
305 55822.1264 90952.3650
306 -933.4333 55822.1264
307 -6260.6167 -933.4333
308 55849.5040 -6260.6167
309 -698.3089 55849.5040
310 -619.9341 -698.3089
311 60022.2974 -619.9341
312 -463.1844 60022.2974
313 2035.8088 -463.1844
314 54553.1757 2035.8088
315 69254.8372 54553.1757
316 53278.9571 69254.8372
317 30262.6125 53278.9571
318 55582.9558 30262.6125
319 53106.8009 55582.9558
320 53185.1757 53106.8009
321 53263.5505 53185.1757
322 53341.9253 53263.5505
323 398.9385 53341.9253
324 53498.6750 398.9385
325 50903.2519 53498.6750
326 634.0630 50903.2519
327 712.4378 634.0630
328 53812.1742 712.4378
329 53890.5490 53812.1742
330 53968.9238 53890.5490
331 54047.2987 53968.9238
332 54125.6735 54047.2987
333 52477.5435 54125.6735
334 55111.2563 52477.5435
335 47153.0532 55111.2563
336 43345.8869 47153.0532
337 54517.5475 43345.8869
338 52429.4205 54517.5475
339 60691.8419 52429.4205
340 28236.5456 60691.8419
341 54831.0468 28236.5456
342 54909.4216 54831.0468
343 54987.7964 54909.4216
344 50159.7441 54987.7964
345 55144.5460 50159.7441
346 45508.0229 55144.5460
347 55301.2957 45508.0229
348 55379.6705 55301.2957
349 19254.3082 55379.6705
350 59492.5143 19254.3082
351 56070.8167 59492.5143
352 55767.3742 56070.8167
353 48970.2762 55767.3742
354 55849.9194 48970.2762
355 55928.2942 55849.9194
356 36949.9413 55928.2942
357 56085.0438 36949.9413
358 56163.4186 56085.0438
359 56241.7934 56163.4186
360 58031.6874 56241.7934
361 56398.5431 58031.6874
362 56476.9179 56398.5431
363 -47304.2524 56476.9179
364 59018.5915 -47304.2524
365 56712.0423 59018.5915
366 56790.4171 56712.0423
367 56868.7920 56790.4171
368 56947.1668 56868.7920
369 59768.6666 56947.1668
370 57103.9164 59768.6666
371 70807.9878 57103.9164
372 46377.8269 70807.9878
373 37095.0566 46377.8269
374 57417.4157 37095.0566
375 57495.7905 57417.4157
376 52666.8413 57495.7905
377 51794.5121 52666.8413
378 2339.7637 51794.5121
379 57809.2897 2339.7637
380 57887.6645 57809.2897
381 46625.4127 57887.6645
382 75947.9325 46625.4127
383 58122.7890 75947.9325
384 61096.1226 58122.7890
385 -88897.9241 61096.1226
386 70681.5594 -88897.9241
387 31901.0989 70681.5594
388 57085.9673 31901.0989
389 44937.4683 57085.9673
390 54448.0084 44937.4683
391 28112.8887 54448.0084
392 -9212.4280 28112.8887
393 80465.3240 -9212.4280
394 -164475.8443 80465.3240
395 20295.3487 -164475.8443
396 -126452.0213 20295.3487
397 26193.5371 -126452.0213
398 22511.8604 26193.5371
399 -9189.2937 22511.8604
400 1051.0932 -9189.2937
401 51968.6029 1051.0932
402 29208.1466 51968.6029
403 -80264.3333 29208.1466
404 143.0929 -80264.3333
405 23533.3418 143.0929
406 36605.7253 23533.3418
407 27744.3089 36605.7253
408 38904.9583 27744.3089
409 -60715.3444 38904.9583
410 24985.3375 -60715.3444
411 -234351.8343 24985.3375
412 -59353.2269 -234351.8343
413 29258.0085 -59353.2269
414 -174529.2171 29258.0085
415 -101353.6737 -174529.2171
416 130756.7812 -101353.6737
417 -46498.0682 130756.7812
418 -21419.1056 -46498.0682
419 -35853.6695 -21419.1056
420 -46844.9274 -35853.6695
421 -91192.4427 -46844.9274
422 -565498.6666 -91192.4427
423 -470534.4341 -565498.6666
424 -28246.5365 -470534.4341
425 -213608.3783 -28246.5365
426 -299955.6190 -213608.3783
427 -1090321.6984 -299955.6190
428 -289716.9918 -1090321.6984
429 -326935.9312 -289716.9918
430 -660177.0992 -326935.9312
> plot(z,main=paste('Residual Lag plot, lowess, and regression line'), ylab='values of Residuals', xlab='lagged values of Residuals')
> lines(lowess(z))
> abline(lm(z))
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/freestat/rcomp/tmp/71d1w1291223383.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> acf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Autocorrelation Function')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/freestat/rcomp/tmp/8b5iz1291223383.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> pacf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Partial Autocorrelation Function')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/freestat/rcomp/tmp/9b5iz1291223383.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> opar <- par(mfrow = c(2,2), oma = c(0, 0, 1.1, 0))
> plot(mylm, las = 1, sub='Residual Diagnostics')
> par(opar)
> dev.off()
null device
1
> if (n > n25) {
+ postscript(file="/var/www/html/freestat/rcomp/tmp/10mwh21291223383.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
+ plot(kp3:nmkm3,gqarr[,2], main='Goldfeld-Quandt test',ylab='2-sided p-value',xlab='breakpoint')
+ grid()
+ dev.off()
+ }
null device
1
>
> #Note: the /var/www/html/freestat/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/www/html/freestat/rcomp/createtable")
>
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Estimated Regression Equation', 1, TRUE)
> a<-table.row.end(a)
> myeq <- colnames(x)[1]
> myeq <- paste(myeq, '[t] = ', sep='')
> for (i in 1:k){
+ if (mysum$coefficients[i,1] > 0) myeq <- paste(myeq, '+', '')
+ myeq <- paste(myeq, mysum$coefficients[i,1], sep=' ')
+ if (rownames(mysum$coefficients)[i] != '(Intercept)') {
+ myeq <- paste(myeq, rownames(mysum$coefficients)[i], sep='')
+ if (rownames(mysum$coefficients)[i] != 't') myeq <- paste(myeq, '[t]', sep='')
+ }
+ }
> myeq <- paste(myeq, ' + e[t]')
> a<-table.row.start(a)
> a<-table.element(a, myeq)
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/www/html/freestat/rcomp/tmp/11qxyq1291223383.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,hyperlink('http://www.xycoon.com/ols1.htm','Multiple Linear Regression - Ordinary Least Squares',''), 6, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Variable',header=TRUE)
> a<-table.element(a,'Parameter',header=TRUE)
> a<-table.element(a,'S.D.',header=TRUE)
> a<-table.element(a,'T-STAT
H0: parameter = 0',header=TRUE)
> a<-table.element(a,'2-tail p-value',header=TRUE)
> a<-table.element(a,'1-tail p-value',header=TRUE)
> a<-table.row.end(a)
> for (i in 1:k){
+ a<-table.row.start(a)
+ a<-table.element(a,rownames(mysum$coefficients)[i],header=TRUE)
+ a<-table.element(a,mysum$coefficients[i,1])
+ a<-table.element(a, round(mysum$coefficients[i,2],6))
+ a<-table.element(a, round(mysum$coefficients[i,3],4))
+ a<-table.element(a, round(mysum$coefficients[i,4],6))
+ a<-table.element(a, round(mysum$coefficients[i,4]/2,6))
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/www/html/freestat/rcomp/tmp/12ioxb1291223383.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Regression Statistics', 2, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple R',1,TRUE)
> a<-table.element(a, sqrt(mysum$r.squared))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'R-squared',1,TRUE)
> a<-table.element(a, mysum$r.squared)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Adjusted R-squared',1,TRUE)
> a<-table.element(a, mysum$adj.r.squared)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (value)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[1])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (DF numerator)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[2])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (DF denominator)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[3])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'p-value',1,TRUE)
> a<-table.element(a, 1-pf(mysum$fstatistic[1],mysum$fstatistic[2],mysum$fstatistic[3]))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Residual Statistics', 2, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Residual Standard Deviation',1,TRUE)
> a<-table.element(a, mysum$sigma)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Sum Squared Residuals',1,TRUE)
> a<-table.element(a, sum(myerror*myerror))
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/www/html/freestat/rcomp/tmp/13p7cn1291223383.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Actuals, Interpolation, and Residuals', 4, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Time or Index', 1, TRUE)
> a<-table.element(a, 'Actuals', 1, TRUE)
> a<-table.element(a, 'Interpolation
Forecast', 1, TRUE)
> a<-table.element(a, 'Residuals
Prediction Error', 1, TRUE)
> a<-table.row.end(a)
> for (i in 1:n) {
+ a<-table.row.start(a)
+ a<-table.element(a,i, 1, TRUE)
+ a<-table.element(a,x[i])
+ a<-table.element(a,x[i]-mysum$resid[i])
+ a<-table.element(a,mysum$resid[i])
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/www/html/freestat/rcomp/tmp/140gtp1291223383.tab")
> if (n > n25) {
+ a<-table.start()
+ a<-table.row.start(a)
+ a<-table.element(a,'Goldfeld-Quandt test for Heteroskedasticity',4,TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'p-values',header=TRUE)
+ a<-table.element(a,'Alternative Hypothesis',3,header=TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'breakpoint index',header=TRUE)
+ a<-table.element(a,'greater',header=TRUE)
+ a<-table.element(a,'2-sided',header=TRUE)
+ a<-table.element(a,'less',header=TRUE)
+ a<-table.row.end(a)
+ for (mypoint in kp3:nmkm3) {
+ a<-table.row.start(a)
+ a<-table.element(a,mypoint,header=TRUE)
+ a<-table.element(a,gqarr[mypoint-kp3+1,1])
+ a<-table.element(a,gqarr[mypoint-kp3+1,2])
+ a<-table.element(a,gqarr[mypoint-kp3+1,3])
+ a<-table.row.end(a)
+ }
+ a<-table.end(a)
+ table.save(a,file="/var/www/html/freestat/rcomp/tmp/15lzsv1291223383.tab")
+ a<-table.start()
+ a<-table.row.start(a)
+ a<-table.element(a,'Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity',4,TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'Description',header=TRUE)
+ a<-table.element(a,'# significant tests',header=TRUE)
+ a<-table.element(a,'% significant tests',header=TRUE)
+ a<-table.element(a,'OK/NOK',header=TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'1% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant1)
+ a<-table.element(a,numsignificant1/numgqtests)
+ if (numsignificant1/numgqtests < 0.01) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'5% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant5)
+ a<-table.element(a,numsignificant5/numgqtests)
+ if (numsignificant5/numgqtests < 0.05) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'10% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant10)
+ a<-table.element(a,numsignificant10/numgqtests)
+ if (numsignificant10/numgqtests < 0.1) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.end(a)
+ table.save(a,file="/var/www/html/freestat/rcomp/tmp/16h9q41291223383.tab")
+ }
>
> try(system("convert tmp/1fdkq1291223383.ps tmp/1fdkq1291223383.png",intern=TRUE))
character(0)
> try(system("convert tmp/2fdkq1291223383.ps tmp/2fdkq1291223383.png",intern=TRUE))
character(0)
> try(system("convert tmp/3fdkq1291223383.ps tmp/3fdkq1291223383.png",intern=TRUE))
character(0)
> try(system("convert tmp/484kt1291223383.ps tmp/484kt1291223383.png",intern=TRUE))
character(0)
> try(system("convert tmp/584kt1291223383.ps tmp/584kt1291223383.png",intern=TRUE))
character(0)
> try(system("convert tmp/684kt1291223383.ps tmp/684kt1291223383.png",intern=TRUE))
character(0)
> try(system("convert tmp/71d1w1291223383.ps tmp/71d1w1291223383.png",intern=TRUE))
character(0)
> try(system("convert tmp/8b5iz1291223383.ps tmp/8b5iz1291223383.png",intern=TRUE))
character(0)
> try(system("convert tmp/9b5iz1291223383.ps tmp/9b5iz1291223383.png",intern=TRUE))
character(0)
> try(system("convert tmp/10mwh21291223383.ps tmp/10mwh21291223383.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
14.417 3.177 14.893