R version 2.9.0 (2009-04-17)
Copyright (C) 2009 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> x <- array(list(1
+ ,162556
+ ,1081
+ ,807
+ ,213118
+ ,6282154
+ ,1
+ ,29790
+ ,309
+ ,444
+ ,81767
+ ,4321023
+ ,1
+ ,87550
+ ,458
+ ,412
+ ,153198
+ ,4111912
+ ,0
+ ,84738
+ ,588
+ ,428
+ ,-26007
+ ,223193
+ ,1
+ ,54660
+ ,302
+ ,315
+ ,126942
+ ,1491348
+ ,1
+ ,42634
+ ,156
+ ,168
+ ,157214
+ ,1629616
+ ,0
+ ,40949
+ ,481
+ ,263
+ ,129352
+ ,1398893
+ ,1
+ ,45187
+ ,353
+ ,267
+ ,234817
+ ,1926517
+ ,1
+ ,37704
+ ,452
+ ,228
+ ,60448
+ ,983660
+ ,1
+ ,16275
+ ,109
+ ,129
+ ,47818
+ ,1443586
+ ,0
+ ,25830
+ ,115
+ ,104
+ ,245546
+ ,1073089
+ ,0
+ ,12679
+ ,110
+ ,122
+ ,48020
+ ,984885
+ ,1
+ ,18014
+ ,239
+ ,393
+ ,-1710
+ ,1405225
+ ,0
+ ,43556
+ ,247
+ ,190
+ ,32648
+ ,227132
+ ,1
+ ,24811
+ ,505
+ ,280
+ ,95350
+ ,929118
+ ,0
+ ,6575
+ ,159
+ ,63
+ ,151352
+ ,1071292
+ ,0
+ ,7123
+ ,109
+ ,102
+ ,288170
+ ,638830
+ ,1
+ ,21950
+ ,519
+ ,265
+ ,114337
+ ,856956
+ ,1
+ ,37597
+ ,248
+ ,234
+ ,37884
+ ,992426
+ ,0
+ ,17821
+ ,373
+ ,277
+ ,122844
+ ,444477
+ ,1
+ ,12988
+ ,119
+ ,73
+ ,82340
+ ,857217
+ ,1
+ ,22330
+ ,84
+ ,67
+ ,79801
+ ,711969
+ ,0
+ ,13326
+ ,102
+ ,103
+ ,165548
+ ,702380
+ ,0
+ ,16189
+ ,295
+ ,290
+ ,116384
+ ,358589
+ ,0
+ ,7146
+ ,105
+ ,83
+ ,134028
+ ,297978
+ ,0
+ ,15824
+ ,64
+ ,56
+ ,63838
+ ,585715
+ ,1
+ ,27664
+ ,282
+ ,236
+ ,74996
+ ,657954
+ ,0
+ ,11920
+ ,182
+ ,73
+ ,31080
+ ,209458
+ ,0
+ ,8568
+ ,37
+ ,34
+ ,32168
+ ,786690
+ ,0
+ ,14416
+ ,361
+ ,139
+ ,49857
+ ,439798
+ ,1
+ ,3369
+ ,28
+ ,26
+ ,87161
+ ,688779
+ ,1
+ ,11819
+ ,85
+ ,70
+ ,106113
+ ,574339
+ ,1
+ ,6984
+ ,45
+ ,40
+ ,80570
+ ,741409
+ ,1
+ ,4519
+ ,49
+ ,42
+ ,102129
+ ,597793
+ ,0
+ ,2220
+ ,22
+ ,12
+ ,301670
+ ,644190
+ ,0
+ ,18562
+ ,155
+ ,211
+ ,102313
+ ,377934
+ ,0
+ ,10327
+ ,91
+ ,74
+ ,88577
+ ,640273
+ ,1
+ ,5336
+ ,81
+ ,80
+ ,112477
+ ,697458
+ ,1
+ ,2365
+ ,79
+ ,83
+ ,191778
+ ,550608
+ ,0
+ ,4069
+ ,145
+ ,131
+ ,79804
+ ,207393
+ ,0
+ ,8636
+ ,855
+ ,203
+ ,128294
+ ,301607
+ ,0
+ ,13718
+ ,61
+ ,56
+ ,96448
+ ,345783
+ ,0
+ ,4525
+ ,226
+ ,89
+ ,93811
+ ,501749
+ ,0
+ ,6869
+ ,105
+ ,88
+ ,117520
+ ,379983
+ ,0
+ ,4628
+ ,62
+ ,39
+ ,69159
+ ,387475
+ ,1
+ ,3689
+ ,25
+ ,25
+ ,101792
+ ,377305
+ ,1
+ ,4891
+ ,217
+ ,49
+ ,210568
+ ,370837
+ ,1
+ ,7489
+ ,322
+ ,149
+ ,136996
+ ,430866
+ ,0
+ ,4901
+ ,84
+ ,58
+ ,121920
+ ,469107
+ ,0
+ ,2284
+ ,33
+ ,41
+ ,76403
+ ,194493
+ ,1
+ ,3160
+ ,108
+ ,90
+ ,108094
+ ,530670
+ ,1
+ ,4150
+ ,150
+ ,136
+ ,134759
+ ,518365
+ ,1
+ ,7285
+ ,115
+ ,97
+ ,188873
+ ,491303
+ ,1
+ ,1134
+ ,162
+ ,63
+ ,146216
+ ,527021
+ ,1
+ ,4658
+ ,158
+ ,114
+ ,156608
+ ,233773
+ ,0
+ ,2384
+ ,97
+ ,77
+ ,61348
+ ,405972
+ ,0
+ ,3748
+ ,9
+ ,6
+ ,50350
+ ,652925
+ ,0
+ ,5371
+ ,66
+ ,47
+ ,87720
+ ,446211
+ ,0
+ ,1285
+ ,107
+ ,51
+ ,99489
+ ,341340
+ ,1
+ ,9327
+ ,101
+ ,85
+ ,87419
+ ,387699
+ ,1
+ ,5565
+ ,47
+ ,43
+ ,94355
+ ,493408
+ ,0
+ ,1528
+ ,38
+ ,32
+ ,60326
+ ,146494
+ ,1
+ ,3122
+ ,34
+ ,25
+ ,94670
+ ,414462
+ ,1
+ ,7561
+ ,87
+ ,77
+ ,82425
+ ,364304
+ ,0
+ ,2675
+ ,79
+ ,54
+ ,59017
+ ,355178
+ ,0
+ ,13253
+ ,947
+ ,251
+ ,90829
+ ,357760
+ ,0
+ ,880
+ ,74
+ ,15
+ ,80791
+ ,261216
+ ,1
+ ,2053
+ ,53
+ ,44
+ ,100423
+ ,397144
+ ,0
+ ,1424
+ ,94
+ ,73
+ ,131116
+ ,374943
+ ,1
+ ,4036
+ ,63
+ ,85
+ ,100269
+ ,424898
+ ,1
+ ,3045
+ ,58
+ ,49
+ ,27330
+ ,202055
+ ,0
+ ,5119
+ ,49
+ ,38
+ ,39039
+ ,378525
+ ,0
+ ,1431
+ ,34
+ ,35
+ ,106885
+ ,310768
+ ,0
+ ,554
+ ,11
+ ,9
+ ,79285
+ ,325738
+ ,0
+ ,1975
+ ,35
+ ,34
+ ,118881
+ ,394510
+ ,1
+ ,1765
+ ,20
+ ,20
+ ,77623
+ ,247060
+ ,0
+ ,1012
+ ,47
+ ,29
+ ,114768
+ ,368078
+ ,0
+ ,810
+ ,43
+ ,11
+ ,74015
+ ,236761
+ ,0
+ ,1280
+ ,117
+ ,52
+ ,69465
+ ,312378
+ ,1
+ ,666
+ ,171
+ ,13
+ ,117869
+ ,339836
+ ,0
+ ,1380
+ ,26
+ ,29
+ ,60982
+ ,347385
+ ,1
+ ,4677
+ ,75
+ ,66
+ ,90131
+ ,426280
+ ,0
+ ,876
+ ,59
+ ,33
+ ,138971
+ ,352850
+ ,0
+ ,814
+ ,18
+ ,15
+ ,39625
+ ,301881
+ ,0
+ ,514
+ ,15
+ ,15
+ ,102725
+ ,377516
+ ,1
+ ,5692
+ ,72
+ ,68
+ ,64239
+ ,357312
+ ,0
+ ,3642
+ ,86
+ ,100
+ ,90262
+ ,458343
+ ,0
+ ,540
+ ,14
+ ,13
+ ,103960
+ ,354228
+ ,0
+ ,2099
+ ,64
+ ,45
+ ,106611
+ ,308636
+ ,0
+ ,567
+ ,11
+ ,14
+ ,103345
+ ,386212
+ ,0
+ ,2001
+ ,52
+ ,36
+ ,95551
+ ,393343
+ ,1
+ ,2949
+ ,41
+ ,40
+ ,82903
+ ,378509
+ ,0
+ ,2253
+ ,99
+ ,68
+ ,63593
+ ,452469
+ ,1
+ ,6533
+ ,75
+ ,29
+ ,126910
+ ,364839
+ ,0
+ ,1889
+ ,45
+ ,43
+ ,37527
+ ,358649
+ ,1
+ ,3055
+ ,43
+ ,30
+ ,60247
+ ,376641
+ ,0
+ ,272
+ ,8
+ ,9
+ ,112995
+ ,429112
+ ,1
+ ,1414
+ ,198
+ ,22
+ ,70184
+ ,330546
+ ,0
+ ,2564
+ ,22
+ ,19
+ ,130140
+ ,403560
+ ,1
+ ,1383
+ ,11
+ ,9
+ ,73221
+ ,317892
+ ,1
+ ,1261
+ ,33
+ ,31
+ ,76114
+ ,307528
+ ,0
+ ,975
+ ,23
+ ,19
+ ,90534
+ ,235133
+ ,0
+ ,3366
+ ,80
+ ,55
+ ,108479
+ ,299243
+ ,0
+ ,576
+ ,18
+ ,8
+ ,113761
+ ,314073
+ ,0
+ ,1686
+ ,40
+ ,28
+ ,68696
+ ,368186
+ ,0
+ ,746
+ ,23
+ ,29
+ ,71561
+ ,269661
+ ,1
+ ,3192
+ ,60
+ ,48
+ ,59831
+ ,125390
+ ,1
+ ,2045
+ ,20
+ ,16
+ ,97890
+ ,510834
+ ,0
+ ,5702
+ ,61
+ ,47
+ ,101481
+ ,321896
+ ,1
+ ,1932
+ ,36
+ ,20
+ ,72954
+ ,249898
+ ,0
+ ,936
+ ,30
+ ,22
+ ,67939
+ ,408881
+ ,1
+ ,3437
+ ,47
+ ,33
+ ,48022
+ ,158492
+ ,0
+ ,5131
+ ,71
+ ,44
+ ,86111
+ ,292154
+ ,1
+ ,2397
+ ,14
+ ,13
+ ,74020
+ ,289513
+ ,1
+ ,1389
+ ,9
+ ,6
+ ,57530
+ ,378049
+ ,0
+ ,1503
+ ,39
+ ,35
+ ,56364
+ ,343466
+ ,0
+ ,402
+ ,26
+ ,8
+ ,84990
+ ,332743
+ ,0
+ ,2239
+ ,21
+ ,17
+ ,88590
+ ,442882
+ ,1
+ ,2234
+ ,16
+ ,11
+ ,77200
+ ,214215
+ ,0
+ ,837
+ ,69
+ ,21
+ ,61262
+ ,315688
+ ,0
+ ,10579
+ ,92
+ ,92
+ ,110309
+ ,375195
+ ,0
+ ,875
+ ,14
+ ,12
+ ,67000
+ ,334280
+ ,0
+ ,1585
+ ,107
+ ,112
+ ,93099
+ ,355864
+ ,1
+ ,1659
+ ,29
+ ,25
+ ,107577
+ ,480382
+ ,1
+ ,2647
+ ,37
+ ,17
+ ,62920
+ ,353058
+ ,1
+ ,3294
+ ,23
+ ,23
+ ,75832
+ ,217193
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,94
+ ,7
+ ,10
+ ,60793
+ ,314533
+ ,0
+ ,422
+ ,28
+ ,23
+ ,57935
+ ,318056
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,34
+ ,8
+ ,7
+ ,60630
+ ,314353
+ ,0
+ ,1558
+ ,63
+ ,25
+ ,55637
+ ,369448
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,43
+ ,3
+ ,20
+ ,60887
+ ,312846
+ ,0
+ ,645
+ ,5
+ ,4
+ ,60720
+ ,312075
+ ,0
+ ,316
+ ,9
+ ,4
+ ,60505
+ ,315009
+ ,0
+ ,115
+ ,13
+ ,10
+ ,60945
+ ,318903
+ ,0
+ ,5
+ ,2
+ ,1
+ ,60720
+ ,314887
+ ,0
+ ,897
+ ,5
+ ,4
+ ,60720
+ ,314913
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,389
+ ,14
+ ,8
+ ,58990
+ ,325506
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1002
+ ,15
+ ,11
+ ,56750
+ ,298568
+ ,0
+ ,36
+ ,3
+ ,4
+ ,60894
+ ,315834
+ ,0
+ ,460
+ ,15
+ ,15
+ ,63346
+ ,329784
+ ,0
+ ,309
+ ,11
+ ,9
+ ,56535
+ ,312878
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,9
+ ,6
+ ,7
+ ,60835
+ ,314987
+ ,0
+ ,271
+ ,2
+ ,2
+ ,60720
+ ,325249
+ ,0
+ ,14
+ ,1
+ ,0
+ ,61016
+ ,315877
+ ,0
+ ,520
+ ,10
+ ,7
+ ,58650
+ ,291650
+ ,0
+ ,1766
+ ,73
+ ,46
+ ,60438
+ ,305959
+ ,0
+ ,0
+ ,0
+ ,5
+ ,60720
+ ,315380
+ ,0
+ ,458
+ ,11
+ ,7
+ ,58625
+ ,297765
+ ,0
+ ,20
+ ,3
+ ,2
+ ,60938
+ ,315245
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,98
+ ,2
+ ,2
+ ,61490
+ ,315236
+ ,0
+ ,405
+ ,7
+ ,5
+ ,60845
+ ,336425
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,483
+ ,27
+ ,7
+ ,60830
+ ,306268
+ ,0
+ ,454
+ ,51
+ ,24
+ ,63261
+ ,302187
+ ,0
+ ,47
+ ,3
+ ,1
+ ,60720
+ ,314882
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,757
+ ,19
+ ,18
+ ,45689
+ ,382712
+ ,0
+ ,4655
+ ,393
+ ,55
+ ,60720
+ ,341570
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,36
+ ,4
+ ,3
+ ,61564
+ ,312412
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,203
+ ,9
+ ,9
+ ,61938
+ ,309596
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,126
+ ,10
+ ,8
+ ,60951
+ ,315547
+ ,0
+ ,400
+ ,152
+ ,113
+ ,60720
+ ,313267
+ ,1
+ ,71
+ ,1
+ ,0
+ ,60745
+ ,316176
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,972
+ ,34
+ ,19
+ ,71642
+ ,359335
+ ,1
+ ,531
+ ,10
+ ,11
+ ,71641
+ ,330068
+ ,0
+ ,2461
+ ,57
+ ,25
+ ,55792
+ ,314289
+ ,1
+ ,378
+ ,52
+ ,16
+ ,71873
+ ,297413
+ ,1
+ ,23
+ ,5
+ ,5
+ ,62555
+ ,314806
+ ,1
+ ,638
+ ,14
+ ,11
+ ,60370
+ ,333210
+ ,1
+ ,2300
+ ,29
+ ,23
+ ,64873
+ ,352108
+ ,0
+ ,149
+ ,5
+ ,6
+ ,62041
+ ,313332
+ ,0
+ ,226
+ ,5
+ ,5
+ ,65745
+ ,291787
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,275
+ ,4
+ ,7
+ ,59500
+ ,318745
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,141
+ ,6
+ ,7
+ ,61630
+ ,315366
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,28
+ ,2
+ ,3
+ ,60890
+ ,315688
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,4980
+ ,91
+ ,89
+ ,113521
+ ,409642
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,472
+ ,20
+ ,19
+ ,80045
+ ,269587
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,203
+ ,27
+ ,12
+ ,50804
+ ,300962
+ ,1
+ ,496
+ ,17
+ ,12
+ ,87390
+ ,325479
+ ,1
+ ,10
+ ,2
+ ,5
+ ,61656
+ ,316155
+ ,1
+ ,63
+ ,4
+ ,2
+ ,65688
+ ,318574
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,1136
+ ,32
+ ,26
+ ,48522
+ ,343613
+ ,1
+ ,265
+ ,31
+ ,3
+ ,60720
+ ,306948
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,267
+ ,32
+ ,11
+ ,57640
+ ,330059
+ ,0
+ ,474
+ ,20
+ ,10
+ ,61977
+ ,288985
+ ,0
+ ,534
+ ,7
+ ,5
+ ,62620
+ ,304485
+ ,0
+ ,0
+ ,0
+ ,2
+ ,60720
+ ,315380
+ ,0
+ ,15
+ ,8
+ ,6
+ ,60831
+ ,315688
+ ,0
+ ,397
+ ,28
+ ,7
+ ,60646
+ ,317736
+ ,0
+ ,0
+ ,0
+ ,2
+ ,60720
+ ,315380
+ ,0
+ ,1866
+ ,29
+ ,28
+ ,56225
+ ,322331
+ ,0
+ ,288
+ ,4
+ ,3
+ ,60510
+ ,296656
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,3
+ ,2
+ ,1
+ ,60698
+ ,315354
+ ,0
+ ,468
+ ,21
+ ,20
+ ,60720
+ ,312161
+ ,0
+ ,20
+ ,2
+ ,1
+ ,60805
+ ,315576
+ ,0
+ ,278
+ ,26
+ ,22
+ ,61404
+ ,314922
+ ,0
+ ,61
+ ,14
+ ,9
+ ,60720
+ ,314551
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,192
+ ,4
+ ,2
+ ,65276
+ ,312339
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,317
+ ,9
+ ,7
+ ,63915
+ ,298700
+ ,0
+ ,738
+ ,10
+ ,9
+ ,60720
+ ,321376
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,368
+ ,17
+ ,13
+ ,61686
+ ,303230
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,2
+ ,1
+ ,0
+ ,60743
+ ,315487
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,53
+ ,6
+ ,6
+ ,60349
+ ,315793
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,94
+ ,3
+ ,3
+ ,61360
+ ,312887
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,24
+ ,8
+ ,7
+ ,59818
+ ,315637
+ ,0
+ ,2332
+ ,4
+ ,2
+ ,72680
+ ,324385
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,131
+ ,11
+ ,15
+ ,61808
+ ,308989
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,206
+ ,9
+ ,9
+ ,53110
+ ,296702
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,167
+ ,2
+ ,1
+ ,64245
+ ,307322
+ ,0
+ ,622
+ ,73
+ ,38
+ ,73007
+ ,304376
+ ,0
+ ,2328
+ ,94
+ ,57
+ ,82732
+ ,253588
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,365
+ ,8
+ ,7
+ ,54820
+ ,309560
+ ,0
+ ,364
+ ,35
+ ,26
+ ,47705
+ ,298466
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,226
+ ,12
+ ,13
+ ,72835
+ ,343929
+ ,0
+ ,307
+ ,15
+ ,10
+ ,58856
+ ,331955
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,188
+ ,11
+ ,9
+ ,77655
+ ,381180
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,138
+ ,6
+ ,26
+ ,69817
+ ,331420
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,125
+ ,12
+ ,19
+ ,60798
+ ,310201
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,282
+ ,30
+ ,12
+ ,62452
+ ,320016
+ ,0
+ ,335
+ ,33
+ ,23
+ ,64175
+ ,320398
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,1324
+ ,117
+ ,29
+ ,67440
+ ,291841
+ ,0
+ ,176
+ ,28
+ ,8
+ ,68136
+ ,310670
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,249
+ ,72
+ ,26
+ ,56726
+ ,313491
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,333
+ ,13
+ ,9
+ ,70811
+ ,331323
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,601
+ ,6
+ ,5
+ ,60720
+ ,319210
+ ,0
+ ,30
+ ,4
+ ,3
+ ,62045
+ ,318098
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,249
+ ,62
+ ,13
+ ,54323
+ ,292754
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,165
+ ,24
+ ,12
+ ,62841
+ ,325176
+ ,0
+ ,453
+ ,21
+ ,19
+ ,81125
+ ,365959
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,53
+ ,14
+ ,10
+ ,59506
+ ,302409
+ ,1
+ ,382
+ ,21
+ ,9
+ ,59365
+ ,340968
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,9
+ ,60720
+ ,315380
+ ,1
+ ,30
+ ,4
+ ,4
+ ,60798
+ ,313164
+ ,0
+ ,290
+ ,2
+ ,1
+ ,58790
+ ,301164
+ ,0
+ ,0
+ ,0
+ ,1
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,366
+ ,53
+ ,14
+ ,61808
+ ,344425
+ ,0
+ ,2
+ ,9
+ ,12
+ ,60735
+ ,315394
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,209
+ ,13
+ ,19
+ ,64016
+ ,316647
+ ,0
+ ,384
+ ,22
+ ,17
+ ,54683
+ ,309836
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,365
+ ,83
+ ,32
+ ,87192
+ ,346611
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,49
+ ,8
+ ,14
+ ,64107
+ ,322031
+ ,0
+ ,3
+ ,4
+ ,8
+ ,60761
+ ,315656
+ ,0
+ ,133
+ ,14
+ ,4
+ ,65990
+ ,339445
+ ,0
+ ,32
+ ,1
+ ,0
+ ,59988
+ ,314964
+ ,0
+ ,368
+ ,17
+ ,20
+ ,61167
+ ,297141
+ ,0
+ ,1
+ ,6
+ ,5
+ ,60719
+ ,315372
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,22
+ ,2
+ ,1
+ ,60722
+ ,312502
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,96
+ ,5
+ ,4
+ ,60379
+ ,313729
+ ,0
+ ,1
+ ,2
+ ,1
+ ,60727
+ ,315388
+ ,0
+ ,314
+ ,5
+ ,4
+ ,60720
+ ,315371
+ ,0
+ ,844
+ ,78
+ ,20
+ ,60925
+ ,296139
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,26
+ ,1
+ ,1
+ ,60896
+ ,313880
+ ,0
+ ,125
+ ,13
+ ,10
+ ,59734
+ ,317698
+ ,0
+ ,304
+ ,15
+ ,12
+ ,62969
+ ,295580
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,621
+ ,48
+ ,13
+ ,60720
+ ,308256
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,119
+ ,6
+ ,3
+ ,59118
+ ,303677
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1595
+ ,17
+ ,10
+ ,60720
+ ,319369
+ ,0
+ ,312
+ ,14
+ ,3
+ ,58598
+ ,318690
+ ,0
+ ,60
+ ,10
+ ,7
+ ,61124
+ ,314049
+ ,0
+ ,587
+ ,12
+ ,10
+ ,59595
+ ,325699
+ ,0
+ ,135
+ ,2
+ ,1
+ ,62065
+ ,314210
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,514
+ ,52
+ ,15
+ ,78780
+ ,322378
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,4
+ ,4
+ ,60722
+ ,315398
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,1763
+ ,24
+ ,28
+ ,61600
+ ,308336
+ ,0
+ ,180
+ ,11
+ ,9
+ ,59635
+ ,316386
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,218
+ ,21
+ ,7
+ ,60720
+ ,315553
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,448
+ ,40
+ ,7
+ ,59781
+ ,323361
+ ,0
+ ,227
+ ,9
+ ,7
+ ,76644
+ ,336639
+ ,0
+ ,174
+ ,1
+ ,3
+ ,64820
+ ,307424
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,121
+ ,24
+ ,11
+ ,56178
+ ,295370
+ ,0
+ ,607
+ ,11
+ ,7
+ ,60436
+ ,322340
+ ,0
+ ,2212
+ ,14
+ ,10
+ ,60720
+ ,319864
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,530
+ ,60
+ ,18
+ ,73433
+ ,317291
+ ,0
+ ,571
+ ,80
+ ,14
+ ,41477
+ ,280398
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,78
+ ,16
+ ,12
+ ,62700
+ ,317330
+ ,0
+ ,2489
+ ,40
+ ,29
+ ,67804
+ ,238125
+ ,0
+ ,131
+ ,6
+ ,3
+ ,59661
+ ,327071
+ ,0
+ ,923
+ ,8
+ ,6
+ ,58620
+ ,309038
+ ,0
+ ,72
+ ,3
+ ,3
+ ,60398
+ ,314210
+ ,0
+ ,572
+ ,16
+ ,8
+ ,58580
+ ,307930
+ ,0
+ ,397
+ ,10
+ ,10
+ ,62710
+ ,322327
+ ,0
+ ,450
+ ,8
+ ,6
+ ,59325
+ ,292136
+ ,0
+ ,622
+ ,7
+ ,8
+ ,60950
+ ,263276
+ ,0
+ ,694
+ ,8
+ ,6
+ ,68060
+ ,367655
+ ,1
+ ,3425
+ ,12
+ ,9
+ ,83620
+ ,283910
+ ,0
+ ,562
+ ,13
+ ,8
+ ,58456
+ ,283587
+ ,0
+ ,4917
+ ,42
+ ,26
+ ,52811
+ ,243650
+ ,1
+ ,1442
+ ,118
+ ,239
+ ,121173
+ ,438493
+ ,0
+ ,529
+ ,9
+ ,7
+ ,63870
+ ,296261
+ ,1
+ ,2126
+ ,138
+ ,41
+ ,21001
+ ,230621
+ ,0
+ ,1061
+ ,5
+ ,3
+ ,70415
+ ,304252
+ ,0
+ ,776
+ ,9
+ ,8
+ ,64230
+ ,333505
+ ,0
+ ,611
+ ,8
+ ,6
+ ,59190
+ ,296919
+ ,1
+ ,1526
+ ,25
+ ,21
+ ,69351
+ ,278990
+ ,0
+ ,592
+ ,7
+ ,7
+ ,64270
+ ,276898
+ ,0
+ ,1182
+ ,13
+ ,11
+ ,70694
+ ,327007
+ ,0
+ ,621
+ ,16
+ ,11
+ ,68005
+ ,317046
+ ,0
+ ,989
+ ,11
+ ,12
+ ,58930
+ ,304555
+ ,0
+ ,438
+ ,11
+ ,9
+ ,58320
+ ,298096
+ ,0
+ ,726
+ ,3
+ ,3
+ ,69980
+ ,231861
+ ,0
+ ,1303
+ ,61
+ ,57
+ ,69863
+ ,309422
+ ,1
+ ,7419
+ ,29
+ ,21
+ ,63255
+ ,286963
+ ,1
+ ,1164
+ ,17
+ ,15
+ ,57320
+ ,269753
+ ,1
+ ,3310
+ ,33
+ ,32
+ ,75230
+ ,448243
+ ,0
+ ,1920
+ ,15
+ ,11
+ ,79420
+ ,165404
+ ,0
+ ,965
+ ,3
+ ,2
+ ,73490
+ ,204325
+ ,0
+ ,3256
+ ,66
+ ,23
+ ,35250
+ ,407159
+ ,1
+ ,1135
+ ,17
+ ,20
+ ,62285
+ ,290476
+ ,0
+ ,1270
+ ,26
+ ,24
+ ,69206
+ ,275311
+ ,0
+ ,661
+ ,3
+ ,1
+ ,65920
+ ,246541
+ ,0
+ ,1013
+ ,2
+ ,1
+ ,69770
+ ,253468
+ ,0
+ ,2844
+ ,67
+ ,74
+ ,72683
+ ,240897
+ ,1
+ ,11528
+ ,70
+ ,68
+ ,-14545
+ ,-83265
+ ,0
+ ,6526
+ ,26
+ ,20
+ ,55830
+ ,-42143
+ ,0
+ ,2264
+ ,24
+ ,20
+ ,55174
+ ,272713
+ ,1
+ ,5109
+ ,97
+ ,82
+ ,67038
+ ,215362
+ ,0
+ ,3999
+ ,30
+ ,21
+ ,51252
+ ,42754
+ ,0
+ ,35624
+ ,223
+ ,244
+ ,157278
+ ,306275
+ ,0
+ ,9252
+ ,48
+ ,32
+ ,79510
+ ,253537
+ ,0
+ ,15236
+ ,90
+ ,86
+ ,77440
+ ,372631
+ ,0
+ ,18073
+ ,180
+ ,69
+ ,27284
+ ,-7170)
+ ,dim=c(6
+ ,431)
+ ,dimnames=list(c('Group'
+ ,'Costs'
+ ,'Trades'
+ ,'Orders'
+ ,'Dividends'
+ ,'Wealth')
+ ,1:431))
> y <- array(NA,dim=c(6,431),dimnames=list(c('Group','Costs','Trades','Orders','Dividends','Wealth'),1:431))
> for (i in 1:dim(x)[1])
+ {
+ for (j in 1:dim(x)[2])
+ {
+ y[i,j] <- as.numeric(x[i,j])
+ }
+ }
> par3 = 'No Linear Trend'
> par2 = 'Do not include Seasonal Dummies'
> par1 = '6'
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Prof. Dr. P. Wessa
> #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #Technical description: Write here your technical program description (don't use hard returns!)
> library(lattice)
> library(lmtest)
Loading required package: zoo
Attaching package: 'zoo'
The following object(s) are masked from package:base :
as.Date.numeric
> n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test
> par1 <- as.numeric(par1)
> x <- t(y)
> k <- length(x[1,])
> n <- length(x[,1])
> x1 <- cbind(x[,par1], x[,1:k!=par1])
> mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1])
> colnames(x1) <- mycolnames #colnames(x)[par1]
> x <- x1
> if (par3 == 'First Differences'){
+ x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep='')))
+ for (i in 1:n-1) {
+ for (j in 1:k) {
+ x2[i,j] <- x[i+1,j] - x[i,j]
+ }
+ }
+ x <- x2
+ }
> if (par2 == 'Include Monthly Dummies'){
+ x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep ='')))
+ for (i in 1:11){
+ x2[seq(i,n,12),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> if (par2 == 'Include Quarterly Dummies'){
+ x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep ='')))
+ for (i in 1:3){
+ x2[seq(i,n,4),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> k <- length(x[1,])
> if (par3 == 'Linear Trend'){
+ x <- cbind(x, c(1:n))
+ colnames(x)[k+1] <- 't'
+ }
> x
Wealth Group Costs Trades Orders Dividends
1 6282154 1 162556 1081 807 213118
2 4321023 1 29790 309 444 81767
3 4111912 1 87550 458 412 153198
4 223193 0 84738 588 428 -26007
5 1491348 1 54660 302 315 126942
6 1629616 1 42634 156 168 157214
7 1398893 0 40949 481 263 129352
8 1926517 1 45187 353 267 234817
9 983660 1 37704 452 228 60448
10 1443586 1 16275 109 129 47818
11 1073089 0 25830 115 104 245546
12 984885 0 12679 110 122 48020
13 1405225 1 18014 239 393 -1710
14 227132 0 43556 247 190 32648
15 929118 1 24811 505 280 95350
16 1071292 0 6575 159 63 151352
17 638830 0 7123 109 102 288170
18 856956 1 21950 519 265 114337
19 992426 1 37597 248 234 37884
20 444477 0 17821 373 277 122844
21 857217 1 12988 119 73 82340
22 711969 1 22330 84 67 79801
23 702380 0 13326 102 103 165548
24 358589 0 16189 295 290 116384
25 297978 0 7146 105 83 134028
26 585715 0 15824 64 56 63838
27 657954 1 27664 282 236 74996
28 209458 0 11920 182 73 31080
29 786690 0 8568 37 34 32168
30 439798 0 14416 361 139 49857
31 688779 1 3369 28 26 87161
32 574339 1 11819 85 70 106113
33 741409 1 6984 45 40 80570
34 597793 1 4519 49 42 102129
35 644190 0 2220 22 12 301670
36 377934 0 18562 155 211 102313
37 640273 0 10327 91 74 88577
38 697458 1 5336 81 80 112477
39 550608 1 2365 79 83 191778
40 207393 0 4069 145 131 79804
41 301607 0 8636 855 203 128294
42 345783 0 13718 61 56 96448
43 501749 0 4525 226 89 93811
44 379983 0 6869 105 88 117520
45 387475 0 4628 62 39 69159
46 377305 1 3689 25 25 101792
47 370837 1 4891 217 49 210568
48 430866 1 7489 322 149 136996
49 469107 0 4901 84 58 121920
50 194493 0 2284 33 41 76403
51 530670 1 3160 108 90 108094
52 518365 1 4150 150 136 134759
53 491303 1 7285 115 97 188873
54 527021 1 1134 162 63 146216
55 233773 1 4658 158 114 156608
56 405972 0 2384 97 77 61348
57 652925 0 3748 9 6 50350
58 446211 0 5371 66 47 87720
59 341340 0 1285 107 51 99489
60 387699 1 9327 101 85 87419
61 493408 1 5565 47 43 94355
62 146494 0 1528 38 32 60326
63 414462 1 3122 34 25 94670
64 364304 1 7561 87 77 82425
65 355178 0 2675 79 54 59017
66 357760 0 13253 947 251 90829
67 261216 0 880 74 15 80791
68 397144 1 2053 53 44 100423
69 374943 0 1424 94 73 131116
70 424898 1 4036 63 85 100269
71 202055 1 3045 58 49 27330
72 378525 0 5119 49 38 39039
73 310768 0 1431 34 35 106885
74 325738 0 554 11 9 79285
75 394510 0 1975 35 34 118881
76 247060 1 1765 20 20 77623
77 368078 0 1012 47 29 114768
78 236761 0 810 43 11 74015
79 312378 0 1280 117 52 69465
80 339836 1 666 171 13 117869
81 347385 0 1380 26 29 60982
82 426280 1 4677 75 66 90131
83 352850 0 876 59 33 138971
84 301881 0 814 18 15 39625
85 377516 0 514 15 15 102725
86 357312 1 5692 72 68 64239
87 458343 0 3642 86 100 90262
88 354228 0 540 14 13 103960
89 308636 0 2099 64 45 106611
90 386212 0 567 11 14 103345
91 393343 0 2001 52 36 95551
92 378509 1 2949 41 40 82903
93 452469 0 2253 99 68 63593
94 364839 1 6533 75 29 126910
95 358649 0 1889 45 43 37527
96 376641 1 3055 43 30 60247
97 429112 0 272 8 9 112995
98 330546 1 1414 198 22 70184
99 403560 0 2564 22 19 130140
100 317892 1 1383 11 9 73221
101 307528 1 1261 33 31 76114
102 235133 0 975 23 19 90534
103 299243 0 3366 80 55 108479
104 314073 0 576 18 8 113761
105 368186 0 1686 40 28 68696
106 269661 0 746 23 29 71561
107 125390 1 3192 60 48 59831
108 510834 1 2045 20 16 97890
109 321896 0 5702 61 47 101481
110 249898 1 1932 36 20 72954
111 408881 0 936 30 22 67939
112 158492 1 3437 47 33 48022
113 292154 0 5131 71 44 86111
114 289513 1 2397 14 13 74020
115 378049 1 1389 9 6 57530
116 343466 0 1503 39 35 56364
117 332743 0 402 26 8 84990
118 442882 0 2239 21 17 88590
119 214215 1 2234 16 11 77200
120 315688 0 837 69 21 61262
121 375195 0 10579 92 92 110309
122 334280 0 875 14 12 67000
123 355864 0 1585 107 112 93099
124 480382 1 1659 29 25 107577
125 353058 1 2647 37 17 62920
126 217193 1 3294 23 23 75832
127 315380 0 0 0 0 60720
128 314533 0 94 7 10 60793
129 318056 0 422 28 23 57935
130 315380 0 0 0 0 60720
131 314353 0 34 8 7 60630
132 369448 0 1558 63 25 55637
133 315380 0 0 0 0 60720
134 312846 0 43 3 20 60887
135 312075 0 645 5 4 60720
136 315009 0 316 9 4 60505
137 318903 0 115 13 10 60945
138 314887 0 5 2 1 60720
139 314913 0 897 5 4 60720
140 315380 0 0 0 0 60720
141 325506 0 389 14 8 58990
142 315380 0 0 0 0 60720
143 298568 0 1002 15 11 56750
144 315834 0 36 3 4 60894
145 329784 0 460 15 15 63346
146 312878 0 309 11 9 56535
147 315380 0 0 0 0 60720
148 314987 0 9 6 7 60835
149 325249 0 271 2 2 60720
150 315877 0 14 1 0 61016
151 291650 0 520 10 7 58650
152 305959 0 1766 73 46 60438
153 315380 0 0 0 5 60720
154 297765 0 458 11 7 58625
155 315245 0 20 3 2 60938
156 315380 0 0 0 0 60720
157 315380 0 0 0 0 60720
158 315236 0 98 2 2 61490
159 336425 0 405 7 5 60845
160 315380 0 0 0 0 60720
161 315380 0 0 0 0 60720
162 315380 0 0 0 0 60720
163 315380 0 0 0 0 60720
164 306268 0 483 27 7 60830
165 302187 0 454 51 24 63261
166 314882 0 47 3 1 60720
167 315380 0 0 0 0 60720
168 382712 0 757 19 18 45689
169 341570 0 4655 393 55 60720
170 315380 0 0 0 0 60720
171 315380 0 0 0 0 60720
172 312412 0 36 4 3 61564
173 315380 0 0 0 0 60720
174 309596 0 203 9 9 61938
175 315380 0 0 0 0 60720
176 315547 0 126 10 8 60951
177 313267 0 400 152 113 60720
178 316176 1 71 1 0 60745
179 315380 0 0 0 0 60720
180 315380 0 0 0 0 60720
181 359335 0 972 34 19 71642
182 330068 1 531 10 11 71641
183 314289 0 2461 57 25 55792
184 297413 1 378 52 16 71873
185 314806 1 23 5 5 62555
186 333210 1 638 14 11 60370
187 352108 1 2300 29 23 64873
188 313332 0 149 5 6 62041
189 291787 0 226 5 5 65745
190 315380 0 0 0 0 60720
191 318745 0 275 4 7 59500
192 315380 0 0 0 0 60720
193 315366 0 141 6 7 61630
194 315380 0 0 0 0 60720
195 315688 0 28 2 3 60890
196 315380 1 0 0 0 60720
197 409642 1 4980 91 89 113521
198 315380 1 0 0 0 60720
199 315380 1 0 0 0 60720
200 269587 1 472 20 19 80045
201 315380 1 0 0 0 60720
202 315380 1 0 0 0 60720
203 315380 1 0 0 0 60720
204 300962 1 203 27 12 50804
205 325479 1 496 17 12 87390
206 316155 1 10 2 5 61656
207 318574 1 63 4 2 65688
208 315380 1 0 0 0 60720
209 343613 1 1136 32 26 48522
210 306948 1 265 31 3 60720
211 315380 0 0 0 0 60720
212 315380 0 0 0 0 60720
213 330059 0 267 32 11 57640
214 288985 0 474 20 10 61977
215 304485 0 534 7 5 62620
216 315380 0 0 0 2 60720
217 315688 0 15 8 6 60831
218 317736 0 397 28 7 60646
219 315380 0 0 0 2 60720
220 322331 0 1866 29 28 56225
221 296656 0 288 4 3 60510
222 315380 0 0 0 0 60720
223 315354 0 3 2 1 60698
224 312161 0 468 21 20 60720
225 315576 0 20 2 1 60805
226 314922 0 278 26 22 61404
227 314551 0 61 14 9 60720
228 315380 0 0 0 0 60720
229 312339 0 192 4 2 65276
230 315380 0 0 0 0 60720
231 298700 0 317 9 7 63915
232 321376 0 738 10 9 60720
233 315380 0 0 0 0 60720
234 303230 0 368 17 13 61686
235 315380 0 0 0 0 60720
236 315487 0 2 1 0 60743
237 315380 0 0 0 0 60720
238 315793 0 53 6 6 60349
239 315380 0 0 0 0 60720
240 315380 0 0 0 0 60720
241 315380 0 0 0 0 60720
242 312887 0 94 3 3 61360
243 315380 0 0 0 0 60720
244 315637 0 24 8 7 59818
245 324385 0 2332 4 2 72680
246 315380 1 0 0 0 60720
247 315380 1 0 0 0 60720
248 308989 0 131 11 15 61808
249 315380 1 0 0 0 60720
250 315380 1 0 0 0 60720
251 296702 0 206 9 9 53110
252 315380 1 0 0 0 60720
253 307322 0 167 2 1 64245
254 304376 0 622 73 38 73007
255 253588 0 2328 94 57 82732
256 315380 0 0 0 0 60720
257 309560 0 365 8 7 54820
258 298466 0 364 35 26 47705
259 315380 1 0 0 0 60720
260 315380 0 0 0 0 60720
261 315380 1 0 0 0 60720
262 315380 1 0 0 0 60720
263 343929 0 226 12 13 72835
264 331955 0 307 15 10 58856
265 315380 1 0 0 0 60720
266 315380 0 0 0 0 60720
267 315380 1 0 0 0 60720
268 381180 0 188 11 9 77655
269 315380 1 0 0 0 60720
270 331420 0 138 6 26 69817
271 315380 1 0 0 0 60720
272 315380 1 0 0 0 60720
273 315380 1 0 0 0 60720
274 310201 0 125 12 19 60798
275 315380 0 0 0 0 60720
276 320016 0 282 30 12 62452
277 320398 0 335 33 23 64175
278 315380 0 0 0 0 60720
279 291841 1 1324 117 29 67440
280 310670 0 176 28 8 68136
281 315380 1 0 0 0 60720
282 315380 0 0 0 0 60720
283 313491 0 249 72 26 56726
284 315380 0 0 0 0 60720
285 331323 0 333 13 9 70811
286 315380 0 0 0 0 60720
287 319210 1 601 6 5 60720
288 318098 0 30 4 3 62045
289 315380 1 0 0 0 60720
290 292754 0 249 62 13 54323
291 315380 1 0 0 0 60720
292 325176 0 165 24 12 62841
293 365959 0 453 21 19 81125
294 315380 0 0 0 0 60720
295 302409 0 53 14 10 59506
296 340968 1 382 21 9 59365
297 315380 0 0 0 0 60720
298 315380 1 0 0 0 60720
299 315380 0 0 0 0 60720
300 315380 0 0 0 9 60720
301 313164 1 30 4 4 60798
302 301164 0 290 2 1 58790
303 315380 0 0 0 1 60720
304 315380 1 0 0 0 60720
305 344425 0 366 53 14 61808
306 315394 0 2 9 12 60735
307 315380 1 0 0 0 60720
308 316647 1 209 13 19 64016
309 309836 0 384 22 17 54683
310 315380 1 0 0 0 60720
311 315380 1 0 0 0 60720
312 346611 0 365 83 32 87192
313 315380 1 0 0 0 60720
314 322031 1 49 8 14 64107
315 315656 0 3 4 8 60761
316 339445 0 133 14 4 65990
317 314964 0 32 1 0 59988
318 297141 0 368 17 20 61167
319 315372 0 1 6 5 60719
320 315380 0 0 0 0 60720
321 315380 0 0 0 0 60720
322 315380 0 0 0 0 60720
323 315380 0 0 0 0 60720
324 315380 1 0 0 0 60720
325 315380 0 0 0 0 60720
326 312502 0 22 2 1 60722
327 315380 1 0 0 0 60720
328 315380 1 0 0 0 60720
329 315380 0 0 0 0 60720
330 315380 0 0 0 0 60720
331 315380 0 0 0 0 60720
332 315380 0 0 0 0 60720
333 315380 0 0 0 0 60720
334 313729 0 96 5 4 60379
335 315388 0 1 2 1 60727
336 315371 0 314 5 4 60720
337 296139 0 844 78 20 60925
338 315380 0 0 0 0 60720
339 313880 0 26 1 1 60896
340 317698 0 125 13 10 59734
341 295580 0 304 15 12 62969
342 315380 0 0 0 0 60720
343 315380 0 0 0 0 60720
344 315380 0 0 0 0 60720
345 308256 0 621 48 13 60720
346 315380 0 0 0 0 60720
347 303677 0 119 6 3 59118
348 315380 0 0 0 0 60720
349 315380 0 0 0 0 60720
350 319369 0 1595 17 10 60720
351 318690 0 312 14 3 58598
352 314049 0 60 10 7 61124
353 325699 0 587 12 10 59595
354 314210 0 135 2 1 62065
355 315380 0 0 0 0 60720
356 315380 0 0 0 0 60720
357 322378 0 514 52 15 78780
358 315380 0 0 0 0 60720
359 315380 0 0 0 0 60720
360 315380 0 0 0 0 60720
361 315398 0 1 4 4 60722
362 315380 0 0 0 0 60720
363 315380 0 0 0 0 60720
364 308336 1 1763 24 28 61600
365 316386 0 180 11 9 59635
366 315380 0 0 0 0 60720
367 315380 0 0 0 0 60720
368 315380 0 0 0 0 60720
369 315380 0 0 0 0 60720
370 315553 0 218 21 7 60720
371 315380 0 0 0 0 60720
372 323361 0 448 40 7 59781
373 336639 0 227 9 7 76644
374 307424 0 174 1 3 64820
375 315380 0 0 0 0 60720
376 315380 0 0 0 0 60720
377 295370 0 121 24 11 56178
378 322340 0 607 11 7 60436
379 319864 0 2212 14 10 60720
380 315380 0 0 0 0 60720
381 315380 0 0 0 0 60720
382 317291 0 530 60 18 73433
383 280398 0 571 80 14 41477
384 315380 0 0 0 0 60720
385 317330 0 78 16 12 62700
386 238125 0 2489 40 29 67804
387 327071 0 131 6 3 59661
388 309038 0 923 8 6 58620
389 314210 0 72 3 3 60398
390 307930 0 572 16 8 58580
391 322327 0 397 10 10 62710
392 292136 0 450 8 6 59325
393 263276 0 622 7 8 60950
394 367655 0 694 8 6 68060
395 283910 1 3425 12 9 83620
396 283587 0 562 13 8 58456
397 243650 0 4917 42 26 52811
398 438493 1 1442 118 239 121173
399 296261 0 529 9 7 63870
400 230621 1 2126 138 41 21001
401 304252 0 1061 5 3 70415
402 333505 0 776 9 8 64230
403 296919 0 611 8 6 59190
404 278990 1 1526 25 21 69351
405 276898 0 592 7 7 64270
406 327007 0 1182 13 11 70694
407 317046 0 621 16 11 68005
408 304555 0 989 11 12 58930
409 298096 0 438 11 9 58320
410 231861 0 726 3 3 69980
411 309422 0 1303 61 57 69863
412 286963 1 7419 29 21 63255
413 269753 1 1164 17 15 57320
414 448243 1 3310 33 32 75230
415 165404 0 1920 15 11 79420
416 204325 0 965 3 2 73490
417 407159 0 3256 66 23 35250
418 290476 1 1135 17 20 62285
419 275311 0 1270 26 24 69206
420 246541 0 661 3 1 65920
421 253468 0 1013 2 1 69770
422 240897 0 2844 67 74 72683
423 -83265 1 11528 70 68 -14545
424 -42143 0 6526 26 20 55830
425 272713 0 2264 24 20 55174
426 215362 1 5109 97 82 67038
427 42754 0 3999 30 21 51252
428 306275 0 35624 223 244 157278
429 253537 0 9252 48 32 79510
430 372631 0 15236 90 86 77440
431 -7170 0 18073 180 69 27284
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) Group Costs Trades Orders Dividends
155435.184 37009.233 21.738 -1157.434 2499.590 1.669
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-2120094 -11235 34448 58624 2592406
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.554e+05 2.986e+04 5.205 3.03e-07 ***
Group 3.701e+04 2.706e+04 1.368 0.172
Costs 2.174e+01 2.079e+00 10.458 < 2e-16 ***
Trades -1.157e+03 2.075e+02 -5.579 4.31e-08 ***
Orders 2.500e+03 4.380e+02 5.707 2.16e-08 ***
Dividends 1.669e+00 3.987e-01 4.185 3.46e-05 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 244500 on 425 degrees of freedom
Multiple R-squared: 0.6813, Adjusted R-squared: 0.6776
F-statistic: 181.7 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 9.974460e-44 4.987230e-44
[2,] 1.0000000 1.331465e-56 6.657326e-57
[3,] 1.0000000 3.386240e-58 1.693120e-58
[4,] 1.0000000 2.143603e-69 1.071801e-69
[5,] 1.0000000 3.132541e-95 1.566271e-95
[6,] 1.0000000 6.821950e-97 3.410975e-97
[7,] 1.0000000 4.695350e-100 2.347675e-100
[8,] 1.0000000 5.381015e-112 2.690507e-112
[9,] 1.0000000 2.505376e-119 1.252688e-119
[10,] 1.0000000 7.021203e-122 3.510601e-122
[11,] 1.0000000 9.838501e-127 4.919250e-127
[12,] 1.0000000 1.197743e-128 5.988716e-129
[13,] 1.0000000 2.719007e-134 1.359503e-134
[14,] 1.0000000 9.146647e-136 4.573323e-136
[15,] 1.0000000 3.136624e-136 1.568312e-136
[16,] 1.0000000 9.139798e-140 4.569899e-140
[17,] 1.0000000 2.180360e-139 1.090180e-139
[18,] 1.0000000 8.671937e-144 4.335969e-144
[19,] 1.0000000 7.624892e-149 3.812446e-149
[20,] 1.0000000 2.918407e-151 1.459204e-151
[21,] 1.0000000 3.739420e-167 1.869710e-167
[22,] 1.0000000 2.372397e-170 1.186199e-170
[23,] 1.0000000 2.406106e-174 1.203053e-174
[24,] 1.0000000 6.504338e-176 3.252169e-176
[25,] 1.0000000 8.757654e-184 4.378827e-184
[26,] 1.0000000 1.027542e-185 5.137711e-186
[27,] 1.0000000 6.846050e-185 3.423025e-185
[28,] 1.0000000 5.178241e-186 2.589120e-186
[29,] 1.0000000 2.884381e-192 1.442190e-192
[30,] 1.0000000 1.097182e-197 5.485912e-198
[31,] 1.0000000 3.130504e-198 1.565252e-198
[32,] 1.0000000 1.165698e-198 5.828489e-199
[33,] 1.0000000 9.045465e-201 4.522733e-201
[34,] 1.0000000 4.955653e-200 2.477827e-200
[35,] 1.0000000 1.752112e-201 8.760562e-202
[36,] 1.0000000 2.088923e-200 1.044461e-200
[37,] 1.0000000 1.119034e-200 5.595171e-201
[38,] 1.0000000 8.642654e-200 4.321327e-200
[39,] 1.0000000 8.808475e-201 4.404237e-201
[40,] 1.0000000 2.117767e-200 1.058883e-200
[41,] 1.0000000 4.560586e-200 2.280293e-200
[42,] 1.0000000 1.971335e-200 9.856676e-201
[43,] 1.0000000 6.608758e-201 3.304379e-201
[44,] 1.0000000 3.746165e-201 1.873082e-201
[45,] 1.0000000 2.295774e-201 1.147887e-201
[46,] 1.0000000 7.235020e-201 3.617510e-201
[47,] 1.0000000 5.084592e-205 2.542296e-205
[48,] 1.0000000 5.164669e-205 2.582334e-205
[49,] 1.0000000 3.872729e-219 1.936364e-219
[50,] 1.0000000 9.255818e-220 4.627909e-220
[51,] 1.0000000 5.048187e-219 2.524093e-219
[52,] 1.0000000 7.279062e-219 3.639531e-219
[53,] 1.0000000 2.317114e-220 1.158557e-220
[54,] 1.0000000 6.229618e-222 3.114809e-222
[55,] 1.0000000 2.342625e-221 1.171313e-221
[56,] 1.0000000 7.841021e-221 3.920511e-221
[57,] 1.0000000 2.544665e-220 1.272332e-220
[58,] 1.0000000 3.602077e-220 1.801039e-220
[59,] 1.0000000 5.512893e-220 2.756447e-220
[60,] 1.0000000 5.420487e-219 2.710244e-219
[61,] 1.0000000 2.974699e-218 1.487350e-218
[62,] 1.0000000 1.072594e-217 5.362970e-218
[63,] 1.0000000 7.797233e-217 3.898617e-217
[64,] 1.0000000 3.990194e-218 1.995097e-218
[65,] 1.0000000 1.195943e-217 5.979713e-218
[66,] 1.0000000 6.821109e-217 3.410555e-217
[67,] 1.0000000 6.917766e-216 3.458883e-216
[68,] 1.0000000 1.573713e-215 7.868564e-216
[69,] 1.0000000 1.125573e-214 5.627867e-215
[70,] 1.0000000 1.098403e-214 5.492016e-215
[71,] 1.0000000 7.561085e-214 3.780543e-214
[72,] 1.0000000 2.294503e-213 1.147251e-213
[73,] 1.0000000 1.045789e-212 5.228945e-213
[74,] 1.0000000 1.574856e-212 7.874280e-213
[75,] 1.0000000 1.314312e-212 6.571559e-213
[76,] 1.0000000 5.045541e-212 2.522770e-212
[77,] 1.0000000 3.980773e-211 1.990387e-211
[78,] 1.0000000 8.482941e-211 4.241471e-211
[79,] 1.0000000 6.368253e-211 3.184127e-211
[80,] 1.0000000 3.933560e-210 1.966780e-210
[81,] 1.0000000 7.409978e-210 3.704989e-210
[82,] 1.0000000 6.023619e-209 3.011809e-209
[83,] 1.0000000 4.661104e-208 2.330552e-208
[84,] 1.0000000 2.957481e-207 1.478741e-207
[85,] 1.0000000 8.978947e-209 4.489473e-209
[86,] 1.0000000 7.403393e-208 3.701696e-208
[87,] 1.0000000 2.091849e-208 1.045925e-208
[88,] 1.0000000 2.033519e-208 1.016760e-208
[89,] 1.0000000 1.524491e-207 7.622456e-208
[90,] 1.0000000 9.484618e-207 4.742309e-207
[91,] 1.0000000 7.294933e-206 3.647466e-206
[92,] 1.0000000 7.830817e-205 3.915408e-205
[93,] 1.0000000 6.983084e-204 3.491542e-204
[94,] 1.0000000 4.854181e-205 2.427091e-205
[95,] 1.0000000 3.645952e-205 1.822976e-205
[96,] 1.0000000 4.109047e-206 2.054524e-206
[97,] 1.0000000 1.530656e-205 7.653278e-206
[98,] 1.0000000 5.007308e-205 2.503654e-205
[99,] 1.0000000 2.458930e-207 1.229465e-207
[100,] 1.0000000 1.357079e-208 6.785396e-209
[101,] 1.0000000 9.411013e-208 4.705506e-208
[102,] 1.0000000 1.782732e-207 8.913662e-208
[103,] 1.0000000 1.030704e-207 5.153522e-208
[104,] 1.0000000 5.776787e-208 2.888393e-208
[105,] 1.0000000 4.026076e-207 2.013038e-207
[106,] 1.0000000 3.190264e-206 1.595132e-206
[107,] 1.0000000 2.355389e-206 1.177695e-206
[108,] 1.0000000 7.820737e-206 3.910368e-206
[109,] 1.0000000 4.127436e-205 2.063718e-205
[110,] 1.0000000 3.031649e-205 1.515825e-205
[111,] 1.0000000 2.678520e-206 1.339260e-206
[112,] 1.0000000 1.987369e-205 9.936844e-206
[113,] 1.0000000 8.449207e-205 4.224604e-205
[114,] 1.0000000 6.549712e-204 3.274856e-204
[115,] 1.0000000 5.085971e-203 2.542986e-203
[116,] 1.0000000 9.430968e-203 4.715484e-203
[117,] 1.0000000 2.928238e-202 1.464119e-202
[118,] 1.0000000 7.892139e-203 3.946069e-203
[119,] 1.0000000 5.908090e-202 2.954045e-202
[120,] 1.0000000 4.937940e-201 2.468970e-201
[121,] 1.0000000 4.008111e-200 2.004056e-200
[122,] 1.0000000 3.098253e-199 1.549126e-199
[123,] 1.0000000 2.551252e-198 1.275626e-198
[124,] 1.0000000 1.883034e-198 9.415169e-199
[125,] 1.0000000 1.511286e-197 7.556430e-198
[126,] 1.0000000 1.459092e-196 7.295462e-197
[127,] 1.0000000 1.244129e-195 6.220647e-196
[128,] 1.0000000 1.025973e-194 5.129864e-195
[129,] 1.0000000 8.684065e-194 4.342033e-194
[130,] 1.0000000 7.199991e-193 3.599995e-193
[131,] 1.0000000 6.026371e-192 3.013186e-192
[132,] 1.0000000 5.009055e-191 2.504527e-191
[133,] 1.0000000 3.584103e-190 1.792052e-190
[134,] 1.0000000 3.003007e-189 1.501503e-189
[135,] 1.0000000 2.759688e-188 1.379844e-188
[136,] 1.0000000 2.376722e-187 1.188361e-187
[137,] 1.0000000 1.965088e-186 9.825441e-187
[138,] 1.0000000 1.600984e-185 8.004922e-186
[139,] 1.0000000 1.351842e-184 6.759209e-185
[140,] 1.0000000 1.188672e-183 5.943361e-184
[141,] 1.0000000 9.092506e-183 4.546253e-183
[142,] 1.0000000 7.693173e-182 3.846587e-182
[143,] 1.0000000 6.892251e-181 3.446126e-181
[144,] 1.0000000 6.657374e-180 3.328687e-180
[145,] 1.0000000 5.807407e-179 2.903703e-179
[146,] 1.0000000 5.215368e-178 2.607684e-178
[147,] 1.0000000 4.442202e-177 2.221101e-177
[148,] 1.0000000 3.745329e-176 1.872664e-176
[149,] 1.0000000 3.156819e-175 1.578409e-175
[150,] 1.0000000 2.708358e-174 1.354179e-174
[151,] 1.0000000 1.697324e-173 8.486618e-174
[152,] 1.0000000 1.428335e-172 7.141675e-173
[153,] 1.0000000 1.200777e-171 6.003887e-172
[154,] 1.0000000 1.008314e-170 5.041571e-171
[155,] 1.0000000 8.456060e-170 4.228030e-170
[156,] 1.0000000 7.225614e-169 3.612807e-169
[157,] 1.0000000 6.065334e-168 3.032667e-168
[158,] 1.0000000 5.071152e-167 2.535576e-167
[159,] 1.0000000 4.222341e-166 2.111170e-166
[160,] 1.0000000 4.007606e-167 2.003803e-167
[161,] 1.0000000 4.989211e-167 2.494606e-167
[162,] 1.0000000 4.257436e-166 2.128718e-166
[163,] 1.0000000 3.624777e-165 1.812388e-165
[164,] 1.0000000 3.156103e-164 1.578051e-164
[165,] 1.0000000 2.673803e-163 1.336902e-163
[166,] 1.0000000 2.368311e-162 1.184155e-162
[167,] 1.0000000 1.995139e-161 9.975693e-162
[168,] 1.0000000 1.705029e-160 8.525144e-161
[169,] 1.0000000 1.502878e-159 7.514390e-160
[170,] 1.0000000 1.337343e-158 6.686715e-159
[171,] 1.0000000 1.110310e-157 5.551551e-158
[172,] 1.0000000 9.192717e-157 4.596359e-157
[173,] 1.0000000 5.503871e-156 2.751935e-156
[174,] 1.0000000 4.889033e-155 2.444516e-155
[175,] 1.0000000 2.485759e-154 1.242879e-154
[176,] 1.0000000 1.434550e-153 7.172752e-154
[177,] 1.0000000 1.257698e-152 6.288492e-153
[178,] 1.0000000 8.765372e-152 4.382686e-152
[179,] 1.0000000 3.445473e-151 1.722737e-151
[180,] 1.0000000 2.881730e-150 1.440865e-150
[181,] 1.0000000 1.943279e-149 9.716393e-150
[182,] 1.0000000 1.563303e-148 7.816515e-149
[183,] 1.0000000 1.205935e-147 6.029674e-148
[184,] 1.0000000 9.638910e-147 4.819455e-147
[185,] 1.0000000 7.863879e-146 3.931939e-146
[186,] 1.0000000 6.243712e-145 3.121856e-145
[187,] 1.0000000 4.979794e-144 2.489897e-144
[188,] 1.0000000 4.150601e-143 2.075301e-143
[189,] 1.0000000 2.513182e-142 1.256591e-142
[190,] 1.0000000 2.075130e-141 1.037565e-141
[191,] 1.0000000 1.705968e-140 8.529840e-141
[192,] 1.0000000 1.459537e-140 7.297686e-141
[193,] 1.0000000 1.204166e-139 6.020828e-140
[194,] 1.0000000 9.891965e-139 4.945982e-139
[195,] 1.0000000 8.090767e-138 4.045383e-138
[196,] 1.0000000 6.161996e-137 3.080998e-137
[197,] 1.0000000 2.329648e-136 1.164824e-136
[198,] 1.0000000 1.899446e-135 9.497231e-136
[199,] 1.0000000 1.512641e-134 7.563207e-135
[200,] 1.0000000 1.215832e-133 6.079160e-134
[201,] 1.0000000 2.010215e-133 1.005107e-133
[202,] 1.0000000 1.597931e-132 7.989654e-133
[203,] 1.0000000 1.208301e-131 6.041506e-132
[204,] 1.0000000 9.105597e-131 4.552798e-131
[205,] 1.0000000 4.911331e-130 2.455665e-130
[206,] 1.0000000 3.507629e-129 1.753814e-129
[207,] 1.0000000 2.710332e-128 1.355166e-128
[208,] 1.0000000 2.028478e-127 1.014239e-127
[209,] 1.0000000 1.512727e-126 7.563634e-127
[210,] 1.0000000 1.061814e-125 5.309071e-126
[211,] 1.0000000 7.859702e-125 3.929851e-125
[212,] 1.0000000 3.491157e-124 1.745578e-124
[213,] 1.0000000 2.597173e-123 1.298586e-123
[214,] 1.0000000 1.893315e-122 9.466576e-123
[215,] 1.0000000 1.374892e-121 6.874462e-122
[216,] 1.0000000 1.014911e-120 5.074555e-121
[217,] 1.0000000 7.309882e-120 3.654941e-120
[218,] 1.0000000 5.350771e-119 2.675385e-119
[219,] 1.0000000 3.852648e-118 1.926324e-118
[220,] 1.0000000 2.742208e-117 1.371104e-117
[221,] 1.0000000 1.979850e-116 9.899251e-117
[222,] 1.0000000 1.398480e-115 6.992399e-116
[223,] 1.0000000 9.665062e-115 4.832531e-115
[224,] 1.0000000 6.226271e-114 3.113136e-114
[225,] 1.0000000 4.347513e-113 2.173757e-113
[226,] 1.0000000 3.133984e-112 1.566992e-112
[227,] 1.0000000 2.170497e-111 1.085248e-111
[228,] 1.0000000 1.495137e-110 7.475685e-111
[229,] 1.0000000 1.027087e-109 5.135435e-110
[230,] 1.0000000 7.003932e-109 3.501966e-109
[231,] 1.0000000 4.771286e-108 2.385643e-108
[232,] 1.0000000 3.236830e-107 1.618415e-107
[233,] 1.0000000 2.186676e-106 1.093338e-106
[234,] 1.0000000 1.497211e-105 7.486054e-106
[235,] 1.0000000 1.002908e-104 5.014542e-105
[236,] 1.0000000 6.620763e-104 3.310382e-104
[237,] 1.0000000 4.536778e-103 2.268389e-103
[238,] 1.0000000 3.113333e-102 1.556667e-102
[239,] 1.0000000 2.126388e-101 1.063194e-101
[240,] 1.0000000 1.445377e-100 7.226885e-101
[241,] 1.0000000 9.777797e-100 4.888899e-100
[242,] 1.0000000 6.582893e-99 3.291447e-99
[243,] 1.0000000 4.223707e-98 2.111854e-98
[244,] 1.0000000 2.817123e-97 1.408562e-97
[245,] 1.0000000 1.828819e-96 9.144094e-97
[246,] 1.0000000 9.954439e-96 4.977220e-96
[247,] 1.0000000 7.311556e-96 3.655778e-96
[248,] 1.0000000 4.678196e-95 2.339098e-95
[249,] 1.0000000 2.698555e-94 1.349278e-94
[250,] 1.0000000 1.411163e-93 7.055814e-94
[251,] 1.0000000 9.244191e-93 4.622095e-93
[252,] 1.0000000 5.811616e-92 2.905808e-92
[253,] 1.0000000 3.770383e-91 1.885192e-91
[254,] 1.0000000 2.433789e-90 1.216895e-90
[255,] 1.0000000 1.544276e-89 7.721379e-90
[256,] 1.0000000 7.471609e-89 3.735804e-89
[257,] 1.0000000 4.760559e-88 2.380280e-88
[258,] 1.0000000 2.912610e-87 1.456305e-87
[259,] 1.0000000 1.837527e-86 9.187635e-87
[260,] 1.0000000 9.263582e-86 4.631791e-86
[261,] 1.0000000 5.792544e-85 2.896272e-85
[262,] 1.0000000 3.618307e-84 1.809153e-84
[263,] 1.0000000 2.239391e-83 1.119696e-83
[264,] 1.0000000 1.378628e-82 6.893139e-83
[265,] 1.0000000 8.441771e-82 4.220886e-82
[266,] 1.0000000 5.127502e-81 2.563751e-81
[267,] 1.0000000 3.006248e-80 1.503124e-80
[268,] 1.0000000 1.754672e-79 8.773361e-80
[269,] 1.0000000 1.044374e-78 5.221868e-79
[270,] 1.0000000 6.035695e-78 3.017848e-78
[271,] 1.0000000 2.551655e-77 1.275827e-77
[272,] 1.0000000 1.404023e-76 7.020114e-77
[273,] 1.0000000 8.301931e-76 4.150966e-76
[274,] 1.0000000 4.707267e-75 2.353633e-75
[275,] 1.0000000 2.614894e-74 1.307447e-74
[276,] 1.0000000 1.468868e-73 7.344341e-74
[277,] 1.0000000 8.481277e-73 4.240638e-73
[278,] 1.0000000 4.717298e-72 2.358649e-72
[279,] 1.0000000 2.621144e-71 1.310572e-71
[280,] 1.0000000 1.452902e-70 7.264509e-71
[281,] 1.0000000 8.241609e-70 4.120804e-70
[282,] 1.0000000 4.548347e-69 2.274173e-69
[283,] 1.0000000 2.552170e-68 1.276085e-68
[284,] 1.0000000 1.363193e-67 6.815967e-68
[285,] 1.0000000 7.507074e-67 3.753537e-67
[286,] 1.0000000 4.010445e-66 2.005222e-66
[287,] 1.0000000 2.203106e-65 1.101553e-65
[288,] 1.0000000 9.412203e-65 4.706101e-65
[289,] 1.0000000 4.960013e-64 2.480006e-64
[290,] 1.0000000 2.675802e-63 1.337901e-63
[291,] 1.0000000 1.395108e-62 6.975539e-63
[292,] 1.0000000 7.247511e-62 3.623756e-62
[293,] 1.0000000 3.867825e-61 1.933913e-61
[294,] 1.0000000 2.035083e-60 1.017542e-60
[295,] 1.0000000 1.038378e-59 5.191889e-60
[296,] 1.0000000 5.413444e-59 2.706722e-59
[297,] 1.0000000 2.245641e-58 1.122820e-58
[298,] 1.0000000 1.136489e-57 5.682447e-58
[299,] 1.0000000 5.833760e-57 2.916880e-57
[300,] 1.0000000 3.003672e-56 1.501836e-56
[301,] 1.0000000 1.363325e-55 6.816627e-56
[302,] 1.0000000 6.877006e-55 3.438503e-55
[303,] 1.0000000 3.442900e-54 1.721450e-54
[304,] 1.0000000 1.275543e-53 6.377713e-54
[305,] 1.0000000 6.312187e-53 3.156093e-53
[306,] 1.0000000 3.122885e-52 1.561443e-52
[307,] 1.0000000 1.498254e-51 7.491270e-52
[308,] 1.0000000 6.826651e-51 3.413325e-51
[309,] 1.0000000 3.201368e-50 1.600684e-50
[310,] 1.0000000 1.563448e-49 7.817241e-50
[311,] 1.0000000 7.346504e-49 3.673252e-49
[312,] 1.0000000 3.413845e-48 1.706922e-48
[313,] 1.0000000 1.576562e-47 7.882809e-48
[314,] 1.0000000 7.235355e-47 3.617677e-47
[315,] 1.0000000 3.299644e-46 1.649822e-46
[316,] 1.0000000 1.522859e-45 7.614296e-46
[317,] 1.0000000 6.859128e-45 3.429564e-45
[318,] 1.0000000 3.109828e-44 1.554914e-44
[319,] 1.0000000 1.402465e-43 7.012323e-44
[320,] 1.0000000 6.245412e-43 3.122706e-43
[321,] 1.0000000 2.742599e-42 1.371300e-42
[322,] 1.0000000 1.196384e-41 5.981920e-42
[323,] 1.0000000 5.183925e-41 2.591963e-41
[324,] 1.0000000 2.230996e-40 1.115498e-40
[325,] 1.0000000 9.535908e-40 4.767954e-40
[326,] 1.0000000 4.065385e-39 2.032692e-39
[327,] 1.0000000 1.717402e-38 8.587008e-39
[328,] 1.0000000 7.129937e-38 3.564968e-38
[329,] 1.0000000 2.953427e-37 1.476714e-37
[330,] 1.0000000 1.219889e-36 6.099447e-37
[331,] 1.0000000 5.047328e-36 2.523664e-36
[332,] 1.0000000 2.021038e-35 1.010519e-35
[333,] 1.0000000 8.471949e-35 4.235974e-35
[334,] 1.0000000 3.402743e-34 1.701372e-34
[335,] 1.0000000 1.356422e-33 6.782111e-34
[336,] 1.0000000 5.365894e-33 2.682947e-33
[337,] 1.0000000 2.193437e-32 1.096718e-32
[338,] 1.0000000 8.551125e-32 4.275562e-32
[339,] 1.0000000 3.404075e-31 1.702038e-31
[340,] 1.0000000 1.306558e-30 6.532789e-31
[341,] 1.0000000 4.974531e-30 2.487266e-30
[342,] 1.0000000 1.711571e-29 8.557853e-30
[343,] 1.0000000 6.135520e-29 3.067760e-29
[344,] 1.0000000 2.326670e-28 1.163335e-28
[345,] 1.0000000 7.694758e-28 3.847379e-28
[346,] 1.0000000 2.869614e-27 1.434807e-27
[347,] 1.0000000 1.040071e-26 5.200357e-27
[348,] 1.0000000 3.736351e-26 1.868176e-26
[349,] 1.0000000 1.187715e-25 5.938573e-26
[350,] 1.0000000 4.208751e-25 2.104375e-25
[351,] 1.0000000 1.477710e-24 7.388550e-25
[352,] 1.0000000 5.139890e-24 2.569945e-24
[353,] 1.0000000 1.779826e-23 8.899129e-24
[354,] 1.0000000 6.073455e-23 3.036728e-23
[355,] 1.0000000 2.052138e-22 1.026069e-22
[356,] 1.0000000 6.862306e-22 3.431153e-22
[357,] 1.0000000 2.247030e-21 1.123515e-21
[358,] 1.0000000 7.382186e-21 3.691093e-21
[359,] 1.0000000 2.399845e-20 1.199923e-20
[360,] 1.0000000 7.718034e-20 3.859017e-20
[361,] 1.0000000 2.455004e-19 1.227502e-19
[362,] 1.0000000 8.005635e-19 4.002818e-19
[363,] 1.0000000 2.492720e-18 1.246360e-18
[364,] 1.0000000 7.878511e-18 3.939255e-18
[365,] 1.0000000 2.606939e-17 1.303470e-17
[366,] 1.0000000 8.430119e-17 4.215060e-17
[367,] 1.0000000 2.519399e-16 1.259699e-16
[368,] 1.0000000 7.435345e-16 3.717673e-16
[369,] 1.0000000 2.307569e-15 1.153784e-15
[370,] 1.0000000 6.295799e-15 3.147899e-15
[371,] 1.0000000 1.545826e-14 7.729130e-15
[372,] 1.0000000 4.333658e-14 2.166829e-14
[373,] 1.0000000 1.197359e-13 5.986794e-14
[374,] 1.0000000 3.415083e-13 1.707542e-13
[375,] 1.0000000 9.950943e-13 4.975472e-13
[376,] 1.0000000 2.663265e-12 1.331632e-12
[377,] 1.0000000 7.381387e-12 3.690694e-12
[378,] 1.0000000 1.816267e-11 9.081333e-12
[379,] 1.0000000 4.355462e-11 2.177731e-11
[380,] 1.0000000 1.071530e-10 5.357651e-11
[381,] 1.0000000 2.693241e-10 1.346621e-10
[382,] 1.0000000 6.772349e-10 3.386174e-10
[383,] 1.0000000 1.632293e-09 8.161463e-10
[384,] 1.0000000 4.215518e-09 2.107759e-09
[385,] 1.0000000 1.140337e-08 5.701687e-09
[386,] 1.0000000 2.093450e-08 1.046725e-08
[387,] 1.0000000 4.762370e-08 2.381185e-08
[388,] 0.9999999 1.191779e-07 5.958896e-08
[389,] 0.9999999 2.814458e-07 1.407229e-07
[390,] 0.9999997 5.818837e-07 2.909419e-07
[391,] 0.9999993 1.409663e-06 7.048315e-07
[392,] 0.9999982 3.501477e-06 1.750739e-06
[393,] 0.9999959 8.283890e-06 4.141945e-06
[394,] 0.9999920 1.593713e-05 7.968567e-06
[395,] 0.9999834 3.318217e-05 1.659109e-05
[396,] 0.9999642 7.165457e-05 3.582729e-05
[397,] 0.9999197 1.605918e-04 8.029588e-05
[398,] 0.9998380 3.240952e-04 1.620476e-04
[399,] 0.9996707 6.585617e-04 3.292808e-04
[400,] 0.9994272 1.145610e-03 5.728050e-04
[401,] 0.9989885 2.022990e-03 1.011495e-03
[402,] 0.9979016 4.196783e-03 2.098392e-03
[403,] 0.9958308 8.338443e-03 4.169221e-03
[404,] 0.9922712 1.545765e-02 7.728824e-03
[405,] 0.9852129 2.957411e-02 1.478706e-02
[406,] 0.9839799 3.204016e-02 1.602008e-02
[407,] 0.9756373 4.872541e-02 2.436271e-02
[408,] 0.9573504 8.529912e-02 4.264956e-02
[409,] 0.9837381 3.252381e-02 1.626191e-02
[410,] 0.9667559 6.648817e-02 3.324409e-02
[411,] 0.9333787 1.332427e-01 6.662133e-02
[412,] 0.8719773 2.560454e-01 1.280227e-01
[413,] 0.7708062 4.583876e-01 2.291938e-01
[414,] 0.6147456 7.705087e-01 3.852544e-01
> postscript(file="/var/www/html/rcomp/tmp/186hp1291392609.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> plot(x[,1], type='l', main='Actuals and Interpolation', ylab='value of Actuals and Interpolation (dots)', xlab='time or index')
> points(x[,1]-mysum$resid)
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/2jfgs1291392609.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> plot(mysum$resid, type='b', pch=19, main='Residuals', ylab='value of Residuals', xlab='time or index')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/3jfgs1291392609.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> hist(mysum$resid, main='Residual Histogram', xlab='values of Residuals')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/4jfgs1291392609.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> densityplot(~mysum$resid,col='black',main='Residual Density Plot', xlab='values of Residuals')
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/5coxd1291392609.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
1.434540e+06 2.592406e+06 1.260986e+06 -2.120094e+06 -5.389179e+05
6 7 8 9 10
8.706658e+03 3.681755e+04 1.011745e+05 -1.759896e+05 6.212848e+05
11 12 13 14 15
-1.804099e+05 2.960787e+05 1.183419e+05 -1.118617e+06 -7.714344e+04
16 17 18 19 20
5.469366e+05 -2.810947e+05 -6.509910e+04 -3.783602e+05 -5.639902e+05
21 22 23 24 25
2.003133e+05 -1.692827e+05 -1.583710e+05 -7.263975e+05 -3.223750e+05
26 27 28 29 30
-8.611960e+04 -5.244867e+05 -2.287674e+05 3.491694e+05 -4.180886e+04
31 32 33 34 35
2.450783e+05 -1.286764e+05 2.148073e+05 8.843046e+04 -6.741604e+04
36 37 38 39 40
-6.997293e+05 3.290698e+04 9.512216e+04 -1.292850e+05 -3.292756e+05
41 42 43 44 45
2.265576e+05 -3.381594e+05 1.305302e+05 -2.193006e+05 -9.686861e+03
46 47 48 49 50
-9.873843e+04 -1.506070e+05 -1.527146e+05 -4.405757e+04 -2.023685e+05
51 52 53 54 55
-1.079644e+04 -1.554851e+05 -2.840178e+05 9.597283e+04 -4.233277e+05
56 57 58 59 60
1.614878e+04 3.274203e+05 -1.344067e+04 -1.167402e+04 -2.489276e+05
61 62 63 64 65
-3.053437e+04 -1.788236e+05 -2.695521e+04 -2.218082e+05 -4.244611e+02
66 67 68 69 70
2.313680e+05 -3.961293e+00 -5.613649e+04 -1.039046e+05 -1.621399e+05
71 72 73 74 75
-1.575332e+05 8.402731e+03 -1.022603e+05 1.619668e+04 -4.670361e+04
76 77 78 79 80
-1.401200e+05 -1.895204e+04 -3.751264e+04 1.864700e+04 1.016586e+05
81 82 83 84 85
1.779951e+04 -9.639347e+04 -6.771926e+04 4.597112e+04 1.936319e+04
86 87 88 89 90
-1.526918e+05 -7.729575e+04 -2.709022e+03 -1.087286e+05 2.374239e+04
91 92 93 94 95
5.171042e+03 -6.890441e+04 8.655845e+04 -1.670658e+05 4.413430e+04
96 97 98 99 100
-7.960861e+03 6.597813e+04 1.644330e+05 -4.679715e+04 -3.656027e+04
101 102 103 104 105
-7.862713e+04 -1.134372e+05 -1.552571e+05 -4.287334e+04 3.778051e+04
106 107 108 109 110
-6.726786e+04 -2.868120e+05 9.374730e+04 -1.737003e+05 -1.146023e+05
111 112 113 114 115
9.946503e+04 -2.168834e+05 -1.463099e+05 -9.484045e+04 5.483301e+04
116 117 118 119 120
1.896172e+04 3.684722e+04 7.276375e+04 -1.645875e+05 6.720513e+04
121 122 123 124 125
-3.177472e+05 3.423384e+04 -1.454836e+05 4.344236e+04 -1.585187e+03
126 127 128 129 130
-2.042616e+05 5.862432e+04 3.871832e+04 3.169185e+04 5.862432e+04
131 132 133 134 135
4.877076e+04 9.773557e+04 5.862432e+04 8.357447e+03 3.708742e+04
136 137 138 139 140
5.216156e+04 4.932280e+04 5.783791e+04 3.444756e+04 5.862432e+04
141 142 143 144 145
5.938854e+04 5.862432e+04 1.652187e+04 5.147936e+04 3.851483e+04
146 147 148 149 150
4.662419e+04 5.862432e+04 4.729126e+04 5.991813e+04 5.948051e+04
151 152 153 154 155
2.112212e+04 -1.920308e+04 4.612637e+04 2.978399e+04 5.616392e+04
156 157 158 159 160
5.862432e+04 5.862432e+04 5.238087e+04 6.626112e+04 5.862432e+04
161 162 163 164 165
5.862432e+04 5.862432e+04 5.862432e+04 5.258312e+04 3.036140e+04
166 167 168 169 170
5.807737e+04 5.862432e+04 1.115811e+05 3.010201e+05 5.862432e+04
171 172 173 174 175
5.862432e+04 5.059639e+04 5.862432e+04 3.431578e+04 5.862432e+04
176 177 178 179 180
4.724455e+04 -5.870740e+04 2.198344e+04 5.862432e+04 5.862432e+04
181 182 183 184 185
5.508597e+04 -9.384035e+03 1.574433e+04 -2.986051e+03 1.076837e+04
186 187 188 189 190
1.486915e+04 -2.250814e+04 4.192277e+04 1.502288e+04 5.862432e+04
191 192 193 194 195
4.517986e+04 5.862432e+04 4.347433e+04 5.862432e+04 5.285609e+04
196 197 198 199 200
2.161509e+04 -1.976193e+05 2.161509e+04 2.161509e+04 -9.102825e+04
201 202 203 204 205
2.161509e+04 2.161509e+04 2.161509e+04 2.058634e+04 -3.388936e+04
206 207 208 209 210
1.042777e+04 1.478032e+04 2.161509e+04 1.755700e+04 3.580431e+04
211 212 213 214 215
5.862432e+04 5.862432e+04 8.218123e+04 1.798101e+04 2.855512e+04
216 217 218 219 220
5.362514e+04 5.268297e+04 6.738501e+04 5.362514e+04 -3.909274e+03
221 222 223 224 225
3.112129e+04 5.862432e+04 5.838509e+04 1.954647e+04 5.805901e+04
226 227 228 229 230
2.608423e+04 5.017709e+04 5.862432e+04 4.343789e+04 5.862432e+04
231 232 233 234 235
2.264195e+04 3.765605e+04 5.862432e+04 2.404469e+04 5.862432e+04
236 237 238 239 240
5.980690e+04 5.862432e+04 5.045136e+04 5.862432e+04 5.862432e+04
241 242 243 244 245
5.862432e+04 4.899359e+04 5.862432e+04 5.162708e+04 -3.389133e+03
246 247 248 249 250
2.161509e+04 2.161509e+04 2.280813e+04 2.161509e+04 2.161509e+04
251 252 253 254 255
3.608742e+04 2.161509e+04 4.086943e+04 3.105108e+03 -1.241809e+05
256 257 258 259 260
5.862432e+04 4.647750e+04 3.103620e+04 2.161509e+04 5.862432e+04
261 262 263 264 265
2.161509e+04 2.161509e+04 4.343946e+04 6.400187e+04 2.161509e+04
266 267 268 269 270
5.862432e+04 2.161509e+04 8.231452e+04 2.161509e+04 -1.559911e+03
271 272 273 274 275
2.161509e+04 2.161509e+04 2.161509e+04 1.699497e+04 5.862432e+04
276 277 278 279 280
5.896816e+04 3.129980e+04 5.862432e+04 2.101390e+04 5.012522e+04
281 282 283 284 285
2.161509e+04 5.862432e+04 7.633316e+04 5.862432e+04 4.304069e+04
286 287 288 289 290
5.862432e+04 6.827479e+03 5.561020e+04 2.161509e+04 8.052626e+04
291 292 293 294 295
2.161509e+04 5.907775e+04 4.212129e+04 5.862432e+04 3.773515e+04
296 297 298 299 300
4.297017e+04 5.862432e+04 2.161509e+04 5.862432e+04 3.612801e+04
301 302 303 304 305
1.324818e+04 4.114021e+04 5.612473e+04 2.161509e+04 1.042476e+05
306 307 308 309 310
3.899164e+04 2.161509e+04 -1.960650e+04 3.777727e+04 2.161509e+04
311 312 313 314 315
2.161509e+04 5.382871e+04 2.161509e+04 -4.185561e+03 4.339971e+04
316 317 318 319 320
7.721015e+04 5.989160e+04 1.324597e+03 5.304290e+04 5.862432e+04
321 322 323 324 325
5.862432e+04 5.862432e+04 5.862432e+04 2.161509e+04 5.862432e+04
326 327 328 329 330
5.508003e+04 2.161509e+04 2.161509e+04 5.862432e+04 5.862432e+04
331 332 333 334 335
5.862432e+04 5.862432e+04 5.862432e+04 5.124433e+04 5.841418e+04
336 337 338 339 340
4.757854e+04 6.098280e+04 5.862432e+04 5.492330e+04 4.992116e+04
341 342 343 344 345
1.582974e+04 5.862432e+04 5.862432e+04 5.862432e+04 6.106346e+04
346 347 348 349 350
5.862432e+04 4.645356e+04 5.862432e+04 5.862432e+04 2.262242e+04
351 352 353 354 355
6.739839e+04 4.939214e+04 4.695392e+04 5.209069e+04 5.862432e+04
356 357 358 359 360
5.862432e+04 4.700609e+04 5.862432e+04 5.862432e+04 5.862432e+04
361 362 363 364 365
5.324862e+04 5.862432e+04 5.862432e+04 -6.743071e+04 4.776351e+04
366 367 368 369 370
5.862432e+04 5.862432e+04 5.862432e+04 5.862432e+04 6.086752e+04
371 372 373 374 375
5.862432e+04 8.723399e+04 4.129707e+04 3.370318e+04 5.862432e+04
376 377 378 379 380
5.862432e+04 4.384601e+04 4.809817e+04 6.233059e+03 5.862432e+04
381 382 383 384 385
5.862432e+04 5.225427e+04 1.009405e+05 5.862432e+04 4.409872e+04
386 387 388 389 390
-1.107469e+05 6.868064e+04 2.998467e+04 5.240005e+04 4.083358e+04
391 392 393 394 395
4.019934e+04 2.218813e+04 -1.927890e+04 7.782750e+04 -1.311252e+05
396 397 398 399 400
1.344257e+04 -1.231689e+05 -4.483172e+05 1.566969e+04 1.416191e+04
401 402 403 404 405
6.543617e+03 4.444421e+04 2.369665e+04 -8.590407e+04 -8.045107e+03
406 407 408 409 410
1.546557e+04 2.565864e+04 1.202447e+04 2.605950e+04 -6.015431e+04
411 412 413 414 415
-6.278733e+04 -1.912285e+05 -6.145847e+04 1.652315e+04 -1.744255e+05
416 417 418 419 420
-9.624296e+04 1.410265e+05 -6.088788e+04 -5.310840e+04 -3.228747e+04
421 422 423 424 425
-4.059382e+04 -2.050639e+05 -5.909810e+05 -4.524967e+05 -4.621551e+04
426 427 428 429 430
-2.926978e+05 -3.028997e+05 -1.237773e+06 -2.601184e+05 -3.540133e+05
431
-5.651288e+05
> postscript(file="/var/www/html/rcomp/tmp/6coxd1291392609.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 1.434540e+06 NA
1 2.592406e+06 1.434540e+06
2 1.260986e+06 2.592406e+06
3 -2.120094e+06 1.260986e+06
4 -5.389179e+05 -2.120094e+06
5 8.706658e+03 -5.389179e+05
6 3.681755e+04 8.706658e+03
7 1.011745e+05 3.681755e+04
8 -1.759896e+05 1.011745e+05
9 6.212848e+05 -1.759896e+05
10 -1.804099e+05 6.212848e+05
11 2.960787e+05 -1.804099e+05
12 1.183419e+05 2.960787e+05
13 -1.118617e+06 1.183419e+05
14 -7.714344e+04 -1.118617e+06
15 5.469366e+05 -7.714344e+04
16 -2.810947e+05 5.469366e+05
17 -6.509910e+04 -2.810947e+05
18 -3.783602e+05 -6.509910e+04
19 -5.639902e+05 -3.783602e+05
20 2.003133e+05 -5.639902e+05
21 -1.692827e+05 2.003133e+05
22 -1.583710e+05 -1.692827e+05
23 -7.263975e+05 -1.583710e+05
24 -3.223750e+05 -7.263975e+05
25 -8.611960e+04 -3.223750e+05
26 -5.244867e+05 -8.611960e+04
27 -2.287674e+05 -5.244867e+05
28 3.491694e+05 -2.287674e+05
29 -4.180886e+04 3.491694e+05
30 2.450783e+05 -4.180886e+04
31 -1.286764e+05 2.450783e+05
32 2.148073e+05 -1.286764e+05
33 8.843046e+04 2.148073e+05
34 -6.741604e+04 8.843046e+04
35 -6.997293e+05 -6.741604e+04
36 3.290698e+04 -6.997293e+05
37 9.512216e+04 3.290698e+04
38 -1.292850e+05 9.512216e+04
39 -3.292756e+05 -1.292850e+05
40 2.265576e+05 -3.292756e+05
41 -3.381594e+05 2.265576e+05
42 1.305302e+05 -3.381594e+05
43 -2.193006e+05 1.305302e+05
44 -9.686861e+03 -2.193006e+05
45 -9.873843e+04 -9.686861e+03
46 -1.506070e+05 -9.873843e+04
47 -1.527146e+05 -1.506070e+05
48 -4.405757e+04 -1.527146e+05
49 -2.023685e+05 -4.405757e+04
50 -1.079644e+04 -2.023685e+05
51 -1.554851e+05 -1.079644e+04
52 -2.840178e+05 -1.554851e+05
53 9.597283e+04 -2.840178e+05
54 -4.233277e+05 9.597283e+04
55 1.614878e+04 -4.233277e+05
56 3.274203e+05 1.614878e+04
57 -1.344067e+04 3.274203e+05
58 -1.167402e+04 -1.344067e+04
59 -2.489276e+05 -1.167402e+04
60 -3.053437e+04 -2.489276e+05
61 -1.788236e+05 -3.053437e+04
62 -2.695521e+04 -1.788236e+05
63 -2.218082e+05 -2.695521e+04
64 -4.244611e+02 -2.218082e+05
65 2.313680e+05 -4.244611e+02
66 -3.961293e+00 2.313680e+05
67 -5.613649e+04 -3.961293e+00
68 -1.039046e+05 -5.613649e+04
69 -1.621399e+05 -1.039046e+05
70 -1.575332e+05 -1.621399e+05
71 8.402731e+03 -1.575332e+05
72 -1.022603e+05 8.402731e+03
73 1.619668e+04 -1.022603e+05
74 -4.670361e+04 1.619668e+04
75 -1.401200e+05 -4.670361e+04
76 -1.895204e+04 -1.401200e+05
77 -3.751264e+04 -1.895204e+04
78 1.864700e+04 -3.751264e+04
79 1.016586e+05 1.864700e+04
80 1.779951e+04 1.016586e+05
81 -9.639347e+04 1.779951e+04
82 -6.771926e+04 -9.639347e+04
83 4.597112e+04 -6.771926e+04
84 1.936319e+04 4.597112e+04
85 -1.526918e+05 1.936319e+04
86 -7.729575e+04 -1.526918e+05
87 -2.709022e+03 -7.729575e+04
88 -1.087286e+05 -2.709022e+03
89 2.374239e+04 -1.087286e+05
90 5.171042e+03 2.374239e+04
91 -6.890441e+04 5.171042e+03
92 8.655845e+04 -6.890441e+04
93 -1.670658e+05 8.655845e+04
94 4.413430e+04 -1.670658e+05
95 -7.960861e+03 4.413430e+04
96 6.597813e+04 -7.960861e+03
97 1.644330e+05 6.597813e+04
98 -4.679715e+04 1.644330e+05
99 -3.656027e+04 -4.679715e+04
100 -7.862713e+04 -3.656027e+04
101 -1.134372e+05 -7.862713e+04
102 -1.552571e+05 -1.134372e+05
103 -4.287334e+04 -1.552571e+05
104 3.778051e+04 -4.287334e+04
105 -6.726786e+04 3.778051e+04
106 -2.868120e+05 -6.726786e+04
107 9.374730e+04 -2.868120e+05
108 -1.737003e+05 9.374730e+04
109 -1.146023e+05 -1.737003e+05
110 9.946503e+04 -1.146023e+05
111 -2.168834e+05 9.946503e+04
112 -1.463099e+05 -2.168834e+05
113 -9.484045e+04 -1.463099e+05
114 5.483301e+04 -9.484045e+04
115 1.896172e+04 5.483301e+04
116 3.684722e+04 1.896172e+04
117 7.276375e+04 3.684722e+04
118 -1.645875e+05 7.276375e+04
119 6.720513e+04 -1.645875e+05
120 -3.177472e+05 6.720513e+04
121 3.423384e+04 -3.177472e+05
122 -1.454836e+05 3.423384e+04
123 4.344236e+04 -1.454836e+05
124 -1.585187e+03 4.344236e+04
125 -2.042616e+05 -1.585187e+03
126 5.862432e+04 -2.042616e+05
127 3.871832e+04 5.862432e+04
128 3.169185e+04 3.871832e+04
129 5.862432e+04 3.169185e+04
130 4.877076e+04 5.862432e+04
131 9.773557e+04 4.877076e+04
132 5.862432e+04 9.773557e+04
133 8.357447e+03 5.862432e+04
134 3.708742e+04 8.357447e+03
135 5.216156e+04 3.708742e+04
136 4.932280e+04 5.216156e+04
137 5.783791e+04 4.932280e+04
138 3.444756e+04 5.783791e+04
139 5.862432e+04 3.444756e+04
140 5.938854e+04 5.862432e+04
141 5.862432e+04 5.938854e+04
142 1.652187e+04 5.862432e+04
143 5.147936e+04 1.652187e+04
144 3.851483e+04 5.147936e+04
145 4.662419e+04 3.851483e+04
146 5.862432e+04 4.662419e+04
147 4.729126e+04 5.862432e+04
148 5.991813e+04 4.729126e+04
149 5.948051e+04 5.991813e+04
150 2.112212e+04 5.948051e+04
151 -1.920308e+04 2.112212e+04
152 4.612637e+04 -1.920308e+04
153 2.978399e+04 4.612637e+04
154 5.616392e+04 2.978399e+04
155 5.862432e+04 5.616392e+04
156 5.862432e+04 5.862432e+04
157 5.238087e+04 5.862432e+04
158 6.626112e+04 5.238087e+04
159 5.862432e+04 6.626112e+04
160 5.862432e+04 5.862432e+04
161 5.862432e+04 5.862432e+04
162 5.862432e+04 5.862432e+04
163 5.258312e+04 5.862432e+04
164 3.036140e+04 5.258312e+04
165 5.807737e+04 3.036140e+04
166 5.862432e+04 5.807737e+04
167 1.115811e+05 5.862432e+04
168 3.010201e+05 1.115811e+05
169 5.862432e+04 3.010201e+05
170 5.862432e+04 5.862432e+04
171 5.059639e+04 5.862432e+04
172 5.862432e+04 5.059639e+04
173 3.431578e+04 5.862432e+04
174 5.862432e+04 3.431578e+04
175 4.724455e+04 5.862432e+04
176 -5.870740e+04 4.724455e+04
177 2.198344e+04 -5.870740e+04
178 5.862432e+04 2.198344e+04
179 5.862432e+04 5.862432e+04
180 5.508597e+04 5.862432e+04
181 -9.384035e+03 5.508597e+04
182 1.574433e+04 -9.384035e+03
183 -2.986051e+03 1.574433e+04
184 1.076837e+04 -2.986051e+03
185 1.486915e+04 1.076837e+04
186 -2.250814e+04 1.486915e+04
187 4.192277e+04 -2.250814e+04
188 1.502288e+04 4.192277e+04
189 5.862432e+04 1.502288e+04
190 4.517986e+04 5.862432e+04
191 5.862432e+04 4.517986e+04
192 4.347433e+04 5.862432e+04
193 5.862432e+04 4.347433e+04
194 5.285609e+04 5.862432e+04
195 2.161509e+04 5.285609e+04
196 -1.976193e+05 2.161509e+04
197 2.161509e+04 -1.976193e+05
198 2.161509e+04 2.161509e+04
199 -9.102825e+04 2.161509e+04
200 2.161509e+04 -9.102825e+04
201 2.161509e+04 2.161509e+04
202 2.161509e+04 2.161509e+04
203 2.058634e+04 2.161509e+04
204 -3.388936e+04 2.058634e+04
205 1.042777e+04 -3.388936e+04
206 1.478032e+04 1.042777e+04
207 2.161509e+04 1.478032e+04
208 1.755700e+04 2.161509e+04
209 3.580431e+04 1.755700e+04
210 5.862432e+04 3.580431e+04
211 5.862432e+04 5.862432e+04
212 8.218123e+04 5.862432e+04
213 1.798101e+04 8.218123e+04
214 2.855512e+04 1.798101e+04
215 5.362514e+04 2.855512e+04
216 5.268297e+04 5.362514e+04
217 6.738501e+04 5.268297e+04
218 5.362514e+04 6.738501e+04
219 -3.909274e+03 5.362514e+04
220 3.112129e+04 -3.909274e+03
221 5.862432e+04 3.112129e+04
222 5.838509e+04 5.862432e+04
223 1.954647e+04 5.838509e+04
224 5.805901e+04 1.954647e+04
225 2.608423e+04 5.805901e+04
226 5.017709e+04 2.608423e+04
227 5.862432e+04 5.017709e+04
228 4.343789e+04 5.862432e+04
229 5.862432e+04 4.343789e+04
230 2.264195e+04 5.862432e+04
231 3.765605e+04 2.264195e+04
232 5.862432e+04 3.765605e+04
233 2.404469e+04 5.862432e+04
234 5.862432e+04 2.404469e+04
235 5.980690e+04 5.862432e+04
236 5.862432e+04 5.980690e+04
237 5.045136e+04 5.862432e+04
238 5.862432e+04 5.045136e+04
239 5.862432e+04 5.862432e+04
240 5.862432e+04 5.862432e+04
241 4.899359e+04 5.862432e+04
242 5.862432e+04 4.899359e+04
243 5.162708e+04 5.862432e+04
244 -3.389133e+03 5.162708e+04
245 2.161509e+04 -3.389133e+03
246 2.161509e+04 2.161509e+04
247 2.280813e+04 2.161509e+04
248 2.161509e+04 2.280813e+04
249 2.161509e+04 2.161509e+04
250 3.608742e+04 2.161509e+04
251 2.161509e+04 3.608742e+04
252 4.086943e+04 2.161509e+04
253 3.105108e+03 4.086943e+04
254 -1.241809e+05 3.105108e+03
255 5.862432e+04 -1.241809e+05
256 4.647750e+04 5.862432e+04
257 3.103620e+04 4.647750e+04
258 2.161509e+04 3.103620e+04
259 5.862432e+04 2.161509e+04
260 2.161509e+04 5.862432e+04
261 2.161509e+04 2.161509e+04
262 4.343946e+04 2.161509e+04
263 6.400187e+04 4.343946e+04
264 2.161509e+04 6.400187e+04
265 5.862432e+04 2.161509e+04
266 2.161509e+04 5.862432e+04
267 8.231452e+04 2.161509e+04
268 2.161509e+04 8.231452e+04
269 -1.559911e+03 2.161509e+04
270 2.161509e+04 -1.559911e+03
271 2.161509e+04 2.161509e+04
272 2.161509e+04 2.161509e+04
273 1.699497e+04 2.161509e+04
274 5.862432e+04 1.699497e+04
275 5.896816e+04 5.862432e+04
276 3.129980e+04 5.896816e+04
277 5.862432e+04 3.129980e+04
278 2.101390e+04 5.862432e+04
279 5.012522e+04 2.101390e+04
280 2.161509e+04 5.012522e+04
281 5.862432e+04 2.161509e+04
282 7.633316e+04 5.862432e+04
283 5.862432e+04 7.633316e+04
284 4.304069e+04 5.862432e+04
285 5.862432e+04 4.304069e+04
286 6.827479e+03 5.862432e+04
287 5.561020e+04 6.827479e+03
288 2.161509e+04 5.561020e+04
289 8.052626e+04 2.161509e+04
290 2.161509e+04 8.052626e+04
291 5.907775e+04 2.161509e+04
292 4.212129e+04 5.907775e+04
293 5.862432e+04 4.212129e+04
294 3.773515e+04 5.862432e+04
295 4.297017e+04 3.773515e+04
296 5.862432e+04 4.297017e+04
297 2.161509e+04 5.862432e+04
298 5.862432e+04 2.161509e+04
299 3.612801e+04 5.862432e+04
300 1.324818e+04 3.612801e+04
301 4.114021e+04 1.324818e+04
302 5.612473e+04 4.114021e+04
303 2.161509e+04 5.612473e+04
304 1.042476e+05 2.161509e+04
305 3.899164e+04 1.042476e+05
306 2.161509e+04 3.899164e+04
307 -1.960650e+04 2.161509e+04
308 3.777727e+04 -1.960650e+04
309 2.161509e+04 3.777727e+04
310 2.161509e+04 2.161509e+04
311 5.382871e+04 2.161509e+04
312 2.161509e+04 5.382871e+04
313 -4.185561e+03 2.161509e+04
314 4.339971e+04 -4.185561e+03
315 7.721015e+04 4.339971e+04
316 5.989160e+04 7.721015e+04
317 1.324597e+03 5.989160e+04
318 5.304290e+04 1.324597e+03
319 5.862432e+04 5.304290e+04
320 5.862432e+04 5.862432e+04
321 5.862432e+04 5.862432e+04
322 5.862432e+04 5.862432e+04
323 2.161509e+04 5.862432e+04
324 5.862432e+04 2.161509e+04
325 5.508003e+04 5.862432e+04
326 2.161509e+04 5.508003e+04
327 2.161509e+04 2.161509e+04
328 5.862432e+04 2.161509e+04
329 5.862432e+04 5.862432e+04
330 5.862432e+04 5.862432e+04
331 5.862432e+04 5.862432e+04
332 5.862432e+04 5.862432e+04
333 5.124433e+04 5.862432e+04
334 5.841418e+04 5.124433e+04
335 4.757854e+04 5.841418e+04
336 6.098280e+04 4.757854e+04
337 5.862432e+04 6.098280e+04
338 5.492330e+04 5.862432e+04
339 4.992116e+04 5.492330e+04
340 1.582974e+04 4.992116e+04
341 5.862432e+04 1.582974e+04
342 5.862432e+04 5.862432e+04
343 5.862432e+04 5.862432e+04
344 6.106346e+04 5.862432e+04
345 5.862432e+04 6.106346e+04
346 4.645356e+04 5.862432e+04
347 5.862432e+04 4.645356e+04
348 5.862432e+04 5.862432e+04
349 2.262242e+04 5.862432e+04
350 6.739839e+04 2.262242e+04
351 4.939214e+04 6.739839e+04
352 4.695392e+04 4.939214e+04
353 5.209069e+04 4.695392e+04
354 5.862432e+04 5.209069e+04
355 5.862432e+04 5.862432e+04
356 4.700609e+04 5.862432e+04
357 5.862432e+04 4.700609e+04
358 5.862432e+04 5.862432e+04
359 5.862432e+04 5.862432e+04
360 5.324862e+04 5.862432e+04
361 5.862432e+04 5.324862e+04
362 5.862432e+04 5.862432e+04
363 -6.743071e+04 5.862432e+04
364 4.776351e+04 -6.743071e+04
365 5.862432e+04 4.776351e+04
366 5.862432e+04 5.862432e+04
367 5.862432e+04 5.862432e+04
368 5.862432e+04 5.862432e+04
369 6.086752e+04 5.862432e+04
370 5.862432e+04 6.086752e+04
371 8.723399e+04 5.862432e+04
372 4.129707e+04 8.723399e+04
373 3.370318e+04 4.129707e+04
374 5.862432e+04 3.370318e+04
375 5.862432e+04 5.862432e+04
376 4.384601e+04 5.862432e+04
377 4.809817e+04 4.384601e+04
378 6.233059e+03 4.809817e+04
379 5.862432e+04 6.233059e+03
380 5.862432e+04 5.862432e+04
381 5.225427e+04 5.862432e+04
382 1.009405e+05 5.225427e+04
383 5.862432e+04 1.009405e+05
384 4.409872e+04 5.862432e+04
385 -1.107469e+05 4.409872e+04
386 6.868064e+04 -1.107469e+05
387 2.998467e+04 6.868064e+04
388 5.240005e+04 2.998467e+04
389 4.083358e+04 5.240005e+04
390 4.019934e+04 4.083358e+04
391 2.218813e+04 4.019934e+04
392 -1.927890e+04 2.218813e+04
393 7.782750e+04 -1.927890e+04
394 -1.311252e+05 7.782750e+04
395 1.344257e+04 -1.311252e+05
396 -1.231689e+05 1.344257e+04
397 -4.483172e+05 -1.231689e+05
398 1.566969e+04 -4.483172e+05
399 1.416191e+04 1.566969e+04
400 6.543617e+03 1.416191e+04
401 4.444421e+04 6.543617e+03
402 2.369665e+04 4.444421e+04
403 -8.590407e+04 2.369665e+04
404 -8.045107e+03 -8.590407e+04
405 1.546557e+04 -8.045107e+03
406 2.565864e+04 1.546557e+04
407 1.202447e+04 2.565864e+04
408 2.605950e+04 1.202447e+04
409 -6.015431e+04 2.605950e+04
410 -6.278733e+04 -6.015431e+04
411 -1.912285e+05 -6.278733e+04
412 -6.145847e+04 -1.912285e+05
413 1.652315e+04 -6.145847e+04
414 -1.744255e+05 1.652315e+04
415 -9.624296e+04 -1.744255e+05
416 1.410265e+05 -9.624296e+04
417 -6.088788e+04 1.410265e+05
418 -5.310840e+04 -6.088788e+04
419 -3.228747e+04 -5.310840e+04
420 -4.059382e+04 -3.228747e+04
421 -2.050639e+05 -4.059382e+04
422 -5.909810e+05 -2.050639e+05
423 -4.524967e+05 -5.909810e+05
424 -4.621551e+04 -4.524967e+05
425 -2.926978e+05 -4.621551e+04
426 -3.028997e+05 -2.926978e+05
427 -1.237773e+06 -3.028997e+05
428 -2.601184e+05 -1.237773e+06
429 -3.540133e+05 -2.601184e+05
430 -5.651288e+05 -3.540133e+05
431 NA -5.651288e+05
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] 2.592406e+06 1.434540e+06
[2,] 1.260986e+06 2.592406e+06
[3,] -2.120094e+06 1.260986e+06
[4,] -5.389179e+05 -2.120094e+06
[5,] 8.706658e+03 -5.389179e+05
[6,] 3.681755e+04 8.706658e+03
[7,] 1.011745e+05 3.681755e+04
[8,] -1.759896e+05 1.011745e+05
[9,] 6.212848e+05 -1.759896e+05
[10,] -1.804099e+05 6.212848e+05
[11,] 2.960787e+05 -1.804099e+05
[12,] 1.183419e+05 2.960787e+05
[13,] -1.118617e+06 1.183419e+05
[14,] -7.714344e+04 -1.118617e+06
[15,] 5.469366e+05 -7.714344e+04
[16,] -2.810947e+05 5.469366e+05
[17,] -6.509910e+04 -2.810947e+05
[18,] -3.783602e+05 -6.509910e+04
[19,] -5.639902e+05 -3.783602e+05
[20,] 2.003133e+05 -5.639902e+05
[21,] -1.692827e+05 2.003133e+05
[22,] -1.583710e+05 -1.692827e+05
[23,] -7.263975e+05 -1.583710e+05
[24,] -3.223750e+05 -7.263975e+05
[25,] -8.611960e+04 -3.223750e+05
[26,] -5.244867e+05 -8.611960e+04
[27,] -2.287674e+05 -5.244867e+05
[28,] 3.491694e+05 -2.287674e+05
[29,] -4.180886e+04 3.491694e+05
[30,] 2.450783e+05 -4.180886e+04
[31,] -1.286764e+05 2.450783e+05
[32,] 2.148073e+05 -1.286764e+05
[33,] 8.843046e+04 2.148073e+05
[34,] -6.741604e+04 8.843046e+04
[35,] -6.997293e+05 -6.741604e+04
[36,] 3.290698e+04 -6.997293e+05
[37,] 9.512216e+04 3.290698e+04
[38,] -1.292850e+05 9.512216e+04
[39,] -3.292756e+05 -1.292850e+05
[40,] 2.265576e+05 -3.292756e+05
[41,] -3.381594e+05 2.265576e+05
[42,] 1.305302e+05 -3.381594e+05
[43,] -2.193006e+05 1.305302e+05
[44,] -9.686861e+03 -2.193006e+05
[45,] -9.873843e+04 -9.686861e+03
[46,] -1.506070e+05 -9.873843e+04
[47,] -1.527146e+05 -1.506070e+05
[48,] -4.405757e+04 -1.527146e+05
[49,] -2.023685e+05 -4.405757e+04
[50,] -1.079644e+04 -2.023685e+05
[51,] -1.554851e+05 -1.079644e+04
[52,] -2.840178e+05 -1.554851e+05
[53,] 9.597283e+04 -2.840178e+05
[54,] -4.233277e+05 9.597283e+04
[55,] 1.614878e+04 -4.233277e+05
[56,] 3.274203e+05 1.614878e+04
[57,] -1.344067e+04 3.274203e+05
[58,] -1.167402e+04 -1.344067e+04
[59,] -2.489276e+05 -1.167402e+04
[60,] -3.053437e+04 -2.489276e+05
[61,] -1.788236e+05 -3.053437e+04
[62,] -2.695521e+04 -1.788236e+05
[63,] -2.218082e+05 -2.695521e+04
[64,] -4.244611e+02 -2.218082e+05
[65,] 2.313680e+05 -4.244611e+02
[66,] -3.961293e+00 2.313680e+05
[67,] -5.613649e+04 -3.961293e+00
[68,] -1.039046e+05 -5.613649e+04
[69,] -1.621399e+05 -1.039046e+05
[70,] -1.575332e+05 -1.621399e+05
[71,] 8.402731e+03 -1.575332e+05
[72,] -1.022603e+05 8.402731e+03
[73,] 1.619668e+04 -1.022603e+05
[74,] -4.670361e+04 1.619668e+04
[75,] -1.401200e+05 -4.670361e+04
[76,] -1.895204e+04 -1.401200e+05
[77,] -3.751264e+04 -1.895204e+04
[78,] 1.864700e+04 -3.751264e+04
[79,] 1.016586e+05 1.864700e+04
[80,] 1.779951e+04 1.016586e+05
[81,] -9.639347e+04 1.779951e+04
[82,] -6.771926e+04 -9.639347e+04
[83,] 4.597112e+04 -6.771926e+04
[84,] 1.936319e+04 4.597112e+04
[85,] -1.526918e+05 1.936319e+04
[86,] -7.729575e+04 -1.526918e+05
[87,] -2.709022e+03 -7.729575e+04
[88,] -1.087286e+05 -2.709022e+03
[89,] 2.374239e+04 -1.087286e+05
[90,] 5.171042e+03 2.374239e+04
[91,] -6.890441e+04 5.171042e+03
[92,] 8.655845e+04 -6.890441e+04
[93,] -1.670658e+05 8.655845e+04
[94,] 4.413430e+04 -1.670658e+05
[95,] -7.960861e+03 4.413430e+04
[96,] 6.597813e+04 -7.960861e+03
[97,] 1.644330e+05 6.597813e+04
[98,] -4.679715e+04 1.644330e+05
[99,] -3.656027e+04 -4.679715e+04
[100,] -7.862713e+04 -3.656027e+04
[101,] -1.134372e+05 -7.862713e+04
[102,] -1.552571e+05 -1.134372e+05
[103,] -4.287334e+04 -1.552571e+05
[104,] 3.778051e+04 -4.287334e+04
[105,] -6.726786e+04 3.778051e+04
[106,] -2.868120e+05 -6.726786e+04
[107,] 9.374730e+04 -2.868120e+05
[108,] -1.737003e+05 9.374730e+04
[109,] -1.146023e+05 -1.737003e+05
[110,] 9.946503e+04 -1.146023e+05
[111,] -2.168834e+05 9.946503e+04
[112,] -1.463099e+05 -2.168834e+05
[113,] -9.484045e+04 -1.463099e+05
[114,] 5.483301e+04 -9.484045e+04
[115,] 1.896172e+04 5.483301e+04
[116,] 3.684722e+04 1.896172e+04
[117,] 7.276375e+04 3.684722e+04
[118,] -1.645875e+05 7.276375e+04
[119,] 6.720513e+04 -1.645875e+05
[120,] -3.177472e+05 6.720513e+04
[121,] 3.423384e+04 -3.177472e+05
[122,] -1.454836e+05 3.423384e+04
[123,] 4.344236e+04 -1.454836e+05
[124,] -1.585187e+03 4.344236e+04
[125,] -2.042616e+05 -1.585187e+03
[126,] 5.862432e+04 -2.042616e+05
[127,] 3.871832e+04 5.862432e+04
[128,] 3.169185e+04 3.871832e+04
[129,] 5.862432e+04 3.169185e+04
[130,] 4.877076e+04 5.862432e+04
[131,] 9.773557e+04 4.877076e+04
[132,] 5.862432e+04 9.773557e+04
[133,] 8.357447e+03 5.862432e+04
[134,] 3.708742e+04 8.357447e+03
[135,] 5.216156e+04 3.708742e+04
[136,] 4.932280e+04 5.216156e+04
[137,] 5.783791e+04 4.932280e+04
[138,] 3.444756e+04 5.783791e+04
[139,] 5.862432e+04 3.444756e+04
[140,] 5.938854e+04 5.862432e+04
[141,] 5.862432e+04 5.938854e+04
[142,] 1.652187e+04 5.862432e+04
[143,] 5.147936e+04 1.652187e+04
[144,] 3.851483e+04 5.147936e+04
[145,] 4.662419e+04 3.851483e+04
[146,] 5.862432e+04 4.662419e+04
[147,] 4.729126e+04 5.862432e+04
[148,] 5.991813e+04 4.729126e+04
[149,] 5.948051e+04 5.991813e+04
[150,] 2.112212e+04 5.948051e+04
[151,] -1.920308e+04 2.112212e+04
[152,] 4.612637e+04 -1.920308e+04
[153,] 2.978399e+04 4.612637e+04
[154,] 5.616392e+04 2.978399e+04
[155,] 5.862432e+04 5.616392e+04
[156,] 5.862432e+04 5.862432e+04
[157,] 5.238087e+04 5.862432e+04
[158,] 6.626112e+04 5.238087e+04
[159,] 5.862432e+04 6.626112e+04
[160,] 5.862432e+04 5.862432e+04
[161,] 5.862432e+04 5.862432e+04
[162,] 5.862432e+04 5.862432e+04
[163,] 5.258312e+04 5.862432e+04
[164,] 3.036140e+04 5.258312e+04
[165,] 5.807737e+04 3.036140e+04
[166,] 5.862432e+04 5.807737e+04
[167,] 1.115811e+05 5.862432e+04
[168,] 3.010201e+05 1.115811e+05
[169,] 5.862432e+04 3.010201e+05
[170,] 5.862432e+04 5.862432e+04
[171,] 5.059639e+04 5.862432e+04
[172,] 5.862432e+04 5.059639e+04
[173,] 3.431578e+04 5.862432e+04
[174,] 5.862432e+04 3.431578e+04
[175,] 4.724455e+04 5.862432e+04
[176,] -5.870740e+04 4.724455e+04
[177,] 2.198344e+04 -5.870740e+04
[178,] 5.862432e+04 2.198344e+04
[179,] 5.862432e+04 5.862432e+04
[180,] 5.508597e+04 5.862432e+04
[181,] -9.384035e+03 5.508597e+04
[182,] 1.574433e+04 -9.384035e+03
[183,] -2.986051e+03 1.574433e+04
[184,] 1.076837e+04 -2.986051e+03
[185,] 1.486915e+04 1.076837e+04
[186,] -2.250814e+04 1.486915e+04
[187,] 4.192277e+04 -2.250814e+04
[188,] 1.502288e+04 4.192277e+04
[189,] 5.862432e+04 1.502288e+04
[190,] 4.517986e+04 5.862432e+04
[191,] 5.862432e+04 4.517986e+04
[192,] 4.347433e+04 5.862432e+04
[193,] 5.862432e+04 4.347433e+04
[194,] 5.285609e+04 5.862432e+04
[195,] 2.161509e+04 5.285609e+04
[196,] -1.976193e+05 2.161509e+04
[197,] 2.161509e+04 -1.976193e+05
[198,] 2.161509e+04 2.161509e+04
[199,] -9.102825e+04 2.161509e+04
[200,] 2.161509e+04 -9.102825e+04
[201,] 2.161509e+04 2.161509e+04
[202,] 2.161509e+04 2.161509e+04
[203,] 2.058634e+04 2.161509e+04
[204,] -3.388936e+04 2.058634e+04
[205,] 1.042777e+04 -3.388936e+04
[206,] 1.478032e+04 1.042777e+04
[207,] 2.161509e+04 1.478032e+04
[208,] 1.755700e+04 2.161509e+04
[209,] 3.580431e+04 1.755700e+04
[210,] 5.862432e+04 3.580431e+04
[211,] 5.862432e+04 5.862432e+04
[212,] 8.218123e+04 5.862432e+04
[213,] 1.798101e+04 8.218123e+04
[214,] 2.855512e+04 1.798101e+04
[215,] 5.362514e+04 2.855512e+04
[216,] 5.268297e+04 5.362514e+04
[217,] 6.738501e+04 5.268297e+04
[218,] 5.362514e+04 6.738501e+04
[219,] -3.909274e+03 5.362514e+04
[220,] 3.112129e+04 -3.909274e+03
[221,] 5.862432e+04 3.112129e+04
[222,] 5.838509e+04 5.862432e+04
[223,] 1.954647e+04 5.838509e+04
[224,] 5.805901e+04 1.954647e+04
[225,] 2.608423e+04 5.805901e+04
[226,] 5.017709e+04 2.608423e+04
[227,] 5.862432e+04 5.017709e+04
[228,] 4.343789e+04 5.862432e+04
[229,] 5.862432e+04 4.343789e+04
[230,] 2.264195e+04 5.862432e+04
[231,] 3.765605e+04 2.264195e+04
[232,] 5.862432e+04 3.765605e+04
[233,] 2.404469e+04 5.862432e+04
[234,] 5.862432e+04 2.404469e+04
[235,] 5.980690e+04 5.862432e+04
[236,] 5.862432e+04 5.980690e+04
[237,] 5.045136e+04 5.862432e+04
[238,] 5.862432e+04 5.045136e+04
[239,] 5.862432e+04 5.862432e+04
[240,] 5.862432e+04 5.862432e+04
[241,] 4.899359e+04 5.862432e+04
[242,] 5.862432e+04 4.899359e+04
[243,] 5.162708e+04 5.862432e+04
[244,] -3.389133e+03 5.162708e+04
[245,] 2.161509e+04 -3.389133e+03
[246,] 2.161509e+04 2.161509e+04
[247,] 2.280813e+04 2.161509e+04
[248,] 2.161509e+04 2.280813e+04
[249,] 2.161509e+04 2.161509e+04
[250,] 3.608742e+04 2.161509e+04
[251,] 2.161509e+04 3.608742e+04
[252,] 4.086943e+04 2.161509e+04
[253,] 3.105108e+03 4.086943e+04
[254,] -1.241809e+05 3.105108e+03
[255,] 5.862432e+04 -1.241809e+05
[256,] 4.647750e+04 5.862432e+04
[257,] 3.103620e+04 4.647750e+04
[258,] 2.161509e+04 3.103620e+04
[259,] 5.862432e+04 2.161509e+04
[260,] 2.161509e+04 5.862432e+04
[261,] 2.161509e+04 2.161509e+04
[262,] 4.343946e+04 2.161509e+04
[263,] 6.400187e+04 4.343946e+04
[264,] 2.161509e+04 6.400187e+04
[265,] 5.862432e+04 2.161509e+04
[266,] 2.161509e+04 5.862432e+04
[267,] 8.231452e+04 2.161509e+04
[268,] 2.161509e+04 8.231452e+04
[269,] -1.559911e+03 2.161509e+04
[270,] 2.161509e+04 -1.559911e+03
[271,] 2.161509e+04 2.161509e+04
[272,] 2.161509e+04 2.161509e+04
[273,] 1.699497e+04 2.161509e+04
[274,] 5.862432e+04 1.699497e+04
[275,] 5.896816e+04 5.862432e+04
[276,] 3.129980e+04 5.896816e+04
[277,] 5.862432e+04 3.129980e+04
[278,] 2.101390e+04 5.862432e+04
[279,] 5.012522e+04 2.101390e+04
[280,] 2.161509e+04 5.012522e+04
[281,] 5.862432e+04 2.161509e+04
[282,] 7.633316e+04 5.862432e+04
[283,] 5.862432e+04 7.633316e+04
[284,] 4.304069e+04 5.862432e+04
[285,] 5.862432e+04 4.304069e+04
[286,] 6.827479e+03 5.862432e+04
[287,] 5.561020e+04 6.827479e+03
[288,] 2.161509e+04 5.561020e+04
[289,] 8.052626e+04 2.161509e+04
[290,] 2.161509e+04 8.052626e+04
[291,] 5.907775e+04 2.161509e+04
[292,] 4.212129e+04 5.907775e+04
[293,] 5.862432e+04 4.212129e+04
[294,] 3.773515e+04 5.862432e+04
[295,] 4.297017e+04 3.773515e+04
[296,] 5.862432e+04 4.297017e+04
[297,] 2.161509e+04 5.862432e+04
[298,] 5.862432e+04 2.161509e+04
[299,] 3.612801e+04 5.862432e+04
[300,] 1.324818e+04 3.612801e+04
[301,] 4.114021e+04 1.324818e+04
[302,] 5.612473e+04 4.114021e+04
[303,] 2.161509e+04 5.612473e+04
[304,] 1.042476e+05 2.161509e+04
[305,] 3.899164e+04 1.042476e+05
[306,] 2.161509e+04 3.899164e+04
[307,] -1.960650e+04 2.161509e+04
[308,] 3.777727e+04 -1.960650e+04
[309,] 2.161509e+04 3.777727e+04
[310,] 2.161509e+04 2.161509e+04
[311,] 5.382871e+04 2.161509e+04
[312,] 2.161509e+04 5.382871e+04
[313,] -4.185561e+03 2.161509e+04
[314,] 4.339971e+04 -4.185561e+03
[315,] 7.721015e+04 4.339971e+04
[316,] 5.989160e+04 7.721015e+04
[317,] 1.324597e+03 5.989160e+04
[318,] 5.304290e+04 1.324597e+03
[319,] 5.862432e+04 5.304290e+04
[320,] 5.862432e+04 5.862432e+04
[321,] 5.862432e+04 5.862432e+04
[322,] 5.862432e+04 5.862432e+04
[323,] 2.161509e+04 5.862432e+04
[324,] 5.862432e+04 2.161509e+04
[325,] 5.508003e+04 5.862432e+04
[326,] 2.161509e+04 5.508003e+04
[327,] 2.161509e+04 2.161509e+04
[328,] 5.862432e+04 2.161509e+04
[329,] 5.862432e+04 5.862432e+04
[330,] 5.862432e+04 5.862432e+04
[331,] 5.862432e+04 5.862432e+04
[332,] 5.862432e+04 5.862432e+04
[333,] 5.124433e+04 5.862432e+04
[334,] 5.841418e+04 5.124433e+04
[335,] 4.757854e+04 5.841418e+04
[336,] 6.098280e+04 4.757854e+04
[337,] 5.862432e+04 6.098280e+04
[338,] 5.492330e+04 5.862432e+04
[339,] 4.992116e+04 5.492330e+04
[340,] 1.582974e+04 4.992116e+04
[341,] 5.862432e+04 1.582974e+04
[342,] 5.862432e+04 5.862432e+04
[343,] 5.862432e+04 5.862432e+04
[344,] 6.106346e+04 5.862432e+04
[345,] 5.862432e+04 6.106346e+04
[346,] 4.645356e+04 5.862432e+04
[347,] 5.862432e+04 4.645356e+04
[348,] 5.862432e+04 5.862432e+04
[349,] 2.262242e+04 5.862432e+04
[350,] 6.739839e+04 2.262242e+04
[351,] 4.939214e+04 6.739839e+04
[352,] 4.695392e+04 4.939214e+04
[353,] 5.209069e+04 4.695392e+04
[354,] 5.862432e+04 5.209069e+04
[355,] 5.862432e+04 5.862432e+04
[356,] 4.700609e+04 5.862432e+04
[357,] 5.862432e+04 4.700609e+04
[358,] 5.862432e+04 5.862432e+04
[359,] 5.862432e+04 5.862432e+04
[360,] 5.324862e+04 5.862432e+04
[361,] 5.862432e+04 5.324862e+04
[362,] 5.862432e+04 5.862432e+04
[363,] -6.743071e+04 5.862432e+04
[364,] 4.776351e+04 -6.743071e+04
[365,] 5.862432e+04 4.776351e+04
[366,] 5.862432e+04 5.862432e+04
[367,] 5.862432e+04 5.862432e+04
[368,] 5.862432e+04 5.862432e+04
[369,] 6.086752e+04 5.862432e+04
[370,] 5.862432e+04 6.086752e+04
[371,] 8.723399e+04 5.862432e+04
[372,] 4.129707e+04 8.723399e+04
[373,] 3.370318e+04 4.129707e+04
[374,] 5.862432e+04 3.370318e+04
[375,] 5.862432e+04 5.862432e+04
[376,] 4.384601e+04 5.862432e+04
[377,] 4.809817e+04 4.384601e+04
[378,] 6.233059e+03 4.809817e+04
[379,] 5.862432e+04 6.233059e+03
[380,] 5.862432e+04 5.862432e+04
[381,] 5.225427e+04 5.862432e+04
[382,] 1.009405e+05 5.225427e+04
[383,] 5.862432e+04 1.009405e+05
[384,] 4.409872e+04 5.862432e+04
[385,] -1.107469e+05 4.409872e+04
[386,] 6.868064e+04 -1.107469e+05
[387,] 2.998467e+04 6.868064e+04
[388,] 5.240005e+04 2.998467e+04
[389,] 4.083358e+04 5.240005e+04
[390,] 4.019934e+04 4.083358e+04
[391,] 2.218813e+04 4.019934e+04
[392,] -1.927890e+04 2.218813e+04
[393,] 7.782750e+04 -1.927890e+04
[394,] -1.311252e+05 7.782750e+04
[395,] 1.344257e+04 -1.311252e+05
[396,] -1.231689e+05 1.344257e+04
[397,] -4.483172e+05 -1.231689e+05
[398,] 1.566969e+04 -4.483172e+05
[399,] 1.416191e+04 1.566969e+04
[400,] 6.543617e+03 1.416191e+04
[401,] 4.444421e+04 6.543617e+03
[402,] 2.369665e+04 4.444421e+04
[403,] -8.590407e+04 2.369665e+04
[404,] -8.045107e+03 -8.590407e+04
[405,] 1.546557e+04 -8.045107e+03
[406,] 2.565864e+04 1.546557e+04
[407,] 1.202447e+04 2.565864e+04
[408,] 2.605950e+04 1.202447e+04
[409,] -6.015431e+04 2.605950e+04
[410,] -6.278733e+04 -6.015431e+04
[411,] -1.912285e+05 -6.278733e+04
[412,] -6.145847e+04 -1.912285e+05
[413,] 1.652315e+04 -6.145847e+04
[414,] -1.744255e+05 1.652315e+04
[415,] -9.624296e+04 -1.744255e+05
[416,] 1.410265e+05 -9.624296e+04
[417,] -6.088788e+04 1.410265e+05
[418,] -5.310840e+04 -6.088788e+04
[419,] -3.228747e+04 -5.310840e+04
[420,] -4.059382e+04 -3.228747e+04
[421,] -2.050639e+05 -4.059382e+04
[422,] -5.909810e+05 -2.050639e+05
[423,] -4.524967e+05 -5.909810e+05
[424,] -4.621551e+04 -4.524967e+05
[425,] -2.926978e+05 -4.621551e+04
[426,] -3.028997e+05 -2.926978e+05
[427,] -1.237773e+06 -3.028997e+05
[428,] -2.601184e+05 -1.237773e+06
[429,] -3.540133e+05 -2.601184e+05
[430,] -5.651288e+05 -3.540133e+05
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 2.592406e+06 1.434540e+06
2 1.260986e+06 2.592406e+06
3 -2.120094e+06 1.260986e+06
4 -5.389179e+05 -2.120094e+06
5 8.706658e+03 -5.389179e+05
6 3.681755e+04 8.706658e+03
7 1.011745e+05 3.681755e+04
8 -1.759896e+05 1.011745e+05
9 6.212848e+05 -1.759896e+05
10 -1.804099e+05 6.212848e+05
11 2.960787e+05 -1.804099e+05
12 1.183419e+05 2.960787e+05
13 -1.118617e+06 1.183419e+05
14 -7.714344e+04 -1.118617e+06
15 5.469366e+05 -7.714344e+04
16 -2.810947e+05 5.469366e+05
17 -6.509910e+04 -2.810947e+05
18 -3.783602e+05 -6.509910e+04
19 -5.639902e+05 -3.783602e+05
20 2.003133e+05 -5.639902e+05
21 -1.692827e+05 2.003133e+05
22 -1.583710e+05 -1.692827e+05
23 -7.263975e+05 -1.583710e+05
24 -3.223750e+05 -7.263975e+05
25 -8.611960e+04 -3.223750e+05
26 -5.244867e+05 -8.611960e+04
27 -2.287674e+05 -5.244867e+05
28 3.491694e+05 -2.287674e+05
29 -4.180886e+04 3.491694e+05
30 2.450783e+05 -4.180886e+04
31 -1.286764e+05 2.450783e+05
32 2.148073e+05 -1.286764e+05
33 8.843046e+04 2.148073e+05
34 -6.741604e+04 8.843046e+04
35 -6.997293e+05 -6.741604e+04
36 3.290698e+04 -6.997293e+05
37 9.512216e+04 3.290698e+04
38 -1.292850e+05 9.512216e+04
39 -3.292756e+05 -1.292850e+05
40 2.265576e+05 -3.292756e+05
41 -3.381594e+05 2.265576e+05
42 1.305302e+05 -3.381594e+05
43 -2.193006e+05 1.305302e+05
44 -9.686861e+03 -2.193006e+05
45 -9.873843e+04 -9.686861e+03
46 -1.506070e+05 -9.873843e+04
47 -1.527146e+05 -1.506070e+05
48 -4.405757e+04 -1.527146e+05
49 -2.023685e+05 -4.405757e+04
50 -1.079644e+04 -2.023685e+05
51 -1.554851e+05 -1.079644e+04
52 -2.840178e+05 -1.554851e+05
53 9.597283e+04 -2.840178e+05
54 -4.233277e+05 9.597283e+04
55 1.614878e+04 -4.233277e+05
56 3.274203e+05 1.614878e+04
57 -1.344067e+04 3.274203e+05
58 -1.167402e+04 -1.344067e+04
59 -2.489276e+05 -1.167402e+04
60 -3.053437e+04 -2.489276e+05
61 -1.788236e+05 -3.053437e+04
62 -2.695521e+04 -1.788236e+05
63 -2.218082e+05 -2.695521e+04
64 -4.244611e+02 -2.218082e+05
65 2.313680e+05 -4.244611e+02
66 -3.961293e+00 2.313680e+05
67 -5.613649e+04 -3.961293e+00
68 -1.039046e+05 -5.613649e+04
69 -1.621399e+05 -1.039046e+05
70 -1.575332e+05 -1.621399e+05
71 8.402731e+03 -1.575332e+05
72 -1.022603e+05 8.402731e+03
73 1.619668e+04 -1.022603e+05
74 -4.670361e+04 1.619668e+04
75 -1.401200e+05 -4.670361e+04
76 -1.895204e+04 -1.401200e+05
77 -3.751264e+04 -1.895204e+04
78 1.864700e+04 -3.751264e+04
79 1.016586e+05 1.864700e+04
80 1.779951e+04 1.016586e+05
81 -9.639347e+04 1.779951e+04
82 -6.771926e+04 -9.639347e+04
83 4.597112e+04 -6.771926e+04
84 1.936319e+04 4.597112e+04
85 -1.526918e+05 1.936319e+04
86 -7.729575e+04 -1.526918e+05
87 -2.709022e+03 -7.729575e+04
88 -1.087286e+05 -2.709022e+03
89 2.374239e+04 -1.087286e+05
90 5.171042e+03 2.374239e+04
91 -6.890441e+04 5.171042e+03
92 8.655845e+04 -6.890441e+04
93 -1.670658e+05 8.655845e+04
94 4.413430e+04 -1.670658e+05
95 -7.960861e+03 4.413430e+04
96 6.597813e+04 -7.960861e+03
97 1.644330e+05 6.597813e+04
98 -4.679715e+04 1.644330e+05
99 -3.656027e+04 -4.679715e+04
100 -7.862713e+04 -3.656027e+04
101 -1.134372e+05 -7.862713e+04
102 -1.552571e+05 -1.134372e+05
103 -4.287334e+04 -1.552571e+05
104 3.778051e+04 -4.287334e+04
105 -6.726786e+04 3.778051e+04
106 -2.868120e+05 -6.726786e+04
107 9.374730e+04 -2.868120e+05
108 -1.737003e+05 9.374730e+04
109 -1.146023e+05 -1.737003e+05
110 9.946503e+04 -1.146023e+05
111 -2.168834e+05 9.946503e+04
112 -1.463099e+05 -2.168834e+05
113 -9.484045e+04 -1.463099e+05
114 5.483301e+04 -9.484045e+04
115 1.896172e+04 5.483301e+04
116 3.684722e+04 1.896172e+04
117 7.276375e+04 3.684722e+04
118 -1.645875e+05 7.276375e+04
119 6.720513e+04 -1.645875e+05
120 -3.177472e+05 6.720513e+04
121 3.423384e+04 -3.177472e+05
122 -1.454836e+05 3.423384e+04
123 4.344236e+04 -1.454836e+05
124 -1.585187e+03 4.344236e+04
125 -2.042616e+05 -1.585187e+03
126 5.862432e+04 -2.042616e+05
127 3.871832e+04 5.862432e+04
128 3.169185e+04 3.871832e+04
129 5.862432e+04 3.169185e+04
130 4.877076e+04 5.862432e+04
131 9.773557e+04 4.877076e+04
132 5.862432e+04 9.773557e+04
133 8.357447e+03 5.862432e+04
134 3.708742e+04 8.357447e+03
135 5.216156e+04 3.708742e+04
136 4.932280e+04 5.216156e+04
137 5.783791e+04 4.932280e+04
138 3.444756e+04 5.783791e+04
139 5.862432e+04 3.444756e+04
140 5.938854e+04 5.862432e+04
141 5.862432e+04 5.938854e+04
142 1.652187e+04 5.862432e+04
143 5.147936e+04 1.652187e+04
144 3.851483e+04 5.147936e+04
145 4.662419e+04 3.851483e+04
146 5.862432e+04 4.662419e+04
147 4.729126e+04 5.862432e+04
148 5.991813e+04 4.729126e+04
149 5.948051e+04 5.991813e+04
150 2.112212e+04 5.948051e+04
151 -1.920308e+04 2.112212e+04
152 4.612637e+04 -1.920308e+04
153 2.978399e+04 4.612637e+04
154 5.616392e+04 2.978399e+04
155 5.862432e+04 5.616392e+04
156 5.862432e+04 5.862432e+04
157 5.238087e+04 5.862432e+04
158 6.626112e+04 5.238087e+04
159 5.862432e+04 6.626112e+04
160 5.862432e+04 5.862432e+04
161 5.862432e+04 5.862432e+04
162 5.862432e+04 5.862432e+04
163 5.258312e+04 5.862432e+04
164 3.036140e+04 5.258312e+04
165 5.807737e+04 3.036140e+04
166 5.862432e+04 5.807737e+04
167 1.115811e+05 5.862432e+04
168 3.010201e+05 1.115811e+05
169 5.862432e+04 3.010201e+05
170 5.862432e+04 5.862432e+04
171 5.059639e+04 5.862432e+04
172 5.862432e+04 5.059639e+04
173 3.431578e+04 5.862432e+04
174 5.862432e+04 3.431578e+04
175 4.724455e+04 5.862432e+04
176 -5.870740e+04 4.724455e+04
177 2.198344e+04 -5.870740e+04
178 5.862432e+04 2.198344e+04
179 5.862432e+04 5.862432e+04
180 5.508597e+04 5.862432e+04
181 -9.384035e+03 5.508597e+04
182 1.574433e+04 -9.384035e+03
183 -2.986051e+03 1.574433e+04
184 1.076837e+04 -2.986051e+03
185 1.486915e+04 1.076837e+04
186 -2.250814e+04 1.486915e+04
187 4.192277e+04 -2.250814e+04
188 1.502288e+04 4.192277e+04
189 5.862432e+04 1.502288e+04
190 4.517986e+04 5.862432e+04
191 5.862432e+04 4.517986e+04
192 4.347433e+04 5.862432e+04
193 5.862432e+04 4.347433e+04
194 5.285609e+04 5.862432e+04
195 2.161509e+04 5.285609e+04
196 -1.976193e+05 2.161509e+04
197 2.161509e+04 -1.976193e+05
198 2.161509e+04 2.161509e+04
199 -9.102825e+04 2.161509e+04
200 2.161509e+04 -9.102825e+04
201 2.161509e+04 2.161509e+04
202 2.161509e+04 2.161509e+04
203 2.058634e+04 2.161509e+04
204 -3.388936e+04 2.058634e+04
205 1.042777e+04 -3.388936e+04
206 1.478032e+04 1.042777e+04
207 2.161509e+04 1.478032e+04
208 1.755700e+04 2.161509e+04
209 3.580431e+04 1.755700e+04
210 5.862432e+04 3.580431e+04
211 5.862432e+04 5.862432e+04
212 8.218123e+04 5.862432e+04
213 1.798101e+04 8.218123e+04
214 2.855512e+04 1.798101e+04
215 5.362514e+04 2.855512e+04
216 5.268297e+04 5.362514e+04
217 6.738501e+04 5.268297e+04
218 5.362514e+04 6.738501e+04
219 -3.909274e+03 5.362514e+04
220 3.112129e+04 -3.909274e+03
221 5.862432e+04 3.112129e+04
222 5.838509e+04 5.862432e+04
223 1.954647e+04 5.838509e+04
224 5.805901e+04 1.954647e+04
225 2.608423e+04 5.805901e+04
226 5.017709e+04 2.608423e+04
227 5.862432e+04 5.017709e+04
228 4.343789e+04 5.862432e+04
229 5.862432e+04 4.343789e+04
230 2.264195e+04 5.862432e+04
231 3.765605e+04 2.264195e+04
232 5.862432e+04 3.765605e+04
233 2.404469e+04 5.862432e+04
234 5.862432e+04 2.404469e+04
235 5.980690e+04 5.862432e+04
236 5.862432e+04 5.980690e+04
237 5.045136e+04 5.862432e+04
238 5.862432e+04 5.045136e+04
239 5.862432e+04 5.862432e+04
240 5.862432e+04 5.862432e+04
241 4.899359e+04 5.862432e+04
242 5.862432e+04 4.899359e+04
243 5.162708e+04 5.862432e+04
244 -3.389133e+03 5.162708e+04
245 2.161509e+04 -3.389133e+03
246 2.161509e+04 2.161509e+04
247 2.280813e+04 2.161509e+04
248 2.161509e+04 2.280813e+04
249 2.161509e+04 2.161509e+04
250 3.608742e+04 2.161509e+04
251 2.161509e+04 3.608742e+04
252 4.086943e+04 2.161509e+04
253 3.105108e+03 4.086943e+04
254 -1.241809e+05 3.105108e+03
255 5.862432e+04 -1.241809e+05
256 4.647750e+04 5.862432e+04
257 3.103620e+04 4.647750e+04
258 2.161509e+04 3.103620e+04
259 5.862432e+04 2.161509e+04
260 2.161509e+04 5.862432e+04
261 2.161509e+04 2.161509e+04
262 4.343946e+04 2.161509e+04
263 6.400187e+04 4.343946e+04
264 2.161509e+04 6.400187e+04
265 5.862432e+04 2.161509e+04
266 2.161509e+04 5.862432e+04
267 8.231452e+04 2.161509e+04
268 2.161509e+04 8.231452e+04
269 -1.559911e+03 2.161509e+04
270 2.161509e+04 -1.559911e+03
271 2.161509e+04 2.161509e+04
272 2.161509e+04 2.161509e+04
273 1.699497e+04 2.161509e+04
274 5.862432e+04 1.699497e+04
275 5.896816e+04 5.862432e+04
276 3.129980e+04 5.896816e+04
277 5.862432e+04 3.129980e+04
278 2.101390e+04 5.862432e+04
279 5.012522e+04 2.101390e+04
280 2.161509e+04 5.012522e+04
281 5.862432e+04 2.161509e+04
282 7.633316e+04 5.862432e+04
283 5.862432e+04 7.633316e+04
284 4.304069e+04 5.862432e+04
285 5.862432e+04 4.304069e+04
286 6.827479e+03 5.862432e+04
287 5.561020e+04 6.827479e+03
288 2.161509e+04 5.561020e+04
289 8.052626e+04 2.161509e+04
290 2.161509e+04 8.052626e+04
291 5.907775e+04 2.161509e+04
292 4.212129e+04 5.907775e+04
293 5.862432e+04 4.212129e+04
294 3.773515e+04 5.862432e+04
295 4.297017e+04 3.773515e+04
296 5.862432e+04 4.297017e+04
297 2.161509e+04 5.862432e+04
298 5.862432e+04 2.161509e+04
299 3.612801e+04 5.862432e+04
300 1.324818e+04 3.612801e+04
301 4.114021e+04 1.324818e+04
302 5.612473e+04 4.114021e+04
303 2.161509e+04 5.612473e+04
304 1.042476e+05 2.161509e+04
305 3.899164e+04 1.042476e+05
306 2.161509e+04 3.899164e+04
307 -1.960650e+04 2.161509e+04
308 3.777727e+04 -1.960650e+04
309 2.161509e+04 3.777727e+04
310 2.161509e+04 2.161509e+04
311 5.382871e+04 2.161509e+04
312 2.161509e+04 5.382871e+04
313 -4.185561e+03 2.161509e+04
314 4.339971e+04 -4.185561e+03
315 7.721015e+04 4.339971e+04
316 5.989160e+04 7.721015e+04
317 1.324597e+03 5.989160e+04
318 5.304290e+04 1.324597e+03
319 5.862432e+04 5.304290e+04
320 5.862432e+04 5.862432e+04
321 5.862432e+04 5.862432e+04
322 5.862432e+04 5.862432e+04
323 2.161509e+04 5.862432e+04
324 5.862432e+04 2.161509e+04
325 5.508003e+04 5.862432e+04
326 2.161509e+04 5.508003e+04
327 2.161509e+04 2.161509e+04
328 5.862432e+04 2.161509e+04
329 5.862432e+04 5.862432e+04
330 5.862432e+04 5.862432e+04
331 5.862432e+04 5.862432e+04
332 5.862432e+04 5.862432e+04
333 5.124433e+04 5.862432e+04
334 5.841418e+04 5.124433e+04
335 4.757854e+04 5.841418e+04
336 6.098280e+04 4.757854e+04
337 5.862432e+04 6.098280e+04
338 5.492330e+04 5.862432e+04
339 4.992116e+04 5.492330e+04
340 1.582974e+04 4.992116e+04
341 5.862432e+04 1.582974e+04
342 5.862432e+04 5.862432e+04
343 5.862432e+04 5.862432e+04
344 6.106346e+04 5.862432e+04
345 5.862432e+04 6.106346e+04
346 4.645356e+04 5.862432e+04
347 5.862432e+04 4.645356e+04
348 5.862432e+04 5.862432e+04
349 2.262242e+04 5.862432e+04
350 6.739839e+04 2.262242e+04
351 4.939214e+04 6.739839e+04
352 4.695392e+04 4.939214e+04
353 5.209069e+04 4.695392e+04
354 5.862432e+04 5.209069e+04
355 5.862432e+04 5.862432e+04
356 4.700609e+04 5.862432e+04
357 5.862432e+04 4.700609e+04
358 5.862432e+04 5.862432e+04
359 5.862432e+04 5.862432e+04
360 5.324862e+04 5.862432e+04
361 5.862432e+04 5.324862e+04
362 5.862432e+04 5.862432e+04
363 -6.743071e+04 5.862432e+04
364 4.776351e+04 -6.743071e+04
365 5.862432e+04 4.776351e+04
366 5.862432e+04 5.862432e+04
367 5.862432e+04 5.862432e+04
368 5.862432e+04 5.862432e+04
369 6.086752e+04 5.862432e+04
370 5.862432e+04 6.086752e+04
371 8.723399e+04 5.862432e+04
372 4.129707e+04 8.723399e+04
373 3.370318e+04 4.129707e+04
374 5.862432e+04 3.370318e+04
375 5.862432e+04 5.862432e+04
376 4.384601e+04 5.862432e+04
377 4.809817e+04 4.384601e+04
378 6.233059e+03 4.809817e+04
379 5.862432e+04 6.233059e+03
380 5.862432e+04 5.862432e+04
381 5.225427e+04 5.862432e+04
382 1.009405e+05 5.225427e+04
383 5.862432e+04 1.009405e+05
384 4.409872e+04 5.862432e+04
385 -1.107469e+05 4.409872e+04
386 6.868064e+04 -1.107469e+05
387 2.998467e+04 6.868064e+04
388 5.240005e+04 2.998467e+04
389 4.083358e+04 5.240005e+04
390 4.019934e+04 4.083358e+04
391 2.218813e+04 4.019934e+04
392 -1.927890e+04 2.218813e+04
393 7.782750e+04 -1.927890e+04
394 -1.311252e+05 7.782750e+04
395 1.344257e+04 -1.311252e+05
396 -1.231689e+05 1.344257e+04
397 -4.483172e+05 -1.231689e+05
398 1.566969e+04 -4.483172e+05
399 1.416191e+04 1.566969e+04
400 6.543617e+03 1.416191e+04
401 4.444421e+04 6.543617e+03
402 2.369665e+04 4.444421e+04
403 -8.590407e+04 2.369665e+04
404 -8.045107e+03 -8.590407e+04
405 1.546557e+04 -8.045107e+03
406 2.565864e+04 1.546557e+04
407 1.202447e+04 2.565864e+04
408 2.605950e+04 1.202447e+04
409 -6.015431e+04 2.605950e+04
410 -6.278733e+04 -6.015431e+04
411 -1.912285e+05 -6.278733e+04
412 -6.145847e+04 -1.912285e+05
413 1.652315e+04 -6.145847e+04
414 -1.744255e+05 1.652315e+04
415 -9.624296e+04 -1.744255e+05
416 1.410265e+05 -9.624296e+04
417 -6.088788e+04 1.410265e+05
418 -5.310840e+04 -6.088788e+04
419 -3.228747e+04 -5.310840e+04
420 -4.059382e+04 -3.228747e+04
421 -2.050639e+05 -4.059382e+04
422 -5.909810e+05 -2.050639e+05
423 -4.524967e+05 -5.909810e+05
424 -4.621551e+04 -4.524967e+05
425 -2.926978e+05 -4.621551e+04
426 -3.028997e+05 -2.926978e+05
427 -1.237773e+06 -3.028997e+05
428 -2.601184e+05 -1.237773e+06
429 -3.540133e+05 -2.601184e+05
430 -5.651288e+05 -3.540133e+05
> plot(z,main=paste('Residual Lag plot, lowess, and regression line'), ylab='values of Residuals', xlab='lagged values of Residuals')
> lines(lowess(z))
> abline(lm(z))
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/7mgwg1291392609.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> acf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Autocorrelation Function')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/8x7d11291392609.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> pacf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Partial Autocorrelation Function')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/9x7d11291392609.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> opar <- par(mfrow = c(2,2), oma = c(0, 0, 1.1, 0))
> plot(mylm, las = 1, sub='Residual Diagnostics')
> par(opar)
> dev.off()
null device
1
> if (n > n25) {
+ postscript(file="/var/www/html/rcomp/tmp/10x7d11291392609.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
+ plot(kp3:nmkm3,gqarr[,2], main='Goldfeld-Quandt test',ylab='2-sided p-value',xlab='breakpoint')
+ grid()
+ dev.off()
+ }
null device
1
>
> #Note: the /var/www/html/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/www/html/rcomp/createtable")
>
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Estimated Regression Equation', 1, TRUE)
> a<-table.row.end(a)
> myeq <- colnames(x)[1]
> myeq <- paste(myeq, '[t] = ', sep='')
> for (i in 1:k){
+ if (mysum$coefficients[i,1] > 0) myeq <- paste(myeq, '+', '')
+ myeq <- paste(myeq, mysum$coefficients[i,1], sep=' ')
+ if (rownames(mysum$coefficients)[i] != '(Intercept)') {
+ myeq <- paste(myeq, rownames(mysum$coefficients)[i], sep='')
+ if (rownames(mysum$coefficients)[i] != 't') myeq <- paste(myeq, '[t]', sep='')
+ }
+ }
> myeq <- paste(myeq, ' + e[t]')
> a<-table.row.start(a)
> a<-table.element(a, myeq)
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/111qcp1291392609.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,hyperlink('http://www.xycoon.com/ols1.htm','Multiple Linear Regression - Ordinary Least Squares',''), 6, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Variable',header=TRUE)
> a<-table.element(a,'Parameter',header=TRUE)
> a<-table.element(a,'S.D.',header=TRUE)
> a<-table.element(a,'T-STAT
H0: parameter = 0',header=TRUE)
> a<-table.element(a,'2-tail p-value',header=TRUE)
> a<-table.element(a,'1-tail p-value',header=TRUE)
> a<-table.row.end(a)
> for (i in 1:k){
+ a<-table.row.start(a)
+ a<-table.element(a,rownames(mysum$coefficients)[i],header=TRUE)
+ a<-table.element(a,mysum$coefficients[i,1])
+ a<-table.element(a, round(mysum$coefficients[i,2],6))
+ a<-table.element(a, round(mysum$coefficients[i,3],4))
+ a<-table.element(a, round(mysum$coefficients[i,4],6))
+ a<-table.element(a, round(mysum$coefficients[i,4]/2,6))
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/12m8tv1291392609.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Regression Statistics', 2, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple R',1,TRUE)
> a<-table.element(a, sqrt(mysum$r.squared))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'R-squared',1,TRUE)
> a<-table.element(a, mysum$r.squared)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Adjusted R-squared',1,TRUE)
> a<-table.element(a, mysum$adj.r.squared)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (value)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[1])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (DF numerator)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[2])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (DF denominator)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[3])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'p-value',1,TRUE)
> a<-table.element(a, 1-pf(mysum$fstatistic[1],mysum$fstatistic[2],mysum$fstatistic[3]))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Residual Statistics', 2, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Residual Standard Deviation',1,TRUE)
> a<-table.element(a, mysum$sigma)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Sum Squared Residuals',1,TRUE)
> a<-table.element(a, sum(myerror*myerror))
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/13t98o1291392609.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Actuals, Interpolation, and Residuals', 4, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Time or Index', 1, TRUE)
> a<-table.element(a, 'Actuals', 1, TRUE)
> a<-table.element(a, 'Interpolation
Forecast', 1, TRUE)
> a<-table.element(a, 'Residuals
Prediction Error', 1, TRUE)
> a<-table.row.end(a)
> for (i in 1:n) {
+ a<-table.row.start(a)
+ a<-table.element(a,i, 1, TRUE)
+ a<-table.element(a,x[i])
+ a<-table.element(a,x[i]-mysum$resid[i])
+ a<-table.element(a,mysum$resid[i])
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/1430p91291392609.tab")
> if (n > n25) {
+ a<-table.start()
+ a<-table.row.start(a)
+ a<-table.element(a,'Goldfeld-Quandt test for Heteroskedasticity',4,TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'p-values',header=TRUE)
+ a<-table.element(a,'Alternative Hypothesis',3,header=TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'breakpoint index',header=TRUE)
+ a<-table.element(a,'greater',header=TRUE)
+ a<-table.element(a,'2-sided',header=TRUE)
+ a<-table.element(a,'less',header=TRUE)
+ a<-table.row.end(a)
+ for (mypoint in kp3:nmkm3) {
+ a<-table.row.start(a)
+ a<-table.element(a,mypoint,header=TRUE)
+ a<-table.element(a,gqarr[mypoint-kp3+1,1])
+ a<-table.element(a,gqarr[mypoint-kp3+1,2])
+ a<-table.element(a,gqarr[mypoint-kp3+1,3])
+ a<-table.row.end(a)
+ }
+ a<-table.end(a)
+ table.save(a,file="/var/www/html/rcomp/tmp/15pj5f1291392609.tab")
+ a<-table.start()
+ a<-table.row.start(a)
+ a<-table.element(a,'Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity',4,TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'Description',header=TRUE)
+ a<-table.element(a,'# significant tests',header=TRUE)
+ a<-table.element(a,'% significant tests',header=TRUE)
+ a<-table.element(a,'OK/NOK',header=TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'1% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant1)
+ a<-table.element(a,numsignificant1/numgqtests)
+ if (numsignificant1/numgqtests < 0.01) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'5% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant5)
+ a<-table.element(a,numsignificant5/numgqtests)
+ if (numsignificant5/numgqtests < 0.05) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'10% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant10)
+ a<-table.element(a,numsignificant10/numgqtests)
+ if (numsignificant10/numgqtests < 0.1) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.end(a)
+ table.save(a,file="/var/www/html/rcomp/tmp/163t361291392609.tab")
+ }
>
> try(system("convert tmp/186hp1291392609.ps tmp/186hp1291392609.png",intern=TRUE))
character(0)
> try(system("convert tmp/2jfgs1291392609.ps tmp/2jfgs1291392609.png",intern=TRUE))
character(0)
> try(system("convert tmp/3jfgs1291392609.ps tmp/3jfgs1291392609.png",intern=TRUE))
character(0)
> try(system("convert tmp/4jfgs1291392609.ps tmp/4jfgs1291392609.png",intern=TRUE))
character(0)
> try(system("convert tmp/5coxd1291392609.ps tmp/5coxd1291392609.png",intern=TRUE))
character(0)
> try(system("convert tmp/6coxd1291392609.ps tmp/6coxd1291392609.png",intern=TRUE))
character(0)
> try(system("convert tmp/7mgwg1291392609.ps tmp/7mgwg1291392609.png",intern=TRUE))
character(0)
> try(system("convert tmp/8x7d11291392609.ps tmp/8x7d11291392609.png",intern=TRUE))
character(0)
> try(system("convert tmp/9x7d11291392609.ps tmp/9x7d11291392609.png",intern=TRUE))
character(0)
> try(system("convert tmp/10x7d11291392609.ps tmp/10x7d11291392609.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
12.522 2.144 28.291