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