R version 2.9.0 (2009-04-17)
Copyright (C) 2009 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> x <- array(list(4321023
+ ,29790
+ ,444
+ ,81767
+ ,1
+ ,1
+ ,4111912
+ ,87550
+ ,412
+ ,153198
+ ,1
+ ,1
+ ,223193
+ ,84738
+ ,428
+ ,-26007
+ ,0
+ ,1
+ ,1491348
+ ,54660
+ ,315
+ ,126942
+ ,1
+ ,1
+ ,1629616
+ ,42634
+ ,168
+ ,157214
+ ,1
+ ,0
+ ,1398893
+ ,40949
+ ,263
+ ,129352
+ ,0
+ ,1
+ ,1926517
+ ,45187
+ ,267
+ ,234817
+ ,1
+ ,0
+ ,983660
+ ,37704
+ ,228
+ ,60448
+ ,1
+ ,1
+ ,1443586
+ ,16275
+ ,129
+ ,47818
+ ,1
+ ,1
+ ,1073089
+ ,25830
+ ,104
+ ,245546
+ ,0
+ ,1
+ ,984885
+ ,12679
+ ,122
+ ,48020
+ ,0
+ ,1
+ ,1405225
+ ,18014
+ ,393
+ ,-1710
+ ,1
+ ,1
+ ,227132
+ ,43556
+ ,190
+ ,32648
+ ,0
+ ,1
+ ,929118
+ ,24811
+ ,280
+ ,95350
+ ,1
+ ,1
+ ,1071292
+ ,6575
+ ,63
+ ,151352
+ ,0
+ ,0
+ ,638830
+ ,7123
+ ,102
+ ,288170
+ ,0
+ ,0
+ ,856956
+ ,21950
+ ,265
+ ,114337
+ ,1
+ ,1
+ ,992426
+ ,37597
+ ,234
+ ,37884
+ ,1
+ ,1
+ ,444477
+ ,17821
+ ,277
+ ,122844
+ ,0
+ ,1
+ ,857217
+ ,12988
+ ,73
+ ,82340
+ ,1
+ ,1
+ ,711969
+ ,22330
+ ,67
+ ,79801
+ ,1
+ ,0
+ ,702380
+ ,13326
+ ,103
+ ,165548
+ ,0
+ ,0
+ ,358589
+ ,16189
+ ,290
+ ,116384
+ ,0
+ ,1
+ ,297978
+ ,7146
+ ,83
+ ,134028
+ ,0
+ ,1
+ ,585715
+ ,15824
+ ,56
+ ,63838
+ ,0
+ ,1
+ ,657954
+ ,27664
+ ,236
+ ,74996
+ ,1
+ ,1
+ ,209458
+ ,11920
+ ,73
+ ,31080
+ ,0
+ ,1
+ ,786690
+ ,8568
+ ,34
+ ,32168
+ ,0
+ ,1
+ ,439798
+ ,14416
+ ,139
+ ,49857
+ ,0
+ ,1
+ ,688779
+ ,3369
+ ,26
+ ,87161
+ ,1
+ ,1
+ ,574339
+ ,11819
+ ,70
+ ,106113
+ ,1
+ ,1
+ ,741409
+ ,6984
+ ,40
+ ,80570
+ ,1
+ ,0
+ ,597793
+ ,4519
+ ,42
+ ,102129
+ ,1
+ ,1
+ ,644190
+ ,2220
+ ,12
+ ,301670
+ ,0
+ ,1
+ ,377934
+ ,18562
+ ,211
+ ,102313
+ ,0
+ ,1
+ ,640273
+ ,10327
+ ,74
+ ,88577
+ ,0
+ ,1
+ ,697458
+ ,5336
+ ,80
+ ,112477
+ ,1
+ ,1
+ ,550608
+ ,2365
+ ,83
+ ,191778
+ ,1
+ ,1
+ ,207393
+ ,4069
+ ,131
+ ,79804
+ ,0
+ ,1
+ ,301607
+ ,8636
+ ,203
+ ,128294
+ ,0
+ ,0
+ ,345783
+ ,13718
+ ,56
+ ,96448
+ ,0
+ ,1
+ ,501749
+ ,4525
+ ,89
+ ,93811
+ ,0
+ ,0
+ ,379983
+ ,6869
+ ,88
+ ,117520
+ ,0
+ ,0
+ ,387475
+ ,4628
+ ,39
+ ,69159
+ ,0
+ ,0
+ ,377305
+ ,3689
+ ,25
+ ,101792
+ ,1
+ ,1
+ ,370837
+ ,4891
+ ,49
+ ,210568
+ ,1
+ ,1
+ ,430866
+ ,7489
+ ,149
+ ,136996
+ ,1
+ ,1
+ ,469107
+ ,4901
+ ,58
+ ,121920
+ ,0
+ ,1
+ ,194493
+ ,2284
+ ,41
+ ,76403
+ ,0
+ ,1
+ ,530670
+ ,3160
+ ,90
+ ,108094
+ ,1
+ ,1
+ ,518365
+ ,4150
+ ,136
+ ,134759
+ ,1
+ ,1
+ ,491303
+ ,7285
+ ,97
+ ,188873
+ ,1
+ ,1
+ ,527021
+ ,1134
+ ,63
+ ,146216
+ ,1
+ ,0
+ ,233773
+ ,4658
+ ,114
+ ,156608
+ ,1
+ ,1
+ ,405972
+ ,2384
+ ,77
+ ,61348
+ ,0
+ ,0
+ ,652925
+ ,3748
+ ,6
+ ,50350
+ ,0
+ ,1
+ ,446211
+ ,5371
+ ,47
+ ,87720
+ ,0
+ ,0
+ ,341340
+ ,1285
+ ,51
+ ,99489
+ ,0
+ ,1
+ ,387699
+ ,9327
+ ,85
+ ,87419
+ ,1
+ ,1
+ ,493408
+ ,5565
+ ,43
+ ,94355
+ ,1
+ ,0
+ ,146494
+ ,1528
+ ,32
+ ,60326
+ ,0
+ ,0
+ ,414462
+ ,3122
+ ,25
+ ,94670
+ ,1
+ ,1
+ ,364304
+ ,7561
+ ,77
+ ,82425
+ ,1
+ ,0
+ ,355178
+ ,2675
+ ,54
+ ,59017
+ ,0
+ ,0
+ ,357760
+ ,13253
+ ,251
+ ,90829
+ ,0
+ ,1
+ ,261216
+ ,880
+ ,15
+ ,80791
+ ,0
+ ,1
+ ,397144
+ ,2053
+ ,44
+ ,100423
+ ,1
+ ,1
+ ,374943
+ ,1424
+ ,73
+ ,131116
+ ,0
+ ,1
+ ,424898
+ ,4036
+ ,85
+ ,100269
+ ,1
+ ,1
+ ,202055
+ ,3045
+ ,49
+ ,27330
+ ,1
+ ,1
+ ,378525
+ ,5119
+ ,38
+ ,39039
+ ,0
+ ,0
+ ,310768
+ ,1431
+ ,35
+ ,106885
+ ,0
+ ,0
+ ,325738
+ ,554
+ ,9
+ ,79285
+ ,0
+ ,1
+ ,394510
+ ,1975
+ ,34
+ ,118881
+ ,0
+ ,1
+ ,247060
+ ,1765
+ ,20
+ ,77623
+ ,1
+ ,1
+ ,368078
+ ,1012
+ ,29
+ ,114768
+ ,0
+ ,0
+ ,236761
+ ,810
+ ,11
+ ,74015
+ ,0
+ ,0
+ ,312378
+ ,1280
+ ,52
+ ,69465
+ ,0
+ ,1
+ ,339836
+ ,666
+ ,13
+ ,117869
+ ,1
+ ,1
+ ,347385
+ ,1380
+ ,29
+ ,60982
+ ,0
+ ,1
+ ,426280
+ ,4677
+ ,66
+ ,90131
+ ,1
+ ,0
+ ,352850
+ ,876
+ ,33
+ ,138971
+ ,0
+ ,1
+ ,301881
+ ,814
+ ,15
+ ,39625
+ ,0
+ ,0
+ ,377516
+ ,514
+ ,15
+ ,102725
+ ,0
+ ,1
+ ,357312
+ ,5692
+ ,68
+ ,64239
+ ,1
+ ,1
+ ,458343
+ ,3642
+ ,100
+ ,90262
+ ,0
+ ,1
+ ,354228
+ ,540
+ ,13
+ ,103960
+ ,0
+ ,0
+ ,308636
+ ,2099
+ ,45
+ ,106611
+ ,0
+ ,0
+ ,386212
+ ,567
+ ,14
+ ,103345
+ ,0
+ ,0
+ ,393343
+ ,2001
+ ,36
+ ,95551
+ ,0
+ ,0
+ ,378509
+ ,2949
+ ,40
+ ,82903
+ ,1
+ ,1
+ ,452469
+ ,2253
+ ,68
+ ,63593
+ ,0
+ ,1
+ ,364839
+ ,6533
+ ,29
+ ,126910
+ ,1
+ ,1
+ ,358649
+ ,1889
+ ,43
+ ,37527
+ ,0
+ ,1
+ ,376641
+ ,3055
+ ,30
+ ,60247
+ ,1
+ ,0
+ ,429112
+ ,272
+ ,9
+ ,112995
+ ,0
+ ,1
+ ,330546
+ ,1414
+ ,22
+ ,70184
+ ,1
+ ,1
+ ,403560
+ ,2564
+ ,19
+ ,130140
+ ,0
+ ,0
+ ,317892
+ ,1383
+ ,9
+ ,73221
+ ,1
+ ,1
+ ,307528
+ ,1261
+ ,31
+ ,76114
+ ,1
+ ,1
+ ,235133
+ ,975
+ ,19
+ ,90534
+ ,0
+ ,1
+ ,299243
+ ,3366
+ ,55
+ ,108479
+ ,0
+ ,1
+ ,314073
+ ,576
+ ,8
+ ,113761
+ ,0
+ ,1
+ ,368186
+ ,1686
+ ,28
+ ,68696
+ ,0
+ ,1
+ ,269661
+ ,746
+ ,29
+ ,71561
+ ,0
+ ,1
+ ,125390
+ ,3192
+ ,48
+ ,59831
+ ,1
+ ,1
+ ,510834
+ ,2045
+ ,16
+ ,97890
+ ,1
+ ,1
+ ,321896
+ ,5702
+ ,47
+ ,101481
+ ,0
+ ,0
+ ,249898
+ ,1932
+ ,20
+ ,72954
+ ,1
+ ,1
+ ,408881
+ ,936
+ ,22
+ ,67939
+ ,0
+ ,1
+ ,158492
+ ,3437
+ ,33
+ ,48022
+ ,1
+ ,1
+ ,292154
+ ,5131
+ ,44
+ ,86111
+ ,0
+ ,0
+ ,289513
+ ,2397
+ ,13
+ ,74020
+ ,1
+ ,1
+ ,378049
+ ,1389
+ ,6
+ ,57530
+ ,1
+ ,1
+ ,343466
+ ,1503
+ ,35
+ ,56364
+ ,0
+ ,0
+ ,332743
+ ,402
+ ,8
+ ,84990
+ ,0
+ ,0
+ ,442882
+ ,2239
+ ,17
+ ,88590
+ ,0
+ ,0
+ ,214215
+ ,2234
+ ,11
+ ,77200
+ ,1
+ ,1
+ ,315688
+ ,837
+ ,21
+ ,61262
+ ,0
+ ,1
+ ,375195
+ ,10579
+ ,92
+ ,110309
+ ,0
+ ,1
+ ,334280
+ ,875
+ ,12
+ ,67000
+ ,0
+ ,0
+ ,355864
+ ,1585
+ ,112
+ ,93099
+ ,0
+ ,1
+ ,480382
+ ,1659
+ ,25
+ ,107577
+ ,1
+ ,1
+ ,353058
+ ,2647
+ ,17
+ ,62920
+ ,1
+ ,1
+ ,217193
+ ,3294
+ ,23
+ ,75832
+ ,1
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,314533
+ ,94
+ ,10
+ ,60793
+ ,0
+ ,0
+ ,318056
+ ,422
+ ,23
+ ,57935
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1
+ ,314353
+ ,34
+ ,7
+ ,60630
+ ,0
+ ,0
+ ,369448
+ ,1558
+ ,25
+ ,55637
+ ,0
+ ,1
+ ,315380
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,1
+ ,312846
+ ,43
+ ,20
+ ,60887
+ ,0
+ ,0
+ ,312075
+ ,645
+ ,4
+ ,60720
+ ,0
+ ,0
+ ,315009
+ ,316
+ ,4
+ ,60505
+ ,0
+ ,0
+ ,318903
+ ,115
+ ,10
+ ,60945
+ ,0
+ ,1
+ ,314887
+ ,5
+ ,1
+ ,60720
+ ,0
+ ,0
+ ,314913
+ ,897
+ ,4
+ ,60720
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325506
+ ,389
+ ,8
+ ,58990
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1
+ ,298568
+ ,1002
+ ,11
+ ,56750
+ ,0
+ ,0
+ ,315834
+ ,36
+ ,4
+ ,60894
+ ,0
+ ,0
+ ,329784
+ ,460
+ ,15
+ ,63346
+ ,0
+ ,1
+ ,312878
+ ,309
+ ,9
+ ,56535
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1
+ ,314987
+ ,9
+ ,7
+ ,60835
+ ,0
+ ,0
+ ,325249
+ ,271
+ ,2
+ ,60720
+ ,0
+ ,1
+ ,315877
+ ,14
+ ,0
+ ,61016
+ ,0
+ ,1
+ ,291650
+ ,520
+ ,7
+ ,58650
+ ,0
+ ,0
+ ,305959
+ ,1766
+ ,46
+ ,60438
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,5
+ ,60720
+ ,0
+ ,0
+ ,297765
+ ,458
+ ,7
+ ,58625
+ ,0
+ ,1
+ ,315245
+ ,20
+ ,2
+ ,60938
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,315236
+ ,98
+ ,2
+ ,61490
+ ,0
+ ,0
+ ,336425
+ ,405
+ ,5
+ ,60845
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,306268
+ ,483
+ ,7
+ ,60830
+ ,0
+ ,0
+ ,302187
+ ,454
+ ,24
+ ,63261
+ ,0
+ ,1
+ ,314882
+ ,47
+ ,1
+ ,60720
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1
+ ,382712
+ ,757
+ ,18
+ ,45689
+ ,0
+ ,1
+ ,341570
+ ,4655
+ ,55
+ ,60720
+ ,0
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,312412
+ ,36
+ ,3
+ ,61564
+ ,0
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1
+ ,309596
+ ,203
+ ,9
+ ,61938
+ ,0
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1
+ ,315547
+ ,126
+ ,8
+ ,60951
+ ,0
+ ,0
+ ,313267
+ ,400
+ ,113
+ ,60720
+ ,0
+ ,0
+ ,316176
+ ,71
+ ,0
+ ,60745
+ ,1
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,359335
+ ,972
+ ,19
+ ,71642
+ ,0
+ ,1
+ ,330068
+ ,531
+ ,11
+ ,71641
+ ,1
+ ,0
+ ,314289
+ ,2461
+ ,25
+ ,55792
+ ,0
+ ,0
+ ,297413
+ ,378
+ ,16
+ ,71873
+ ,1
+ ,0
+ ,314806
+ ,23
+ ,5
+ ,62555
+ ,1
+ ,1
+ ,333210
+ ,638
+ ,11
+ ,60370
+ ,1
+ ,1
+ ,352108
+ ,2300
+ ,23
+ ,64873
+ ,1
+ ,0
+ ,313332
+ ,149
+ ,6
+ ,62041
+ ,0
+ ,0
+ ,291787
+ ,226
+ ,5
+ ,65745
+ ,0
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,318745
+ ,275
+ ,7
+ ,59500
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,315366
+ ,141
+ ,7
+ ,61630
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1
+ ,315688
+ ,28
+ ,3
+ ,60890
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,1
+ ,409642
+ ,4980
+ ,89
+ ,113521
+ ,1
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,1
+ ,269587
+ ,472
+ ,19
+ ,80045
+ ,1
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,1
+ ,300962
+ ,203
+ ,12
+ ,50804
+ ,1
+ ,0
+ ,325479
+ ,496
+ ,12
+ ,87390
+ ,1
+ ,1
+ ,316155
+ ,10
+ ,5
+ ,61656
+ ,1
+ ,1
+ ,318574
+ ,63
+ ,2
+ ,65688
+ ,1
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,1
+ ,343613
+ ,1136
+ ,26
+ ,48522
+ ,1
+ ,0
+ ,306948
+ ,265
+ ,3
+ ,60720
+ ,1
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1
+ ,330059
+ ,267
+ ,11
+ ,57640
+ ,0
+ ,1
+ ,288985
+ ,474
+ ,10
+ ,61977
+ ,0
+ ,1
+ ,304485
+ ,534
+ ,5
+ ,62620
+ ,0
+ ,1
+ ,315380
+ ,0
+ ,2
+ ,60720
+ ,0
+ ,1
+ ,315688
+ ,15
+ ,6
+ ,60831
+ ,0
+ ,0
+ ,317736
+ ,397
+ ,7
+ ,60646
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,2
+ ,60720
+ ,0
+ ,0
+ ,322331
+ ,1866
+ ,28
+ ,56225
+ ,0
+ ,0
+ ,296656
+ ,288
+ ,3
+ ,60510
+ ,0
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1
+ ,315354
+ ,3
+ ,1
+ ,60698
+ ,0
+ ,1
+ ,312161
+ ,468
+ ,20
+ ,60720
+ ,0
+ ,0
+ ,315576
+ ,20
+ ,1
+ ,60805
+ ,0
+ ,0
+ ,314922
+ ,278
+ ,22
+ ,61404
+ ,0
+ ,0
+ ,314551
+ ,61
+ ,9
+ ,60720
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1
+ ,312339
+ ,192
+ ,2
+ ,65276
+ ,0
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,298700
+ ,317
+ ,7
+ ,63915
+ ,0
+ ,1
+ ,321376
+ ,738
+ ,9
+ ,60720
+ ,0
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1
+ ,303230
+ ,368
+ ,13
+ ,61686
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,315487
+ ,2
+ ,0
+ ,60743
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,315793
+ ,53
+ ,6
+ ,60349
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,312887
+ ,94
+ ,3
+ ,61360
+ ,0
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,315637
+ ,24
+ ,7
+ ,59818
+ ,0
+ ,1
+ ,324385
+ ,2332
+ ,2
+ ,72680
+ ,0
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,1
+ ,308989
+ ,131
+ ,15
+ ,61808
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,296702
+ ,206
+ ,9
+ ,53110
+ ,0
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,1
+ ,307322
+ ,167
+ ,1
+ ,64245
+ ,0
+ ,1
+ ,304376
+ ,622
+ ,38
+ ,73007
+ ,0
+ ,0
+ ,253588
+ ,2328
+ ,57
+ ,82732
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,309560
+ ,365
+ ,7
+ ,54820
+ ,0
+ ,0
+ ,298466
+ ,364
+ ,26
+ ,47705
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,1
+ ,343929
+ ,226
+ ,13
+ ,72835
+ ,0
+ ,0
+ ,331955
+ ,307
+ ,10
+ ,58856
+ ,0
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,381180
+ ,188
+ ,9
+ ,77655
+ ,0
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,1
+ ,331420
+ ,138
+ ,26
+ ,69817
+ ,0
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,1
+ ,310201
+ ,125
+ ,19
+ ,60798
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,320016
+ ,282
+ ,12
+ ,62452
+ ,0
+ ,1
+ ,320398
+ ,335
+ ,23
+ ,64175
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1
+ ,291841
+ ,1324
+ ,29
+ ,67440
+ ,1
+ ,1
+ ,310670
+ ,176
+ ,8
+ ,68136
+ ,0
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,313491
+ ,249
+ ,26
+ ,56726
+ ,0
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,331323
+ ,333
+ ,9
+ ,70811
+ ,0
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1
+ ,319210
+ ,601
+ ,5
+ ,60720
+ ,1
+ ,1
+ ,318098
+ ,30
+ ,3
+ ,62045
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,292754
+ ,249
+ ,13
+ ,54323
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,1
+ ,325176
+ ,165
+ ,12
+ ,62841
+ ,0
+ ,0
+ ,365959
+ ,453
+ ,19
+ ,81125
+ ,0
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,302409
+ ,53
+ ,10
+ ,59506
+ ,0
+ ,1
+ ,340968
+ ,382
+ ,9
+ ,59365
+ ,1
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1
+ ,315380
+ ,0
+ ,9
+ ,60720
+ ,0
+ ,0
+ ,313164
+ ,30
+ ,4
+ ,60798
+ ,1
+ ,1
+ ,301164
+ ,290
+ ,1
+ ,58790
+ ,0
+ ,1
+ ,315380
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,1
+ ,344425
+ ,366
+ ,14
+ ,61808
+ ,0
+ ,0
+ ,315394
+ ,2
+ ,12
+ ,60735
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,1
+ ,316647
+ ,209
+ ,19
+ ,64016
+ ,1
+ ,1
+ ,309836
+ ,384
+ ,17
+ ,54683
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,1
+ ,346611
+ ,365
+ ,32
+ ,87192
+ ,0
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,1
+ ,322031
+ ,49
+ ,14
+ ,64107
+ ,1
+ ,0
+ ,315656
+ ,3
+ ,8
+ ,60761
+ ,0
+ ,0
+ ,339445
+ ,133
+ ,4
+ ,65990
+ ,0
+ ,0
+ ,314964
+ ,32
+ ,0
+ ,59988
+ ,0
+ ,0
+ ,297141
+ ,368
+ ,20
+ ,61167
+ ,0
+ ,0
+ ,315372
+ ,1
+ ,5
+ ,60719
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1
+ ,312502
+ ,22
+ ,1
+ ,60722
+ ,0
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1
+ ,313729
+ ,96
+ ,4
+ ,60379
+ ,0
+ ,1
+ ,315388
+ ,1
+ ,1
+ ,60727
+ ,0
+ ,0
+ ,315371
+ ,314
+ ,4
+ ,60720
+ ,0
+ ,0
+ ,296139
+ ,844
+ ,20
+ ,60925
+ ,0
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,313880
+ ,26
+ ,1
+ ,60896
+ ,0
+ ,0
+ ,317698
+ ,125
+ ,10
+ ,59734
+ ,0
+ ,1
+ ,295580
+ ,304
+ ,12
+ ,62969
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1
+ ,308256
+ ,621
+ ,13
+ ,60720
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,303677
+ ,119
+ ,3
+ ,59118
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1
+ ,319369
+ ,1595
+ ,10
+ ,60720
+ ,0
+ ,0
+ ,318690
+ ,312
+ ,3
+ ,58598
+ ,0
+ ,1
+ ,314049
+ ,60
+ ,7
+ ,61124
+ ,0
+ ,0
+ ,325699
+ ,587
+ ,10
+ ,59595
+ ,0
+ ,1
+ ,314210
+ ,135
+ ,1
+ ,62065
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,322378
+ ,514
+ ,15
+ ,78780
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1
+ ,315398
+ ,1
+ ,4
+ ,60722
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,308336
+ ,1763
+ ,28
+ ,61600
+ ,1
+ ,1
+ ,316386
+ ,180
+ ,9
+ ,59635
+ ,0
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1
+ ,315553
+ ,218
+ ,7
+ ,60720
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1
+ ,323361
+ ,448
+ ,7
+ ,59781
+ ,0
+ ,1
+ ,336639
+ ,227
+ ,7
+ ,76644
+ ,0
+ ,0
+ ,307424
+ ,174
+ ,3
+ ,64820
+ ,0
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1
+ ,295370
+ ,121
+ ,11
+ ,56178
+ ,0
+ ,1
+ ,322340
+ ,607
+ ,7
+ ,60436
+ ,0
+ ,0
+ ,319864
+ ,2212
+ ,10
+ ,60720
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1
+ ,317291
+ ,530
+ ,18
+ ,73433
+ ,0
+ ,0
+ ,280398
+ ,571
+ ,14
+ ,41477
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1
+ ,317330
+ ,78
+ ,12
+ ,62700
+ ,0
+ ,0
+ ,238125
+ ,2489
+ ,29
+ ,67804
+ ,0
+ ,1
+ ,327071
+ ,131
+ ,3
+ ,59661
+ ,0
+ ,0
+ ,309038
+ ,923
+ ,6
+ ,58620
+ ,0
+ ,1
+ ,314210
+ ,72
+ ,3
+ ,60398
+ ,0
+ ,0
+ ,307930
+ ,572
+ ,8
+ ,58580
+ ,0
+ ,1
+ ,322327
+ ,397
+ ,10
+ ,62710
+ ,0
+ ,1
+ ,292136
+ ,450
+ ,6
+ ,59325
+ ,0
+ ,0
+ ,263276
+ ,622
+ ,8
+ ,60950
+ ,0
+ ,1
+ ,367655
+ ,694
+ ,6
+ ,68060
+ ,0
+ ,0
+ ,283910
+ ,3425
+ ,9
+ ,83620
+ ,1
+ ,0
+ ,283587
+ ,562
+ ,8
+ ,58456
+ ,0
+ ,1
+ ,243650
+ ,4917
+ ,26
+ ,52811
+ ,0
+ ,1
+ ,438493
+ ,1442
+ ,239
+ ,121173
+ ,1
+ ,1
+ ,296261
+ ,529
+ ,7
+ ,63870
+ ,0
+ ,1
+ ,230621
+ ,2126
+ ,41
+ ,21001
+ ,1
+ ,1
+ ,304252
+ ,1061
+ ,3
+ ,70415
+ ,0
+ ,0
+ ,333505
+ ,776
+ ,8
+ ,64230
+ ,0
+ ,0
+ ,296919
+ ,611
+ ,6
+ ,59190
+ ,0
+ ,0
+ ,278990
+ ,1526
+ ,21
+ ,69351
+ ,1
+ ,0
+ ,276898
+ ,592
+ ,7
+ ,64270
+ ,0
+ ,0
+ ,327007
+ ,1182
+ ,11
+ ,70694
+ ,0
+ ,0
+ ,317046
+ ,621
+ ,11
+ ,68005
+ ,0
+ ,1
+ ,304555
+ ,989
+ ,12
+ ,58930
+ ,0
+ ,1
+ ,298096
+ ,438
+ ,9
+ ,58320
+ ,0
+ ,0
+ ,231861
+ ,726
+ ,3
+ ,69980
+ ,0
+ ,1
+ ,309422
+ ,1303
+ ,57
+ ,69863
+ ,0
+ ,1
+ ,286963
+ ,7419
+ ,21
+ ,63255
+ ,1
+ ,1
+ ,269753
+ ,1164
+ ,15
+ ,57320
+ ,1
+ ,1
+ ,448243
+ ,3310
+ ,32
+ ,75230
+ ,1
+ ,1
+ ,165404
+ ,1920
+ ,11
+ ,79420
+ ,0
+ ,1
+ ,204325
+ ,965
+ ,2
+ ,73490
+ ,0
+ ,0
+ ,407159
+ ,3256
+ ,23
+ ,35250
+ ,0
+ ,1
+ ,290476
+ ,1135
+ ,20
+ ,62285
+ ,1
+ ,0
+ ,275311
+ ,1270
+ ,24
+ ,69206
+ ,0
+ ,0
+ ,246541
+ ,661
+ ,1
+ ,65920
+ ,0
+ ,0
+ ,253468
+ ,1013
+ ,1
+ ,69770
+ ,0
+ ,0
+ ,240897
+ ,2844
+ ,74
+ ,72683
+ ,0
+ ,1
+ ,-83265
+ ,11528
+ ,68
+ ,-14545
+ ,1
+ ,1
+ ,-42143
+ ,6526
+ ,20
+ ,55830
+ ,0
+ ,1
+ ,272713
+ ,2264
+ ,20
+ ,55174
+ ,0
+ ,1
+ ,215362
+ ,5109
+ ,82
+ ,67038
+ ,1
+ ,1
+ ,42754
+ ,3999
+ ,21
+ ,51252
+ ,0
+ ,1
+ ,306275
+ ,35624
+ ,244
+ ,157278
+ ,0
+ ,1
+ ,253537
+ ,9252
+ ,32
+ ,79510
+ ,0
+ ,1
+ ,372631
+ ,15236
+ ,86
+ ,77440
+ ,0
+ ,1
+ ,-7170
+ ,18073
+ ,69
+ ,27284
+ ,0
+ ,1)
+ ,dim=c(6
+ ,430)
+ ,dimnames=list(c('Wealth'
+ ,'Costs'
+ ,'Orders'
+ ,'Dividends'
+ ,'Group'
+ ,'Gender
')
+ ,1:430))
> y <- array(NA,dim=c(6,430),dimnames=list(c('Wealth','Costs','Orders','Dividends','Group','Gender
'),1:430))
> for (i in 1:dim(x)[1])
+ {
+ for (j in 1:dim(x)[2])
+ {
+ y[i,j] <- as.numeric(x[i,j])
+ }
+ }
> par3 = 'No Linear Trend'
> par2 = 'Do not include Seasonal Dummies'
> par1 = '1'
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Prof. Dr. P. Wessa
> #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #Technical description: Write here your technical program description (don't use hard returns!)
> library(lattice)
> library(lmtest)
Loading required package: zoo
Attaching package: 'zoo'
The following object(s) are masked from package:base :
as.Date.numeric
> n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test
> par1 <- as.numeric(par1)
> x <- t(y)
> k <- length(x[1,])
> n <- length(x[,1])
> x1 <- cbind(x[,par1], x[,1:k!=par1])
> mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1])
> colnames(x1) <- mycolnames #colnames(x)[par1]
> x <- x1
> if (par3 == 'First Differences'){
+ x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep='')))
+ for (i in 1:n-1) {
+ for (j in 1:k) {
+ x2[i,j] <- x[i+1,j] - x[i,j]
+ }
+ }
+ x <- x2
+ }
> if (par2 == 'Include Monthly Dummies'){
+ x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep ='')))
+ for (i in 1:11){
+ x2[seq(i,n,12),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> if (par2 == 'Include Quarterly Dummies'){
+ x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep ='')))
+ for (i in 1:3){
+ x2[seq(i,n,4),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> k <- length(x[1,])
> if (par3 == 'Linear Trend'){
+ x <- cbind(x, c(1:n))
+ colnames(x)[k+1] <- 't'
+ }
> x
Wealth Costs Orders Dividends Group Gender\r
1 4321023 29790 444 81767 1 1
2 4111912 87550 412 153198 1 1
3 223193 84738 428 -26007 0 1
4 1491348 54660 315 126942 1 1
5 1629616 42634 168 157214 1 0
6 1398893 40949 263 129352 0 1
7 1926517 45187 267 234817 1 0
8 983660 37704 228 60448 1 1
9 1443586 16275 129 47818 1 1
10 1073089 25830 104 245546 0 1
11 984885 12679 122 48020 0 1
12 1405225 18014 393 -1710 1 1
13 227132 43556 190 32648 0 1
14 929118 24811 280 95350 1 1
15 1071292 6575 63 151352 0 0
16 638830 7123 102 288170 0 0
17 856956 21950 265 114337 1 1
18 992426 37597 234 37884 1 1
19 444477 17821 277 122844 0 1
20 857217 12988 73 82340 1 1
21 711969 22330 67 79801 1 0
22 702380 13326 103 165548 0 0
23 358589 16189 290 116384 0 1
24 297978 7146 83 134028 0 1
25 585715 15824 56 63838 0 1
26 657954 27664 236 74996 1 1
27 209458 11920 73 31080 0 1
28 786690 8568 34 32168 0 1
29 439798 14416 139 49857 0 1
30 688779 3369 26 87161 1 1
31 574339 11819 70 106113 1 1
32 741409 6984 40 80570 1 0
33 597793 4519 42 102129 1 1
34 644190 2220 12 301670 0 1
35 377934 18562 211 102313 0 1
36 640273 10327 74 88577 0 1
37 697458 5336 80 112477 1 1
38 550608 2365 83 191778 1 1
39 207393 4069 131 79804 0 1
40 301607 8636 203 128294 0 0
41 345783 13718 56 96448 0 1
42 501749 4525 89 93811 0 0
43 379983 6869 88 117520 0 0
44 387475 4628 39 69159 0 0
45 377305 3689 25 101792 1 1
46 370837 4891 49 210568 1 1
47 430866 7489 149 136996 1 1
48 469107 4901 58 121920 0 1
49 194493 2284 41 76403 0 1
50 530670 3160 90 108094 1 1
51 518365 4150 136 134759 1 1
52 491303 7285 97 188873 1 1
53 527021 1134 63 146216 1 0
54 233773 4658 114 156608 1 1
55 405972 2384 77 61348 0 0
56 652925 3748 6 50350 0 1
57 446211 5371 47 87720 0 0
58 341340 1285 51 99489 0 1
59 387699 9327 85 87419 1 1
60 493408 5565 43 94355 1 0
61 146494 1528 32 60326 0 0
62 414462 3122 25 94670 1 1
63 364304 7561 77 82425 1 0
64 355178 2675 54 59017 0 0
65 357760 13253 251 90829 0 1
66 261216 880 15 80791 0 1
67 397144 2053 44 100423 1 1
68 374943 1424 73 131116 0 1
69 424898 4036 85 100269 1 1
70 202055 3045 49 27330 1 1
71 378525 5119 38 39039 0 0
72 310768 1431 35 106885 0 0
73 325738 554 9 79285 0 1
74 394510 1975 34 118881 0 1
75 247060 1765 20 77623 1 1
76 368078 1012 29 114768 0 0
77 236761 810 11 74015 0 0
78 312378 1280 52 69465 0 1
79 339836 666 13 117869 1 1
80 347385 1380 29 60982 0 1
81 426280 4677 66 90131 1 0
82 352850 876 33 138971 0 1
83 301881 814 15 39625 0 0
84 377516 514 15 102725 0 1
85 357312 5692 68 64239 1 1
86 458343 3642 100 90262 0 1
87 354228 540 13 103960 0 0
88 308636 2099 45 106611 0 0
89 386212 567 14 103345 0 0
90 393343 2001 36 95551 0 0
91 378509 2949 40 82903 1 1
92 452469 2253 68 63593 0 1
93 364839 6533 29 126910 1 1
94 358649 1889 43 37527 0 1
95 376641 3055 30 60247 1 0
96 429112 272 9 112995 0 1
97 330546 1414 22 70184 1 1
98 403560 2564 19 130140 0 0
99 317892 1383 9 73221 1 1
100 307528 1261 31 76114 1 1
101 235133 975 19 90534 0 1
102 299243 3366 55 108479 0 1
103 314073 576 8 113761 0 1
104 368186 1686 28 68696 0 1
105 269661 746 29 71561 0 1
106 125390 3192 48 59831 1 1
107 510834 2045 16 97890 1 1
108 321896 5702 47 101481 0 0
109 249898 1932 20 72954 1 1
110 408881 936 22 67939 0 1
111 158492 3437 33 48022 1 1
112 292154 5131 44 86111 0 0
113 289513 2397 13 74020 1 1
114 378049 1389 6 57530 1 1
115 343466 1503 35 56364 0 0
116 332743 402 8 84990 0 0
117 442882 2239 17 88590 0 0
118 214215 2234 11 77200 1 1
119 315688 837 21 61262 0 1
120 375195 10579 92 110309 0 1
121 334280 875 12 67000 0 0
122 355864 1585 112 93099 0 1
123 480382 1659 25 107577 1 1
124 353058 2647 17 62920 1 1
125 217193 3294 23 75832 1 0
126 315380 0 0 60720 0 0
127 314533 94 10 60793 0 0
128 318056 422 23 57935 0 0
129 315380 0 0 60720 0 1
130 314353 34 7 60630 0 0
131 369448 1558 25 55637 0 1
132 315380 0 1 60720 0 1
133 312846 43 20 60887 0 0
134 312075 645 4 60720 0 0
135 315009 316 4 60505 0 0
136 318903 115 10 60945 0 1
137 314887 5 1 60720 0 0
138 314913 897 4 60720 0 0
139 315380 0 0 60720 0 0
140 325506 389 8 58990 0 0
141 315380 0 0 60720 0 1
142 298568 1002 11 56750 0 0
143 315834 36 4 60894 0 0
144 329784 460 15 63346 0 1
145 312878 309 9 56535 0 0
146 315380 0 0 60720 0 1
147 314987 9 7 60835 0 0
148 325249 271 2 60720 0 1
149 315877 14 0 61016 0 1
150 291650 520 7 58650 0 0
151 305959 1766 46 60438 0 0
152 315380 0 5 60720 0 0
153 297765 458 7 58625 0 1
154 315245 20 2 60938 0 0
155 315380 0 0 60720 0 0
156 315380 0 0 60720 0 0
157 315236 98 2 61490 0 0
158 336425 405 5 60845 0 0
159 315380 0 0 60720 0 0
160 315380 0 0 60720 0 0
161 315380 0 0 60720 0 1
162 315380 0 0 60720 0 0
163 306268 483 7 60830 0 0
164 302187 454 24 63261 0 1
165 314882 47 1 60720 0 0
166 315380 0 0 60720 0 1
167 382712 757 18 45689 0 1
168 341570 4655 55 60720 0 1
169 315380 0 0 60720 0 0
170 315380 0 0 60720 0 0
171 312412 36 3 61564 0 1
172 315380 0 0 60720 0 1
173 309596 203 9 61938 0 1
174 315380 0 0 60720 0 1
175 315547 126 8 60951 0 0
176 313267 400 113 60720 0 0
177 316176 71 0 60745 1 1
178 315380 0 0 60720 0 0
179 315380 0 0 60720 0 0
180 359335 972 19 71642 0 1
181 330068 531 11 71641 1 0
182 314289 2461 25 55792 0 0
183 297413 378 16 71873 1 0
184 314806 23 5 62555 1 1
185 333210 638 11 60370 1 1
186 352108 2300 23 64873 1 0
187 313332 149 6 62041 0 0
188 291787 226 5 65745 0 1
189 315380 0 0 60720 0 0
190 318745 275 7 59500 0 0
191 315380 0 0 60720 0 0
192 315366 141 7 61630 0 0
193 315380 0 0 60720 0 1
194 315688 28 3 60890 0 0
195 315380 0 0 60720 1 1
196 409642 4980 89 113521 1 1
197 315380 0 0 60720 1 1
198 315380 0 0 60720 1 1
199 269587 472 19 80045 1 1
200 315380 0 0 60720 1 1
201 315380 0 0 60720 1 1
202 315380 0 0 60720 1 1
203 300962 203 12 50804 1 0
204 325479 496 12 87390 1 1
205 316155 10 5 61656 1 1
206 318574 63 2 65688 1 1
207 315380 0 0 60720 1 1
208 343613 1136 26 48522 1 0
209 306948 265 3 60720 1 0
210 315380 0 0 60720 0 0
211 315380 0 0 60720 0 1
212 330059 267 11 57640 0 1
213 288985 474 10 61977 0 1
214 304485 534 5 62620 0 1
215 315380 0 2 60720 0 1
216 315688 15 6 60831 0 0
217 317736 397 7 60646 0 0
218 315380 0 2 60720 0 0
219 322331 1866 28 56225 0 0
220 296656 288 3 60510 0 1
221 315380 0 0 60720 0 1
222 315354 3 1 60698 0 1
223 312161 468 20 60720 0 0
224 315576 20 1 60805 0 0
225 314922 278 22 61404 0 0
226 314551 61 9 60720 0 0
227 315380 0 0 60720 0 1
228 312339 192 2 65276 0 1
229 315380 0 0 60720 0 0
230 298700 317 7 63915 0 1
231 321376 738 9 60720 0 1
232 315380 0 0 60720 0 1
233 303230 368 13 61686 0 0
234 315380 0 0 60720 0 0
235 315487 2 0 60743 0 0
236 315380 0 0 60720 0 0
237 315793 53 6 60349 0 0
238 315380 0 0 60720 0 1
239 315380 0 0 60720 0 1
240 315380 0 0 60720 0 0
241 312887 94 3 61360 0 1
242 315380 0 0 60720 0 0
243 315637 24 7 59818 0 1
244 324385 2332 2 72680 0 1
245 315380 0 0 60720 1 1
246 315380 0 0 60720 1 1
247 308989 131 15 61808 0 0
248 315380 0 0 60720 1 1
249 315380 0 0 60720 1 0
250 296702 206 9 53110 0 1
251 315380 0 0 60720 1 1
252 307322 167 1 64245 0 1
253 304376 622 38 73007 0 0
254 253588 2328 57 82732 0 0
255 315380 0 0 60720 0 0
256 309560 365 7 54820 0 0
257 298466 364 26 47705 0 0
258 315380 0 0 60720 1 1
259 315380 0 0 60720 0 1
260 315380 0 0 60720 1 1
261 315380 0 0 60720 1 1
262 343929 226 13 72835 0 0
263 331955 307 10 58856 0 1
264 315380 0 0 60720 1 1
265 315380 0 0 60720 0 1
266 315380 0 0 60720 1 0
267 381180 188 9 77655 0 1
268 315380 0 0 60720 1 1
269 331420 138 26 69817 0 1
270 315380 0 0 60720 1 1
271 315380 0 0 60720 1 1
272 315380 0 0 60720 1 1
273 310201 125 19 60798 0 0
274 315380 0 0 60720 0 0
275 320016 282 12 62452 0 1
276 320398 335 23 64175 0 0
277 315380 0 0 60720 0 1
278 291841 1324 29 67440 1 1
279 310670 176 8 68136 0 1
280 315380 0 0 60720 1 1
281 315380 0 0 60720 0 0
282 313491 249 26 56726 0 1
283 315380 0 0 60720 0 0
284 331323 333 9 70811 0 1
285 315380 0 0 60720 0 1
286 319210 601 5 60720 1 1
287 318098 30 3 62045 0 0
288 315380 0 0 60720 1 0
289 292754 249 13 54323 0 0
290 315380 0 0 60720 1 1
291 325176 165 12 62841 0 0
292 365959 453 19 81125 0 1
293 315380 0 0 60720 0 0
294 302409 53 10 59506 0 1
295 340968 382 9 59365 1 1
296 315380 0 0 60720 0 0
297 315380 0 0 60720 1 0
298 315380 0 0 60720 0 1
299 315380 0 9 60720 0 0
300 313164 30 4 60798 1 1
301 301164 290 1 58790 0 1
302 315380 0 1 60720 0 1
303 315380 0 0 60720 1 1
304 344425 366 14 61808 0 0
305 315394 2 12 60735 0 0
306 315380 0 0 60720 1 1
307 316647 209 19 64016 1 1
308 309836 384 17 54683 0 0
309 315380 0 0 60720 1 1
310 315380 0 0 60720 1 1
311 346611 365 32 87192 0 1
312 315380 0 0 60720 1 1
313 322031 49 14 64107 1 0
314 315656 3 8 60761 0 0
315 339445 133 4 65990 0 0
316 314964 32 0 59988 0 0
317 297141 368 20 61167 0 0
318 315372 1 5 60719 0 0
319 315380 0 0 60720 0 0
320 315380 0 0 60720 0 1
321 315380 0 0 60720 0 1
322 315380 0 0 60720 0 1
323 315380 0 0 60720 1 0
324 315380 0 0 60720 0 1
325 312502 22 1 60722 0 1
326 315380 0 0 60720 1 0
327 315380 0 0 60720 1 1
328 315380 0 0 60720 0 1
329 315380 0 0 60720 0 1
330 315380 0 0 60720 0 1
331 315380 0 0 60720 0 0
332 315380 0 0 60720 0 1
333 313729 96 4 60379 0 1
334 315388 1 1 60727 0 0
335 315371 314 4 60720 0 0
336 296139 844 20 60925 0 1
337 315380 0 0 60720 0 0
338 313880 26 1 60896 0 0
339 317698 125 10 59734 0 1
340 295580 304 12 62969 0 0
341 315380 0 0 60720 0 0
342 315380 0 0 60720 0 0
343 315380 0 0 60720 0 1
344 308256 621 13 60720 0 0
345 315380 0 0 60720 0 0
346 303677 119 3 59118 0 0
347 315380 0 0 60720 0 0
348 315380 0 0 60720 0 1
349 319369 1595 10 60720 0 0
350 318690 312 3 58598 0 1
351 314049 60 7 61124 0 0
352 325699 587 10 59595 0 1
353 314210 135 1 62065 0 0
354 315380 0 0 60720 0 0
355 315380 0 0 60720 0 0
356 322378 514 15 78780 0 0
357 315380 0 0 60720 0 0
358 315380 0 0 60720 0 0
359 315380 0 0 60720 0 1
360 315398 1 4 60722 0 0
361 315380 0 0 60720 0 1
362 315380 0 0 60720 0 0
363 308336 1763 28 61600 1 1
364 316386 180 9 59635 0 1
365 315380 0 0 60720 0 1
366 315380 0 0 60720 0 1
367 315380 0 0 60720 0 0
368 315380 0 0 60720 0 1
369 315553 218 7 60720 0 0
370 315380 0 0 60720 0 1
371 323361 448 7 59781 0 1
372 336639 227 7 76644 0 0
373 307424 174 3 64820 0 1
374 315380 0 0 60720 0 1
375 315380 0 0 60720 0 1
376 295370 121 11 56178 0 1
377 322340 607 7 60436 0 0
378 319864 2212 10 60720 0 0
379 315380 0 0 60720 0 1
380 315380 0 0 60720 0 1
381 317291 530 18 73433 0 0
382 280398 571 14 41477 0 0
383 315380 0 0 60720 0 1
384 317330 78 12 62700 0 0
385 238125 2489 29 67804 0 1
386 327071 131 3 59661 0 0
387 309038 923 6 58620 0 1
388 314210 72 3 60398 0 0
389 307930 572 8 58580 0 1
390 322327 397 10 62710 0 1
391 292136 450 6 59325 0 0
392 263276 622 8 60950 0 1
393 367655 694 6 68060 0 0
394 283910 3425 9 83620 1 0
395 283587 562 8 58456 0 1
396 243650 4917 26 52811 0 1
397 438493 1442 239 121173 1 1
398 296261 529 7 63870 0 1
399 230621 2126 41 21001 1 1
400 304252 1061 3 70415 0 0
401 333505 776 8 64230 0 0
402 296919 611 6 59190 0 0
403 278990 1526 21 69351 1 0
404 276898 592 7 64270 0 0
405 327007 1182 11 70694 0 0
406 317046 621 11 68005 0 1
407 304555 989 12 58930 0 1
408 298096 438 9 58320 0 0
409 231861 726 3 69980 0 1
410 309422 1303 57 69863 0 1
411 286963 7419 21 63255 1 1
412 269753 1164 15 57320 1 1
413 448243 3310 32 75230 1 1
414 165404 1920 11 79420 0 1
415 204325 965 2 73490 0 0
416 407159 3256 23 35250 0 1
417 290476 1135 20 62285 1 0
418 275311 1270 24 69206 0 0
419 246541 661 1 65920 0 0
420 253468 1013 1 69770 0 0
421 240897 2844 74 72683 0 1
422 -83265 11528 68 -14545 1 1
423 -42143 6526 20 55830 0 1
424 272713 2264 20 55174 0 1
425 215362 5109 82 67038 1 1
426 42754 3999 21 51252 0 1
427 306275 35624 244 157278 0 1
428 253537 9252 32 79510 0 1
429 372631 15236 86 77440 0 1
430 -7170 18073 69 27284 0 1
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) Costs Orders Dividends Group `Gender\r`
197015.803 12.618 1501.735 1.356 66885.361 -42165.685
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-1608369 -36822 20648 51704 2945709
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.970e+05 3.124e+04 6.307 7.14e-10 ***
Costs 1.262e+01 2.322e+00 5.434 9.32e-08 ***
Orders 1.502e+03 3.337e+02 4.500 8.78e-06 ***
Dividends 1.356e+00 3.839e-01 3.533 0.000456 ***
Group 6.689e+04 2.674e+04 2.501 0.012745 *
`Gender\r` -4.217e+04 2.403e+04 -1.755 0.080049 .
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 235600 on 424 degrees of freedom
Multiple R-squared: 0.4759, Adjusted R-squared: 0.4697
F-statistic: 77.01 on 5 and 424 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 3.119064e-53 1.559532e-53
[2,] 1.0000000 9.391202e-54 4.695601e-54
[3,] 1.0000000 2.804233e-65 1.402117e-65
[4,] 1.0000000 2.178794e-89 1.089397e-89
[5,] 1.0000000 4.181877e-91 2.090939e-91
[6,] 1.0000000 2.087251e-106 1.043626e-106
[7,] 1.0000000 2.655190e-117 1.327595e-117
[8,] 1.0000000 5.038473e-122 2.519236e-122
[9,] 1.0000000 2.112852e-132 1.056426e-132
[10,] 1.0000000 2.896492e-137 1.448246e-137
[11,] 1.0000000 1.732373e-140 8.661863e-141
[12,] 1.0000000 2.895087e-145 1.447543e-145
[13,] 1.0000000 1.463764e-146 7.318819e-147
[14,] 1.0000000 3.203861e-147 1.601930e-147
[15,] 1.0000000 1.626336e-149 8.131679e-150
[16,] 1.0000000 5.983008e-149 2.991504e-149
[17,] 1.0000000 1.149654e-152 5.748272e-153
[18,] 1.0000000 2.928685e-158 1.464343e-158
[19,] 1.0000000 6.008239e-158 3.004119e-158
[20,] 1.0000000 1.739595e-172 8.697973e-173
[21,] 1.0000000 1.161595e-172 5.807976e-173
[22,] 1.0000000 1.261577e-176 6.307883e-177
[23,] 1.0000000 1.086674e-178 5.433368e-179
[24,] 1.0000000 3.225017e-186 1.612509e-186
[25,] 1.0000000 2.193136e-188 1.096568e-188
[26,] 1.0000000 2.002572e-187 1.001286e-187
[27,] 1.0000000 3.010327e-187 1.505164e-187
[28,] 1.0000000 1.676785e-193 8.383923e-194
[29,] 1.0000000 2.675121e-199 1.337560e-199
[30,] 1.0000000 4.075137e-200 2.037568e-200
[31,] 1.0000000 2.007825e-200 1.003913e-200
[32,] 1.0000000 2.315420e-201 1.157710e-201
[33,] 1.0000000 1.145022e-200 5.725110e-201
[34,] 1.0000000 3.235417e-201 1.617709e-201
[35,] 1.0000000 3.952575e-200 1.976288e-200
[36,] 1.0000000 5.801512e-200 2.900756e-200
[37,] 1.0000000 2.549246e-199 1.274623e-199
[38,] 1.0000000 1.921825e-201 9.609124e-202
[39,] 1.0000000 3.638143e-202 1.819071e-202
[40,] 1.0000000 6.835358e-202 3.417679e-202
[41,] 1.0000000 3.434304e-202 1.717152e-202
[42,] 1.0000000 8.509666e-203 4.254833e-203
[43,] 1.0000000 4.624594e-203 2.312297e-203
[44,] 1.0000000 5.660809e-203 2.830405e-203
[45,] 1.0000000 9.425505e-203 4.712752e-203
[46,] 1.0000000 5.852011e-207 2.926006e-207
[47,] 1.0000000 1.494129e-206 7.470643e-207
[48,] 1.0000000 1.149188e-220 5.745941e-221
[49,] 1.0000000 5.270960e-221 2.635480e-221
[50,] 1.0000000 3.863063e-220 1.931532e-220
[51,] 1.0000000 4.912123e-220 2.456062e-220
[52,] 1.0000000 1.562514e-221 7.812569e-222
[53,] 1.0000000 4.198600e-223 2.099300e-223
[54,] 1.0000000 1.543670e-222 7.718348e-223
[55,] 1.0000000 2.851558e-222 1.425779e-222
[56,] 1.0000000 1.712246e-221 8.561229e-222
[57,] 1.0000000 3.167759e-221 1.583880e-221
[58,] 1.0000000 6.530479e-221 3.265240e-221
[59,] 1.0000000 6.028345e-220 3.014173e-220
[60,] 1.0000000 3.164635e-219 1.582317e-219
[61,] 1.0000000 1.180140e-218 5.900698e-219
[62,] 1.0000000 7.196561e-218 3.598280e-218
[63,] 1.0000000 6.924825e-219 3.462413e-219
[64,] 1.0000000 2.197226e-218 1.098613e-218
[65,] 1.0000000 1.189270e-217 5.946349e-218
[66,] 1.0000000 1.008801e-216 5.044007e-217
[67,] 1.0000000 2.120505e-216 1.060252e-216
[68,] 1.0000000 1.795763e-215 8.978817e-216
[69,] 1.0000000 2.357403e-215 1.178702e-215
[70,] 1.0000000 1.985358e-214 9.926792e-215
[71,] 1.0000000 7.346452e-214 3.673226e-214
[72,] 1.0000000 3.003100e-213 1.501550e-213
[73,] 1.0000000 3.503251e-213 1.751626e-213
[74,] 1.0000000 3.078370e-213 1.539185e-213
[75,] 1.0000000 1.976286e-212 9.881428e-213
[76,] 1.0000000 1.354311e-211 6.771557e-212
[77,] 1.0000000 2.780740e-211 1.390370e-211
[78,] 1.0000000 1.956036e-211 9.780180e-212
[79,] 1.0000000 1.403677e-210 7.018385e-211
[80,] 1.0000000 2.735503e-210 1.367752e-210
[81,] 1.0000000 2.592666e-209 1.296333e-209
[82,] 1.0000000 2.335285e-208 1.167643e-208
[83,] 1.0000000 1.449452e-207 7.247261e-208
[84,] 1.0000000 3.665838e-209 1.832919e-209
[85,] 1.0000000 2.861349e-208 1.430674e-208
[86,] 1.0000000 7.671012e-209 3.835506e-209
[87,] 1.0000000 8.183865e-209 4.091932e-209
[88,] 1.0000000 5.041356e-208 2.520678e-208
[89,] 1.0000000 5.299331e-207 2.649666e-207
[90,] 1.0000000 4.175145e-206 2.087573e-206
[91,] 1.0000000 4.538039e-205 2.269019e-205
[92,] 1.0000000 3.976922e-204 1.988461e-204
[93,] 1.0000000 3.148935e-205 1.574467e-205
[94,] 1.0000000 3.000841e-205 1.500421e-205
[95,] 1.0000000 4.306849e-206 2.153424e-206
[96,] 1.0000000 1.524646e-205 7.623232e-206
[97,] 1.0000000 4.739711e-205 2.369855e-205
[98,] 1.0000000 2.011683e-207 1.005841e-207
[99,] 1.0000000 1.005795e-208 5.028975e-209
[100,] 1.0000000 7.115215e-208 3.557608e-208
[101,] 1.0000000 1.389550e-207 6.947748e-208
[102,] 1.0000000 7.618880e-208 3.809440e-208
[103,] 1.0000000 3.668680e-208 1.834340e-208
[104,] 1.0000000 2.631100e-207 1.315550e-207
[105,] 1.0000000 2.076715e-206 1.038358e-206
[106,] 1.0000000 1.737609e-206 8.688046e-207
[107,] 1.0000000 7.196635e-206 3.598318e-206
[108,] 1.0000000 4.943050e-205 2.471525e-205
[109,] 1.0000000 4.028912e-205 2.014456e-205
[110,] 1.0000000 3.825979e-206 1.912989e-206
[111,] 1.0000000 3.229987e-205 1.614994e-205
[112,] 1.0000000 1.984586e-204 9.922928e-205
[113,] 1.0000000 1.832957e-203 9.164785e-204
[114,] 1.0000000 1.575958e-202 7.879789e-203
[115,] 1.0000000 2.464870e-202 1.232435e-202
[116,] 1.0000000 7.421702e-202 3.710851e-202
[117,] 1.0000000 1.415200e-202 7.076000e-203
[118,] 1.0000000 1.322493e-201 6.612466e-202
[119,] 1.0000000 1.308347e-200 6.541737e-201
[120,] 1.0000000 1.258909e-199 6.294546e-200
[121,] 1.0000000 9.268900e-199 4.634450e-199
[122,] 1.0000000 9.083554e-198 4.541777e-198
[123,] 1.0000000 5.931070e-198 2.965535e-198
[124,] 1.0000000 4.582489e-197 2.291245e-197
[125,] 1.0000000 4.730272e-196 2.365136e-196
[126,] 1.0000000 4.648341e-195 2.324170e-195
[127,] 1.0000000 4.509954e-194 2.254977e-194
[128,] 1.0000000 3.629916e-193 1.814958e-193
[129,] 1.0000000 3.501883e-192 1.750941e-192
[130,] 1.0000000 3.325466e-191 1.662733e-191
[131,] 1.0000000 3.188987e-190 1.594494e-190
[132,] 1.0000000 2.673568e-189 1.336784e-189
[133,] 1.0000000 2.085916e-188 1.042958e-188
[134,] 1.0000000 2.094039e-187 1.047019e-187
[135,] 1.0000000 2.019585e-186 1.009792e-186
[136,] 1.0000000 1.540432e-185 7.702161e-186
[137,] 1.0000000 1.415270e-184 7.076352e-185
[138,] 1.0000000 1.114503e-183 5.572515e-184
[139,] 1.0000000 1.075315e-182 5.376575e-183
[140,] 1.0000000 7.617526e-182 3.808763e-182
[141,] 1.0000000 6.075973e-181 3.037986e-181
[142,] 1.0000000 5.731942e-180 2.865971e-180
[143,] 1.0000000 5.553901e-179 2.776951e-179
[144,] 1.0000000 5.228394e-178 2.614197e-178
[145,] 1.0000000 4.500132e-177 2.250066e-177
[146,] 1.0000000 4.186493e-176 2.093246e-176
[147,] 1.0000000 3.857105e-175 1.928553e-175
[148,] 1.0000000 3.543484e-174 1.771742e-174
[149,] 1.0000000 3.276616e-173 1.638308e-173
[150,] 1.0000000 2.293475e-172 1.146737e-172
[151,] 1.0000000 2.088166e-171 1.044083e-171
[152,] 1.0000000 1.894897e-170 9.474486e-171
[153,] 1.0000000 1.476460e-169 7.382300e-170
[154,] 1.0000000 1.330815e-168 6.654075e-169
[155,] 1.0000000 1.229603e-167 6.148016e-168
[156,] 1.0000000 1.045383e-166 5.226914e-167
[157,] 1.0000000 9.340998e-166 4.670499e-166
[158,] 1.0000000 7.230059e-165 3.615029e-165
[159,] 1.0000000 6.864815e-166 3.432408e-166
[160,] 1.0000000 1.050247e-165 5.251235e-166
[161,] 1.0000000 9.443733e-165 4.721866e-165
[162,] 1.0000000 8.459355e-164 4.229678e-164
[163,] 1.0000000 6.882020e-163 3.441010e-163
[164,] 1.0000000 5.403056e-162 2.701528e-162
[165,] 1.0000000 4.533571e-161 2.266785e-161
[166,] 1.0000000 3.562971e-160 1.781486e-160
[167,] 1.0000000 3.161754e-159 1.580877e-159
[168,] 1.0000000 2.490193e-158 1.245096e-158
[169,] 1.0000000 2.216461e-157 1.108231e-157
[170,] 1.0000000 1.918934e-156 9.594671e-157
[171,] 1.0000000 1.654776e-155 8.273881e-156
[172,] 1.0000000 9.287793e-155 4.643897e-155
[173,] 1.0000000 7.685498e-154 3.842749e-154
[174,] 1.0000000 4.331669e-153 2.165834e-153
[175,] 1.0000000 1.951747e-152 9.758736e-153
[176,] 1.0000000 1.707014e-151 8.535070e-152
[177,] 1.0000000 1.163927e-150 5.819633e-151
[178,] 1.0000000 4.752783e-150 2.376391e-150
[179,] 1.0000000 4.055371e-149 2.027686e-149
[180,] 1.0000000 2.714794e-148 1.357397e-148
[181,] 1.0000000 2.255632e-147 1.127816e-147
[182,] 1.0000000 1.785160e-146 8.925800e-147
[183,] 1.0000000 1.471050e-145 7.355251e-146
[184,] 1.0000000 1.223294e-144 6.116472e-145
[185,] 1.0000000 9.133870e-144 4.566935e-144
[186,] 1.0000000 7.461621e-143 3.730811e-143
[187,] 1.0000000 6.190266e-142 3.095133e-142
[188,] 1.0000000 4.034319e-141 2.017159e-141
[189,] 1.0000000 3.314928e-140 1.657464e-140
[190,] 1.0000000 2.710295e-139 1.355147e-139
[191,] 1.0000000 2.282429e-139 1.141215e-139
[192,] 1.0000000 1.874873e-138 9.374367e-139
[193,] 1.0000000 1.532662e-137 7.663309e-138
[194,] 1.0000000 1.246856e-136 6.234280e-137
[195,] 1.0000000 9.342345e-136 4.671172e-136
[196,] 1.0000000 3.479348e-135 1.739674e-135
[197,] 1.0000000 2.812736e-134 1.406368e-134
[198,] 1.0000000 2.230367e-133 1.115183e-133
[199,] 1.0000000 1.781902e-132 8.909511e-133
[200,] 1.0000000 3.024369e-132 1.512184e-132
[201,] 1.0000000 2.350495e-131 1.175248e-131
[202,] 1.0000000 1.811408e-130 9.057038e-131
[203,] 1.0000000 1.294973e-129 6.474867e-130
[204,] 1.0000000 6.814186e-129 3.407093e-129
[205,] 1.0000000 4.776280e-128 2.388140e-128
[206,] 1.0000000 3.562670e-127 1.781335e-127
[207,] 1.0000000 2.540584e-126 1.270292e-126
[208,] 1.0000000 1.923010e-125 9.615050e-126
[209,] 1.0000000 1.403194e-124 7.015970e-125
[210,] 1.0000000 1.049727e-123 5.248633e-124
[211,] 1.0000000 4.658856e-123 2.329428e-123
[212,] 1.0000000 3.354149e-122 1.677075e-122
[213,] 1.0000000 2.341114e-121 1.170557e-121
[214,] 1.0000000 1.632779e-120 8.163896e-121
[215,] 1.0000000 1.202557e-119 6.012784e-120
[216,] 1.0000000 8.760590e-119 4.380295e-119
[217,] 1.0000000 6.409275e-118 3.204637e-118
[218,] 1.0000000 4.658212e-117 2.329106e-117
[219,] 1.0000000 3.182352e-116 1.591176e-116
[220,] 1.0000000 2.222739e-115 1.111370e-115
[221,] 1.0000000 1.583412e-114 7.917060e-115
[222,] 1.0000000 1.066319e-113 5.331593e-114
[223,] 1.0000000 6.614849e-113 3.307424e-113
[224,] 1.0000000 4.453289e-112 2.226645e-112
[225,] 1.0000000 3.197515e-111 1.598758e-111
[226,] 1.0000000 2.230605e-110 1.115302e-110
[227,] 1.0000000 1.548521e-109 7.742604e-110
[228,] 1.0000000 1.070353e-108 5.351767e-109
[229,] 1.0000000 7.320488e-108 3.660244e-108
[230,] 1.0000000 4.810673e-107 2.405337e-107
[231,] 1.0000000 3.151298e-106 1.575649e-106
[232,] 1.0000000 2.139560e-105 1.069780e-105
[233,] 1.0000000 1.418807e-104 7.094035e-105
[234,] 1.0000000 9.545832e-104 4.772916e-104
[235,] 1.0000000 6.111366e-103 3.055683e-103
[236,] 1.0000000 4.070336e-102 2.035168e-102
[237,] 1.0000000 2.769436e-101 1.384718e-101
[238,] 1.0000000 1.875004e-100 9.375022e-101
[239,] 1.0000000 1.264913e-99 6.324564e-100
[240,] 1.0000000 8.480102e-99 4.240051e-99
[241,] 1.0000000 5.616045e-98 2.808022e-98
[242,] 1.0000000 3.508506e-97 1.754253e-97
[243,] 1.0000000 2.318360e-96 1.159180e-96
[244,] 1.0000000 1.467934e-95 7.339668e-96
[245,] 1.0000000 7.787919e-95 3.893960e-95
[246,] 1.0000000 5.870566e-95 2.935283e-95
[247,] 1.0000000 3.786195e-94 1.893098e-94
[248,] 1.0000000 2.229166e-93 1.114583e-93
[249,] 1.0000000 1.181037e-92 5.905187e-93
[250,] 1.0000000 7.659204e-92 3.829602e-92
[251,] 1.0000000 4.678480e-91 2.339240e-91
[252,] 1.0000000 3.004381e-90 1.502190e-90
[253,] 1.0000000 1.919313e-89 9.596566e-90
[254,] 1.0000000 1.211357e-88 6.056786e-89
[255,] 1.0000000 5.657238e-88 2.828619e-88
[256,] 1.0000000 3.566104e-87 1.783052e-87
[257,] 1.0000000 2.120619e-86 1.060310e-86
[258,] 1.0000000 1.319610e-85 6.598048e-86
[259,] 1.0000000 6.234355e-85 3.117177e-85
[260,] 1.0000000 3.856098e-84 1.928049e-84
[261,] 1.0000000 2.356679e-83 1.178340e-83
[262,] 1.0000000 1.442532e-82 7.212658e-83
[263,] 1.0000000 8.781296e-82 4.390648e-82
[264,] 1.0000000 5.315850e-81 2.657925e-81
[265,] 1.0000000 3.213512e-80 1.606756e-80
[266,] 1.0000000 1.896571e-79 9.482857e-80
[267,] 1.0000000 1.076193e-78 5.380964e-79
[268,] 1.0000000 6.368686e-78 3.184343e-78
[269,] 1.0000000 3.581063e-77 1.790531e-77
[270,] 1.0000000 1.803531e-76 9.017655e-77
[271,] 1.0000000 1.014375e-75 5.071875e-76
[272,] 1.0000000 5.917060e-75 2.958530e-75
[273,] 1.0000000 3.376103e-74 1.688052e-74
[274,] 1.0000000 1.753947e-73 8.769736e-74
[275,] 1.0000000 9.910888e-73 4.955444e-73
[276,] 1.0000000 5.585390e-72 2.792695e-72
[277,] 1.0000000 3.030375e-71 1.515188e-71
[278,] 1.0000000 1.652394e-70 8.261968e-71
[279,] 1.0000000 9.183571e-70 4.591786e-70
[280,] 1.0000000 5.137785e-69 2.568892e-69
[281,] 1.0000000 2.856848e-68 1.428424e-68
[282,] 1.0000000 1.580626e-67 7.903132e-68
[283,] 1.0000000 8.394063e-67 4.197031e-67
[284,] 1.0000000 4.450574e-66 2.225287e-66
[285,] 1.0000000 2.387072e-65 1.193536e-65
[286,] 1.0000000 1.285730e-64 6.428651e-65
[287,] 1.0000000 5.203901e-64 2.601950e-64
[288,] 1.0000000 2.753913e-63 1.376956e-63
[289,] 1.0000000 1.467843e-62 7.339214e-63
[290,] 1.0000000 7.467417e-62 3.733709e-62
[291,] 1.0000000 3.901597e-61 1.950799e-61
[292,] 1.0000000 2.049672e-60 1.024836e-60
[293,] 1.0000000 1.056153e-59 5.280763e-60
[294,] 1.0000000 5.263601e-59 2.631801e-59
[295,] 1.0000000 2.700465e-58 1.350232e-58
[296,] 1.0000000 1.068365e-57 5.341826e-58
[297,] 1.0000000 5.405417e-57 2.702708e-57
[298,] 1.0000000 2.732371e-56 1.366186e-56
[299,] 1.0000000 1.385621e-55 6.928106e-56
[300,] 1.0000000 6.338311e-55 3.169156e-55
[301,] 1.0000000 3.147763e-54 1.573881e-54
[302,] 1.0000000 1.550330e-53 7.751648e-54
[303,] 1.0000000 6.913106e-53 3.456553e-53
[304,] 1.0000000 3.362089e-52 1.681044e-52
[305,] 1.0000000 1.631869e-51 8.159347e-52
[306,] 1.0000000 7.834841e-51 3.917420e-51
[307,] 1.0000000 3.519905e-50 1.759952e-50
[308,] 1.0000000 1.652140e-49 8.260698e-50
[309,] 1.0000000 7.913803e-49 3.956902e-49
[310,] 1.0000000 3.703998e-48 1.851999e-48
[311,] 1.0000000 1.718858e-47 8.594288e-48
[312,] 1.0000000 7.789872e-47 3.894936e-47
[313,] 1.0000000 3.508059e-46 1.754030e-46
[314,] 1.0000000 1.569654e-45 7.848270e-46
[315,] 1.0000000 7.154481e-45 3.577240e-45
[316,] 1.0000000 3.162775e-44 1.581388e-44
[317,] 1.0000000 1.405488e-43 7.027441e-44
[318,] 1.0000000 6.272830e-43 3.136415e-43
[319,] 1.0000000 2.748366e-42 1.374183e-42
[320,] 1.0000000 1.183991e-41 5.919953e-42
[321,] 1.0000000 5.064765e-41 2.532382e-41
[322,] 1.0000000 2.151038e-40 1.075519e-40
[323,] 1.0000000 9.281042e-40 4.640521e-40
[324,] 1.0000000 3.889104e-39 1.944552e-39
[325,] 1.0000000 1.622848e-38 8.114239e-39
[326,] 1.0000000 6.870085e-38 3.435043e-38
[327,] 1.0000000 2.862942e-37 1.431471e-37
[328,] 1.0000000 1.222146e-36 6.110729e-37
[329,] 1.0000000 5.059980e-36 2.529990e-36
[330,] 1.0000000 2.092263e-35 1.046132e-35
[331,] 1.0000000 8.173019e-35 4.086509e-35
[332,] 1.0000000 3.352694e-34 1.676347e-34
[333,] 1.0000000 1.350493e-33 6.752466e-34
[334,] 1.0000000 5.397024e-33 2.698512e-33
[335,] 1.0000000 2.097327e-32 1.048664e-32
[336,] 1.0000000 8.345836e-32 4.172918e-32
[337,] 1.0000000 3.259795e-31 1.629897e-31
[338,] 1.0000000 1.287244e-30 6.436219e-31
[339,] 1.0000000 4.945802e-30 2.472901e-30
[340,] 1.0000000 1.854176e-29 9.270879e-30
[341,] 1.0000000 6.385035e-29 3.192517e-29
[342,] 1.0000000 2.217132e-28 1.108566e-28
[343,] 1.0000000 8.333414e-28 4.166707e-28
[344,] 1.0000000 2.719755e-27 1.359877e-27
[345,] 1.0000000 1.010475e-26 5.052376e-27
[346,] 1.0000000 3.672330e-26 1.836165e-26
[347,] 1.0000000 1.321805e-25 6.609023e-26
[348,] 1.0000000 4.695365e-25 2.347683e-25
[349,] 1.0000000 1.665203e-24 8.326014e-25
[350,] 1.0000000 5.845746e-24 2.922873e-24
[351,] 1.0000000 2.012322e-23 1.006161e-23
[352,] 1.0000000 6.948032e-23 3.474016e-23
[353,] 1.0000000 2.349383e-22 1.174691e-22
[354,] 1.0000000 7.948465e-22 3.974233e-22
[355,] 1.0000000 2.632301e-21 1.316150e-21
[356,] 1.0000000 8.492560e-21 4.246280e-21
[357,] 1.0000000 2.762021e-20 1.381010e-20
[358,] 1.0000000 8.880549e-20 4.440274e-20
[359,] 1.0000000 2.880110e-19 1.440055e-19
[360,] 1.0000000 9.075136e-19 4.537568e-19
[361,] 1.0000000 2.871064e-18 1.435532e-18
[362,] 1.0000000 8.862034e-18 4.431017e-18
[363,] 1.0000000 2.504471e-17 1.252235e-17
[364,] 1.0000000 8.090228e-17 4.045114e-17
[365,] 1.0000000 2.556178e-16 1.278089e-16
[366,] 1.0000000 7.517475e-16 3.758738e-16
[367,] 1.0000000 2.177930e-15 1.088965e-15
[368,] 1.0000000 6.462148e-15 3.231074e-15
[369,] 1.0000000 1.804762e-14 9.023808e-15
[370,] 1.0000000 4.659063e-14 2.329531e-14
[371,] 1.0000000 1.291315e-13 6.456577e-14
[372,] 1.0000000 3.515439e-13 1.757719e-13
[373,] 1.0000000 1.041513e-12 5.207565e-13
[374,] 1.0000000 2.588360e-12 1.294180e-12
[375,] 1.0000000 6.816873e-12 3.408437e-12
[376,] 1.0000000 1.859086e-11 9.295429e-12
[377,] 1.0000000 4.744356e-11 2.372178e-11
[378,] 1.0000000 1.142402e-10 5.712012e-11
[379,] 1.0000000 2.817406e-10 1.408703e-10
[380,] 1.0000000 7.112249e-10 3.556124e-10
[381,] 1.0000000 1.732707e-09 8.663534e-10
[382,] 1.0000000 4.091044e-09 2.045522e-09
[383,] 1.0000000 1.051183e-08 5.255917e-09
[384,] 1.0000000 2.757690e-08 1.378845e-08
[385,] 1.0000000 5.064842e-08 2.532421e-08
[386,] 0.9999999 1.137411e-07 5.687056e-08
[387,] 0.9999999 2.775260e-07 1.387630e-07
[388,] 0.9999997 6.440629e-07 3.220315e-07
[389,] 0.9999994 1.206373e-06 6.031866e-07
[390,] 0.9999986 2.823280e-06 1.411640e-06
[391,] 0.9999971 5.876508e-06 2.938254e-06
[392,] 0.9999931 1.373613e-05 6.868066e-06
[393,] 0.9999865 2.694751e-05 1.347375e-05
[394,] 0.9999719 5.613207e-05 2.806604e-05
[395,] 0.9999393 1.214554e-04 6.072768e-05
[396,] 0.9998668 2.664745e-04 1.332373e-04
[397,] 0.9997393 5.214706e-04 2.607353e-04
[398,] 0.9994818 1.036438e-03 5.182188e-04
[399,] 0.9990712 1.857575e-03 9.287875e-04
[400,] 0.9983556 3.288838e-03 1.644419e-03
[401,] 0.9967621 6.475853e-03 3.237926e-03
[402,] 0.9938543 1.229134e-02 6.145668e-03
[403,] 0.9887839 2.243216e-02 1.121608e-02
[404,] 0.9793192 4.136155e-02 2.068078e-02
[405,] 0.9810969 3.780613e-02 1.890307e-02
[406,] 0.9692122 6.157567e-02 3.078784e-02
[407,] 0.9489260 1.021480e-01 5.107402e-02
[408,] 0.9843747 3.125068e-02 1.562534e-02
[409,] 0.9668674 6.626526e-02 3.313263e-02
[410,] 0.9315350 1.369300e-01 6.846500e-02
[411,] 0.8669904 2.660191e-01 1.330096e-01
[412,] 0.7595129 4.809743e-01 2.404871e-01
[413,] 0.6009843 7.980313e-01 3.990157e-01
> postscript(file="/var/www/html/rcomp/tmp/13jx21292932180.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> plot(x[,1], type='l', main='Actuals and Interpolation', ylab='value of Actuals and Interpolation (dots)', xlab='time or index')
> points(x[,1]-mysum$resid)
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/23jx21292932180.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> plot(mysum$resid, type='b', pch=19, main='Residuals', ylab='value of Residuals', xlab='time or index')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/33jx21292932180.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> hist(mysum$resid, main='Residual Histogram', xlab='values of Residuals')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/4esxn1292932180.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> densityplot(~mysum$resid,col='black',main='Residual Density Plot', xlab='values of Residuals')
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/5esxn1292932180.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 = 430
Frequency = 1
1 2 3 4 5
2.945709e+06 1.958932e+06 -1.608369e+06 -6.533468e+04 3.622083e+05
6 7 8 9 10
1.569258e+05 3.729608e+05 -1.382227e+05 7.579032e+05 1.030653e+05
11 12 13 14 15
4.217010e+05 3.683217e+05 -8.069330e+05 -1.555096e+05 4.914049e+05
16 17 18 19 20
-1.921223e+05 -1.947991e+05 -1.065107e+05 -5.178519e+05 2.502813e+05
21 22 23 24 25
-4.255778e+04 -4.201810e+04 -5.939069e+05 -2.534845e+05 6.050526e+04
26 27 28 29 30
-3.689886e+05 -2.475860e+05 4.290343e+05 -1.733252e+05 2.672605e+05
31 32 33 34 35
-4.558695e+04 2.200257e+05 1.174329e+05 3.411520e+04 -4.667818e+05
36 37 38 39 40
1.238381e+05 1.356866e+05 -8.574527e+04 -3.037759e+05 -4.832529e+05
41 42 43 44 45
-1.970856e+05 -1.326601e+04 -1.952668e+05 -2.031444e+04 -6.659535e+04
46 47 48 49 50
-2.718181e+05 -2.949503e+05 -6.060906e+01 -1.543829e+05 -1.271625e+04
51 52 53 54 55
-1.427621e+05 -2.242160e+05 -4.412895e+04 -4.303624e+05 -1.997302e+04
56 57 58 59 60
3.734755e+05 -8.144209e+03 -4.126198e+04 -1.979514e+05 -3.327342e+04
61 62 63 64 65
-1.996854e+05 -1.262337e+04 -2.224404e+05 -3.673721e+04 -4.644577e+05
66 67 68 69 70
-3.685077e+04 -5.278891e+04 -8.535061e+04 -1.114192e+05 -1.687591e+05
71 72 73 74 75
6.897109e+03 -1.018463e+05 4.283796e+04 2.427217e+03 -1.322708e+05
76 77 78 79 80
-4.093145e+04 -8.739013e+04 -3.093746e+04 -6.970571e+04 4.885414e+04
81 82 83 84 85
-1.180068e+05 -5.111407e+04 1.831975e+04 5.431579e+04 -1.254997e+05
86 87 88 89 90
-1.506956e+04 -1.013768e+04 -1.270529e+05 2.083809e+04 -1.259173e+04
91 92 93 94 95
-5.295845e+04 8.081311e+04 -1.550252e+05 6.448599e+04 -5.258125e+04
96 97 98 99 100
1.040454e+05 -3.726888e+04 -3.086653e+04 -3.412860e+04 -7.991547e+04
101 102 103 104 105
-8.335505e+04 -1.278187e+05 -1.436685e+04 5.683226e+04 -3.521946e+04
106 107 108 109 110
-2.898622e+05 1.064864e+05 -1.553016e+05 -1.252069e+05 1.170282e+05
111 112 113 114 115
-2.213077e+05 -1.524851e+05 -8.239323e+04 5.174149e+04 -1.529056e+03
116 117 118 119 120
3.358610e+03 7.191914e+04 -1.569444e+05 3.564294e+04 -2.009287e+05
121 122 123 124 125
1.732224e+04 -1.134618e+05 5.424982e+04 -1.295347e+04 -2.256727e+05
126 127 128 129 130
3.600235e+04 1.885287e+04 2.591166e+03 7.816804e+04 2.415626e+04
131 132 133 134 135
8.192810e+04 7.666630e+04 2.664539e+03 1.855164e+04 2.592867e+04
136 137 138 139 140
6.491739e+04 3.394453e+04 1.820984e+04 3.600235e+04 3.155258e+04
141 142 143 144 145
7.816804e+04 -4.587239e+03 2.975914e+04 6.067965e+04 2.176231e+04
146 147 148 149 150
7.816804e+04 2.482765e+04 8.161402e+04 7.808688e+04 -1.993497e+03
151 152 153 154 155
-6.439980e+04 2.849368e+04 4.710343e+04 3.231582e+04 3.600235e+04
156 157 158 159 160
3.600235e+04 3.057385e+04 4.425873e+04 3.600235e+04 3.600235e+04
161 162 163 164 165
7.816804e+04 3.600235e+04 1.013438e+04 1.975804e+04 3.340956e+04
166 167 168 169 170
7.816804e+04 1.293051e+05 -3.697539e+04 3.600235e+04 3.600235e+04
171 172 173 174 175
6.909576e+04 7.816804e+04 5.465479e+04 7.816804e+04 2.225224e+04
176 177 178 179 180
-1.408540e+05 1.114887e+04 3.600235e+04 3.600235e+04 6.651030e+04
181 182 183 184 185
-5.422786e+04 -2.700112e+04 -9.277563e+04 4.207503e+02 5.017889e+03
186 187 188 189 190
-6.335012e+04 2.127199e+04 3.739862e+04 3.600235e+04 2.704002e+04
191 192 193 194 195
3.600235e+04 2.246269e+04 7.816804e+04 3.122124e+04 1.128268e+04
196 197 198 199 200
-1.625690e+05 1.128268e+04 1.128268e+04 -9.521193e+04 1.128268e+04
201 202 203 204 205
1.128268e+04 1.128268e+04 -5.243307e+04 -3.907353e+04 3.153210e+03
206 207 208 209 210
3.939561e+03 1.128268e+04 -3.948385e+04 -4.716405e+04 3.600235e+04
211 212 213 214 215
7.816804e+04 7.713665e+04 2.906961e+04 5.044902e+04 7.516457e+04
216 217 218 219 220
2.696010e+04 2.293713e+04 3.299888e+04 -1.654379e+04 5.158962e+04
221 222 223 224 225
7.816804e+04 7.663229e+04 -3.156698e+03 3.432896e+04 -1.929492e+03
226 227 228 229 230
2.088802e+04 7.816804e+04 6.352101e+04 3.600235e+04 4.264214e+04
231 232 233 234 235
6.133615e+04 7.816804e+04 -1.624027e+03 3.600235e+04 3.605292e+04
236 237 238 239 240
3.600235e+04 2.723941e+04 7.816804e+04 7.816804e+04 3.600235e+04
241 242 243 244 245
6.911561e+04 3.600235e+04 6.883354e+04 3.852101e+04 1.128268e+04
246 247 248 249 250
1.128268e+04 3.956546e+03 1.128268e+04 -3.088301e+04 5.369742e+04
251 252 253 254 255
1.128268e+04 6.171967e+04 -5.658248e+04 -1.706214e+05 3.600235e+04
256 257 258 259 260
2.306742e+04 -6.896001e+03 1.128268e+04 7.816804e+04 1.128268e+04
261 262 263 264 265
1.128268e+04 2.574404e+04 7.838025e+04 1.128268e+04 7.816804e+04
266 267 268 269 270
-3.088301e+04 1.051092e+05 1.128268e+04 4.108224e+04 1.128268e+04
271 272 273 274 275
1.128268e+04 1.128268e+04 6.072985e+02 3.600235e+04 5.887555e+04
276 277 278 279 280
-2.433108e+03 7.816804e+04 -8.162836e+04 4.916413e+04 1.128268e+04
281 282 283 284 285
3.600235e+04 3.950952e+04 3.600235e+04 6.270590e+04 7.816804e+04
286 287 288 289 290
2.043001e+01 3.203934e+04 -3.088301e+04 -6.111311e+02 1.128268e+04
291 292 293 294 295
2.281855e+04 6.682024e+04 3.600235e+04 5.115761e+04 2.037284e+04
296 297 298 299 300
3.600235e+04 -3.088301e+04 7.816804e+04 2.248674e+04 2.575388e+03
301 302 303 304 305
6.140890e+04 7.666630e+04 1.128268e+04 3.792899e+04 1.794995e+04
306 307 308 309 310
1.128268e+04 -2.309127e+04 8.272152e+03 1.128268e+04 1.128268e+04
311 312 313 314 315
2.083068e+04 1.128268e+04 -5.046879e+04 2.417100e+04 4.523385e+04
316 317 318 319 320
3.617547e+04 -1.752119e+04 2.847441e+04 3.600235e+04 7.816804e+04
321 322 323 324 325
7.816804e+04 7.816804e+04 -3.088301e+04 7.816804e+04 7.350799e+04
326 327 328 329 330
-3.088301e+04 1.128268e+04 7.816804e+04 7.816804e+04 7.816804e+04
331 332 333 334 335
3.600235e+04 7.816804e+04 6.976128e+04 3.448651e+04 2.602428e+04
336 337 338 339 340
1.796446e+04 3.600235e+04 3.243381e+04 6.522883e+04 -8.705010e+03
341 342 343 344 345
3.600235e+04 3.600235e+04 7.816804e+04 1.519857e+03 3.600235e+04
346 347 348 349 350
2.046556e+04 3.600235e+04 7.816804e+04 4.847882e+03 7.591426e+04
351 352 353 354 355
2.285412e+04 6.758874e+04 2.980277e+04 3.600235e+04 3.600235e+04
356 357 358 359 360
-1.050841e+04 3.600235e+04 3.600235e+04 7.816804e+04 2.999808e+04
361 362 363 364 365
7.816804e+04 3.600235e+04 -6.124955e+04 6.485885e+04 7.816804e+04
366 367 368 369 370
7.816804e+04 3.600235e+04 7.816804e+04 2.291243e+04 7.816804e+04
371 372 373 374 375
7.125759e+04 2.228522e+04 5.794993e+04 7.816804e+04 7.816804e+04
376 377 378 379 380
4.627300e+04 2.517615e+04 -2.442581e+03 7.816804e+04 7.816804e+04
381 382 383 384 385
-1.304973e+04 -1.107369e+03 7.816804e+04 1.626159e+04 -8.365301e+04
386 387 388 389 390
4.297161e+04 5.401746e+04 2.985540e+04 5.438925e+04 6.238896e+04
391 392 393 394 395
-3.806724e+01 5.889624e+03 6.055375e+04 -1.501482e+05 3.034063e+04
396 397 398 399 400
-8.392311e+04 -3.247143e+05 3.758911e+04 -1.079982e+05 -6.169317e+03
401 402 403 404 405
2.756067e+04 2.896510e+03 -1.297722e+05 -2.527709e+04 2.666550e+03
406 407 408 409 410
4.559749e+04 3.927075e+04 2.931348e+03 -3.157747e+04 -4.223222e+04
411 412 413 414 415
-1.457241e+05 -6.694617e+04 3.464206e+04 -1.379192e+05 -1.075542e+05
416 417 418 419 420
1.288701e+05 -1.022662e+05 -6.764406e+04 -4.973243e+04 -5.246928e+04
421 422 423 424 425
-1.595566e+05 -5.328526e+05 -3.851035e+05 -1.557869e+04 -2.849142e+05
426 427 428 429 430
-2.636122e+05 -8.778464e+05 -1.739617e+05 -2.086613e+05 -5.306982e+05
> postscript(file="/var/www/html/rcomp/tmp/6esxn1292932180.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 = 430
Frequency = 1
lag(myerror, k = 1) myerror
0 2.945709e+06 NA
1 1.958932e+06 2.945709e+06
2 -1.608369e+06 1.958932e+06
3 -6.533468e+04 -1.608369e+06
4 3.622083e+05 -6.533468e+04
5 1.569258e+05 3.622083e+05
6 3.729608e+05 1.569258e+05
7 -1.382227e+05 3.729608e+05
8 7.579032e+05 -1.382227e+05
9 1.030653e+05 7.579032e+05
10 4.217010e+05 1.030653e+05
11 3.683217e+05 4.217010e+05
12 -8.069330e+05 3.683217e+05
13 -1.555096e+05 -8.069330e+05
14 4.914049e+05 -1.555096e+05
15 -1.921223e+05 4.914049e+05
16 -1.947991e+05 -1.921223e+05
17 -1.065107e+05 -1.947991e+05
18 -5.178519e+05 -1.065107e+05
19 2.502813e+05 -5.178519e+05
20 -4.255778e+04 2.502813e+05
21 -4.201810e+04 -4.255778e+04
22 -5.939069e+05 -4.201810e+04
23 -2.534845e+05 -5.939069e+05
24 6.050526e+04 -2.534845e+05
25 -3.689886e+05 6.050526e+04
26 -2.475860e+05 -3.689886e+05
27 4.290343e+05 -2.475860e+05
28 -1.733252e+05 4.290343e+05
29 2.672605e+05 -1.733252e+05
30 -4.558695e+04 2.672605e+05
31 2.200257e+05 -4.558695e+04
32 1.174329e+05 2.200257e+05
33 3.411520e+04 1.174329e+05
34 -4.667818e+05 3.411520e+04
35 1.238381e+05 -4.667818e+05
36 1.356866e+05 1.238381e+05
37 -8.574527e+04 1.356866e+05
38 -3.037759e+05 -8.574527e+04
39 -4.832529e+05 -3.037759e+05
40 -1.970856e+05 -4.832529e+05
41 -1.326601e+04 -1.970856e+05
42 -1.952668e+05 -1.326601e+04
43 -2.031444e+04 -1.952668e+05
44 -6.659535e+04 -2.031444e+04
45 -2.718181e+05 -6.659535e+04
46 -2.949503e+05 -2.718181e+05
47 -6.060906e+01 -2.949503e+05
48 -1.543829e+05 -6.060906e+01
49 -1.271625e+04 -1.543829e+05
50 -1.427621e+05 -1.271625e+04
51 -2.242160e+05 -1.427621e+05
52 -4.412895e+04 -2.242160e+05
53 -4.303624e+05 -4.412895e+04
54 -1.997302e+04 -4.303624e+05
55 3.734755e+05 -1.997302e+04
56 -8.144209e+03 3.734755e+05
57 -4.126198e+04 -8.144209e+03
58 -1.979514e+05 -4.126198e+04
59 -3.327342e+04 -1.979514e+05
60 -1.996854e+05 -3.327342e+04
61 -1.262337e+04 -1.996854e+05
62 -2.224404e+05 -1.262337e+04
63 -3.673721e+04 -2.224404e+05
64 -4.644577e+05 -3.673721e+04
65 -3.685077e+04 -4.644577e+05
66 -5.278891e+04 -3.685077e+04
67 -8.535061e+04 -5.278891e+04
68 -1.114192e+05 -8.535061e+04
69 -1.687591e+05 -1.114192e+05
70 6.897109e+03 -1.687591e+05
71 -1.018463e+05 6.897109e+03
72 4.283796e+04 -1.018463e+05
73 2.427217e+03 4.283796e+04
74 -1.322708e+05 2.427217e+03
75 -4.093145e+04 -1.322708e+05
76 -8.739013e+04 -4.093145e+04
77 -3.093746e+04 -8.739013e+04
78 -6.970571e+04 -3.093746e+04
79 4.885414e+04 -6.970571e+04
80 -1.180068e+05 4.885414e+04
81 -5.111407e+04 -1.180068e+05
82 1.831975e+04 -5.111407e+04
83 5.431579e+04 1.831975e+04
84 -1.254997e+05 5.431579e+04
85 -1.506956e+04 -1.254997e+05
86 -1.013768e+04 -1.506956e+04
87 -1.270529e+05 -1.013768e+04
88 2.083809e+04 -1.270529e+05
89 -1.259173e+04 2.083809e+04
90 -5.295845e+04 -1.259173e+04
91 8.081311e+04 -5.295845e+04
92 -1.550252e+05 8.081311e+04
93 6.448599e+04 -1.550252e+05
94 -5.258125e+04 6.448599e+04
95 1.040454e+05 -5.258125e+04
96 -3.726888e+04 1.040454e+05
97 -3.086653e+04 -3.726888e+04
98 -3.412860e+04 -3.086653e+04
99 -7.991547e+04 -3.412860e+04
100 -8.335505e+04 -7.991547e+04
101 -1.278187e+05 -8.335505e+04
102 -1.436685e+04 -1.278187e+05
103 5.683226e+04 -1.436685e+04
104 -3.521946e+04 5.683226e+04
105 -2.898622e+05 -3.521946e+04
106 1.064864e+05 -2.898622e+05
107 -1.553016e+05 1.064864e+05
108 -1.252069e+05 -1.553016e+05
109 1.170282e+05 -1.252069e+05
110 -2.213077e+05 1.170282e+05
111 -1.524851e+05 -2.213077e+05
112 -8.239323e+04 -1.524851e+05
113 5.174149e+04 -8.239323e+04
114 -1.529056e+03 5.174149e+04
115 3.358610e+03 -1.529056e+03
116 7.191914e+04 3.358610e+03
117 -1.569444e+05 7.191914e+04
118 3.564294e+04 -1.569444e+05
119 -2.009287e+05 3.564294e+04
120 1.732224e+04 -2.009287e+05
121 -1.134618e+05 1.732224e+04
122 5.424982e+04 -1.134618e+05
123 -1.295347e+04 5.424982e+04
124 -2.256727e+05 -1.295347e+04
125 3.600235e+04 -2.256727e+05
126 1.885287e+04 3.600235e+04
127 2.591166e+03 1.885287e+04
128 7.816804e+04 2.591166e+03
129 2.415626e+04 7.816804e+04
130 8.192810e+04 2.415626e+04
131 7.666630e+04 8.192810e+04
132 2.664539e+03 7.666630e+04
133 1.855164e+04 2.664539e+03
134 2.592867e+04 1.855164e+04
135 6.491739e+04 2.592867e+04
136 3.394453e+04 6.491739e+04
137 1.820984e+04 3.394453e+04
138 3.600235e+04 1.820984e+04
139 3.155258e+04 3.600235e+04
140 7.816804e+04 3.155258e+04
141 -4.587239e+03 7.816804e+04
142 2.975914e+04 -4.587239e+03
143 6.067965e+04 2.975914e+04
144 2.176231e+04 6.067965e+04
145 7.816804e+04 2.176231e+04
146 2.482765e+04 7.816804e+04
147 8.161402e+04 2.482765e+04
148 7.808688e+04 8.161402e+04
149 -1.993497e+03 7.808688e+04
150 -6.439980e+04 -1.993497e+03
151 2.849368e+04 -6.439980e+04
152 4.710343e+04 2.849368e+04
153 3.231582e+04 4.710343e+04
154 3.600235e+04 3.231582e+04
155 3.600235e+04 3.600235e+04
156 3.057385e+04 3.600235e+04
157 4.425873e+04 3.057385e+04
158 3.600235e+04 4.425873e+04
159 3.600235e+04 3.600235e+04
160 7.816804e+04 3.600235e+04
161 3.600235e+04 7.816804e+04
162 1.013438e+04 3.600235e+04
163 1.975804e+04 1.013438e+04
164 3.340956e+04 1.975804e+04
165 7.816804e+04 3.340956e+04
166 1.293051e+05 7.816804e+04
167 -3.697539e+04 1.293051e+05
168 3.600235e+04 -3.697539e+04
169 3.600235e+04 3.600235e+04
170 6.909576e+04 3.600235e+04
171 7.816804e+04 6.909576e+04
172 5.465479e+04 7.816804e+04
173 7.816804e+04 5.465479e+04
174 2.225224e+04 7.816804e+04
175 -1.408540e+05 2.225224e+04
176 1.114887e+04 -1.408540e+05
177 3.600235e+04 1.114887e+04
178 3.600235e+04 3.600235e+04
179 6.651030e+04 3.600235e+04
180 -5.422786e+04 6.651030e+04
181 -2.700112e+04 -5.422786e+04
182 -9.277563e+04 -2.700112e+04
183 4.207503e+02 -9.277563e+04
184 5.017889e+03 4.207503e+02
185 -6.335012e+04 5.017889e+03
186 2.127199e+04 -6.335012e+04
187 3.739862e+04 2.127199e+04
188 3.600235e+04 3.739862e+04
189 2.704002e+04 3.600235e+04
190 3.600235e+04 2.704002e+04
191 2.246269e+04 3.600235e+04
192 7.816804e+04 2.246269e+04
193 3.122124e+04 7.816804e+04
194 1.128268e+04 3.122124e+04
195 -1.625690e+05 1.128268e+04
196 1.128268e+04 -1.625690e+05
197 1.128268e+04 1.128268e+04
198 -9.521193e+04 1.128268e+04
199 1.128268e+04 -9.521193e+04
200 1.128268e+04 1.128268e+04
201 1.128268e+04 1.128268e+04
202 -5.243307e+04 1.128268e+04
203 -3.907353e+04 -5.243307e+04
204 3.153210e+03 -3.907353e+04
205 3.939561e+03 3.153210e+03
206 1.128268e+04 3.939561e+03
207 -3.948385e+04 1.128268e+04
208 -4.716405e+04 -3.948385e+04
209 3.600235e+04 -4.716405e+04
210 7.816804e+04 3.600235e+04
211 7.713665e+04 7.816804e+04
212 2.906961e+04 7.713665e+04
213 5.044902e+04 2.906961e+04
214 7.516457e+04 5.044902e+04
215 2.696010e+04 7.516457e+04
216 2.293713e+04 2.696010e+04
217 3.299888e+04 2.293713e+04
218 -1.654379e+04 3.299888e+04
219 5.158962e+04 -1.654379e+04
220 7.816804e+04 5.158962e+04
221 7.663229e+04 7.816804e+04
222 -3.156698e+03 7.663229e+04
223 3.432896e+04 -3.156698e+03
224 -1.929492e+03 3.432896e+04
225 2.088802e+04 -1.929492e+03
226 7.816804e+04 2.088802e+04
227 6.352101e+04 7.816804e+04
228 3.600235e+04 6.352101e+04
229 4.264214e+04 3.600235e+04
230 6.133615e+04 4.264214e+04
231 7.816804e+04 6.133615e+04
232 -1.624027e+03 7.816804e+04
233 3.600235e+04 -1.624027e+03
234 3.605292e+04 3.600235e+04
235 3.600235e+04 3.605292e+04
236 2.723941e+04 3.600235e+04
237 7.816804e+04 2.723941e+04
238 7.816804e+04 7.816804e+04
239 3.600235e+04 7.816804e+04
240 6.911561e+04 3.600235e+04
241 3.600235e+04 6.911561e+04
242 6.883354e+04 3.600235e+04
243 3.852101e+04 6.883354e+04
244 1.128268e+04 3.852101e+04
245 1.128268e+04 1.128268e+04
246 3.956546e+03 1.128268e+04
247 1.128268e+04 3.956546e+03
248 -3.088301e+04 1.128268e+04
249 5.369742e+04 -3.088301e+04
250 1.128268e+04 5.369742e+04
251 6.171967e+04 1.128268e+04
252 -5.658248e+04 6.171967e+04
253 -1.706214e+05 -5.658248e+04
254 3.600235e+04 -1.706214e+05
255 2.306742e+04 3.600235e+04
256 -6.896001e+03 2.306742e+04
257 1.128268e+04 -6.896001e+03
258 7.816804e+04 1.128268e+04
259 1.128268e+04 7.816804e+04
260 1.128268e+04 1.128268e+04
261 2.574404e+04 1.128268e+04
262 7.838025e+04 2.574404e+04
263 1.128268e+04 7.838025e+04
264 7.816804e+04 1.128268e+04
265 -3.088301e+04 7.816804e+04
266 1.051092e+05 -3.088301e+04
267 1.128268e+04 1.051092e+05
268 4.108224e+04 1.128268e+04
269 1.128268e+04 4.108224e+04
270 1.128268e+04 1.128268e+04
271 1.128268e+04 1.128268e+04
272 6.072985e+02 1.128268e+04
273 3.600235e+04 6.072985e+02
274 5.887555e+04 3.600235e+04
275 -2.433108e+03 5.887555e+04
276 7.816804e+04 -2.433108e+03
277 -8.162836e+04 7.816804e+04
278 4.916413e+04 -8.162836e+04
279 1.128268e+04 4.916413e+04
280 3.600235e+04 1.128268e+04
281 3.950952e+04 3.600235e+04
282 3.600235e+04 3.950952e+04
283 6.270590e+04 3.600235e+04
284 7.816804e+04 6.270590e+04
285 2.043001e+01 7.816804e+04
286 3.203934e+04 2.043001e+01
287 -3.088301e+04 3.203934e+04
288 -6.111311e+02 -3.088301e+04
289 1.128268e+04 -6.111311e+02
290 2.281855e+04 1.128268e+04
291 6.682024e+04 2.281855e+04
292 3.600235e+04 6.682024e+04
293 5.115761e+04 3.600235e+04
294 2.037284e+04 5.115761e+04
295 3.600235e+04 2.037284e+04
296 -3.088301e+04 3.600235e+04
297 7.816804e+04 -3.088301e+04
298 2.248674e+04 7.816804e+04
299 2.575388e+03 2.248674e+04
300 6.140890e+04 2.575388e+03
301 7.666630e+04 6.140890e+04
302 1.128268e+04 7.666630e+04
303 3.792899e+04 1.128268e+04
304 1.794995e+04 3.792899e+04
305 1.128268e+04 1.794995e+04
306 -2.309127e+04 1.128268e+04
307 8.272152e+03 -2.309127e+04
308 1.128268e+04 8.272152e+03
309 1.128268e+04 1.128268e+04
310 2.083068e+04 1.128268e+04
311 1.128268e+04 2.083068e+04
312 -5.046879e+04 1.128268e+04
313 2.417100e+04 -5.046879e+04
314 4.523385e+04 2.417100e+04
315 3.617547e+04 4.523385e+04
316 -1.752119e+04 3.617547e+04
317 2.847441e+04 -1.752119e+04
318 3.600235e+04 2.847441e+04
319 7.816804e+04 3.600235e+04
320 7.816804e+04 7.816804e+04
321 7.816804e+04 7.816804e+04
322 -3.088301e+04 7.816804e+04
323 7.816804e+04 -3.088301e+04
324 7.350799e+04 7.816804e+04
325 -3.088301e+04 7.350799e+04
326 1.128268e+04 -3.088301e+04
327 7.816804e+04 1.128268e+04
328 7.816804e+04 7.816804e+04
329 7.816804e+04 7.816804e+04
330 3.600235e+04 7.816804e+04
331 7.816804e+04 3.600235e+04
332 6.976128e+04 7.816804e+04
333 3.448651e+04 6.976128e+04
334 2.602428e+04 3.448651e+04
335 1.796446e+04 2.602428e+04
336 3.600235e+04 1.796446e+04
337 3.243381e+04 3.600235e+04
338 6.522883e+04 3.243381e+04
339 -8.705010e+03 6.522883e+04
340 3.600235e+04 -8.705010e+03
341 3.600235e+04 3.600235e+04
342 7.816804e+04 3.600235e+04
343 1.519857e+03 7.816804e+04
344 3.600235e+04 1.519857e+03
345 2.046556e+04 3.600235e+04
346 3.600235e+04 2.046556e+04
347 7.816804e+04 3.600235e+04
348 4.847882e+03 7.816804e+04
349 7.591426e+04 4.847882e+03
350 2.285412e+04 7.591426e+04
351 6.758874e+04 2.285412e+04
352 2.980277e+04 6.758874e+04
353 3.600235e+04 2.980277e+04
354 3.600235e+04 3.600235e+04
355 -1.050841e+04 3.600235e+04
356 3.600235e+04 -1.050841e+04
357 3.600235e+04 3.600235e+04
358 7.816804e+04 3.600235e+04
359 2.999808e+04 7.816804e+04
360 7.816804e+04 2.999808e+04
361 3.600235e+04 7.816804e+04
362 -6.124955e+04 3.600235e+04
363 6.485885e+04 -6.124955e+04
364 7.816804e+04 6.485885e+04
365 7.816804e+04 7.816804e+04
366 3.600235e+04 7.816804e+04
367 7.816804e+04 3.600235e+04
368 2.291243e+04 7.816804e+04
369 7.816804e+04 2.291243e+04
370 7.125759e+04 7.816804e+04
371 2.228522e+04 7.125759e+04
372 5.794993e+04 2.228522e+04
373 7.816804e+04 5.794993e+04
374 7.816804e+04 7.816804e+04
375 4.627300e+04 7.816804e+04
376 2.517615e+04 4.627300e+04
377 -2.442581e+03 2.517615e+04
378 7.816804e+04 -2.442581e+03
379 7.816804e+04 7.816804e+04
380 -1.304973e+04 7.816804e+04
381 -1.107369e+03 -1.304973e+04
382 7.816804e+04 -1.107369e+03
383 1.626159e+04 7.816804e+04
384 -8.365301e+04 1.626159e+04
385 4.297161e+04 -8.365301e+04
386 5.401746e+04 4.297161e+04
387 2.985540e+04 5.401746e+04
388 5.438925e+04 2.985540e+04
389 6.238896e+04 5.438925e+04
390 -3.806724e+01 6.238896e+04
391 5.889624e+03 -3.806724e+01
392 6.055375e+04 5.889624e+03
393 -1.501482e+05 6.055375e+04
394 3.034063e+04 -1.501482e+05
395 -8.392311e+04 3.034063e+04
396 -3.247143e+05 -8.392311e+04
397 3.758911e+04 -3.247143e+05
398 -1.079982e+05 3.758911e+04
399 -6.169317e+03 -1.079982e+05
400 2.756067e+04 -6.169317e+03
401 2.896510e+03 2.756067e+04
402 -1.297722e+05 2.896510e+03
403 -2.527709e+04 -1.297722e+05
404 2.666550e+03 -2.527709e+04
405 4.559749e+04 2.666550e+03
406 3.927075e+04 4.559749e+04
407 2.931348e+03 3.927075e+04
408 -3.157747e+04 2.931348e+03
409 -4.223222e+04 -3.157747e+04
410 -1.457241e+05 -4.223222e+04
411 -6.694617e+04 -1.457241e+05
412 3.464206e+04 -6.694617e+04
413 -1.379192e+05 3.464206e+04
414 -1.075542e+05 -1.379192e+05
415 1.288701e+05 -1.075542e+05
416 -1.022662e+05 1.288701e+05
417 -6.764406e+04 -1.022662e+05
418 -4.973243e+04 -6.764406e+04
419 -5.246928e+04 -4.973243e+04
420 -1.595566e+05 -5.246928e+04
421 -5.328526e+05 -1.595566e+05
422 -3.851035e+05 -5.328526e+05
423 -1.557869e+04 -3.851035e+05
424 -2.849142e+05 -1.557869e+04
425 -2.636122e+05 -2.849142e+05
426 -8.778464e+05 -2.636122e+05
427 -1.739617e+05 -8.778464e+05
428 -2.086613e+05 -1.739617e+05
429 -5.306982e+05 -2.086613e+05
430 NA -5.306982e+05
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] 1.958932e+06 2.945709e+06
[2,] -1.608369e+06 1.958932e+06
[3,] -6.533468e+04 -1.608369e+06
[4,] 3.622083e+05 -6.533468e+04
[5,] 1.569258e+05 3.622083e+05
[6,] 3.729608e+05 1.569258e+05
[7,] -1.382227e+05 3.729608e+05
[8,] 7.579032e+05 -1.382227e+05
[9,] 1.030653e+05 7.579032e+05
[10,] 4.217010e+05 1.030653e+05
[11,] 3.683217e+05 4.217010e+05
[12,] -8.069330e+05 3.683217e+05
[13,] -1.555096e+05 -8.069330e+05
[14,] 4.914049e+05 -1.555096e+05
[15,] -1.921223e+05 4.914049e+05
[16,] -1.947991e+05 -1.921223e+05
[17,] -1.065107e+05 -1.947991e+05
[18,] -5.178519e+05 -1.065107e+05
[19,] 2.502813e+05 -5.178519e+05
[20,] -4.255778e+04 2.502813e+05
[21,] -4.201810e+04 -4.255778e+04
[22,] -5.939069e+05 -4.201810e+04
[23,] -2.534845e+05 -5.939069e+05
[24,] 6.050526e+04 -2.534845e+05
[25,] -3.689886e+05 6.050526e+04
[26,] -2.475860e+05 -3.689886e+05
[27,] 4.290343e+05 -2.475860e+05
[28,] -1.733252e+05 4.290343e+05
[29,] 2.672605e+05 -1.733252e+05
[30,] -4.558695e+04 2.672605e+05
[31,] 2.200257e+05 -4.558695e+04
[32,] 1.174329e+05 2.200257e+05
[33,] 3.411520e+04 1.174329e+05
[34,] -4.667818e+05 3.411520e+04
[35,] 1.238381e+05 -4.667818e+05
[36,] 1.356866e+05 1.238381e+05
[37,] -8.574527e+04 1.356866e+05
[38,] -3.037759e+05 -8.574527e+04
[39,] -4.832529e+05 -3.037759e+05
[40,] -1.970856e+05 -4.832529e+05
[41,] -1.326601e+04 -1.970856e+05
[42,] -1.952668e+05 -1.326601e+04
[43,] -2.031444e+04 -1.952668e+05
[44,] -6.659535e+04 -2.031444e+04
[45,] -2.718181e+05 -6.659535e+04
[46,] -2.949503e+05 -2.718181e+05
[47,] -6.060906e+01 -2.949503e+05
[48,] -1.543829e+05 -6.060906e+01
[49,] -1.271625e+04 -1.543829e+05
[50,] -1.427621e+05 -1.271625e+04
[51,] -2.242160e+05 -1.427621e+05
[52,] -4.412895e+04 -2.242160e+05
[53,] -4.303624e+05 -4.412895e+04
[54,] -1.997302e+04 -4.303624e+05
[55,] 3.734755e+05 -1.997302e+04
[56,] -8.144209e+03 3.734755e+05
[57,] -4.126198e+04 -8.144209e+03
[58,] -1.979514e+05 -4.126198e+04
[59,] -3.327342e+04 -1.979514e+05
[60,] -1.996854e+05 -3.327342e+04
[61,] -1.262337e+04 -1.996854e+05
[62,] -2.224404e+05 -1.262337e+04
[63,] -3.673721e+04 -2.224404e+05
[64,] -4.644577e+05 -3.673721e+04
[65,] -3.685077e+04 -4.644577e+05
[66,] -5.278891e+04 -3.685077e+04
[67,] -8.535061e+04 -5.278891e+04
[68,] -1.114192e+05 -8.535061e+04
[69,] -1.687591e+05 -1.114192e+05
[70,] 6.897109e+03 -1.687591e+05
[71,] -1.018463e+05 6.897109e+03
[72,] 4.283796e+04 -1.018463e+05
[73,] 2.427217e+03 4.283796e+04
[74,] -1.322708e+05 2.427217e+03
[75,] -4.093145e+04 -1.322708e+05
[76,] -8.739013e+04 -4.093145e+04
[77,] -3.093746e+04 -8.739013e+04
[78,] -6.970571e+04 -3.093746e+04
[79,] 4.885414e+04 -6.970571e+04
[80,] -1.180068e+05 4.885414e+04
[81,] -5.111407e+04 -1.180068e+05
[82,] 1.831975e+04 -5.111407e+04
[83,] 5.431579e+04 1.831975e+04
[84,] -1.254997e+05 5.431579e+04
[85,] -1.506956e+04 -1.254997e+05
[86,] -1.013768e+04 -1.506956e+04
[87,] -1.270529e+05 -1.013768e+04
[88,] 2.083809e+04 -1.270529e+05
[89,] -1.259173e+04 2.083809e+04
[90,] -5.295845e+04 -1.259173e+04
[91,] 8.081311e+04 -5.295845e+04
[92,] -1.550252e+05 8.081311e+04
[93,] 6.448599e+04 -1.550252e+05
[94,] -5.258125e+04 6.448599e+04
[95,] 1.040454e+05 -5.258125e+04
[96,] -3.726888e+04 1.040454e+05
[97,] -3.086653e+04 -3.726888e+04
[98,] -3.412860e+04 -3.086653e+04
[99,] -7.991547e+04 -3.412860e+04
[100,] -8.335505e+04 -7.991547e+04
[101,] -1.278187e+05 -8.335505e+04
[102,] -1.436685e+04 -1.278187e+05
[103,] 5.683226e+04 -1.436685e+04
[104,] -3.521946e+04 5.683226e+04
[105,] -2.898622e+05 -3.521946e+04
[106,] 1.064864e+05 -2.898622e+05
[107,] -1.553016e+05 1.064864e+05
[108,] -1.252069e+05 -1.553016e+05
[109,] 1.170282e+05 -1.252069e+05
[110,] -2.213077e+05 1.170282e+05
[111,] -1.524851e+05 -2.213077e+05
[112,] -8.239323e+04 -1.524851e+05
[113,] 5.174149e+04 -8.239323e+04
[114,] -1.529056e+03 5.174149e+04
[115,] 3.358610e+03 -1.529056e+03
[116,] 7.191914e+04 3.358610e+03
[117,] -1.569444e+05 7.191914e+04
[118,] 3.564294e+04 -1.569444e+05
[119,] -2.009287e+05 3.564294e+04
[120,] 1.732224e+04 -2.009287e+05
[121,] -1.134618e+05 1.732224e+04
[122,] 5.424982e+04 -1.134618e+05
[123,] -1.295347e+04 5.424982e+04
[124,] -2.256727e+05 -1.295347e+04
[125,] 3.600235e+04 -2.256727e+05
[126,] 1.885287e+04 3.600235e+04
[127,] 2.591166e+03 1.885287e+04
[128,] 7.816804e+04 2.591166e+03
[129,] 2.415626e+04 7.816804e+04
[130,] 8.192810e+04 2.415626e+04
[131,] 7.666630e+04 8.192810e+04
[132,] 2.664539e+03 7.666630e+04
[133,] 1.855164e+04 2.664539e+03
[134,] 2.592867e+04 1.855164e+04
[135,] 6.491739e+04 2.592867e+04
[136,] 3.394453e+04 6.491739e+04
[137,] 1.820984e+04 3.394453e+04
[138,] 3.600235e+04 1.820984e+04
[139,] 3.155258e+04 3.600235e+04
[140,] 7.816804e+04 3.155258e+04
[141,] -4.587239e+03 7.816804e+04
[142,] 2.975914e+04 -4.587239e+03
[143,] 6.067965e+04 2.975914e+04
[144,] 2.176231e+04 6.067965e+04
[145,] 7.816804e+04 2.176231e+04
[146,] 2.482765e+04 7.816804e+04
[147,] 8.161402e+04 2.482765e+04
[148,] 7.808688e+04 8.161402e+04
[149,] -1.993497e+03 7.808688e+04
[150,] -6.439980e+04 -1.993497e+03
[151,] 2.849368e+04 -6.439980e+04
[152,] 4.710343e+04 2.849368e+04
[153,] 3.231582e+04 4.710343e+04
[154,] 3.600235e+04 3.231582e+04
[155,] 3.600235e+04 3.600235e+04
[156,] 3.057385e+04 3.600235e+04
[157,] 4.425873e+04 3.057385e+04
[158,] 3.600235e+04 4.425873e+04
[159,] 3.600235e+04 3.600235e+04
[160,] 7.816804e+04 3.600235e+04
[161,] 3.600235e+04 7.816804e+04
[162,] 1.013438e+04 3.600235e+04
[163,] 1.975804e+04 1.013438e+04
[164,] 3.340956e+04 1.975804e+04
[165,] 7.816804e+04 3.340956e+04
[166,] 1.293051e+05 7.816804e+04
[167,] -3.697539e+04 1.293051e+05
[168,] 3.600235e+04 -3.697539e+04
[169,] 3.600235e+04 3.600235e+04
[170,] 6.909576e+04 3.600235e+04
[171,] 7.816804e+04 6.909576e+04
[172,] 5.465479e+04 7.816804e+04
[173,] 7.816804e+04 5.465479e+04
[174,] 2.225224e+04 7.816804e+04
[175,] -1.408540e+05 2.225224e+04
[176,] 1.114887e+04 -1.408540e+05
[177,] 3.600235e+04 1.114887e+04
[178,] 3.600235e+04 3.600235e+04
[179,] 6.651030e+04 3.600235e+04
[180,] -5.422786e+04 6.651030e+04
[181,] -2.700112e+04 -5.422786e+04
[182,] -9.277563e+04 -2.700112e+04
[183,] 4.207503e+02 -9.277563e+04
[184,] 5.017889e+03 4.207503e+02
[185,] -6.335012e+04 5.017889e+03
[186,] 2.127199e+04 -6.335012e+04
[187,] 3.739862e+04 2.127199e+04
[188,] 3.600235e+04 3.739862e+04
[189,] 2.704002e+04 3.600235e+04
[190,] 3.600235e+04 2.704002e+04
[191,] 2.246269e+04 3.600235e+04
[192,] 7.816804e+04 2.246269e+04
[193,] 3.122124e+04 7.816804e+04
[194,] 1.128268e+04 3.122124e+04
[195,] -1.625690e+05 1.128268e+04
[196,] 1.128268e+04 -1.625690e+05
[197,] 1.128268e+04 1.128268e+04
[198,] -9.521193e+04 1.128268e+04
[199,] 1.128268e+04 -9.521193e+04
[200,] 1.128268e+04 1.128268e+04
[201,] 1.128268e+04 1.128268e+04
[202,] -5.243307e+04 1.128268e+04
[203,] -3.907353e+04 -5.243307e+04
[204,] 3.153210e+03 -3.907353e+04
[205,] 3.939561e+03 3.153210e+03
[206,] 1.128268e+04 3.939561e+03
[207,] -3.948385e+04 1.128268e+04
[208,] -4.716405e+04 -3.948385e+04
[209,] 3.600235e+04 -4.716405e+04
[210,] 7.816804e+04 3.600235e+04
[211,] 7.713665e+04 7.816804e+04
[212,] 2.906961e+04 7.713665e+04
[213,] 5.044902e+04 2.906961e+04
[214,] 7.516457e+04 5.044902e+04
[215,] 2.696010e+04 7.516457e+04
[216,] 2.293713e+04 2.696010e+04
[217,] 3.299888e+04 2.293713e+04
[218,] -1.654379e+04 3.299888e+04
[219,] 5.158962e+04 -1.654379e+04
[220,] 7.816804e+04 5.158962e+04
[221,] 7.663229e+04 7.816804e+04
[222,] -3.156698e+03 7.663229e+04
[223,] 3.432896e+04 -3.156698e+03
[224,] -1.929492e+03 3.432896e+04
[225,] 2.088802e+04 -1.929492e+03
[226,] 7.816804e+04 2.088802e+04
[227,] 6.352101e+04 7.816804e+04
[228,] 3.600235e+04 6.352101e+04
[229,] 4.264214e+04 3.600235e+04
[230,] 6.133615e+04 4.264214e+04
[231,] 7.816804e+04 6.133615e+04
[232,] -1.624027e+03 7.816804e+04
[233,] 3.600235e+04 -1.624027e+03
[234,] 3.605292e+04 3.600235e+04
[235,] 3.600235e+04 3.605292e+04
[236,] 2.723941e+04 3.600235e+04
[237,] 7.816804e+04 2.723941e+04
[238,] 7.816804e+04 7.816804e+04
[239,] 3.600235e+04 7.816804e+04
[240,] 6.911561e+04 3.600235e+04
[241,] 3.600235e+04 6.911561e+04
[242,] 6.883354e+04 3.600235e+04
[243,] 3.852101e+04 6.883354e+04
[244,] 1.128268e+04 3.852101e+04
[245,] 1.128268e+04 1.128268e+04
[246,] 3.956546e+03 1.128268e+04
[247,] 1.128268e+04 3.956546e+03
[248,] -3.088301e+04 1.128268e+04
[249,] 5.369742e+04 -3.088301e+04
[250,] 1.128268e+04 5.369742e+04
[251,] 6.171967e+04 1.128268e+04
[252,] -5.658248e+04 6.171967e+04
[253,] -1.706214e+05 -5.658248e+04
[254,] 3.600235e+04 -1.706214e+05
[255,] 2.306742e+04 3.600235e+04
[256,] -6.896001e+03 2.306742e+04
[257,] 1.128268e+04 -6.896001e+03
[258,] 7.816804e+04 1.128268e+04
[259,] 1.128268e+04 7.816804e+04
[260,] 1.128268e+04 1.128268e+04
[261,] 2.574404e+04 1.128268e+04
[262,] 7.838025e+04 2.574404e+04
[263,] 1.128268e+04 7.838025e+04
[264,] 7.816804e+04 1.128268e+04
[265,] -3.088301e+04 7.816804e+04
[266,] 1.051092e+05 -3.088301e+04
[267,] 1.128268e+04 1.051092e+05
[268,] 4.108224e+04 1.128268e+04
[269,] 1.128268e+04 4.108224e+04
[270,] 1.128268e+04 1.128268e+04
[271,] 1.128268e+04 1.128268e+04
[272,] 6.072985e+02 1.128268e+04
[273,] 3.600235e+04 6.072985e+02
[274,] 5.887555e+04 3.600235e+04
[275,] -2.433108e+03 5.887555e+04
[276,] 7.816804e+04 -2.433108e+03
[277,] -8.162836e+04 7.816804e+04
[278,] 4.916413e+04 -8.162836e+04
[279,] 1.128268e+04 4.916413e+04
[280,] 3.600235e+04 1.128268e+04
[281,] 3.950952e+04 3.600235e+04
[282,] 3.600235e+04 3.950952e+04
[283,] 6.270590e+04 3.600235e+04
[284,] 7.816804e+04 6.270590e+04
[285,] 2.043001e+01 7.816804e+04
[286,] 3.203934e+04 2.043001e+01
[287,] -3.088301e+04 3.203934e+04
[288,] -6.111311e+02 -3.088301e+04
[289,] 1.128268e+04 -6.111311e+02
[290,] 2.281855e+04 1.128268e+04
[291,] 6.682024e+04 2.281855e+04
[292,] 3.600235e+04 6.682024e+04
[293,] 5.115761e+04 3.600235e+04
[294,] 2.037284e+04 5.115761e+04
[295,] 3.600235e+04 2.037284e+04
[296,] -3.088301e+04 3.600235e+04
[297,] 7.816804e+04 -3.088301e+04
[298,] 2.248674e+04 7.816804e+04
[299,] 2.575388e+03 2.248674e+04
[300,] 6.140890e+04 2.575388e+03
[301,] 7.666630e+04 6.140890e+04
[302,] 1.128268e+04 7.666630e+04
[303,] 3.792899e+04 1.128268e+04
[304,] 1.794995e+04 3.792899e+04
[305,] 1.128268e+04 1.794995e+04
[306,] -2.309127e+04 1.128268e+04
[307,] 8.272152e+03 -2.309127e+04
[308,] 1.128268e+04 8.272152e+03
[309,] 1.128268e+04 1.128268e+04
[310,] 2.083068e+04 1.128268e+04
[311,] 1.128268e+04 2.083068e+04
[312,] -5.046879e+04 1.128268e+04
[313,] 2.417100e+04 -5.046879e+04
[314,] 4.523385e+04 2.417100e+04
[315,] 3.617547e+04 4.523385e+04
[316,] -1.752119e+04 3.617547e+04
[317,] 2.847441e+04 -1.752119e+04
[318,] 3.600235e+04 2.847441e+04
[319,] 7.816804e+04 3.600235e+04
[320,] 7.816804e+04 7.816804e+04
[321,] 7.816804e+04 7.816804e+04
[322,] -3.088301e+04 7.816804e+04
[323,] 7.816804e+04 -3.088301e+04
[324,] 7.350799e+04 7.816804e+04
[325,] -3.088301e+04 7.350799e+04
[326,] 1.128268e+04 -3.088301e+04
[327,] 7.816804e+04 1.128268e+04
[328,] 7.816804e+04 7.816804e+04
[329,] 7.816804e+04 7.816804e+04
[330,] 3.600235e+04 7.816804e+04
[331,] 7.816804e+04 3.600235e+04
[332,] 6.976128e+04 7.816804e+04
[333,] 3.448651e+04 6.976128e+04
[334,] 2.602428e+04 3.448651e+04
[335,] 1.796446e+04 2.602428e+04
[336,] 3.600235e+04 1.796446e+04
[337,] 3.243381e+04 3.600235e+04
[338,] 6.522883e+04 3.243381e+04
[339,] -8.705010e+03 6.522883e+04
[340,] 3.600235e+04 -8.705010e+03
[341,] 3.600235e+04 3.600235e+04
[342,] 7.816804e+04 3.600235e+04
[343,] 1.519857e+03 7.816804e+04
[344,] 3.600235e+04 1.519857e+03
[345,] 2.046556e+04 3.600235e+04
[346,] 3.600235e+04 2.046556e+04
[347,] 7.816804e+04 3.600235e+04
[348,] 4.847882e+03 7.816804e+04
[349,] 7.591426e+04 4.847882e+03
[350,] 2.285412e+04 7.591426e+04
[351,] 6.758874e+04 2.285412e+04
[352,] 2.980277e+04 6.758874e+04
[353,] 3.600235e+04 2.980277e+04
[354,] 3.600235e+04 3.600235e+04
[355,] -1.050841e+04 3.600235e+04
[356,] 3.600235e+04 -1.050841e+04
[357,] 3.600235e+04 3.600235e+04
[358,] 7.816804e+04 3.600235e+04
[359,] 2.999808e+04 7.816804e+04
[360,] 7.816804e+04 2.999808e+04
[361,] 3.600235e+04 7.816804e+04
[362,] -6.124955e+04 3.600235e+04
[363,] 6.485885e+04 -6.124955e+04
[364,] 7.816804e+04 6.485885e+04
[365,] 7.816804e+04 7.816804e+04
[366,] 3.600235e+04 7.816804e+04
[367,] 7.816804e+04 3.600235e+04
[368,] 2.291243e+04 7.816804e+04
[369,] 7.816804e+04 2.291243e+04
[370,] 7.125759e+04 7.816804e+04
[371,] 2.228522e+04 7.125759e+04
[372,] 5.794993e+04 2.228522e+04
[373,] 7.816804e+04 5.794993e+04
[374,] 7.816804e+04 7.816804e+04
[375,] 4.627300e+04 7.816804e+04
[376,] 2.517615e+04 4.627300e+04
[377,] -2.442581e+03 2.517615e+04
[378,] 7.816804e+04 -2.442581e+03
[379,] 7.816804e+04 7.816804e+04
[380,] -1.304973e+04 7.816804e+04
[381,] -1.107369e+03 -1.304973e+04
[382,] 7.816804e+04 -1.107369e+03
[383,] 1.626159e+04 7.816804e+04
[384,] -8.365301e+04 1.626159e+04
[385,] 4.297161e+04 -8.365301e+04
[386,] 5.401746e+04 4.297161e+04
[387,] 2.985540e+04 5.401746e+04
[388,] 5.438925e+04 2.985540e+04
[389,] 6.238896e+04 5.438925e+04
[390,] -3.806724e+01 6.238896e+04
[391,] 5.889624e+03 -3.806724e+01
[392,] 6.055375e+04 5.889624e+03
[393,] -1.501482e+05 6.055375e+04
[394,] 3.034063e+04 -1.501482e+05
[395,] -8.392311e+04 3.034063e+04
[396,] -3.247143e+05 -8.392311e+04
[397,] 3.758911e+04 -3.247143e+05
[398,] -1.079982e+05 3.758911e+04
[399,] -6.169317e+03 -1.079982e+05
[400,] 2.756067e+04 -6.169317e+03
[401,] 2.896510e+03 2.756067e+04
[402,] -1.297722e+05 2.896510e+03
[403,] -2.527709e+04 -1.297722e+05
[404,] 2.666550e+03 -2.527709e+04
[405,] 4.559749e+04 2.666550e+03
[406,] 3.927075e+04 4.559749e+04
[407,] 2.931348e+03 3.927075e+04
[408,] -3.157747e+04 2.931348e+03
[409,] -4.223222e+04 -3.157747e+04
[410,] -1.457241e+05 -4.223222e+04
[411,] -6.694617e+04 -1.457241e+05
[412,] 3.464206e+04 -6.694617e+04
[413,] -1.379192e+05 3.464206e+04
[414,] -1.075542e+05 -1.379192e+05
[415,] 1.288701e+05 -1.075542e+05
[416,] -1.022662e+05 1.288701e+05
[417,] -6.764406e+04 -1.022662e+05
[418,] -4.973243e+04 -6.764406e+04
[419,] -5.246928e+04 -4.973243e+04
[420,] -1.595566e+05 -5.246928e+04
[421,] -5.328526e+05 -1.595566e+05
[422,] -3.851035e+05 -5.328526e+05
[423,] -1.557869e+04 -3.851035e+05
[424,] -2.849142e+05 -1.557869e+04
[425,] -2.636122e+05 -2.849142e+05
[426,] -8.778464e+05 -2.636122e+05
[427,] -1.739617e+05 -8.778464e+05
[428,] -2.086613e+05 -1.739617e+05
[429,] -5.306982e+05 -2.086613e+05
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 1.958932e+06 2.945709e+06
2 -1.608369e+06 1.958932e+06
3 -6.533468e+04 -1.608369e+06
4 3.622083e+05 -6.533468e+04
5 1.569258e+05 3.622083e+05
6 3.729608e+05 1.569258e+05
7 -1.382227e+05 3.729608e+05
8 7.579032e+05 -1.382227e+05
9 1.030653e+05 7.579032e+05
10 4.217010e+05 1.030653e+05
11 3.683217e+05 4.217010e+05
12 -8.069330e+05 3.683217e+05
13 -1.555096e+05 -8.069330e+05
14 4.914049e+05 -1.555096e+05
15 -1.921223e+05 4.914049e+05
16 -1.947991e+05 -1.921223e+05
17 -1.065107e+05 -1.947991e+05
18 -5.178519e+05 -1.065107e+05
19 2.502813e+05 -5.178519e+05
20 -4.255778e+04 2.502813e+05
21 -4.201810e+04 -4.255778e+04
22 -5.939069e+05 -4.201810e+04
23 -2.534845e+05 -5.939069e+05
24 6.050526e+04 -2.534845e+05
25 -3.689886e+05 6.050526e+04
26 -2.475860e+05 -3.689886e+05
27 4.290343e+05 -2.475860e+05
28 -1.733252e+05 4.290343e+05
29 2.672605e+05 -1.733252e+05
30 -4.558695e+04 2.672605e+05
31 2.200257e+05 -4.558695e+04
32 1.174329e+05 2.200257e+05
33 3.411520e+04 1.174329e+05
34 -4.667818e+05 3.411520e+04
35 1.238381e+05 -4.667818e+05
36 1.356866e+05 1.238381e+05
37 -8.574527e+04 1.356866e+05
38 -3.037759e+05 -8.574527e+04
39 -4.832529e+05 -3.037759e+05
40 -1.970856e+05 -4.832529e+05
41 -1.326601e+04 -1.970856e+05
42 -1.952668e+05 -1.326601e+04
43 -2.031444e+04 -1.952668e+05
44 -6.659535e+04 -2.031444e+04
45 -2.718181e+05 -6.659535e+04
46 -2.949503e+05 -2.718181e+05
47 -6.060906e+01 -2.949503e+05
48 -1.543829e+05 -6.060906e+01
49 -1.271625e+04 -1.543829e+05
50 -1.427621e+05 -1.271625e+04
51 -2.242160e+05 -1.427621e+05
52 -4.412895e+04 -2.242160e+05
53 -4.303624e+05 -4.412895e+04
54 -1.997302e+04 -4.303624e+05
55 3.734755e+05 -1.997302e+04
56 -8.144209e+03 3.734755e+05
57 -4.126198e+04 -8.144209e+03
58 -1.979514e+05 -4.126198e+04
59 -3.327342e+04 -1.979514e+05
60 -1.996854e+05 -3.327342e+04
61 -1.262337e+04 -1.996854e+05
62 -2.224404e+05 -1.262337e+04
63 -3.673721e+04 -2.224404e+05
64 -4.644577e+05 -3.673721e+04
65 -3.685077e+04 -4.644577e+05
66 -5.278891e+04 -3.685077e+04
67 -8.535061e+04 -5.278891e+04
68 -1.114192e+05 -8.535061e+04
69 -1.687591e+05 -1.114192e+05
70 6.897109e+03 -1.687591e+05
71 -1.018463e+05 6.897109e+03
72 4.283796e+04 -1.018463e+05
73 2.427217e+03 4.283796e+04
74 -1.322708e+05 2.427217e+03
75 -4.093145e+04 -1.322708e+05
76 -8.739013e+04 -4.093145e+04
77 -3.093746e+04 -8.739013e+04
78 -6.970571e+04 -3.093746e+04
79 4.885414e+04 -6.970571e+04
80 -1.180068e+05 4.885414e+04
81 -5.111407e+04 -1.180068e+05
82 1.831975e+04 -5.111407e+04
83 5.431579e+04 1.831975e+04
84 -1.254997e+05 5.431579e+04
85 -1.506956e+04 -1.254997e+05
86 -1.013768e+04 -1.506956e+04
87 -1.270529e+05 -1.013768e+04
88 2.083809e+04 -1.270529e+05
89 -1.259173e+04 2.083809e+04
90 -5.295845e+04 -1.259173e+04
91 8.081311e+04 -5.295845e+04
92 -1.550252e+05 8.081311e+04
93 6.448599e+04 -1.550252e+05
94 -5.258125e+04 6.448599e+04
95 1.040454e+05 -5.258125e+04
96 -3.726888e+04 1.040454e+05
97 -3.086653e+04 -3.726888e+04
98 -3.412860e+04 -3.086653e+04
99 -7.991547e+04 -3.412860e+04
100 -8.335505e+04 -7.991547e+04
101 -1.278187e+05 -8.335505e+04
102 -1.436685e+04 -1.278187e+05
103 5.683226e+04 -1.436685e+04
104 -3.521946e+04 5.683226e+04
105 -2.898622e+05 -3.521946e+04
106 1.064864e+05 -2.898622e+05
107 -1.553016e+05 1.064864e+05
108 -1.252069e+05 -1.553016e+05
109 1.170282e+05 -1.252069e+05
110 -2.213077e+05 1.170282e+05
111 -1.524851e+05 -2.213077e+05
112 -8.239323e+04 -1.524851e+05
113 5.174149e+04 -8.239323e+04
114 -1.529056e+03 5.174149e+04
115 3.358610e+03 -1.529056e+03
116 7.191914e+04 3.358610e+03
117 -1.569444e+05 7.191914e+04
118 3.564294e+04 -1.569444e+05
119 -2.009287e+05 3.564294e+04
120 1.732224e+04 -2.009287e+05
121 -1.134618e+05 1.732224e+04
122 5.424982e+04 -1.134618e+05
123 -1.295347e+04 5.424982e+04
124 -2.256727e+05 -1.295347e+04
125 3.600235e+04 -2.256727e+05
126 1.885287e+04 3.600235e+04
127 2.591166e+03 1.885287e+04
128 7.816804e+04 2.591166e+03
129 2.415626e+04 7.816804e+04
130 8.192810e+04 2.415626e+04
131 7.666630e+04 8.192810e+04
132 2.664539e+03 7.666630e+04
133 1.855164e+04 2.664539e+03
134 2.592867e+04 1.855164e+04
135 6.491739e+04 2.592867e+04
136 3.394453e+04 6.491739e+04
137 1.820984e+04 3.394453e+04
138 3.600235e+04 1.820984e+04
139 3.155258e+04 3.600235e+04
140 7.816804e+04 3.155258e+04
141 -4.587239e+03 7.816804e+04
142 2.975914e+04 -4.587239e+03
143 6.067965e+04 2.975914e+04
144 2.176231e+04 6.067965e+04
145 7.816804e+04 2.176231e+04
146 2.482765e+04 7.816804e+04
147 8.161402e+04 2.482765e+04
148 7.808688e+04 8.161402e+04
149 -1.993497e+03 7.808688e+04
150 -6.439980e+04 -1.993497e+03
151 2.849368e+04 -6.439980e+04
152 4.710343e+04 2.849368e+04
153 3.231582e+04 4.710343e+04
154 3.600235e+04 3.231582e+04
155 3.600235e+04 3.600235e+04
156 3.057385e+04 3.600235e+04
157 4.425873e+04 3.057385e+04
158 3.600235e+04 4.425873e+04
159 3.600235e+04 3.600235e+04
160 7.816804e+04 3.600235e+04
161 3.600235e+04 7.816804e+04
162 1.013438e+04 3.600235e+04
163 1.975804e+04 1.013438e+04
164 3.340956e+04 1.975804e+04
165 7.816804e+04 3.340956e+04
166 1.293051e+05 7.816804e+04
167 -3.697539e+04 1.293051e+05
168 3.600235e+04 -3.697539e+04
169 3.600235e+04 3.600235e+04
170 6.909576e+04 3.600235e+04
171 7.816804e+04 6.909576e+04
172 5.465479e+04 7.816804e+04
173 7.816804e+04 5.465479e+04
174 2.225224e+04 7.816804e+04
175 -1.408540e+05 2.225224e+04
176 1.114887e+04 -1.408540e+05
177 3.600235e+04 1.114887e+04
178 3.600235e+04 3.600235e+04
179 6.651030e+04 3.600235e+04
180 -5.422786e+04 6.651030e+04
181 -2.700112e+04 -5.422786e+04
182 -9.277563e+04 -2.700112e+04
183 4.207503e+02 -9.277563e+04
184 5.017889e+03 4.207503e+02
185 -6.335012e+04 5.017889e+03
186 2.127199e+04 -6.335012e+04
187 3.739862e+04 2.127199e+04
188 3.600235e+04 3.739862e+04
189 2.704002e+04 3.600235e+04
190 3.600235e+04 2.704002e+04
191 2.246269e+04 3.600235e+04
192 7.816804e+04 2.246269e+04
193 3.122124e+04 7.816804e+04
194 1.128268e+04 3.122124e+04
195 -1.625690e+05 1.128268e+04
196 1.128268e+04 -1.625690e+05
197 1.128268e+04 1.128268e+04
198 -9.521193e+04 1.128268e+04
199 1.128268e+04 -9.521193e+04
200 1.128268e+04 1.128268e+04
201 1.128268e+04 1.128268e+04
202 -5.243307e+04 1.128268e+04
203 -3.907353e+04 -5.243307e+04
204 3.153210e+03 -3.907353e+04
205 3.939561e+03 3.153210e+03
206 1.128268e+04 3.939561e+03
207 -3.948385e+04 1.128268e+04
208 -4.716405e+04 -3.948385e+04
209 3.600235e+04 -4.716405e+04
210 7.816804e+04 3.600235e+04
211 7.713665e+04 7.816804e+04
212 2.906961e+04 7.713665e+04
213 5.044902e+04 2.906961e+04
214 7.516457e+04 5.044902e+04
215 2.696010e+04 7.516457e+04
216 2.293713e+04 2.696010e+04
217 3.299888e+04 2.293713e+04
218 -1.654379e+04 3.299888e+04
219 5.158962e+04 -1.654379e+04
220 7.816804e+04 5.158962e+04
221 7.663229e+04 7.816804e+04
222 -3.156698e+03 7.663229e+04
223 3.432896e+04 -3.156698e+03
224 -1.929492e+03 3.432896e+04
225 2.088802e+04 -1.929492e+03
226 7.816804e+04 2.088802e+04
227 6.352101e+04 7.816804e+04
228 3.600235e+04 6.352101e+04
229 4.264214e+04 3.600235e+04
230 6.133615e+04 4.264214e+04
231 7.816804e+04 6.133615e+04
232 -1.624027e+03 7.816804e+04
233 3.600235e+04 -1.624027e+03
234 3.605292e+04 3.600235e+04
235 3.600235e+04 3.605292e+04
236 2.723941e+04 3.600235e+04
237 7.816804e+04 2.723941e+04
238 7.816804e+04 7.816804e+04
239 3.600235e+04 7.816804e+04
240 6.911561e+04 3.600235e+04
241 3.600235e+04 6.911561e+04
242 6.883354e+04 3.600235e+04
243 3.852101e+04 6.883354e+04
244 1.128268e+04 3.852101e+04
245 1.128268e+04 1.128268e+04
246 3.956546e+03 1.128268e+04
247 1.128268e+04 3.956546e+03
248 -3.088301e+04 1.128268e+04
249 5.369742e+04 -3.088301e+04
250 1.128268e+04 5.369742e+04
251 6.171967e+04 1.128268e+04
252 -5.658248e+04 6.171967e+04
253 -1.706214e+05 -5.658248e+04
254 3.600235e+04 -1.706214e+05
255 2.306742e+04 3.600235e+04
256 -6.896001e+03 2.306742e+04
257 1.128268e+04 -6.896001e+03
258 7.816804e+04 1.128268e+04
259 1.128268e+04 7.816804e+04
260 1.128268e+04 1.128268e+04
261 2.574404e+04 1.128268e+04
262 7.838025e+04 2.574404e+04
263 1.128268e+04 7.838025e+04
264 7.816804e+04 1.128268e+04
265 -3.088301e+04 7.816804e+04
266 1.051092e+05 -3.088301e+04
267 1.128268e+04 1.051092e+05
268 4.108224e+04 1.128268e+04
269 1.128268e+04 4.108224e+04
270 1.128268e+04 1.128268e+04
271 1.128268e+04 1.128268e+04
272 6.072985e+02 1.128268e+04
273 3.600235e+04 6.072985e+02
274 5.887555e+04 3.600235e+04
275 -2.433108e+03 5.887555e+04
276 7.816804e+04 -2.433108e+03
277 -8.162836e+04 7.816804e+04
278 4.916413e+04 -8.162836e+04
279 1.128268e+04 4.916413e+04
280 3.600235e+04 1.128268e+04
281 3.950952e+04 3.600235e+04
282 3.600235e+04 3.950952e+04
283 6.270590e+04 3.600235e+04
284 7.816804e+04 6.270590e+04
285 2.043001e+01 7.816804e+04
286 3.203934e+04 2.043001e+01
287 -3.088301e+04 3.203934e+04
288 -6.111311e+02 -3.088301e+04
289 1.128268e+04 -6.111311e+02
290 2.281855e+04 1.128268e+04
291 6.682024e+04 2.281855e+04
292 3.600235e+04 6.682024e+04
293 5.115761e+04 3.600235e+04
294 2.037284e+04 5.115761e+04
295 3.600235e+04 2.037284e+04
296 -3.088301e+04 3.600235e+04
297 7.816804e+04 -3.088301e+04
298 2.248674e+04 7.816804e+04
299 2.575388e+03 2.248674e+04
300 6.140890e+04 2.575388e+03
301 7.666630e+04 6.140890e+04
302 1.128268e+04 7.666630e+04
303 3.792899e+04 1.128268e+04
304 1.794995e+04 3.792899e+04
305 1.128268e+04 1.794995e+04
306 -2.309127e+04 1.128268e+04
307 8.272152e+03 -2.309127e+04
308 1.128268e+04 8.272152e+03
309 1.128268e+04 1.128268e+04
310 2.083068e+04 1.128268e+04
311 1.128268e+04 2.083068e+04
312 -5.046879e+04 1.128268e+04
313 2.417100e+04 -5.046879e+04
314 4.523385e+04 2.417100e+04
315 3.617547e+04 4.523385e+04
316 -1.752119e+04 3.617547e+04
317 2.847441e+04 -1.752119e+04
318 3.600235e+04 2.847441e+04
319 7.816804e+04 3.600235e+04
320 7.816804e+04 7.816804e+04
321 7.816804e+04 7.816804e+04
322 -3.088301e+04 7.816804e+04
323 7.816804e+04 -3.088301e+04
324 7.350799e+04 7.816804e+04
325 -3.088301e+04 7.350799e+04
326 1.128268e+04 -3.088301e+04
327 7.816804e+04 1.128268e+04
328 7.816804e+04 7.816804e+04
329 7.816804e+04 7.816804e+04
330 3.600235e+04 7.816804e+04
331 7.816804e+04 3.600235e+04
332 6.976128e+04 7.816804e+04
333 3.448651e+04 6.976128e+04
334 2.602428e+04 3.448651e+04
335 1.796446e+04 2.602428e+04
336 3.600235e+04 1.796446e+04
337 3.243381e+04 3.600235e+04
338 6.522883e+04 3.243381e+04
339 -8.705010e+03 6.522883e+04
340 3.600235e+04 -8.705010e+03
341 3.600235e+04 3.600235e+04
342 7.816804e+04 3.600235e+04
343 1.519857e+03 7.816804e+04
344 3.600235e+04 1.519857e+03
345 2.046556e+04 3.600235e+04
346 3.600235e+04 2.046556e+04
347 7.816804e+04 3.600235e+04
348 4.847882e+03 7.816804e+04
349 7.591426e+04 4.847882e+03
350 2.285412e+04 7.591426e+04
351 6.758874e+04 2.285412e+04
352 2.980277e+04 6.758874e+04
353 3.600235e+04 2.980277e+04
354 3.600235e+04 3.600235e+04
355 -1.050841e+04 3.600235e+04
356 3.600235e+04 -1.050841e+04
357 3.600235e+04 3.600235e+04
358 7.816804e+04 3.600235e+04
359 2.999808e+04 7.816804e+04
360 7.816804e+04 2.999808e+04
361 3.600235e+04 7.816804e+04
362 -6.124955e+04 3.600235e+04
363 6.485885e+04 -6.124955e+04
364 7.816804e+04 6.485885e+04
365 7.816804e+04 7.816804e+04
366 3.600235e+04 7.816804e+04
367 7.816804e+04 3.600235e+04
368 2.291243e+04 7.816804e+04
369 7.816804e+04 2.291243e+04
370 7.125759e+04 7.816804e+04
371 2.228522e+04 7.125759e+04
372 5.794993e+04 2.228522e+04
373 7.816804e+04 5.794993e+04
374 7.816804e+04 7.816804e+04
375 4.627300e+04 7.816804e+04
376 2.517615e+04 4.627300e+04
377 -2.442581e+03 2.517615e+04
378 7.816804e+04 -2.442581e+03
379 7.816804e+04 7.816804e+04
380 -1.304973e+04 7.816804e+04
381 -1.107369e+03 -1.304973e+04
382 7.816804e+04 -1.107369e+03
383 1.626159e+04 7.816804e+04
384 -8.365301e+04 1.626159e+04
385 4.297161e+04 -8.365301e+04
386 5.401746e+04 4.297161e+04
387 2.985540e+04 5.401746e+04
388 5.438925e+04 2.985540e+04
389 6.238896e+04 5.438925e+04
390 -3.806724e+01 6.238896e+04
391 5.889624e+03 -3.806724e+01
392 6.055375e+04 5.889624e+03
393 -1.501482e+05 6.055375e+04
394 3.034063e+04 -1.501482e+05
395 -8.392311e+04 3.034063e+04
396 -3.247143e+05 -8.392311e+04
397 3.758911e+04 -3.247143e+05
398 -1.079982e+05 3.758911e+04
399 -6.169317e+03 -1.079982e+05
400 2.756067e+04 -6.169317e+03
401 2.896510e+03 2.756067e+04
402 -1.297722e+05 2.896510e+03
403 -2.527709e+04 -1.297722e+05
404 2.666550e+03 -2.527709e+04
405 4.559749e+04 2.666550e+03
406 3.927075e+04 4.559749e+04
407 2.931348e+03 3.927075e+04
408 -3.157747e+04 2.931348e+03
409 -4.223222e+04 -3.157747e+04
410 -1.457241e+05 -4.223222e+04
411 -6.694617e+04 -1.457241e+05
412 3.464206e+04 -6.694617e+04
413 -1.379192e+05 3.464206e+04
414 -1.075542e+05 -1.379192e+05
415 1.288701e+05 -1.075542e+05
416 -1.022662e+05 1.288701e+05
417 -6.764406e+04 -1.022662e+05
418 -4.973243e+04 -6.764406e+04
419 -5.246928e+04 -4.973243e+04
420 -1.595566e+05 -5.246928e+04
421 -5.328526e+05 -1.595566e+05
422 -3.851035e+05 -5.328526e+05
423 -1.557869e+04 -3.851035e+05
424 -2.849142e+05 -1.557869e+04
425 -2.636122e+05 -2.849142e+05
426 -8.778464e+05 -2.636122e+05
427 -1.739617e+05 -8.778464e+05
428 -2.086613e+05 -1.739617e+05
429 -5.306982e+05 -2.086613e+05
> plot(z,main=paste('Residual Lag plot, lowess, and regression line'), ylab='values of Residuals', xlab='lagged values of Residuals')
> lines(lowess(z))
> abline(lm(z))
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/7p2wq1292932180.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> acf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Autocorrelation Function')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/8zbds1292932180.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> pacf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Partial Autocorrelation Function')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/9zbds1292932180.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> opar <- par(mfrow = c(2,2), oma = c(0, 0, 1.1, 0))
> plot(mylm, las = 1, sub='Residual Diagnostics')
> par(opar)
> dev.off()
null device
1
> if (n > n25) {
+ postscript(file="/var/www/html/rcomp/tmp/10zbds1292932180.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
+ plot(kp3:nmkm3,gqarr[,2], main='Goldfeld-Quandt test',ylab='2-sided p-value',xlab='breakpoint')
+ grid()
+ dev.off()
+ }
null device
1
>
> #Note: the /var/www/html/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/www/html/rcomp/createtable")
>
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Estimated Regression Equation', 1, TRUE)
> a<-table.row.end(a)
> myeq <- colnames(x)[1]
> myeq <- paste(myeq, '[t] = ', sep='')
> for (i in 1:k){
+ if (mysum$coefficients[i,1] > 0) myeq <- paste(myeq, '+', '')
+ myeq <- paste(myeq, mysum$coefficients[i,1], sep=' ')
+ if (rownames(mysum$coefficients)[i] != '(Intercept)') {
+ myeq <- paste(myeq, rownames(mysum$coefficients)[i], sep='')
+ if (rownames(mysum$coefficients)[i] != 't') myeq <- paste(myeq, '[t]', sep='')
+ }
+ }
> myeq <- paste(myeq, ' + e[t]')
> a<-table.row.start(a)
> a<-table.element(a, myeq)
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/11lccg1292932180.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,hyperlink('http://www.xycoon.com/ols1.htm','Multiple Linear Regression - Ordinary Least Squares',''), 6, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Variable',header=TRUE)
> a<-table.element(a,'Parameter',header=TRUE)
> a<-table.element(a,'S.D.',header=TRUE)
> a<-table.element(a,'T-STAT
H0: parameter = 0',header=TRUE)
> a<-table.element(a,'2-tail p-value',header=TRUE)
> a<-table.element(a,'1-tail p-value',header=TRUE)
> a<-table.row.end(a)
> for (i in 1:k){
+ a<-table.row.start(a)
+ a<-table.element(a,rownames(mysum$coefficients)[i],header=TRUE)
+ a<-table.element(a,mysum$coefficients[i,1])
+ a<-table.element(a, round(mysum$coefficients[i,2],6))
+ a<-table.element(a, round(mysum$coefficients[i,3],4))
+ a<-table.element(a, round(mysum$coefficients[i,4],6))
+ a<-table.element(a, round(mysum$coefficients[i,4]/2,6))
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/12oca41292932180.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Regression Statistics', 2, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple R',1,TRUE)
> a<-table.element(a, sqrt(mysum$r.squared))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'R-squared',1,TRUE)
> a<-table.element(a, mysum$r.squared)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Adjusted R-squared',1,TRUE)
> a<-table.element(a, mysum$adj.r.squared)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (value)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[1])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (DF numerator)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[2])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (DF denominator)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[3])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'p-value',1,TRUE)
> a<-table.element(a, 1-pf(mysum$fstatistic[1],mysum$fstatistic[2],mysum$fstatistic[3]))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Residual Statistics', 2, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Residual Standard Deviation',1,TRUE)
> a<-table.element(a, mysum$sigma)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Sum Squared Residuals',1,TRUE)
> a<-table.element(a, sum(myerror*myerror))
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/13dvpg1292932180.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Actuals, Interpolation, and Residuals', 4, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Time or Index', 1, TRUE)
> a<-table.element(a, 'Actuals', 1, TRUE)
> a<-table.element(a, 'Interpolation
Forecast', 1, TRUE)
> a<-table.element(a, 'Residuals
Prediction Error', 1, TRUE)
> a<-table.row.end(a)
> for (i in 1:n) {
+ a<-table.row.start(a)
+ a<-table.element(a,i, 1, TRUE)
+ a<-table.element(a,x[i])
+ a<-table.element(a,x[i]-mysum$resid[i])
+ a<-table.element(a,mysum$resid[i])
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/1464o11292932180.tab")
> if (n > n25) {
+ a<-table.start()
+ a<-table.row.start(a)
+ a<-table.element(a,'Goldfeld-Quandt test for Heteroskedasticity',4,TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'p-values',header=TRUE)
+ a<-table.element(a,'Alternative Hypothesis',3,header=TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'breakpoint index',header=TRUE)
+ a<-table.element(a,'greater',header=TRUE)
+ a<-table.element(a,'2-sided',header=TRUE)
+ a<-table.element(a,'less',header=TRUE)
+ a<-table.row.end(a)
+ for (mypoint in kp3:nmkm3) {
+ a<-table.row.start(a)
+ a<-table.element(a,mypoint,header=TRUE)
+ a<-table.element(a,gqarr[mypoint-kp3+1,1])
+ a<-table.element(a,gqarr[mypoint-kp3+1,2])
+ a<-table.element(a,gqarr[mypoint-kp3+1,3])
+ a<-table.row.end(a)
+ }
+ a<-table.end(a)
+ table.save(a,file="/var/www/html/rcomp/tmp/15kema1292932180.tab")
+ a<-table.start()
+ a<-table.row.start(a)
+ a<-table.element(a,'Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity',4,TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'Description',header=TRUE)
+ a<-table.element(a,'# significant tests',header=TRUE)
+ a<-table.element(a,'% significant tests',header=TRUE)
+ a<-table.element(a,'OK/NOK',header=TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'1% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant1)
+ a<-table.element(a,numsignificant1/numgqtests)
+ if (numsignificant1/numgqtests < 0.01) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'5% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant5)
+ a<-table.element(a,numsignificant5/numgqtests)
+ if (numsignificant5/numgqtests < 0.05) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'10% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant10)
+ a<-table.element(a,numsignificant10/numgqtests)
+ if (numsignificant10/numgqtests < 0.1) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.end(a)
+ table.save(a,file="/var/www/html/rcomp/tmp/16ne3x1292932180.tab")
+ }
>
> try(system("convert tmp/13jx21292932180.ps tmp/13jx21292932180.png",intern=TRUE))
character(0)
> try(system("convert tmp/23jx21292932180.ps tmp/23jx21292932180.png",intern=TRUE))
character(0)
> try(system("convert tmp/33jx21292932180.ps tmp/33jx21292932180.png",intern=TRUE))
character(0)
> try(system("convert tmp/4esxn1292932180.ps tmp/4esxn1292932180.png",intern=TRUE))
character(0)
> try(system("convert tmp/5esxn1292932180.ps tmp/5esxn1292932180.png",intern=TRUE))
character(0)
> try(system("convert tmp/6esxn1292932180.ps tmp/6esxn1292932180.png",intern=TRUE))
character(0)
> try(system("convert tmp/7p2wq1292932180.ps tmp/7p2wq1292932180.png",intern=TRUE))
character(0)
> try(system("convert tmp/8zbds1292932180.ps tmp/8zbds1292932180.png",intern=TRUE))
character(0)
> try(system("convert tmp/9zbds1292932180.ps tmp/9zbds1292932180.png",intern=TRUE))
character(0)
> try(system("convert tmp/10zbds1292932180.ps tmp/10zbds1292932180.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
12.987 2.234 31.619