R version 2.9.0 (2009-04-17)
Copyright (C) 2009 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> x <- array(list(1
+ ,162556
+ ,162556
+ ,807
+ ,807
+ ,213118
+ ,213118
+ ,6282154
+ ,1
+ ,29790
+ ,29790
+ ,444
+ ,444
+ ,81767
+ ,81767
+ ,4321023
+ ,1
+ ,87550
+ ,87550
+ ,412
+ ,412
+ ,153198
+ ,153198
+ ,4111912
+ ,0
+ ,84738
+ ,0
+ ,428
+ ,0
+ ,-26007
+ ,0
+ ,223193
+ ,1
+ ,54660
+ ,54660
+ ,315
+ ,315
+ ,126942
+ ,126942
+ ,1491348
+ ,1
+ ,42634
+ ,42634
+ ,168
+ ,168
+ ,157214
+ ,157214
+ ,1629616
+ ,0
+ ,40949
+ ,0
+ ,263
+ ,0
+ ,129352
+ ,0
+ ,1398893
+ ,1
+ ,45187
+ ,45187
+ ,267
+ ,267
+ ,234817
+ ,234817
+ ,1926517
+ ,1
+ ,37704
+ ,37704
+ ,228
+ ,228
+ ,60448
+ ,60448
+ ,983660
+ ,1
+ ,16275
+ ,16275
+ ,129
+ ,129
+ ,47818
+ ,47818
+ ,1443586
+ ,0
+ ,25830
+ ,0
+ ,104
+ ,0
+ ,245546
+ ,0
+ ,1073089
+ ,0
+ ,12679
+ ,0
+ ,122
+ ,0
+ ,48020
+ ,0
+ ,984885
+ ,1
+ ,18014
+ ,18014
+ ,393
+ ,393
+ ,-1710
+ ,-1710
+ ,1405225
+ ,0
+ ,43556
+ ,0
+ ,190
+ ,0
+ ,32648
+ ,0
+ ,227132
+ ,1
+ ,24811
+ ,24811
+ ,280
+ ,280
+ ,95350
+ ,95350
+ ,929118
+ ,0
+ ,6575
+ ,0
+ ,63
+ ,0
+ ,151352
+ ,0
+ ,1071292
+ ,0
+ ,7123
+ ,0
+ ,102
+ ,0
+ ,288170
+ ,0
+ ,638830
+ ,1
+ ,21950
+ ,21950
+ ,265
+ ,265
+ ,114337
+ ,114337
+ ,856956
+ ,1
+ ,37597
+ ,37597
+ ,234
+ ,234
+ ,37884
+ ,37884
+ ,992426
+ ,0
+ ,17821
+ ,0
+ ,277
+ ,0
+ ,122844
+ ,0
+ ,444477
+ ,1
+ ,12988
+ ,12988
+ ,73
+ ,73
+ ,82340
+ ,82340
+ ,857217
+ ,1
+ ,22330
+ ,22330
+ ,67
+ ,67
+ ,79801
+ ,79801
+ ,711969
+ ,0
+ ,13326
+ ,0
+ ,103
+ ,0
+ ,165548
+ ,0
+ ,702380
+ ,0
+ ,16189
+ ,0
+ ,290
+ ,0
+ ,116384
+ ,0
+ ,358589
+ ,0
+ ,7146
+ ,0
+ ,83
+ ,0
+ ,134028
+ ,0
+ ,297978
+ ,0
+ ,15824
+ ,0
+ ,56
+ ,0
+ ,63838
+ ,0
+ ,585715
+ ,1
+ ,27664
+ ,27664
+ ,236
+ ,236
+ ,74996
+ ,74996
+ ,657954
+ ,0
+ ,11920
+ ,0
+ ,73
+ ,0
+ ,31080
+ ,0
+ ,209458
+ ,0
+ ,8568
+ ,0
+ ,34
+ ,0
+ ,32168
+ ,0
+ ,786690
+ ,0
+ ,14416
+ ,0
+ ,139
+ ,0
+ ,49857
+ ,0
+ ,439798
+ ,1
+ ,3369
+ ,3369
+ ,26
+ ,26
+ ,87161
+ ,87161
+ ,688779
+ ,1
+ ,11819
+ ,11819
+ ,70
+ ,70
+ ,106113
+ ,106113
+ ,574339
+ ,1
+ ,6984
+ ,6984
+ ,40
+ ,40
+ ,80570
+ ,80570
+ ,741409
+ ,1
+ ,4519
+ ,4519
+ ,42
+ ,42
+ ,102129
+ ,102129
+ ,597793
+ ,0
+ ,2220
+ ,0
+ ,12
+ ,0
+ ,301670
+ ,0
+ ,644190
+ ,0
+ ,18562
+ ,0
+ ,211
+ ,0
+ ,102313
+ ,0
+ ,377934
+ ,0
+ ,10327
+ ,0
+ ,74
+ ,0
+ ,88577
+ ,0
+ ,640273
+ ,1
+ ,5336
+ ,5336
+ ,80
+ ,80
+ ,112477
+ ,112477
+ ,697458
+ ,1
+ ,2365
+ ,2365
+ ,83
+ ,83
+ ,191778
+ ,191778
+ ,550608
+ ,0
+ ,4069
+ ,0
+ ,131
+ ,0
+ ,79804
+ ,0
+ ,207393
+ ,0
+ ,8636
+ ,0
+ ,203
+ ,0
+ ,128294
+ ,0
+ ,301607
+ ,0
+ ,13718
+ ,0
+ ,56
+ ,0
+ ,96448
+ ,0
+ ,345783
+ ,0
+ ,4525
+ ,0
+ ,89
+ ,0
+ ,93811
+ ,0
+ ,501749
+ ,0
+ ,6869
+ ,0
+ ,88
+ ,0
+ ,117520
+ ,0
+ ,379983
+ ,0
+ ,4628
+ ,0
+ ,39
+ ,0
+ ,69159
+ ,0
+ ,387475
+ ,1
+ ,3689
+ ,3689
+ ,25
+ ,25
+ ,101792
+ ,101792
+ ,377305
+ ,1
+ ,4891
+ ,4891
+ ,49
+ ,49
+ ,210568
+ ,210568
+ ,370837
+ ,1
+ ,7489
+ ,7489
+ ,149
+ ,149
+ ,136996
+ ,136996
+ ,430866
+ ,0
+ ,4901
+ ,0
+ ,58
+ ,0
+ ,121920
+ ,0
+ ,469107
+ ,0
+ ,2284
+ ,0
+ ,41
+ ,0
+ ,76403
+ ,0
+ ,194493
+ ,1
+ ,3160
+ ,3160
+ ,90
+ ,90
+ ,108094
+ ,108094
+ ,530670
+ ,1
+ ,4150
+ ,4150
+ ,136
+ ,136
+ ,134759
+ ,134759
+ ,518365
+ ,1
+ ,7285
+ ,7285
+ ,97
+ ,97
+ ,188873
+ ,188873
+ ,491303
+ ,1
+ ,1134
+ ,1134
+ ,63
+ ,63
+ ,146216
+ ,146216
+ ,527021
+ ,1
+ ,4658
+ ,4658
+ ,114
+ ,114
+ ,156608
+ ,156608
+ ,233773
+ ,0
+ ,2384
+ ,0
+ ,77
+ ,0
+ ,61348
+ ,0
+ ,405972
+ ,0
+ ,3748
+ ,0
+ ,6
+ ,0
+ ,50350
+ ,0
+ ,652925
+ ,0
+ ,5371
+ ,0
+ ,47
+ ,0
+ ,87720
+ ,0
+ ,446211
+ ,0
+ ,1285
+ ,0
+ ,51
+ ,0
+ ,99489
+ ,0
+ ,341340
+ ,1
+ ,9327
+ ,9327
+ ,85
+ ,85
+ ,87419
+ ,87419
+ ,387699
+ ,1
+ ,5565
+ ,5565
+ ,43
+ ,43
+ ,94355
+ ,94355
+ ,493408
+ ,0
+ ,1528
+ ,0
+ ,32
+ ,0
+ ,60326
+ ,0
+ ,146494
+ ,1
+ ,3122
+ ,3122
+ ,25
+ ,25
+ ,94670
+ ,94670
+ ,414462
+ ,1
+ ,7561
+ ,7561
+ ,77
+ ,77
+ ,82425
+ ,82425
+ ,364304
+ ,0
+ ,2675
+ ,0
+ ,54
+ ,0
+ ,59017
+ ,0
+ ,355178
+ ,0
+ ,13253
+ ,0
+ ,251
+ ,0
+ ,90829
+ ,0
+ ,357760
+ ,0
+ ,880
+ ,0
+ ,15
+ ,0
+ ,80791
+ ,0
+ ,261216
+ ,1
+ ,2053
+ ,2053
+ ,44
+ ,44
+ ,100423
+ ,100423
+ ,397144
+ ,0
+ ,1424
+ ,0
+ ,73
+ ,0
+ ,131116
+ ,0
+ ,374943
+ ,1
+ ,4036
+ ,4036
+ ,85
+ ,85
+ ,100269
+ ,100269
+ ,424898
+ ,1
+ ,3045
+ ,3045
+ ,49
+ ,49
+ ,27330
+ ,27330
+ ,202055
+ ,0
+ ,5119
+ ,0
+ ,38
+ ,0
+ ,39039
+ ,0
+ ,378525
+ ,0
+ ,1431
+ ,0
+ ,35
+ ,0
+ ,106885
+ ,0
+ ,310768
+ ,0
+ ,554
+ ,0
+ ,9
+ ,0
+ ,79285
+ ,0
+ ,325738
+ ,0
+ ,1975
+ ,0
+ ,34
+ ,0
+ ,118881
+ ,0
+ ,394510
+ ,1
+ ,1765
+ ,1765
+ ,20
+ ,20
+ ,77623
+ ,77623
+ ,247060
+ ,0
+ ,1012
+ ,0
+ ,29
+ ,0
+ ,114768
+ ,0
+ ,368078
+ ,0
+ ,810
+ ,0
+ ,11
+ ,0
+ ,74015
+ ,0
+ ,236761
+ ,0
+ ,1280
+ ,0
+ ,52
+ ,0
+ ,69465
+ ,0
+ ,312378
+ ,1
+ ,666
+ ,666
+ ,13
+ ,13
+ ,117869
+ ,117869
+ ,339836
+ ,0
+ ,1380
+ ,0
+ ,29
+ ,0
+ ,60982
+ ,0
+ ,347385
+ ,1
+ ,4677
+ ,4677
+ ,66
+ ,66
+ ,90131
+ ,90131
+ ,426280
+ ,0
+ ,876
+ ,0
+ ,33
+ ,0
+ ,138971
+ ,0
+ ,352850
+ ,0
+ ,814
+ ,0
+ ,15
+ ,0
+ ,39625
+ ,0
+ ,301881
+ ,0
+ ,514
+ ,0
+ ,15
+ ,0
+ ,102725
+ ,0
+ ,377516
+ ,1
+ ,5692
+ ,5692
+ ,68
+ ,68
+ ,64239
+ ,64239
+ ,357312
+ ,0
+ ,3642
+ ,0
+ ,100
+ ,0
+ ,90262
+ ,0
+ ,458343
+ ,0
+ ,540
+ ,0
+ ,13
+ ,0
+ ,103960
+ ,0
+ ,354228
+ ,0
+ ,2099
+ ,0
+ ,45
+ ,0
+ ,106611
+ ,0
+ ,308636
+ ,0
+ ,567
+ ,0
+ ,14
+ ,0
+ ,103345
+ ,0
+ ,386212
+ ,0
+ ,2001
+ ,0
+ ,36
+ ,0
+ ,95551
+ ,0
+ ,393343
+ ,1
+ ,2949
+ ,2949
+ ,40
+ ,40
+ ,82903
+ ,82903
+ ,378509
+ ,0
+ ,2253
+ ,0
+ ,68
+ ,0
+ ,63593
+ ,0
+ ,452469
+ ,1
+ ,6533
+ ,6533
+ ,29
+ ,29
+ ,126910
+ ,126910
+ ,364839
+ ,0
+ ,1889
+ ,0
+ ,43
+ ,0
+ ,37527
+ ,0
+ ,358649
+ ,1
+ ,3055
+ ,3055
+ ,30
+ ,30
+ ,60247
+ ,60247
+ ,376641
+ ,0
+ ,272
+ ,0
+ ,9
+ ,0
+ ,112995
+ ,0
+ ,429112
+ ,1
+ ,1414
+ ,1414
+ ,22
+ ,22
+ ,70184
+ ,70184
+ ,330546
+ ,0
+ ,2564
+ ,0
+ ,19
+ ,0
+ ,130140
+ ,0
+ ,403560
+ ,1
+ ,1383
+ ,1383
+ ,9
+ ,9
+ ,73221
+ ,73221
+ ,317892
+ ,1
+ ,1261
+ ,1261
+ ,31
+ ,31
+ ,76114
+ ,76114
+ ,307528
+ ,0
+ ,975
+ ,0
+ ,19
+ ,0
+ ,90534
+ ,0
+ ,235133
+ ,0
+ ,3366
+ ,0
+ ,55
+ ,0
+ ,108479
+ ,0
+ ,299243
+ ,0
+ ,576
+ ,0
+ ,8
+ ,0
+ ,113761
+ ,0
+ ,314073
+ ,0
+ ,1686
+ ,0
+ ,28
+ ,0
+ ,68696
+ ,0
+ ,368186
+ ,0
+ ,746
+ ,0
+ ,29
+ ,0
+ ,71561
+ ,0
+ ,269661
+ ,1
+ ,3192
+ ,3192
+ ,48
+ ,48
+ ,59831
+ ,59831
+ ,125390
+ ,1
+ ,2045
+ ,2045
+ ,16
+ ,16
+ ,97890
+ ,97890
+ ,510834
+ ,0
+ ,5702
+ ,0
+ ,47
+ ,0
+ ,101481
+ ,0
+ ,321896
+ ,1
+ ,1932
+ ,1932
+ ,20
+ ,20
+ ,72954
+ ,72954
+ ,249898
+ ,0
+ ,936
+ ,0
+ ,22
+ ,0
+ ,67939
+ ,0
+ ,408881
+ ,1
+ ,3437
+ ,3437
+ ,33
+ ,33
+ ,48022
+ ,48022
+ ,158492
+ ,0
+ ,5131
+ ,0
+ ,44
+ ,0
+ ,86111
+ ,0
+ ,292154
+ ,1
+ ,2397
+ ,2397
+ ,13
+ ,13
+ ,74020
+ ,74020
+ ,289513
+ ,1
+ ,1389
+ ,1389
+ ,6
+ ,6
+ ,57530
+ ,57530
+ ,378049
+ ,0
+ ,1503
+ ,0
+ ,35
+ ,0
+ ,56364
+ ,0
+ ,343466
+ ,0
+ ,402
+ ,0
+ ,8
+ ,0
+ ,84990
+ ,0
+ ,332743
+ ,0
+ ,2239
+ ,0
+ ,17
+ ,0
+ ,88590
+ ,0
+ ,442882
+ ,1
+ ,2234
+ ,2234
+ ,11
+ ,11
+ ,77200
+ ,77200
+ ,214215
+ ,0
+ ,837
+ ,0
+ ,21
+ ,0
+ ,61262
+ ,0
+ ,315688
+ ,0
+ ,10579
+ ,0
+ ,92
+ ,0
+ ,110309
+ ,0
+ ,375195
+ ,0
+ ,875
+ ,0
+ ,12
+ ,0
+ ,67000
+ ,0
+ ,334280
+ ,0
+ ,1585
+ ,0
+ ,112
+ ,0
+ ,93099
+ ,0
+ ,355864
+ ,1
+ ,1659
+ ,1659
+ ,25
+ ,25
+ ,107577
+ ,107577
+ ,480382
+ ,1
+ ,2647
+ ,2647
+ ,17
+ ,17
+ ,62920
+ ,62920
+ ,353058
+ ,1
+ ,3294
+ ,3294
+ ,23
+ ,23
+ ,75832
+ ,75832
+ ,217193
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,94
+ ,0
+ ,10
+ ,0
+ ,60793
+ ,0
+ ,314533
+ ,0
+ ,422
+ ,0
+ ,23
+ ,0
+ ,57935
+ ,0
+ ,318056
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,34
+ ,0
+ ,7
+ ,0
+ ,60630
+ ,0
+ ,314353
+ ,0
+ ,1558
+ ,0
+ ,25
+ ,0
+ ,55637
+ ,0
+ ,369448
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,43
+ ,0
+ ,20
+ ,0
+ ,60887
+ ,0
+ ,312846
+ ,0
+ ,645
+ ,0
+ ,4
+ ,0
+ ,60720
+ ,0
+ ,312075
+ ,0
+ ,316
+ ,0
+ ,4
+ ,0
+ ,60505
+ ,0
+ ,315009
+ ,0
+ ,115
+ ,0
+ ,10
+ ,0
+ ,60945
+ ,0
+ ,318903
+ ,0
+ ,5
+ ,0
+ ,1
+ ,0
+ ,60720
+ ,0
+ ,314887
+ ,0
+ ,897
+ ,0
+ ,4
+ ,0
+ ,60720
+ ,0
+ ,314913
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,389
+ ,0
+ ,8
+ ,0
+ ,58990
+ ,0
+ ,325506
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,1002
+ ,0
+ ,11
+ ,0
+ ,56750
+ ,0
+ ,298568
+ ,0
+ ,36
+ ,0
+ ,4
+ ,0
+ ,60894
+ ,0
+ ,315834
+ ,0
+ ,460
+ ,0
+ ,15
+ ,0
+ ,63346
+ ,0
+ ,329784
+ ,0
+ ,309
+ ,0
+ ,9
+ ,0
+ ,56535
+ ,0
+ ,312878
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,9
+ ,0
+ ,7
+ ,0
+ ,60835
+ ,0
+ ,314987
+ ,0
+ ,271
+ ,0
+ ,2
+ ,0
+ ,60720
+ ,0
+ ,325249
+ ,0
+ ,14
+ ,0
+ ,0
+ ,0
+ ,61016
+ ,0
+ ,315877
+ ,0
+ ,520
+ ,0
+ ,7
+ ,0
+ ,58650
+ ,0
+ ,291650
+ ,0
+ ,1766
+ ,0
+ ,46
+ ,0
+ ,60438
+ ,0
+ ,305959
+ ,0
+ ,0
+ ,0
+ ,5
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,458
+ ,0
+ ,7
+ ,0
+ ,58625
+ ,0
+ ,297765
+ ,0
+ ,20
+ ,0
+ ,2
+ ,0
+ ,60938
+ ,0
+ ,315245
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,98
+ ,0
+ ,2
+ ,0
+ ,61490
+ ,0
+ ,315236
+ ,0
+ ,405
+ ,0
+ ,5
+ ,0
+ ,60845
+ ,0
+ ,336425
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,483
+ ,0
+ ,7
+ ,0
+ ,60830
+ ,0
+ ,306268
+ ,0
+ ,454
+ ,0
+ ,24
+ ,0
+ ,63261
+ ,0
+ ,302187
+ ,0
+ ,47
+ ,0
+ ,1
+ ,0
+ ,60720
+ ,0
+ ,314882
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,757
+ ,0
+ ,18
+ ,0
+ ,45689
+ ,0
+ ,382712
+ ,0
+ ,4655
+ ,0
+ ,55
+ ,0
+ ,60720
+ ,0
+ ,341570
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,36
+ ,0
+ ,3
+ ,0
+ ,61564
+ ,0
+ ,312412
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,203
+ ,0
+ ,9
+ ,0
+ ,61938
+ ,0
+ ,309596
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,126
+ ,0
+ ,8
+ ,0
+ ,60951
+ ,0
+ ,315547
+ ,0
+ ,400
+ ,0
+ ,113
+ ,0
+ ,60720
+ ,0
+ ,313267
+ ,1
+ ,71
+ ,71
+ ,0
+ ,0
+ ,60745
+ ,60745
+ ,316176
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,972
+ ,0
+ ,19
+ ,0
+ ,71642
+ ,0
+ ,359335
+ ,1
+ ,531
+ ,531
+ ,11
+ ,11
+ ,71641
+ ,71641
+ ,330068
+ ,0
+ ,2461
+ ,0
+ ,25
+ ,0
+ ,55792
+ ,0
+ ,314289
+ ,1
+ ,378
+ ,378
+ ,16
+ ,16
+ ,71873
+ ,71873
+ ,297413
+ ,1
+ ,23
+ ,23
+ ,5
+ ,5
+ ,62555
+ ,62555
+ ,314806
+ ,1
+ ,638
+ ,638
+ ,11
+ ,11
+ ,60370
+ ,60370
+ ,333210
+ ,1
+ ,2300
+ ,2300
+ ,23
+ ,23
+ ,64873
+ ,64873
+ ,352108
+ ,0
+ ,149
+ ,0
+ ,6
+ ,0
+ ,62041
+ ,0
+ ,313332
+ ,0
+ ,226
+ ,0
+ ,5
+ ,0
+ ,65745
+ ,0
+ ,291787
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,275
+ ,0
+ ,7
+ ,0
+ ,59500
+ ,0
+ ,318745
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,141
+ ,0
+ ,7
+ ,0
+ ,61630
+ ,0
+ ,315366
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,28
+ ,0
+ ,3
+ ,0
+ ,60890
+ ,0
+ ,315688
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,4980
+ ,4980
+ ,89
+ ,89
+ ,113521
+ ,113521
+ ,409642
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,472
+ ,472
+ ,19
+ ,19
+ ,80045
+ ,80045
+ ,269587
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,203
+ ,203
+ ,12
+ ,12
+ ,50804
+ ,50804
+ ,300962
+ ,1
+ ,496
+ ,496
+ ,12
+ ,12
+ ,87390
+ ,87390
+ ,325479
+ ,1
+ ,10
+ ,10
+ ,5
+ ,5
+ ,61656
+ ,61656
+ ,316155
+ ,1
+ ,63
+ ,63
+ ,2
+ ,2
+ ,65688
+ ,65688
+ ,318574
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,1136
+ ,1136
+ ,26
+ ,26
+ ,48522
+ ,48522
+ ,343613
+ ,1
+ ,265
+ ,265
+ ,3
+ ,3
+ ,60720
+ ,60720
+ ,306948
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,267
+ ,0
+ ,11
+ ,0
+ ,57640
+ ,0
+ ,330059
+ ,0
+ ,474
+ ,0
+ ,10
+ ,0
+ ,61977
+ ,0
+ ,288985
+ ,0
+ ,534
+ ,0
+ ,5
+ ,0
+ ,62620
+ ,0
+ ,304485
+ ,0
+ ,0
+ ,0
+ ,2
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,15
+ ,0
+ ,6
+ ,0
+ ,60831
+ ,0
+ ,315688
+ ,0
+ ,397
+ ,0
+ ,7
+ ,0
+ ,60646
+ ,0
+ ,317736
+ ,0
+ ,0
+ ,0
+ ,2
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,1866
+ ,0
+ ,28
+ ,0
+ ,56225
+ ,0
+ ,322331
+ ,0
+ ,288
+ ,0
+ ,3
+ ,0
+ ,60510
+ ,0
+ ,296656
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,3
+ ,0
+ ,1
+ ,0
+ ,60698
+ ,0
+ ,315354
+ ,0
+ ,468
+ ,0
+ ,20
+ ,0
+ ,60720
+ ,0
+ ,312161
+ ,0
+ ,20
+ ,0
+ ,1
+ ,0
+ ,60805
+ ,0
+ ,315576
+ ,0
+ ,278
+ ,0
+ ,22
+ ,0
+ ,61404
+ ,0
+ ,314922
+ ,0
+ ,61
+ ,0
+ ,9
+ ,0
+ ,60720
+ ,0
+ ,314551
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,192
+ ,0
+ ,2
+ ,0
+ ,65276
+ ,0
+ ,312339
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,317
+ ,0
+ ,7
+ ,0
+ ,63915
+ ,0
+ ,298700
+ ,0
+ ,738
+ ,0
+ ,9
+ ,0
+ ,60720
+ ,0
+ ,321376
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,368
+ ,0
+ ,13
+ ,0
+ ,61686
+ ,0
+ ,303230
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,2
+ ,0
+ ,0
+ ,0
+ ,60743
+ ,0
+ ,315487
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,53
+ ,0
+ ,6
+ ,0
+ ,60349
+ ,0
+ ,315793
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,94
+ ,0
+ ,3
+ ,0
+ ,61360
+ ,0
+ ,312887
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,24
+ ,0
+ ,7
+ ,0
+ ,59818
+ ,0
+ ,315637
+ ,0
+ ,2332
+ ,0
+ ,2
+ ,0
+ ,72680
+ ,0
+ ,324385
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,131
+ ,0
+ ,15
+ ,0
+ ,61808
+ ,0
+ ,308989
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,206
+ ,0
+ ,9
+ ,0
+ ,53110
+ ,0
+ ,296702
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,167
+ ,0
+ ,1
+ ,0
+ ,64245
+ ,0
+ ,307322
+ ,0
+ ,622
+ ,0
+ ,38
+ ,0
+ ,73007
+ ,0
+ ,304376
+ ,0
+ ,2328
+ ,0
+ ,57
+ ,0
+ ,82732
+ ,0
+ ,253588
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,365
+ ,0
+ ,7
+ ,0
+ ,54820
+ ,0
+ ,309560
+ ,0
+ ,364
+ ,0
+ ,26
+ ,0
+ ,47705
+ ,0
+ ,298466
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,226
+ ,0
+ ,13
+ ,0
+ ,72835
+ ,0
+ ,343929
+ ,0
+ ,307
+ ,0
+ ,10
+ ,0
+ ,58856
+ ,0
+ ,331955
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,188
+ ,0
+ ,9
+ ,0
+ ,77655
+ ,0
+ ,381180
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,138
+ ,0
+ ,26
+ ,0
+ ,69817
+ ,0
+ ,331420
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,125
+ ,0
+ ,19
+ ,0
+ ,60798
+ ,0
+ ,310201
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,282
+ ,0
+ ,12
+ ,0
+ ,62452
+ ,0
+ ,320016
+ ,0
+ ,335
+ ,0
+ ,23
+ ,0
+ ,64175
+ ,0
+ ,320398
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,1324
+ ,1324
+ ,29
+ ,29
+ ,67440
+ ,67440
+ ,291841
+ ,0
+ ,176
+ ,0
+ ,8
+ ,0
+ ,68136
+ ,0
+ ,310670
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,249
+ ,0
+ ,26
+ ,0
+ ,56726
+ ,0
+ ,313491
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,333
+ ,0
+ ,9
+ ,0
+ ,70811
+ ,0
+ ,331323
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,601
+ ,601
+ ,5
+ ,5
+ ,60720
+ ,60720
+ ,319210
+ ,0
+ ,30
+ ,0
+ ,3
+ ,0
+ ,62045
+ ,0
+ ,318098
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,249
+ ,0
+ ,13
+ ,0
+ ,54323
+ ,0
+ ,292754
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,165
+ ,0
+ ,12
+ ,0
+ ,62841
+ ,0
+ ,325176
+ ,0
+ ,453
+ ,0
+ ,19
+ ,0
+ ,81125
+ ,0
+ ,365959
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,53
+ ,0
+ ,10
+ ,0
+ ,59506
+ ,0
+ ,302409
+ ,1
+ ,382
+ ,382
+ ,9
+ ,9
+ ,59365
+ ,59365
+ ,340968
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,9
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,30
+ ,30
+ ,4
+ ,4
+ ,60798
+ ,60798
+ ,313164
+ ,0
+ ,290
+ ,0
+ ,1
+ ,0
+ ,58790
+ ,0
+ ,301164
+ ,0
+ ,0
+ ,0
+ ,1
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,366
+ ,0
+ ,14
+ ,0
+ ,61808
+ ,0
+ ,344425
+ ,0
+ ,2
+ ,0
+ ,12
+ ,0
+ ,60735
+ ,0
+ ,315394
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,209
+ ,209
+ ,19
+ ,19
+ ,64016
+ ,64016
+ ,316647
+ ,0
+ ,384
+ ,0
+ ,17
+ ,0
+ ,54683
+ ,0
+ ,309836
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,365
+ ,0
+ ,32
+ ,0
+ ,87192
+ ,0
+ ,346611
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,49
+ ,49
+ ,14
+ ,14
+ ,64107
+ ,64107
+ ,322031
+ ,0
+ ,3
+ ,0
+ ,8
+ ,0
+ ,60761
+ ,0
+ ,315656
+ ,0
+ ,133
+ ,0
+ ,4
+ ,0
+ ,65990
+ ,0
+ ,339445
+ ,0
+ ,32
+ ,0
+ ,0
+ ,0
+ ,59988
+ ,0
+ ,314964
+ ,0
+ ,368
+ ,0
+ ,20
+ ,0
+ ,61167
+ ,0
+ ,297141
+ ,0
+ ,1
+ ,0
+ ,5
+ ,0
+ ,60719
+ ,0
+ ,315372
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,22
+ ,0
+ ,1
+ ,0
+ ,60722
+ ,0
+ ,312502
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,96
+ ,0
+ ,4
+ ,0
+ ,60379
+ ,0
+ ,313729
+ ,0
+ ,1
+ ,0
+ ,1
+ ,0
+ ,60727
+ ,0
+ ,315388
+ ,0
+ ,314
+ ,0
+ ,4
+ ,0
+ ,60720
+ ,0
+ ,315371
+ ,0
+ ,844
+ ,0
+ ,20
+ ,0
+ ,60925
+ ,0
+ ,296139
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,26
+ ,0
+ ,1
+ ,0
+ ,60896
+ ,0
+ ,313880
+ ,0
+ ,125
+ ,0
+ ,10
+ ,0
+ ,59734
+ ,0
+ ,317698
+ ,0
+ ,304
+ ,0
+ ,12
+ ,0
+ ,62969
+ ,0
+ ,295580
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,621
+ ,0
+ ,13
+ ,0
+ ,60720
+ ,0
+ ,308256
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,119
+ ,0
+ ,3
+ ,0
+ ,59118
+ ,0
+ ,303677
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,1595
+ ,0
+ ,10
+ ,0
+ ,60720
+ ,0
+ ,319369
+ ,0
+ ,312
+ ,0
+ ,3
+ ,0
+ ,58598
+ ,0
+ ,318690
+ ,0
+ ,60
+ ,0
+ ,7
+ ,0
+ ,61124
+ ,0
+ ,314049
+ ,0
+ ,587
+ ,0
+ ,10
+ ,0
+ ,59595
+ ,0
+ ,325699
+ ,0
+ ,135
+ ,0
+ ,1
+ ,0
+ ,62065
+ ,0
+ ,314210
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,514
+ ,0
+ ,15
+ ,0
+ ,78780
+ ,0
+ ,322378
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,1
+ ,0
+ ,4
+ ,0
+ ,60722
+ ,0
+ ,315398
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,1763
+ ,1763
+ ,28
+ ,28
+ ,61600
+ ,61600
+ ,308336
+ ,0
+ ,180
+ ,0
+ ,9
+ ,0
+ ,59635
+ ,0
+ ,316386
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,218
+ ,0
+ ,7
+ ,0
+ ,60720
+ ,0
+ ,315553
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,448
+ ,0
+ ,7
+ ,0
+ ,59781
+ ,0
+ ,323361
+ ,0
+ ,227
+ ,0
+ ,7
+ ,0
+ ,76644
+ ,0
+ ,336639
+ ,0
+ ,174
+ ,0
+ ,3
+ ,0
+ ,64820
+ ,0
+ ,307424
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,121
+ ,0
+ ,11
+ ,0
+ ,56178
+ ,0
+ ,295370
+ ,0
+ ,607
+ ,0
+ ,7
+ ,0
+ ,60436
+ ,0
+ ,322340
+ ,0
+ ,2212
+ ,0
+ ,10
+ ,0
+ ,60720
+ ,0
+ ,319864
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,530
+ ,0
+ ,18
+ ,0
+ ,73433
+ ,0
+ ,317291
+ ,0
+ ,571
+ ,0
+ ,14
+ ,0
+ ,41477
+ ,0
+ ,280398
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,78
+ ,0
+ ,12
+ ,0
+ ,62700
+ ,0
+ ,317330
+ ,0
+ ,2489
+ ,0
+ ,29
+ ,0
+ ,67804
+ ,0
+ ,238125
+ ,0
+ ,131
+ ,0
+ ,3
+ ,0
+ ,59661
+ ,0
+ ,327071
+ ,0
+ ,923
+ ,0
+ ,6
+ ,0
+ ,58620
+ ,0
+ ,309038
+ ,0
+ ,72
+ ,0
+ ,3
+ ,0
+ ,60398
+ ,0
+ ,314210
+ ,0
+ ,572
+ ,0
+ ,8
+ ,0
+ ,58580
+ ,0
+ ,307930
+ ,0
+ ,397
+ ,0
+ ,10
+ ,0
+ ,62710
+ ,0
+ ,322327
+ ,0
+ ,450
+ ,0
+ ,6
+ ,0
+ ,59325
+ ,0
+ ,292136
+ ,0
+ ,622
+ ,0
+ ,8
+ ,0
+ ,60950
+ ,0
+ ,263276
+ ,0
+ ,694
+ ,0
+ ,6
+ ,0
+ ,68060
+ ,0
+ ,367655
+ ,1
+ ,3425
+ ,3425
+ ,9
+ ,9
+ ,83620
+ ,83620
+ ,283910
+ ,0
+ ,562
+ ,0
+ ,8
+ ,0
+ ,58456
+ ,0
+ ,283587
+ ,0
+ ,4917
+ ,0
+ ,26
+ ,0
+ ,52811
+ ,0
+ ,243650
+ ,1
+ ,1442
+ ,1442
+ ,239
+ ,239
+ ,121173
+ ,121173
+ ,438493
+ ,0
+ ,529
+ ,0
+ ,7
+ ,0
+ ,63870
+ ,0
+ ,296261
+ ,1
+ ,2126
+ ,2126
+ ,41
+ ,41
+ ,21001
+ ,21001
+ ,230621
+ ,0
+ ,1061
+ ,0
+ ,3
+ ,0
+ ,70415
+ ,0
+ ,304252
+ ,0
+ ,776
+ ,0
+ ,8
+ ,0
+ ,64230
+ ,0
+ ,333505
+ ,0
+ ,611
+ ,0
+ ,6
+ ,0
+ ,59190
+ ,0
+ ,296919
+ ,1
+ ,1526
+ ,1526
+ ,21
+ ,21
+ ,69351
+ ,69351
+ ,278990
+ ,0
+ ,592
+ ,0
+ ,7
+ ,0
+ ,64270
+ ,0
+ ,276898
+ ,0
+ ,1182
+ ,0
+ ,11
+ ,0
+ ,70694
+ ,0
+ ,327007
+ ,0
+ ,621
+ ,0
+ ,11
+ ,0
+ ,68005
+ ,0
+ ,317046
+ ,0
+ ,989
+ ,0
+ ,12
+ ,0
+ ,58930
+ ,0
+ ,304555
+ ,0
+ ,438
+ ,0
+ ,9
+ ,0
+ ,58320
+ ,0
+ ,298096
+ ,0
+ ,726
+ ,0
+ ,3
+ ,0
+ ,69980
+ ,0
+ ,231861
+ ,0
+ ,1303
+ ,0
+ ,57
+ ,0
+ ,69863
+ ,0
+ ,309422
+ ,1
+ ,7419
+ ,7419
+ ,21
+ ,21
+ ,63255
+ ,63255
+ ,286963
+ ,1
+ ,1164
+ ,1164
+ ,15
+ ,15
+ ,57320
+ ,57320
+ ,269753
+ ,1
+ ,3310
+ ,3310
+ ,32
+ ,32
+ ,75230
+ ,75230
+ ,448243
+ ,0
+ ,1920
+ ,0
+ ,11
+ ,0
+ ,79420
+ ,0
+ ,165404
+ ,0
+ ,965
+ ,0
+ ,2
+ ,0
+ ,73490
+ ,0
+ ,204325
+ ,0
+ ,3256
+ ,0
+ ,23
+ ,0
+ ,35250
+ ,0
+ ,407159
+ ,1
+ ,1135
+ ,1135
+ ,20
+ ,20
+ ,62285
+ ,62285
+ ,290476
+ ,0
+ ,1270
+ ,0
+ ,24
+ ,0
+ ,69206
+ ,0
+ ,275311
+ ,0
+ ,661
+ ,0
+ ,1
+ ,0
+ ,65920
+ ,0
+ ,246541
+ ,0
+ ,1013
+ ,0
+ ,1
+ ,0
+ ,69770
+ ,0
+ ,253468
+ ,0
+ ,2844
+ ,0
+ ,74
+ ,0
+ ,72683
+ ,0
+ ,240897
+ ,1
+ ,11528
+ ,11528
+ ,68
+ ,68
+ ,-14545
+ ,-14545
+ ,-83265
+ ,0
+ ,6526
+ ,0
+ ,20
+ ,0
+ ,55830
+ ,0
+ ,-42143
+ ,0
+ ,2264
+ ,0
+ ,20
+ ,0
+ ,55174
+ ,0
+ ,272713
+ ,1
+ ,5109
+ ,5109
+ ,82
+ ,82
+ ,67038
+ ,67038
+ ,215362
+ ,0
+ ,3999
+ ,0
+ ,21
+ ,0
+ ,51252
+ ,0
+ ,42754
+ ,0
+ ,35624
+ ,0
+ ,244
+ ,0
+ ,157278
+ ,0
+ ,306275
+ ,0
+ ,9252
+ ,0
+ ,32
+ ,0
+ ,79510
+ ,0
+ ,253537
+ ,0
+ ,15236
+ ,0
+ ,86
+ ,0
+ ,77440
+ ,0
+ ,372631
+ ,0
+ ,18073
+ ,0
+ ,69
+ ,0
+ ,27284
+ ,0
+ ,-7170)
+ ,dim=c(8
+ ,431)
+ ,dimnames=list(c('Group'
+ ,'Costs'
+ ,'Costs_g'
+ ,'Orders'
+ ,'Orders_g'
+ ,'Dividends'
+ ,'Dividends_g'
+ ,'Wealth')
+ ,1:431))
> y <- array(NA,dim=c(8,431),dimnames=list(c('Group','Costs','Costs_g','Orders','Orders_g','Dividends','Dividends_g','Wealth'),1:431))
> for (i in 1:dim(x)[1])
+ {
+ for (j in 1:dim(x)[2])
+ {
+ y[i,j] <- as.numeric(x[i,j])
+ }
+ }
> par3 = 'No Linear Trend'
> par2 = 'Do not include Seasonal Dummies'
> par1 = '8'
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Prof. Dr. P. Wessa
> #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #Technical description: Write here your technical program description (don't use hard returns!)
> library(lattice)
> library(lmtest)
Loading required package: zoo
Attaching package: 'zoo'
The following object(s) are masked from package:base :
as.Date.numeric
> n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test
> par1 <- as.numeric(par1)
> x <- t(y)
> k <- length(x[1,])
> n <- length(x[,1])
> x1 <- cbind(x[,par1], x[,1:k!=par1])
> mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1])
> colnames(x1) <- mycolnames #colnames(x)[par1]
> x <- x1
> if (par3 == 'First Differences'){
+ x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep='')))
+ for (i in 1:n-1) {
+ for (j in 1:k) {
+ x2[i,j] <- x[i+1,j] - x[i,j]
+ }
+ }
+ x <- x2
+ }
> if (par2 == 'Include Monthly Dummies'){
+ x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep ='')))
+ for (i in 1:11){
+ x2[seq(i,n,12),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> if (par2 == 'Include Quarterly Dummies'){
+ x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep ='')))
+ for (i in 1:3){
+ x2[seq(i,n,4),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> k <- length(x[1,])
> if (par3 == 'Linear Trend'){
+ x <- cbind(x, c(1:n))
+ colnames(x)[k+1] <- 't'
+ }
> x
Wealth Group Costs Costs_g Orders Orders_g Dividends Dividends_g
1 6282154 1 162556 162556 807 807 213118 213118
2 4321023 1 29790 29790 444 444 81767 81767
3 4111912 1 87550 87550 412 412 153198 153198
4 223193 0 84738 0 428 0 -26007 0
5 1491348 1 54660 54660 315 315 126942 126942
6 1629616 1 42634 42634 168 168 157214 157214
7 1398893 0 40949 0 263 0 129352 0
8 1926517 1 45187 45187 267 267 234817 234817
9 983660 1 37704 37704 228 228 60448 60448
10 1443586 1 16275 16275 129 129 47818 47818
11 1073089 0 25830 0 104 0 245546 0
12 984885 0 12679 0 122 0 48020 0
13 1405225 1 18014 18014 393 393 -1710 -1710
14 227132 0 43556 0 190 0 32648 0
15 929118 1 24811 24811 280 280 95350 95350
16 1071292 0 6575 0 63 0 151352 0
17 638830 0 7123 0 102 0 288170 0
18 856956 1 21950 21950 265 265 114337 114337
19 992426 1 37597 37597 234 234 37884 37884
20 444477 0 17821 0 277 0 122844 0
21 857217 1 12988 12988 73 73 82340 82340
22 711969 1 22330 22330 67 67 79801 79801
23 702380 0 13326 0 103 0 165548 0
24 358589 0 16189 0 290 0 116384 0
25 297978 0 7146 0 83 0 134028 0
26 585715 0 15824 0 56 0 63838 0
27 657954 1 27664 27664 236 236 74996 74996
28 209458 0 11920 0 73 0 31080 0
29 786690 0 8568 0 34 0 32168 0
30 439798 0 14416 0 139 0 49857 0
31 688779 1 3369 3369 26 26 87161 87161
32 574339 1 11819 11819 70 70 106113 106113
33 741409 1 6984 6984 40 40 80570 80570
34 597793 1 4519 4519 42 42 102129 102129
35 644190 0 2220 0 12 0 301670 0
36 377934 0 18562 0 211 0 102313 0
37 640273 0 10327 0 74 0 88577 0
38 697458 1 5336 5336 80 80 112477 112477
39 550608 1 2365 2365 83 83 191778 191778
40 207393 0 4069 0 131 0 79804 0
41 301607 0 8636 0 203 0 128294 0
42 345783 0 13718 0 56 0 96448 0
43 501749 0 4525 0 89 0 93811 0
44 379983 0 6869 0 88 0 117520 0
45 387475 0 4628 0 39 0 69159 0
46 377305 1 3689 3689 25 25 101792 101792
47 370837 1 4891 4891 49 49 210568 210568
48 430866 1 7489 7489 149 149 136996 136996
49 469107 0 4901 0 58 0 121920 0
50 194493 0 2284 0 41 0 76403 0
51 530670 1 3160 3160 90 90 108094 108094
52 518365 1 4150 4150 136 136 134759 134759
53 491303 1 7285 7285 97 97 188873 188873
54 527021 1 1134 1134 63 63 146216 146216
55 233773 1 4658 4658 114 114 156608 156608
56 405972 0 2384 0 77 0 61348 0
57 652925 0 3748 0 6 0 50350 0
58 446211 0 5371 0 47 0 87720 0
59 341340 0 1285 0 51 0 99489 0
60 387699 1 9327 9327 85 85 87419 87419
61 493408 1 5565 5565 43 43 94355 94355
62 146494 0 1528 0 32 0 60326 0
63 414462 1 3122 3122 25 25 94670 94670
64 364304 1 7561 7561 77 77 82425 82425
65 355178 0 2675 0 54 0 59017 0
66 357760 0 13253 0 251 0 90829 0
67 261216 0 880 0 15 0 80791 0
68 397144 1 2053 2053 44 44 100423 100423
69 374943 0 1424 0 73 0 131116 0
70 424898 1 4036 4036 85 85 100269 100269
71 202055 1 3045 3045 49 49 27330 27330
72 378525 0 5119 0 38 0 39039 0
73 310768 0 1431 0 35 0 106885 0
74 325738 0 554 0 9 0 79285 0
75 394510 0 1975 0 34 0 118881 0
76 247060 1 1765 1765 20 20 77623 77623
77 368078 0 1012 0 29 0 114768 0
78 236761 0 810 0 11 0 74015 0
79 312378 0 1280 0 52 0 69465 0
80 339836 1 666 666 13 13 117869 117869
81 347385 0 1380 0 29 0 60982 0
82 426280 1 4677 4677 66 66 90131 90131
83 352850 0 876 0 33 0 138971 0
84 301881 0 814 0 15 0 39625 0
85 377516 0 514 0 15 0 102725 0
86 357312 1 5692 5692 68 68 64239 64239
87 458343 0 3642 0 100 0 90262 0
88 354228 0 540 0 13 0 103960 0
89 308636 0 2099 0 45 0 106611 0
90 386212 0 567 0 14 0 103345 0
91 393343 0 2001 0 36 0 95551 0
92 378509 1 2949 2949 40 40 82903 82903
93 452469 0 2253 0 68 0 63593 0
94 364839 1 6533 6533 29 29 126910 126910
95 358649 0 1889 0 43 0 37527 0
96 376641 1 3055 3055 30 30 60247 60247
97 429112 0 272 0 9 0 112995 0
98 330546 1 1414 1414 22 22 70184 70184
99 403560 0 2564 0 19 0 130140 0
100 317892 1 1383 1383 9 9 73221 73221
101 307528 1 1261 1261 31 31 76114 76114
102 235133 0 975 0 19 0 90534 0
103 299243 0 3366 0 55 0 108479 0
104 314073 0 576 0 8 0 113761 0
105 368186 0 1686 0 28 0 68696 0
106 269661 0 746 0 29 0 71561 0
107 125390 1 3192 3192 48 48 59831 59831
108 510834 1 2045 2045 16 16 97890 97890
109 321896 0 5702 0 47 0 101481 0
110 249898 1 1932 1932 20 20 72954 72954
111 408881 0 936 0 22 0 67939 0
112 158492 1 3437 3437 33 33 48022 48022
113 292154 0 5131 0 44 0 86111 0
114 289513 1 2397 2397 13 13 74020 74020
115 378049 1 1389 1389 6 6 57530 57530
116 343466 0 1503 0 35 0 56364 0
117 332743 0 402 0 8 0 84990 0
118 442882 0 2239 0 17 0 88590 0
119 214215 1 2234 2234 11 11 77200 77200
120 315688 0 837 0 21 0 61262 0
121 375195 0 10579 0 92 0 110309 0
122 334280 0 875 0 12 0 67000 0
123 355864 0 1585 0 112 0 93099 0
124 480382 1 1659 1659 25 25 107577 107577
125 353058 1 2647 2647 17 17 62920 62920
126 217193 1 3294 3294 23 23 75832 75832
127 315380 0 0 0 0 0 60720 0
128 314533 0 94 0 10 0 60793 0
129 318056 0 422 0 23 0 57935 0
130 315380 0 0 0 0 0 60720 0
131 314353 0 34 0 7 0 60630 0
132 369448 0 1558 0 25 0 55637 0
133 315380 0 0 0 0 0 60720 0
134 312846 0 43 0 20 0 60887 0
135 312075 0 645 0 4 0 60720 0
136 315009 0 316 0 4 0 60505 0
137 318903 0 115 0 10 0 60945 0
138 314887 0 5 0 1 0 60720 0
139 314913 0 897 0 4 0 60720 0
140 315380 0 0 0 0 0 60720 0
141 325506 0 389 0 8 0 58990 0
142 315380 0 0 0 0 0 60720 0
143 298568 0 1002 0 11 0 56750 0
144 315834 0 36 0 4 0 60894 0
145 329784 0 460 0 15 0 63346 0
146 312878 0 309 0 9 0 56535 0
147 315380 0 0 0 0 0 60720 0
148 314987 0 9 0 7 0 60835 0
149 325249 0 271 0 2 0 60720 0
150 315877 0 14 0 0 0 61016 0
151 291650 0 520 0 7 0 58650 0
152 305959 0 1766 0 46 0 60438 0
153 315380 0 0 0 5 0 60720 0
154 297765 0 458 0 7 0 58625 0
155 315245 0 20 0 2 0 60938 0
156 315380 0 0 0 0 0 60720 0
157 315380 0 0 0 0 0 60720 0
158 315236 0 98 0 2 0 61490 0
159 336425 0 405 0 5 0 60845 0
160 315380 0 0 0 0 0 60720 0
161 315380 0 0 0 0 0 60720 0
162 315380 0 0 0 0 0 60720 0
163 315380 0 0 0 0 0 60720 0
164 306268 0 483 0 7 0 60830 0
165 302187 0 454 0 24 0 63261 0
166 314882 0 47 0 1 0 60720 0
167 315380 0 0 0 0 0 60720 0
168 382712 0 757 0 18 0 45689 0
169 341570 0 4655 0 55 0 60720 0
170 315380 0 0 0 0 0 60720 0
171 315380 0 0 0 0 0 60720 0
172 312412 0 36 0 3 0 61564 0
173 315380 0 0 0 0 0 60720 0
174 309596 0 203 0 9 0 61938 0
175 315380 0 0 0 0 0 60720 0
176 315547 0 126 0 8 0 60951 0
177 313267 0 400 0 113 0 60720 0
178 316176 1 71 71 0 0 60745 60745
179 315380 0 0 0 0 0 60720 0
180 315380 0 0 0 0 0 60720 0
181 359335 0 972 0 19 0 71642 0
182 330068 1 531 531 11 11 71641 71641
183 314289 0 2461 0 25 0 55792 0
184 297413 1 378 378 16 16 71873 71873
185 314806 1 23 23 5 5 62555 62555
186 333210 1 638 638 11 11 60370 60370
187 352108 1 2300 2300 23 23 64873 64873
188 313332 0 149 0 6 0 62041 0
189 291787 0 226 0 5 0 65745 0
190 315380 0 0 0 0 0 60720 0
191 318745 0 275 0 7 0 59500 0
192 315380 0 0 0 0 0 60720 0
193 315366 0 141 0 7 0 61630 0
194 315380 0 0 0 0 0 60720 0
195 315688 0 28 0 3 0 60890 0
196 315380 1 0 0 0 0 60720 60720
197 409642 1 4980 4980 89 89 113521 113521
198 315380 1 0 0 0 0 60720 60720
199 315380 1 0 0 0 0 60720 60720
200 269587 1 472 472 19 19 80045 80045
201 315380 1 0 0 0 0 60720 60720
202 315380 1 0 0 0 0 60720 60720
203 315380 1 0 0 0 0 60720 60720
204 300962 1 203 203 12 12 50804 50804
205 325479 1 496 496 12 12 87390 87390
206 316155 1 10 10 5 5 61656 61656
207 318574 1 63 63 2 2 65688 65688
208 315380 1 0 0 0 0 60720 60720
209 343613 1 1136 1136 26 26 48522 48522
210 306948 1 265 265 3 3 60720 60720
211 315380 0 0 0 0 0 60720 0
212 315380 0 0 0 0 0 60720 0
213 330059 0 267 0 11 0 57640 0
214 288985 0 474 0 10 0 61977 0
215 304485 0 534 0 5 0 62620 0
216 315380 0 0 0 2 0 60720 0
217 315688 0 15 0 6 0 60831 0
218 317736 0 397 0 7 0 60646 0
219 315380 0 0 0 2 0 60720 0
220 322331 0 1866 0 28 0 56225 0
221 296656 0 288 0 3 0 60510 0
222 315380 0 0 0 0 0 60720 0
223 315354 0 3 0 1 0 60698 0
224 312161 0 468 0 20 0 60720 0
225 315576 0 20 0 1 0 60805 0
226 314922 0 278 0 22 0 61404 0
227 314551 0 61 0 9 0 60720 0
228 315380 0 0 0 0 0 60720 0
229 312339 0 192 0 2 0 65276 0
230 315380 0 0 0 0 0 60720 0
231 298700 0 317 0 7 0 63915 0
232 321376 0 738 0 9 0 60720 0
233 315380 0 0 0 0 0 60720 0
234 303230 0 368 0 13 0 61686 0
235 315380 0 0 0 0 0 60720 0
236 315487 0 2 0 0 0 60743 0
237 315380 0 0 0 0 0 60720 0
238 315793 0 53 0 6 0 60349 0
239 315380 0 0 0 0 0 60720 0
240 315380 0 0 0 0 0 60720 0
241 315380 0 0 0 0 0 60720 0
242 312887 0 94 0 3 0 61360 0
243 315380 0 0 0 0 0 60720 0
244 315637 0 24 0 7 0 59818 0
245 324385 0 2332 0 2 0 72680 0
246 315380 1 0 0 0 0 60720 60720
247 315380 1 0 0 0 0 60720 60720
248 308989 0 131 0 15 0 61808 0
249 315380 1 0 0 0 0 60720 60720
250 315380 1 0 0 0 0 60720 60720
251 296702 0 206 0 9 0 53110 0
252 315380 1 0 0 0 0 60720 60720
253 307322 0 167 0 1 0 64245 0
254 304376 0 622 0 38 0 73007 0
255 253588 0 2328 0 57 0 82732 0
256 315380 0 0 0 0 0 60720 0
257 309560 0 365 0 7 0 54820 0
258 298466 0 364 0 26 0 47705 0
259 315380 1 0 0 0 0 60720 60720
260 315380 0 0 0 0 0 60720 0
261 315380 1 0 0 0 0 60720 60720
262 315380 1 0 0 0 0 60720 60720
263 343929 0 226 0 13 0 72835 0
264 331955 0 307 0 10 0 58856 0
265 315380 1 0 0 0 0 60720 60720
266 315380 0 0 0 0 0 60720 0
267 315380 1 0 0 0 0 60720 60720
268 381180 0 188 0 9 0 77655 0
269 315380 1 0 0 0 0 60720 60720
270 331420 0 138 0 26 0 69817 0
271 315380 1 0 0 0 0 60720 60720
272 315380 1 0 0 0 0 60720 60720
273 315380 1 0 0 0 0 60720 60720
274 310201 0 125 0 19 0 60798 0
275 315380 0 0 0 0 0 60720 0
276 320016 0 282 0 12 0 62452 0
277 320398 0 335 0 23 0 64175 0
278 315380 0 0 0 0 0 60720 0
279 291841 1 1324 1324 29 29 67440 67440
280 310670 0 176 0 8 0 68136 0
281 315380 1 0 0 0 0 60720 60720
282 315380 0 0 0 0 0 60720 0
283 313491 0 249 0 26 0 56726 0
284 315380 0 0 0 0 0 60720 0
285 331323 0 333 0 9 0 70811 0
286 315380 0 0 0 0 0 60720 0
287 319210 1 601 601 5 5 60720 60720
288 318098 0 30 0 3 0 62045 0
289 315380 1 0 0 0 0 60720 60720
290 292754 0 249 0 13 0 54323 0
291 315380 1 0 0 0 0 60720 60720
292 325176 0 165 0 12 0 62841 0
293 365959 0 453 0 19 0 81125 0
294 315380 0 0 0 0 0 60720 0
295 302409 0 53 0 10 0 59506 0
296 340968 1 382 382 9 9 59365 59365
297 315380 0 0 0 0 0 60720 0
298 315380 1 0 0 0 0 60720 60720
299 315380 0 0 0 0 0 60720 0
300 315380 0 0 0 9 0 60720 0
301 313164 1 30 30 4 4 60798 60798
302 301164 0 290 0 1 0 58790 0
303 315380 0 0 0 1 0 60720 0
304 315380 1 0 0 0 0 60720 60720
305 344425 0 366 0 14 0 61808 0
306 315394 0 2 0 12 0 60735 0
307 315380 1 0 0 0 0 60720 60720
308 316647 1 209 209 19 19 64016 64016
309 309836 0 384 0 17 0 54683 0
310 315380 1 0 0 0 0 60720 60720
311 315380 1 0 0 0 0 60720 60720
312 346611 0 365 0 32 0 87192 0
313 315380 1 0 0 0 0 60720 60720
314 322031 1 49 49 14 14 64107 64107
315 315656 0 3 0 8 0 60761 0
316 339445 0 133 0 4 0 65990 0
317 314964 0 32 0 0 0 59988 0
318 297141 0 368 0 20 0 61167 0
319 315372 0 1 0 5 0 60719 0
320 315380 0 0 0 0 0 60720 0
321 315380 0 0 0 0 0 60720 0
322 315380 0 0 0 0 0 60720 0
323 315380 0 0 0 0 0 60720 0
324 315380 1 0 0 0 0 60720 60720
325 315380 0 0 0 0 0 60720 0
326 312502 0 22 0 1 0 60722 0
327 315380 1 0 0 0 0 60720 60720
328 315380 1 0 0 0 0 60720 60720
329 315380 0 0 0 0 0 60720 0
330 315380 0 0 0 0 0 60720 0
331 315380 0 0 0 0 0 60720 0
332 315380 0 0 0 0 0 60720 0
333 315380 0 0 0 0 0 60720 0
334 313729 0 96 0 4 0 60379 0
335 315388 0 1 0 1 0 60727 0
336 315371 0 314 0 4 0 60720 0
337 296139 0 844 0 20 0 60925 0
338 315380 0 0 0 0 0 60720 0
339 313880 0 26 0 1 0 60896 0
340 317698 0 125 0 10 0 59734 0
341 295580 0 304 0 12 0 62969 0
342 315380 0 0 0 0 0 60720 0
343 315380 0 0 0 0 0 60720 0
344 315380 0 0 0 0 0 60720 0
345 308256 0 621 0 13 0 60720 0
346 315380 0 0 0 0 0 60720 0
347 303677 0 119 0 3 0 59118 0
348 315380 0 0 0 0 0 60720 0
349 315380 0 0 0 0 0 60720 0
350 319369 0 1595 0 10 0 60720 0
351 318690 0 312 0 3 0 58598 0
352 314049 0 60 0 7 0 61124 0
353 325699 0 587 0 10 0 59595 0
354 314210 0 135 0 1 0 62065 0
355 315380 0 0 0 0 0 60720 0
356 315380 0 0 0 0 0 60720 0
357 322378 0 514 0 15 0 78780 0
358 315380 0 0 0 0 0 60720 0
359 315380 0 0 0 0 0 60720 0
360 315380 0 0 0 0 0 60720 0
361 315398 0 1 0 4 0 60722 0
362 315380 0 0 0 0 0 60720 0
363 315380 0 0 0 0 0 60720 0
364 308336 1 1763 1763 28 28 61600 61600
365 316386 0 180 0 9 0 59635 0
366 315380 0 0 0 0 0 60720 0
367 315380 0 0 0 0 0 60720 0
368 315380 0 0 0 0 0 60720 0
369 315380 0 0 0 0 0 60720 0
370 315553 0 218 0 7 0 60720 0
371 315380 0 0 0 0 0 60720 0
372 323361 0 448 0 7 0 59781 0
373 336639 0 227 0 7 0 76644 0
374 307424 0 174 0 3 0 64820 0
375 315380 0 0 0 0 0 60720 0
376 315380 0 0 0 0 0 60720 0
377 295370 0 121 0 11 0 56178 0
378 322340 0 607 0 7 0 60436 0
379 319864 0 2212 0 10 0 60720 0
380 315380 0 0 0 0 0 60720 0
381 315380 0 0 0 0 0 60720 0
382 317291 0 530 0 18 0 73433 0
383 280398 0 571 0 14 0 41477 0
384 315380 0 0 0 0 0 60720 0
385 317330 0 78 0 12 0 62700 0
386 238125 0 2489 0 29 0 67804 0
387 327071 0 131 0 3 0 59661 0
388 309038 0 923 0 6 0 58620 0
389 314210 0 72 0 3 0 60398 0
390 307930 0 572 0 8 0 58580 0
391 322327 0 397 0 10 0 62710 0
392 292136 0 450 0 6 0 59325 0
393 263276 0 622 0 8 0 60950 0
394 367655 0 694 0 6 0 68060 0
395 283910 1 3425 3425 9 9 83620 83620
396 283587 0 562 0 8 0 58456 0
397 243650 0 4917 0 26 0 52811 0
398 438493 1 1442 1442 239 239 121173 121173
399 296261 0 529 0 7 0 63870 0
400 230621 1 2126 2126 41 41 21001 21001
401 304252 0 1061 0 3 0 70415 0
402 333505 0 776 0 8 0 64230 0
403 296919 0 611 0 6 0 59190 0
404 278990 1 1526 1526 21 21 69351 69351
405 276898 0 592 0 7 0 64270 0
406 327007 0 1182 0 11 0 70694 0
407 317046 0 621 0 11 0 68005 0
408 304555 0 989 0 12 0 58930 0
409 298096 0 438 0 9 0 58320 0
410 231861 0 726 0 3 0 69980 0
411 309422 0 1303 0 57 0 69863 0
412 286963 1 7419 7419 21 21 63255 63255
413 269753 1 1164 1164 15 15 57320 57320
414 448243 1 3310 3310 32 32 75230 75230
415 165404 0 1920 0 11 0 79420 0
416 204325 0 965 0 2 0 73490 0
417 407159 0 3256 0 23 0 35250 0
418 290476 1 1135 1135 20 20 62285 62285
419 275311 0 1270 0 24 0 69206 0
420 246541 0 661 0 1 0 65920 0
421 253468 0 1013 0 1 0 69770 0
422 240897 0 2844 0 74 0 72683 0
423 -83265 1 11528 11528 68 68 -14545 -14545
424 -42143 0 6526 0 20 0 55830 0
425 272713 0 2264 0 20 0 55174 0
426 215362 1 5109 5109 82 82 67038 67038
427 42754 0 3999 0 21 0 51252 0
428 306275 0 35624 0 244 0 157278 0
429 253537 0 9252 0 32 0 79510 0
430 372631 0 15236 0 86 0 77440 0
431 -7170 0 18073 0 69 0 27284 0
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) Group Costs Costs_g Orders Orders_g
195390.250 13727.221 5.051 18.168 -186.494 2735.481
Dividends Dividends_g
1.899 -1.742
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-809715 -16639 4696 20133 2275663
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.954e+05 2.778e+04 7.033 8.19e-12 ***
Group 1.373e+04 4.962e+04 0.277 0.78218
Costs 5.051e+00 2.855e+00 1.769 0.07756 .
Costs_g 1.817e+01 3.481e+00 5.219 2.83e-07 ***
Orders -1.865e+02 4.150e+02 -0.449 0.65339
Orders_g 2.735e+03 5.325e+02 5.137 4.26e-07 ***
Dividends 1.899e+00 3.894e-01 4.876 1.54e-06 ***
Dividends_g -1.742e+00 6.362e-01 -2.738 0.00644 **
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 186000 on 423 degrees of freedom
Multiple R-squared: 0.8165, Adjusted R-squared: 0.8134
F-statistic: 268.8 on 7 and 423 DF, p-value: < 2.2e-16
> if (n > n25) {
+ kp3 <- k + 3
+ nmkm3 <- n - k - 3
+ gqarr <- array(NA, dim=c(nmkm3-kp3+1,3))
+ numgqtests <- 0
+ numsignificant1 <- 0
+ numsignificant5 <- 0
+ numsignificant10 <- 0
+ for (mypoint in kp3:nmkm3) {
+ j <- 0
+ numgqtests <- numgqtests + 1
+ for (myalt in c('greater', 'two.sided', 'less')) {
+ j <- j + 1
+ gqarr[mypoint-kp3+1,j] <- gqtest(mylm, point=mypoint, alternative=myalt)$p.value
+ }
+ if (gqarr[mypoint-kp3+1,2] < 0.01) numsignificant1 <- numsignificant1 + 1
+ if (gqarr[mypoint-kp3+1,2] < 0.05) numsignificant5 <- numsignificant5 + 1
+ if (gqarr[mypoint-kp3+1,2] < 0.10) numsignificant10 <- numsignificant10 + 1
+ }
+ gqarr
+ }
[,1] [,2] [,3]
[1,] 1.0000000 1.859270e-53 9.296352e-54
[2,] 1.0000000 1.002634e-57 5.013170e-58
[3,] 1.0000000 1.106020e-88 5.530099e-89
[4,] 1.0000000 8.162639e-88 4.081319e-88
[5,] 1.0000000 1.057709e-96 5.288546e-97
[6,] 1.0000000 5.741315e-105 2.870658e-105
[7,] 1.0000000 1.546406e-111 7.732028e-112
[8,] 1.0000000 5.957946e-117 2.978973e-117
[9,] 1.0000000 3.571888e-118 1.785944e-118
[10,] 1.0000000 4.811510e-119 2.405755e-119
[11,] 1.0000000 4.323433e-125 2.161716e-125
[12,] 1.0000000 1.338208e-125 6.691038e-126
[13,] 1.0000000 3.738779e-127 1.869390e-127
[14,] 1.0000000 1.765264e-126 8.826319e-127
[15,] 1.0000000 1.119176e-127 5.595880e-128
[16,] 1.0000000 3.845026e-129 1.922513e-129
[17,] 1.0000000 1.192014e-133 5.960070e-134
[18,] 1.0000000 1.646223e-133 8.231116e-134
[19,] 1.0000000 4.150782e-143 2.075391e-143
[20,] 1.0000000 3.023552e-143 1.511776e-143
[21,] 1.0000000 7.806372e-150 3.903186e-150
[22,] 1.0000000 2.756734e-150 1.378367e-150
[23,] 1.0000000 1.602737e-159 8.013686e-160
[24,] 1.0000000 7.728812e-162 3.864406e-162
[25,] 1.0000000 6.560252e-162 3.280126e-162
[26,] 1.0000000 3.388429e-161 1.694214e-161
[27,] 1.0000000 1.880456e-166 9.402279e-167
[28,] 1.0000000 8.574646e-172 4.287323e-172
[29,] 1.0000000 3.998124e-171 1.999062e-171
[30,] 1.0000000 4.767691e-172 2.383845e-172
[31,] 1.0000000 4.083583e-172 2.041791e-172
[32,] 1.0000000 5.363968e-172 2.681984e-172
[33,] 1.0000000 2.233137e-172 1.116568e-172
[34,] 1.0000000 1.304259e-171 6.521295e-172
[35,] 1.0000000 3.583430e-171 1.791715e-171
[36,] 1.0000000 2.672162e-170 1.336081e-170
[37,] 1.0000000 1.521574e-172 7.607869e-173
[38,] 1.0000000 4.469866e-173 2.234933e-173
[39,] 1.0000000 5.129913e-173 2.564956e-173
[40,] 1.0000000 7.752304e-174 3.876152e-174
[41,] 1.0000000 3.050100e-174 1.525050e-174
[42,] 1.0000000 3.621719e-174 1.810860e-174
[43,] 1.0000000 1.515167e-173 7.575833e-174
[44,] 1.0000000 8.020644e-173 4.010322e-173
[45,] 1.0000000 1.118021e-178 5.590103e-179
[46,] 1.0000000 4.265351e-178 2.132676e-178
[47,] 1.0000000 3.039374e-190 1.519687e-190
[48,] 1.0000000 6.006216e-191 3.003108e-191
[49,] 1.0000000 4.021584e-190 2.010792e-190
[50,] 1.0000000 7.905830e-190 3.952915e-190
[51,] 1.0000000 3.546664e-191 1.773332e-191
[52,] 1.0000000 1.171620e-193 5.858101e-194
[53,] 1.0000000 3.564972e-193 1.782486e-193
[54,] 1.0000000 6.993212e-193 3.496606e-193
[55,] 1.0000000 5.066030e-192 2.533015e-192
[56,] 1.0000000 4.541343e-191 2.270672e-191
[57,] 1.0000000 8.524810e-191 4.262405e-191
[58,] 1.0000000 6.712301e-190 3.356151e-190
[59,] 1.0000000 3.767170e-189 1.883585e-189
[60,] 1.0000000 1.544637e-188 7.723184e-189
[61,] 1.0000000 1.513599e-187 7.567996e-188
[62,] 1.0000000 5.227659e-188 2.613829e-188
[63,] 1.0000000 1.500655e-187 7.503275e-188
[64,] 1.0000000 1.260969e-186 6.304844e-187
[65,] 1.0000000 9.862948e-186 4.931474e-186
[66,] 1.0000000 1.174082e-185 5.870411e-186
[67,] 1.0000000 9.393846e-185 4.696923e-185
[68,] 1.0000000 9.110194e-185 4.555097e-185
[69,] 1.0000000 6.653659e-184 3.326830e-184
[70,] 1.0000000 1.162629e-184 5.813146e-185
[71,] 1.0000000 8.523358e-184 4.261679e-184
[72,] 1.0000000 8.743694e-184 4.371847e-184
[73,] 1.0000000 2.143346e-183 1.071673e-183
[74,] 1.0000000 1.858758e-182 9.293790e-183
[75,] 1.0000000 1.588619e-181 7.943096e-182
[76,] 1.0000000 4.934589e-182 2.467294e-182
[77,] 1.0000000 6.133660e-182 3.066830e-182
[78,] 1.0000000 4.970145e-181 2.485072e-181
[79,] 1.0000000 1.150320e-180 5.751601e-181
[80,] 1.0000000 9.290865e-180 4.645433e-180
[81,] 1.0000000 5.684153e-179 2.842077e-179
[82,] 1.0000000 2.290821e-178 1.145411e-178
[83,] 1.0000000 2.400047e-179 1.200023e-179
[84,] 1.0000000 1.377684e-178 6.888419e-179
[85,] 1.0000000 2.030862e-178 1.015431e-178
[86,] 1.0000000 4.713138e-179 2.356569e-179
[87,] 1.0000000 2.226320e-178 1.113160e-178
[88,] 1.0000000 1.384287e-177 6.921436e-178
[89,] 1.0000000 9.860386e-177 4.930193e-177
[90,] 1.0000000 5.385393e-176 2.692696e-176
[91,] 1.0000000 3.414243e-175 1.707121e-175
[92,] 1.0000000 5.073707e-176 2.536854e-176
[93,] 1.0000000 7.426381e-176 3.713191e-176
[94,] 1.0000000 7.198109e-176 3.599054e-176
[95,] 1.0000000 3.065713e-175 1.532857e-175
[96,] 1.0000000 7.386030e-175 3.693015e-175
[97,] 1.0000000 1.174400e-176 5.872002e-177
[98,] 1.0000000 1.926782e-178 9.633909e-179
[99,] 1.0000000 1.175430e-177 5.877149e-178
[100,] 1.0000000 2.289161e-177 1.144580e-177
[101,] 1.0000000 1.741724e-177 8.708619e-178
[102,] 1.0000000 1.886852e-177 9.434261e-178
[103,] 1.0000000 1.023256e-176 5.116282e-177
[104,] 1.0000000 6.174353e-176 3.087177e-176
[105,] 1.0000000 1.907022e-176 9.535111e-177
[106,] 1.0000000 8.889662e-176 4.444831e-176
[107,] 1.0000000 6.432145e-175 3.216072e-175
[108,] 1.0000000 2.463997e-175 1.231998e-175
[109,] 1.0000000 2.196189e-176 1.098094e-176
[110,] 1.0000000 1.951773e-175 9.758867e-176
[111,] 1.0000000 7.608143e-175 3.804071e-175
[112,] 1.0000000 6.038643e-174 3.019321e-174
[113,] 1.0000000 3.227808e-173 1.613904e-173
[114,] 1.0000000 3.111016e-173 1.555508e-173
[115,] 1.0000000 4.216492e-173 2.108246e-173
[116,] 1.0000000 1.741588e-173 8.707942e-174
[117,] 1.0000000 1.522130e-172 7.610649e-173
[118,] 1.0000000 1.342996e-171 6.714982e-172
[119,] 1.0000000 1.158537e-170 5.792685e-171
[120,] 1.0000000 9.955929e-170 4.977965e-170
[121,] 1.0000000 8.620859e-169 4.310429e-169
[122,] 1.0000000 1.013482e-168 5.067411e-169
[123,] 1.0000000 8.614615e-168 4.307308e-168
[124,] 1.0000000 7.369994e-167 3.684997e-167
[125,] 1.0000000 6.130761e-166 3.065381e-166
[126,] 1.0000000 5.085850e-165 2.542925e-165
[127,] 1.0000000 4.237830e-164 2.118915e-164
[128,] 1.0000000 3.518238e-163 1.759119e-163
[129,] 1.0000000 2.750148e-162 1.375074e-162
[130,] 1.0000000 2.249791e-161 1.124896e-161
[131,] 1.0000000 1.608721e-160 8.043607e-161
[132,] 1.0000000 1.301862e-159 6.509309e-160
[133,] 1.0000000 1.057447e-158 5.287235e-159
[134,] 1.0000000 8.516373e-158 4.258186e-158
[135,] 1.0000000 6.398076e-157 3.199038e-157
[136,] 1.0000000 4.941617e-156 2.470809e-156
[137,] 1.0000000 3.893192e-155 1.946596e-155
[138,] 1.0000000 3.092629e-154 1.546315e-154
[139,] 1.0000000 2.174795e-153 1.087398e-153
[140,] 1.0000000 1.685108e-152 8.425541e-153
[141,] 1.0000000 1.260904e-151 6.304518e-152
[142,] 1.0000000 9.802149e-151 4.901075e-151
[143,] 1.0000000 7.557074e-150 3.778537e-150
[144,] 1.0000000 5.760767e-149 2.880383e-149
[145,] 1.0000000 4.377791e-148 2.188896e-148
[146,] 1.0000000 3.291658e-147 1.645829e-147
[147,] 1.0000000 2.461976e-146 1.230988e-146
[148,] 1.0000000 1.842507e-145 9.212537e-146
[149,] 1.0000000 1.029251e-144 5.146256e-145
[150,] 1.0000000 7.581243e-144 3.790622e-144
[151,] 1.0000000 5.555053e-143 2.777526e-143
[152,] 1.0000000 4.049198e-142 2.024599e-142
[153,] 1.0000000 2.936212e-141 1.468106e-141
[154,] 1.0000000 2.153285e-140 1.076642e-140
[155,] 1.0000000 1.444617e-139 7.223084e-140
[156,] 1.0000000 1.034609e-138 5.173043e-139
[157,] 1.0000000 7.355165e-138 3.677582e-138
[158,] 1.0000000 1.576302e-138 7.881511e-139
[159,] 1.0000000 3.480937e-138 1.740469e-138
[160,] 1.0000000 2.480282e-137 1.240141e-137
[161,] 1.0000000 1.758249e-136 8.791243e-137
[162,] 1.0000000 1.262759e-135 6.313797e-136
[163,] 1.0000000 8.861637e-135 4.430818e-135
[164,] 1.0000000 6.311552e-134 3.155776e-134
[165,] 1.0000000 4.385392e-133 2.192696e-133
[166,] 1.0000000 3.058349e-132 1.529174e-132
[167,] 1.0000000 1.363451e-131 6.817257e-132
[168,] 1.0000000 6.326054e-131 3.163027e-131
[169,] 1.0000000 4.306306e-130 2.153153e-130
[170,] 1.0000000 2.916511e-129 1.458256e-129
[171,] 1.0000000 1.434010e-128 7.170048e-129
[172,] 1.0000000 7.659527e-128 3.829763e-128
[173,] 1.0000000 3.898096e-127 1.949048e-127
[174,] 1.0000000 1.373586e-126 6.867931e-127
[175,] 1.0000000 7.052713e-126 3.526356e-126
[176,] 1.0000000 3.125534e-125 1.562767e-125
[177,] 1.0000000 7.533314e-125 3.766657e-125
[178,] 1.0000000 5.089730e-124 2.544865e-124
[179,] 1.0000000 2.823128e-123 1.411564e-123
[180,] 1.0000000 1.848844e-122 9.244221e-123
[181,] 1.0000000 1.161433e-121 5.807166e-122
[182,] 1.0000000 7.530845e-121 3.765423e-121
[183,] 1.0000000 4.945832e-120 2.472916e-120
[184,] 1.0000000 3.175377e-119 1.587689e-119
[185,] 1.0000000 2.042847e-118 1.021424e-118
[186,] 1.0000000 1.013895e-117 5.069476e-118
[187,] 1.0000000 3.938831e-117 1.969415e-117
[188,] 1.0000000 1.977995e-116 9.889977e-117
[189,] 1.0000000 1.007420e-115 5.037100e-116
[190,] 1.0000000 6.375272e-116 3.187636e-116
[191,] 1.0000000 3.314475e-115 1.657238e-115
[192,] 1.0000000 1.741010e-114 8.705049e-115
[193,] 1.0000000 9.227527e-114 4.613763e-114
[194,] 1.0000000 5.265745e-113 2.632873e-113
[195,] 1.0000000 1.187410e-112 5.937050e-113
[196,] 1.0000000 6.599407e-112 3.299704e-112
[197,] 1.0000000 3.442296e-111 1.721148e-111
[198,] 1.0000000 1.867705e-110 9.338525e-111
[199,] 1.0000000 1.851183e-110 9.255914e-111
[200,] 1.0000000 1.052821e-109 5.264107e-110
[201,] 1.0000000 6.465497e-109 3.232749e-109
[202,] 1.0000000 3.950933e-108 1.975467e-108
[203,] 1.0000000 1.971419e-107 9.857097e-108
[204,] 1.0000000 1.091789e-106 5.458946e-107
[205,] 1.0000000 6.739202e-106 3.369601e-106
[206,] 1.0000000 4.070813e-105 2.035406e-105
[207,] 1.0000000 2.467736e-104 1.233868e-104
[208,] 1.0000000 1.437769e-103 7.188847e-104
[209,] 1.0000000 8.559181e-103 4.279590e-103
[210,] 1.0000000 3.890513e-102 1.945257e-102
[211,] 1.0000000 2.307902e-101 1.153951e-101
[212,] 1.0000000 1.349082e-100 6.745409e-101
[213,] 1.0000000 7.866339e-100 3.933170e-100
[214,] 1.0000000 4.680174e-99 2.340087e-99
[215,] 1.0000000 2.699914e-98 1.349957e-98
[216,] 1.0000000 1.593521e-97 7.967605e-98
[217,] 1.0000000 9.277306e-97 4.638653e-97
[218,] 1.0000000 5.267942e-96 2.633971e-96
[219,] 1.0000000 3.056520e-95 1.528260e-95
[220,] 1.0000000 1.719012e-94 8.595062e-95
[221,] 1.0000000 9.423447e-94 4.711724e-94
[222,] 1.0000000 4.864129e-93 2.432065e-93
[223,] 1.0000000 2.698320e-92 1.349160e-92
[224,] 1.0000000 1.515895e-91 7.579474e-92
[225,] 1.0000000 8.330060e-91 4.165030e-91
[226,] 1.0000000 4.552405e-90 2.276202e-90
[227,] 1.0000000 2.476545e-89 1.238272e-89
[228,] 1.0000000 1.350140e-88 6.750698e-89
[229,] 1.0000000 7.271984e-88 3.635992e-88
[230,] 1.0000000 3.896967e-87 1.948484e-87
[231,] 1.0000000 2.077760e-86 1.038880e-86
[232,] 1.0000000 1.123832e-85 5.619158e-86
[233,] 1.0000000 5.932133e-85 2.966066e-85
[234,] 1.0000000 3.138790e-84 1.569395e-84
[235,] 1.0000000 1.498153e-83 7.490766e-84
[236,] 1.0000000 7.200266e-83 3.600133e-83
[237,] 1.0000000 3.465835e-82 1.732918e-82
[238,] 1.0000000 1.840094e-81 9.200468e-82
[239,] 1.0000000 8.828882e-81 4.414441e-81
[240,] 1.0000000 4.238824e-80 2.119412e-80
[241,] 1.0000000 2.215489e-79 1.107744e-79
[242,] 1.0000000 1.059322e-78 5.296609e-79
[243,] 1.0000000 5.498331e-78 2.749165e-78
[244,] 1.0000000 1.987154e-77 9.935772e-78
[245,] 1.0000000 5.802622e-78 2.901311e-78
[246,] 1.0000000 2.928944e-77 1.464472e-77
[247,] 1.0000000 1.427303e-76 7.136516e-77
[248,] 1.0000000 7.388195e-76 3.694097e-76
[249,] 1.0000000 3.511975e-75 1.755987e-75
[250,] 1.0000000 1.738928e-74 8.694640e-75
[251,] 1.0000000 8.225467e-74 4.112734e-74
[252,] 1.0000000 3.887067e-73 1.943533e-73
[253,] 1.0000000 1.936282e-72 9.681411e-73
[254,] 1.0000000 8.101872e-72 4.050936e-72
[255,] 1.0000000 3.794102e-71 1.897051e-71
[256,] 1.0000000 1.826726e-70 9.133631e-71
[257,] 1.0000000 8.503974e-70 4.251987e-70
[258,] 1.0000000 2.881226e-69 1.440613e-69
[259,] 1.0000000 1.336075e-68 6.680375e-69
[260,] 1.0000000 6.434815e-68 3.217408e-68
[261,] 1.0000000 2.964194e-67 1.482097e-67
[262,] 1.0000000 1.362323e-66 6.811614e-67
[263,] 1.0000000 6.245535e-66 3.122767e-66
[264,] 1.0000000 2.973588e-65 1.486794e-65
[265,] 1.0000000 1.372955e-64 6.864776e-65
[266,] 1.0000000 6.483862e-64 3.241931e-64
[267,] 1.0000000 3.065062e-63 1.532531e-63
[268,] 1.0000000 1.393436e-62 6.967182e-63
[269,] 1.0000000 5.613624e-62 2.806812e-62
[270,] 1.0000000 2.563572e-61 1.281786e-61
[271,] 1.0000000 1.133035e-60 5.665176e-61
[272,] 1.0000000 5.057098e-60 2.528549e-60
[273,] 1.0000000 2.331364e-59 1.165682e-59
[274,] 1.0000000 1.029533e-58 5.147665e-59
[275,] 1.0000000 4.661273e-58 2.330636e-58
[276,] 1.0000000 2.035881e-57 1.017940e-57
[277,] 1.0000000 8.859554e-57 4.429777e-57
[278,] 1.0000000 3.874733e-56 1.937366e-56
[279,] 1.0000000 1.653626e-55 8.268128e-56
[280,] 1.0000000 7.306698e-55 3.653349e-55
[281,] 1.0000000 3.091235e-54 1.545617e-54
[282,] 1.0000000 1.332340e-53 6.661700e-54
[283,] 1.0000000 5.603748e-53 2.801874e-53
[284,] 1.0000000 2.347421e-52 1.173710e-52
[285,] 1.0000000 1.011522e-51 5.057611e-52
[286,] 1.0000000 3.470761e-51 1.735380e-51
[287,] 1.0000000 1.433157e-50 7.165787e-51
[288,] 1.0000000 5.892636e-50 2.946318e-50
[289,] 1.0000000 2.406769e-49 1.203385e-49
[290,] 1.0000000 1.010638e-48 5.053188e-49
[291,] 1.0000000 4.131765e-48 2.065883e-48
[292,] 1.0000000 1.709760e-47 8.548800e-48
[293,] 1.0000000 6.862904e-47 3.431452e-47
[294,] 1.0000000 2.739420e-46 1.369710e-46
[295,] 1.0000000 9.243566e-46 4.621783e-46
[296,] 1.0000000 3.773930e-45 1.886965e-45
[297,] 1.0000000 1.485752e-44 7.428762e-45
[298,] 1.0000000 5.953781e-44 2.976891e-44
[299,] 1.0000000 2.355398e-43 1.177699e-43
[300,] 1.0000000 9.123191e-43 4.561596e-43
[301,] 1.0000000 3.515290e-42 1.757645e-42
[302,] 1.0000000 1.202977e-41 6.014883e-42
[303,] 1.0000000 4.589251e-41 2.294626e-41
[304,] 1.0000000 1.769288e-40 8.846439e-41
[305,] 1.0000000 6.816804e-40 3.408402e-40
[306,] 1.0000000 2.343019e-39 1.171510e-39
[307,] 1.0000000 8.608296e-39 4.304148e-39
[308,] 1.0000000 3.114822e-38 1.557411e-38
[309,] 1.0000000 1.163686e-37 5.818430e-38
[310,] 1.0000000 4.230651e-37 2.115326e-37
[311,] 1.0000000 1.527991e-36 7.639955e-37
[312,] 1.0000000 5.482204e-36 2.741102e-36
[313,] 1.0000000 1.953839e-35 9.769197e-36
[314,] 1.0000000 6.994808e-35 3.497404e-35
[315,] 1.0000000 2.461911e-34 1.230956e-34
[316,] 1.0000000 8.743055e-34 4.371528e-34
[317,] 1.0000000 3.067151e-33 1.533576e-33
[318,] 1.0000000 1.063840e-32 5.319201e-33
[319,] 1.0000000 3.650664e-32 1.825332e-32
[320,] 1.0000000 1.244026e-31 6.220131e-32
[321,] 1.0000000 4.209433e-31 2.104717e-31
[322,] 1.0000000 1.414255e-30 7.071274e-31
[323,] 1.0000000 4.717527e-30 2.358763e-30
[324,] 1.0000000 1.590148e-29 7.950739e-30
[325,] 1.0000000 5.255731e-29 2.627866e-29
[326,] 1.0000000 1.716936e-28 8.584679e-29
[327,] 1.0000000 5.720107e-28 2.860054e-28
[328,] 1.0000000 1.842462e-27 9.212309e-28
[329,] 1.0000000 5.958614e-27 2.979307e-27
[330,] 1.0000000 1.933843e-26 9.669214e-27
[331,] 1.0000000 6.223739e-26 3.111869e-26
[332,] 1.0000000 1.949054e-25 9.745268e-26
[333,] 1.0000000 6.056215e-25 3.028108e-25
[334,] 1.0000000 1.867015e-24 9.335076e-25
[335,] 1.0000000 5.958287e-24 2.979144e-24
[336,] 1.0000000 1.808832e-23 9.044158e-24
[337,] 1.0000000 5.651498e-23 2.825749e-23
[338,] 1.0000000 1.688598e-22 8.442991e-23
[339,] 1.0000000 5.003528e-22 2.501764e-22
[340,] 1.0000000 1.345220e-21 6.726102e-22
[341,] 1.0000000 3.750855e-21 1.875428e-21
[342,] 1.0000000 1.124412e-20 5.622058e-21
[343,] 1.0000000 3.069099e-20 1.534550e-20
[344,] 1.0000000 8.841589e-20 4.420794e-20
[345,] 1.0000000 2.491506e-19 1.245753e-19
[346,] 1.0000000 6.957594e-19 3.478797e-19
[347,] 1.0000000 2.010850e-18 1.005425e-18
[348,] 1.0000000 5.528211e-18 2.764105e-18
[349,] 1.0000000 1.505568e-17 7.527842e-18
[350,] 1.0000000 4.061301e-17 2.030651e-17
[351,] 1.0000000 1.108940e-16 5.544701e-17
[352,] 1.0000000 2.934615e-16 1.467308e-16
[353,] 1.0000000 7.688602e-16 3.844301e-16
[354,] 1.0000000 2.114519e-15 1.057260e-15
[355,] 1.0000000 5.581670e-15 2.790835e-15
[356,] 1.0000000 1.422068e-14 7.110342e-15
[357,] 1.0000000 3.584977e-14 1.792488e-14
[358,] 1.0000000 8.940791e-14 4.470396e-14
[359,] 1.0000000 2.205456e-13 1.102728e-13
[360,] 1.0000000 5.527686e-13 2.763843e-13
[361,] 1.0000000 1.334095e-12 6.670476e-13
[362,] 1.0000000 3.072227e-12 1.536114e-12
[363,] 1.0000000 7.622041e-12 3.811021e-12
[364,] 1.0000000 1.890643e-11 9.453213e-12
[365,] 1.0000000 4.358722e-11 2.179361e-11
[366,] 1.0000000 9.921869e-11 4.960935e-11
[367,] 1.0000000 2.427350e-10 1.213675e-10
[368,] 1.0000000 5.225316e-10 2.612658e-10
[369,] 1.0000000 9.975512e-10 4.987756e-10
[370,] 1.0000000 2.160301e-09 1.080151e-09
[371,] 1.0000000 4.611232e-09 2.305616e-09
[372,] 1.0000000 1.084087e-08 5.420436e-09
[373,] 1.0000000 2.416476e-08 1.208238e-08
[374,] 1.0000000 4.973695e-08 2.486848e-08
[375,] 0.9999999 1.086228e-07 5.431139e-08
[376,] 0.9999999 2.072455e-07 1.036227e-07
[377,] 0.9999998 3.905439e-07 1.952720e-07
[378,] 0.9999996 7.636398e-07 3.818199e-07
[379,] 0.9999992 1.508221e-06 7.541106e-07
[380,] 0.9999985 2.981606e-06 1.490803e-06
[381,] 0.9999971 5.748432e-06 2.874216e-06
[382,] 0.9999942 1.162291e-05 5.811454e-06
[383,] 0.9999880 2.408334e-05 1.204167e-05
[384,] 0.9999838 3.237006e-05 1.618503e-05
[385,] 0.9999746 5.085836e-05 2.542918e-05
[386,] 0.9999498 1.004109e-04 5.020546e-05
[387,] 0.9999033 1.934453e-04 9.672267e-05
[388,] 0.9998421 3.158941e-04 1.579471e-04
[389,] 0.9997056 5.888583e-04 2.944291e-04
[390,] 0.9996335 7.330176e-04 3.665088e-04
[391,] 0.9993539 1.292101e-03 6.460506e-04
[392,] 0.9990169 1.966253e-03 9.831267e-04
[393,] 0.9983338 3.332432e-03 1.666216e-03
[394,] 0.9971950 5.610010e-03 2.805005e-03
[395,] 0.9950441 9.911888e-03 4.955944e-03
[396,] 0.9925095 1.498099e-02 7.490493e-03
[397,] 0.9883788 2.324236e-02 1.162118e-02
[398,] 0.9825882 3.482356e-02 1.741178e-02
[399,] 0.9739988 5.200233e-02 2.600117e-02
[400,] 0.9571630 8.567394e-02 4.283697e-02
[401,] 0.9313077 1.373846e-01 6.869229e-02
[402,] 0.9152946 1.694107e-01 8.470536e-02
[403,] 0.8692038 2.615924e-01 1.307962e-01
[404,] 0.8056023 3.887953e-01 1.943977e-01
[405,] 0.7418097 5.163806e-01 2.581903e-01
[406,] 0.6467826 7.064347e-01 3.532174e-01
[407,] 0.7987507 4.024987e-01 2.012493e-01
[408,] 0.6871032 6.257936e-01 3.128968e-01
[409,] 0.5494810 9.010379e-01 4.505190e-01
[410,] 0.3868022 7.736045e-01 6.131978e-01
> postscript(file="/var/www/html/rcomp/tmp/1ph0j1291225443.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/2ph0j1291225443.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/3ph0j1291225443.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/40qhm1291225443.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/50qhm1291225443.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> qqnorm(mysum$resid, main='Residual Normal Q-Q Plot')
> qqline(mysum$resid)
> grid()
> dev.off()
null device
1
> (myerror <- as.ts(mysum$resid))
Time Series:
Start = 1
End = 431
Frequency = 1
1 2 3 4 5
208308.99237 2275662.88080 795828.19898 -270987.19826 -809715.40781
6 7 8 9 10
-22266.82622 800116.23761 -49155.83317 -691527.85593 520274.61695
11 12 13 14 15
300395.15083 657029.06624 -223632.43021 -214806.51231 -584732.95539
16 17 18 19 20
567057.94638 -120685.92990 -555208.39746 -692034.64645 -22517.65388
21 22 23 24 25
147556.65102 -198905.77113 144552.96763 -85472.26929 -172515.67738
26 27 28 29 30
199630.89796 -806793.29430 -91537.29089 493285.69665 102851.24420
31 32 33 34 35
321503.16649 -104258.33967 255545.92790 160686.08391 -132980.09112
36 37 38 39 40
-66128.03087 238335.88236 142898.18103 44953.18110 -135648.29772
41 42 43 44 45
-143145.37863 -91583.43595 121975.57705 -56834.07852 44665.36706
46 47 48 49 50
2854.95073 -109747.17134 -353406.30019 28280.34861 -149859.43355
51 52 53 54 55
1830.56568 -154893.59582 -163816.82868 108069.30970 -398627.52899
56 57 58 59 60
96414.42623 344119.88167 65897.34530 -39937.02809 -268342.71677
61 62 63 64 65
30684.17243 -165191.81686 54293.11167 -229559.34975 44287.34309
66 67 68 69 70
-30222.19008 -89225.87105 12463.09493 -62985.93646 -110309.30452
71 72 73 74 75
-206946.64271 90240.21711 -88273.72819 -21316.72675 -30243.60472
76 77 78 79 80
-66184.57835 -44934.51284 -101207.15920 -11678.29645 63643.05378
81 82 83 84 85
34641.61231 -73790.49159 -104685.81508 29937.67397 -12725.14358
86 87 88 89 90
-167364.74107 91819.60560 -38862.44579 -91394.42511 -5660.57195
91 92 93 94 95
13129.61159 -14033.50196 137631.86438 -89777.21321 90481.59227
96 97 98 99 100
10678.39226 19473.69877 21519.12241 -48344.28494 42245.68748
101 102 103 104 105
-21816.82929 -133542.53336 -108868.71359 -98741.69172 39063.39389
106 107 108 109 110
-59967.40542 -289570.08055 198107.49533 -86218.57294 -66492.21629
111 112 113 114 115
83864.87293 -222070.80804 -84451.81816 -19997.95631 112369.79821
116 117 118 119 120
39988.90315 -24563.04824 71140.54934 -86911.82797 3663.56051
121 122 123 124 125
-65922.30668 9489.98057 -3418.89374 152158.54615 29286.34957
126 127 128 129 130
-138918.74718 4695.80258 5100.36378 14817.85613 4695.80258
131 132 133 134 135
4973.42744 65208.58145 4695.80258 5357.41047 -1120.95879
136 137 138 139 140
3882.97829 9075.68325 4364.04313 444.25026 4695.80258
141 142 143 144 145
17633.90268 4695.80258 -7587.46063 5383.56487 14587.67173
146 147 148 149 150
10257.96358 4695.80258 5344.44639 13569.03592 4560.05297
151 152 153 154 155
-16424.65558 -4530.63985 5628.27428 -9949.03909 4418.84197
156 157 158 159 160
4695.80258 4695.80258 2967.75607 24390.36985 4695.80258
161 162 163 164 165
4695.80258 4695.80258 4695.80258 -5759.11887 -11139.17833
166 167 168 169 170
4146.91131 4695.80258 100101.84568 17631.71406 4695.80258
171 172 173 174 175
4695.80258 502.88764 4695.80258 -2747.76674 4695.80258
176 177 178 179 180
5279.74355 21636.35992 95888.64101 4695.80258 4695.80258
181 182 183 184 185
26546.37872 69353.44777 5194.43625 27469.56632 82604.48615
186 187 188 189 190
71777.73429 20793.09613 505.91678 -28647.56861 4695.80258
191 192 193 194 195
10293.81670 4695.80258 3547.21593 4695.80258 5099.07171
196 197 198 199 200
96745.06614 -159756.58648 96745.06614 96745.06614 -11466.83054
201 202 203 204 205
96745.06614 96745.06614 96745.06614 48580.16067 60559.55258
206 207 208 209 210
84396.23774 92599.63230 96745.06614 34240.26905 74513.23008
211 212 213 214 215
4695.80258 4695.80258 25925.93172 -24615.08018 -11571.51316
216 217 218 219 220
5068.79126 5836.24263 6492.62189 5068.79126 15978.94819
221 222 223 224 225
-14524.58874 4695.80258 4882.91782 2842.93476 4815.88543
226 227 228 229 230
5637.80160 5237.15542 4695.80258 -7592.79784 4695.80258
231 232 233 234 235
-18346.43070 8642.79244 4695.80258 -8722.67173 4695.80258
236 237 238 239 240
4749.02911 4695.80258 6664.52600 4695.80258 4695.80258
241 242 243 244 245
4695.80258 1072.29490 4695.80258 7849.74472 -20413.98966
246 247 248 249 250
96745.06614 96745.06614 -1625.30469 96745.06614 96745.06614
251 252 253 254 255
1105.51459 96745.06614 -10712.38062 -25693.29890 -100020.13876
256 257 258 259 260
4695.80258 9540.54090 15504.80685 96745.06614 4695.80258
261 262 263 264 265
96745.06614 96745.06614 11524.03289 25124.49007 96745.06614
266 267 268 269 270
4695.80258 96745.06614 39068.86267 96745.06614 7614.44532
271 272 273 274 275
96745.06614 96745.06614 96745.06614 2280.74545 4695.80258
276 277 278 279 280
6856.73328 5750.88219 4695.80258 -32509.02284 -13492.53170
281 282 283 284 285
96745.06614 4695.80258 13981.74577 4695.80258 1474.75604
286 287 288 289 290
4695.80258 73875.87396 5305.87880 96745.06614 -4616.91129
291 292 293 294 295
96745.06614 11869.04635 17785.58720 4695.80258 -4372.82474
296 297 298 299 300
90735.14542 4695.80258 96745.06614 4695.80258 6374.25164
301 302 303 304 305
83624.34209 -7133.77660 4882.29692 96745.06614 32437.27291
306 307 308 309 310
6909.15144 96745.06614 44212.05015 11845.65282 96745.06614
311 312 313 314 315
96745.06614 -10213.42054 96745.06614 66041.66334 6370.75503
316 317 318 319 320
18828.45632 5508.08560 -12520.74431 5617.12230 4695.80258
321 322 323 324 325
4695.80258 4695.80258 4695.80258 96745.06614 4695.80258
326 327 328 329 330
1889.38269 96745.06614 96745.06614 4695.80258 4695.80258
331 332 333 334 335
4695.80258 4695.80258 4695.80258 3953.39132 4871.95470
336 337 338 339 340
3846.84201 -15467.40013 4695.80258 2916.79186 10119.59876
341 342 343 344 345
-18672.05287 4695.80258 4695.80258 4695.80258 -3140.29155
346 347 348 349 350
4695.80258 -4006.90833 4695.80258 4695.80258 2493.78740
351 352 353 354 355
11018.66125 3600.11020 16051.07917 476.58495 4695.80258
356 357 358 359 360
4695.80258 -22396.84628 4695.80258 4695.80258 4695.80258
361 362 363 364 365
5450.93162 4695.80258 4695.80258 -22742.53307 8531.29201
366 367 368 369 370
4695.80258 4695.80258 4695.80258 4695.80258 5073.19778
371 372 373 374 375
4695.80258 13502.47832 -4122.43855 -11364.54599 4695.80258
376 377 378 379 380
4695.80258 -5249.64067 10434.70666 -127.53011 4695.80258
381 382 383 384 385
4695.80258 -16852.39645 5978.97241 4695.80258 4730.19031
386 387 388 389 390
-93173.15808 18295.44481 -1201.64187 4333.03834 -87.88601
391 392 393 394 395
7724.02211 -17053.27366 -49494.53350 40647.49443 -40778.26494
396 397 398 399 400
-24144.92966 -72002.46594 -432306.25565 -21770.74622 -135659.00978
401 402 403 404 405
-29640.24458 13728.65001 -12827.11030 -29957.78483 -42211.45613
406 407 408 409 410
-6534.19128 -8555.89563 -5487.64777 -8564.90722 -99513.27125
411 412 413 414 415
-14573.70679 -157855.32051 -13610.02368 68913.23945 -188433.40848
416 417 418 419 420
-135107.61599 132680.84269 -5736.85380 -53424.84627 -77168.91215
421 422 423 424 425
-79330.08351 -93066.08147 -731095.45915 -372773.51927 -35145.58987
426 427 428 429 430
-331902.94857 -266234.14336 -322175.21212 -133587.07324 -30715.63764
431
-332780.81006
> postscript(file="/var/www/html/rcomp/tmp/6shy71291225443.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> dum <- cbind(lag(myerror,k=1),myerror)
> dum
Time Series:
Start = 0
End = 431
Frequency = 1
lag(myerror, k = 1) myerror
0 208308.99237 NA
1 2275662.88080 208308.99237
2 795828.19898 2275662.88080
3 -270987.19826 795828.19898
4 -809715.40781 -270987.19826
5 -22266.82622 -809715.40781
6 800116.23761 -22266.82622
7 -49155.83317 800116.23761
8 -691527.85593 -49155.83317
9 520274.61695 -691527.85593
10 300395.15083 520274.61695
11 657029.06624 300395.15083
12 -223632.43021 657029.06624
13 -214806.51231 -223632.43021
14 -584732.95539 -214806.51231
15 567057.94638 -584732.95539
16 -120685.92990 567057.94638
17 -555208.39746 -120685.92990
18 -692034.64645 -555208.39746
19 -22517.65388 -692034.64645
20 147556.65102 -22517.65388
21 -198905.77113 147556.65102
22 144552.96763 -198905.77113
23 -85472.26929 144552.96763
24 -172515.67738 -85472.26929
25 199630.89796 -172515.67738
26 -806793.29430 199630.89796
27 -91537.29089 -806793.29430
28 493285.69665 -91537.29089
29 102851.24420 493285.69665
30 321503.16649 102851.24420
31 -104258.33967 321503.16649
32 255545.92790 -104258.33967
33 160686.08391 255545.92790
34 -132980.09112 160686.08391
35 -66128.03087 -132980.09112
36 238335.88236 -66128.03087
37 142898.18103 238335.88236
38 44953.18110 142898.18103
39 -135648.29772 44953.18110
40 -143145.37863 -135648.29772
41 -91583.43595 -143145.37863
42 121975.57705 -91583.43595
43 -56834.07852 121975.57705
44 44665.36706 -56834.07852
45 2854.95073 44665.36706
46 -109747.17134 2854.95073
47 -353406.30019 -109747.17134
48 28280.34861 -353406.30019
49 -149859.43355 28280.34861
50 1830.56568 -149859.43355
51 -154893.59582 1830.56568
52 -163816.82868 -154893.59582
53 108069.30970 -163816.82868
54 -398627.52899 108069.30970
55 96414.42623 -398627.52899
56 344119.88167 96414.42623
57 65897.34530 344119.88167
58 -39937.02809 65897.34530
59 -268342.71677 -39937.02809
60 30684.17243 -268342.71677
61 -165191.81686 30684.17243
62 54293.11167 -165191.81686
63 -229559.34975 54293.11167
64 44287.34309 -229559.34975
65 -30222.19008 44287.34309
66 -89225.87105 -30222.19008
67 12463.09493 -89225.87105
68 -62985.93646 12463.09493
69 -110309.30452 -62985.93646
70 -206946.64271 -110309.30452
71 90240.21711 -206946.64271
72 -88273.72819 90240.21711
73 -21316.72675 -88273.72819
74 -30243.60472 -21316.72675
75 -66184.57835 -30243.60472
76 -44934.51284 -66184.57835
77 -101207.15920 -44934.51284
78 -11678.29645 -101207.15920
79 63643.05378 -11678.29645
80 34641.61231 63643.05378
81 -73790.49159 34641.61231
82 -104685.81508 -73790.49159
83 29937.67397 -104685.81508
84 -12725.14358 29937.67397
85 -167364.74107 -12725.14358
86 91819.60560 -167364.74107
87 -38862.44579 91819.60560
88 -91394.42511 -38862.44579
89 -5660.57195 -91394.42511
90 13129.61159 -5660.57195
91 -14033.50196 13129.61159
92 137631.86438 -14033.50196
93 -89777.21321 137631.86438
94 90481.59227 -89777.21321
95 10678.39226 90481.59227
96 19473.69877 10678.39226
97 21519.12241 19473.69877
98 -48344.28494 21519.12241
99 42245.68748 -48344.28494
100 -21816.82929 42245.68748
101 -133542.53336 -21816.82929
102 -108868.71359 -133542.53336
103 -98741.69172 -108868.71359
104 39063.39389 -98741.69172
105 -59967.40542 39063.39389
106 -289570.08055 -59967.40542
107 198107.49533 -289570.08055
108 -86218.57294 198107.49533
109 -66492.21629 -86218.57294
110 83864.87293 -66492.21629
111 -222070.80804 83864.87293
112 -84451.81816 -222070.80804
113 -19997.95631 -84451.81816
114 112369.79821 -19997.95631
115 39988.90315 112369.79821
116 -24563.04824 39988.90315
117 71140.54934 -24563.04824
118 -86911.82797 71140.54934
119 3663.56051 -86911.82797
120 -65922.30668 3663.56051
121 9489.98057 -65922.30668
122 -3418.89374 9489.98057
123 152158.54615 -3418.89374
124 29286.34957 152158.54615
125 -138918.74718 29286.34957
126 4695.80258 -138918.74718
127 5100.36378 4695.80258
128 14817.85613 5100.36378
129 4695.80258 14817.85613
130 4973.42744 4695.80258
131 65208.58145 4973.42744
132 4695.80258 65208.58145
133 5357.41047 4695.80258
134 -1120.95879 5357.41047
135 3882.97829 -1120.95879
136 9075.68325 3882.97829
137 4364.04313 9075.68325
138 444.25026 4364.04313
139 4695.80258 444.25026
140 17633.90268 4695.80258
141 4695.80258 17633.90268
142 -7587.46063 4695.80258
143 5383.56487 -7587.46063
144 14587.67173 5383.56487
145 10257.96358 14587.67173
146 4695.80258 10257.96358
147 5344.44639 4695.80258
148 13569.03592 5344.44639
149 4560.05297 13569.03592
150 -16424.65558 4560.05297
151 -4530.63985 -16424.65558
152 5628.27428 -4530.63985
153 -9949.03909 5628.27428
154 4418.84197 -9949.03909
155 4695.80258 4418.84197
156 4695.80258 4695.80258
157 2967.75607 4695.80258
158 24390.36985 2967.75607
159 4695.80258 24390.36985
160 4695.80258 4695.80258
161 4695.80258 4695.80258
162 4695.80258 4695.80258
163 -5759.11887 4695.80258
164 -11139.17833 -5759.11887
165 4146.91131 -11139.17833
166 4695.80258 4146.91131
167 100101.84568 4695.80258
168 17631.71406 100101.84568
169 4695.80258 17631.71406
170 4695.80258 4695.80258
171 502.88764 4695.80258
172 4695.80258 502.88764
173 -2747.76674 4695.80258
174 4695.80258 -2747.76674
175 5279.74355 4695.80258
176 21636.35992 5279.74355
177 95888.64101 21636.35992
178 4695.80258 95888.64101
179 4695.80258 4695.80258
180 26546.37872 4695.80258
181 69353.44777 26546.37872
182 5194.43625 69353.44777
183 27469.56632 5194.43625
184 82604.48615 27469.56632
185 71777.73429 82604.48615
186 20793.09613 71777.73429
187 505.91678 20793.09613
188 -28647.56861 505.91678
189 4695.80258 -28647.56861
190 10293.81670 4695.80258
191 4695.80258 10293.81670
192 3547.21593 4695.80258
193 4695.80258 3547.21593
194 5099.07171 4695.80258
195 96745.06614 5099.07171
196 -159756.58648 96745.06614
197 96745.06614 -159756.58648
198 96745.06614 96745.06614
199 -11466.83054 96745.06614
200 96745.06614 -11466.83054
201 96745.06614 96745.06614
202 96745.06614 96745.06614
203 48580.16067 96745.06614
204 60559.55258 48580.16067
205 84396.23774 60559.55258
206 92599.63230 84396.23774
207 96745.06614 92599.63230
208 34240.26905 96745.06614
209 74513.23008 34240.26905
210 4695.80258 74513.23008
211 4695.80258 4695.80258
212 25925.93172 4695.80258
213 -24615.08018 25925.93172
214 -11571.51316 -24615.08018
215 5068.79126 -11571.51316
216 5836.24263 5068.79126
217 6492.62189 5836.24263
218 5068.79126 6492.62189
219 15978.94819 5068.79126
220 -14524.58874 15978.94819
221 4695.80258 -14524.58874
222 4882.91782 4695.80258
223 2842.93476 4882.91782
224 4815.88543 2842.93476
225 5637.80160 4815.88543
226 5237.15542 5637.80160
227 4695.80258 5237.15542
228 -7592.79784 4695.80258
229 4695.80258 -7592.79784
230 -18346.43070 4695.80258
231 8642.79244 -18346.43070
232 4695.80258 8642.79244
233 -8722.67173 4695.80258
234 4695.80258 -8722.67173
235 4749.02911 4695.80258
236 4695.80258 4749.02911
237 6664.52600 4695.80258
238 4695.80258 6664.52600
239 4695.80258 4695.80258
240 4695.80258 4695.80258
241 1072.29490 4695.80258
242 4695.80258 1072.29490
243 7849.74472 4695.80258
244 -20413.98966 7849.74472
245 96745.06614 -20413.98966
246 96745.06614 96745.06614
247 -1625.30469 96745.06614
248 96745.06614 -1625.30469
249 96745.06614 96745.06614
250 1105.51459 96745.06614
251 96745.06614 1105.51459
252 -10712.38062 96745.06614
253 -25693.29890 -10712.38062
254 -100020.13876 -25693.29890
255 4695.80258 -100020.13876
256 9540.54090 4695.80258
257 15504.80685 9540.54090
258 96745.06614 15504.80685
259 4695.80258 96745.06614
260 96745.06614 4695.80258
261 96745.06614 96745.06614
262 11524.03289 96745.06614
263 25124.49007 11524.03289
264 96745.06614 25124.49007
265 4695.80258 96745.06614
266 96745.06614 4695.80258
267 39068.86267 96745.06614
268 96745.06614 39068.86267
269 7614.44532 96745.06614
270 96745.06614 7614.44532
271 96745.06614 96745.06614
272 96745.06614 96745.06614
273 2280.74545 96745.06614
274 4695.80258 2280.74545
275 6856.73328 4695.80258
276 5750.88219 6856.73328
277 4695.80258 5750.88219
278 -32509.02284 4695.80258
279 -13492.53170 -32509.02284
280 96745.06614 -13492.53170
281 4695.80258 96745.06614
282 13981.74577 4695.80258
283 4695.80258 13981.74577
284 1474.75604 4695.80258
285 4695.80258 1474.75604
286 73875.87396 4695.80258
287 5305.87880 73875.87396
288 96745.06614 5305.87880
289 -4616.91129 96745.06614
290 96745.06614 -4616.91129
291 11869.04635 96745.06614
292 17785.58720 11869.04635
293 4695.80258 17785.58720
294 -4372.82474 4695.80258
295 90735.14542 -4372.82474
296 4695.80258 90735.14542
297 96745.06614 4695.80258
298 4695.80258 96745.06614
299 6374.25164 4695.80258
300 83624.34209 6374.25164
301 -7133.77660 83624.34209
302 4882.29692 -7133.77660
303 96745.06614 4882.29692
304 32437.27291 96745.06614
305 6909.15144 32437.27291
306 96745.06614 6909.15144
307 44212.05015 96745.06614
308 11845.65282 44212.05015
309 96745.06614 11845.65282
310 96745.06614 96745.06614
311 -10213.42054 96745.06614
312 96745.06614 -10213.42054
313 66041.66334 96745.06614
314 6370.75503 66041.66334
315 18828.45632 6370.75503
316 5508.08560 18828.45632
317 -12520.74431 5508.08560
318 5617.12230 -12520.74431
319 4695.80258 5617.12230
320 4695.80258 4695.80258
321 4695.80258 4695.80258
322 4695.80258 4695.80258
323 96745.06614 4695.80258
324 4695.80258 96745.06614
325 1889.38269 4695.80258
326 96745.06614 1889.38269
327 96745.06614 96745.06614
328 4695.80258 96745.06614
329 4695.80258 4695.80258
330 4695.80258 4695.80258
331 4695.80258 4695.80258
332 4695.80258 4695.80258
333 3953.39132 4695.80258
334 4871.95470 3953.39132
335 3846.84201 4871.95470
336 -15467.40013 3846.84201
337 4695.80258 -15467.40013
338 2916.79186 4695.80258
339 10119.59876 2916.79186
340 -18672.05287 10119.59876
341 4695.80258 -18672.05287
342 4695.80258 4695.80258
343 4695.80258 4695.80258
344 -3140.29155 4695.80258
345 4695.80258 -3140.29155
346 -4006.90833 4695.80258
347 4695.80258 -4006.90833
348 4695.80258 4695.80258
349 2493.78740 4695.80258
350 11018.66125 2493.78740
351 3600.11020 11018.66125
352 16051.07917 3600.11020
353 476.58495 16051.07917
354 4695.80258 476.58495
355 4695.80258 4695.80258
356 -22396.84628 4695.80258
357 4695.80258 -22396.84628
358 4695.80258 4695.80258
359 4695.80258 4695.80258
360 5450.93162 4695.80258
361 4695.80258 5450.93162
362 4695.80258 4695.80258
363 -22742.53307 4695.80258
364 8531.29201 -22742.53307
365 4695.80258 8531.29201
366 4695.80258 4695.80258
367 4695.80258 4695.80258
368 4695.80258 4695.80258
369 5073.19778 4695.80258
370 4695.80258 5073.19778
371 13502.47832 4695.80258
372 -4122.43855 13502.47832
373 -11364.54599 -4122.43855
374 4695.80258 -11364.54599
375 4695.80258 4695.80258
376 -5249.64067 4695.80258
377 10434.70666 -5249.64067
378 -127.53011 10434.70666
379 4695.80258 -127.53011
380 4695.80258 4695.80258
381 -16852.39645 4695.80258
382 5978.97241 -16852.39645
383 4695.80258 5978.97241
384 4730.19031 4695.80258
385 -93173.15808 4730.19031
386 18295.44481 -93173.15808
387 -1201.64187 18295.44481
388 4333.03834 -1201.64187
389 -87.88601 4333.03834
390 7724.02211 -87.88601
391 -17053.27366 7724.02211
392 -49494.53350 -17053.27366
393 40647.49443 -49494.53350
394 -40778.26494 40647.49443
395 -24144.92966 -40778.26494
396 -72002.46594 -24144.92966
397 -432306.25565 -72002.46594
398 -21770.74622 -432306.25565
399 -135659.00978 -21770.74622
400 -29640.24458 -135659.00978
401 13728.65001 -29640.24458
402 -12827.11030 13728.65001
403 -29957.78483 -12827.11030
404 -42211.45613 -29957.78483
405 -6534.19128 -42211.45613
406 -8555.89563 -6534.19128
407 -5487.64777 -8555.89563
408 -8564.90722 -5487.64777
409 -99513.27125 -8564.90722
410 -14573.70679 -99513.27125
411 -157855.32051 -14573.70679
412 -13610.02368 -157855.32051
413 68913.23945 -13610.02368
414 -188433.40848 68913.23945
415 -135107.61599 -188433.40848
416 132680.84269 -135107.61599
417 -5736.85380 132680.84269
418 -53424.84627 -5736.85380
419 -77168.91215 -53424.84627
420 -79330.08351 -77168.91215
421 -93066.08147 -79330.08351
422 -731095.45915 -93066.08147
423 -372773.51927 -731095.45915
424 -35145.58987 -372773.51927
425 -331902.94857 -35145.58987
426 -266234.14336 -331902.94857
427 -322175.21212 -266234.14336
428 -133587.07324 -322175.21212
429 -30715.63764 -133587.07324
430 -332780.81006 -30715.63764
431 NA -332780.81006
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] 2275662.88080 208308.99237
[2,] 795828.19898 2275662.88080
[3,] -270987.19826 795828.19898
[4,] -809715.40781 -270987.19826
[5,] -22266.82622 -809715.40781
[6,] 800116.23761 -22266.82622
[7,] -49155.83317 800116.23761
[8,] -691527.85593 -49155.83317
[9,] 520274.61695 -691527.85593
[10,] 300395.15083 520274.61695
[11,] 657029.06624 300395.15083
[12,] -223632.43021 657029.06624
[13,] -214806.51231 -223632.43021
[14,] -584732.95539 -214806.51231
[15,] 567057.94638 -584732.95539
[16,] -120685.92990 567057.94638
[17,] -555208.39746 -120685.92990
[18,] -692034.64645 -555208.39746
[19,] -22517.65388 -692034.64645
[20,] 147556.65102 -22517.65388
[21,] -198905.77113 147556.65102
[22,] 144552.96763 -198905.77113
[23,] -85472.26929 144552.96763
[24,] -172515.67738 -85472.26929
[25,] 199630.89796 -172515.67738
[26,] -806793.29430 199630.89796
[27,] -91537.29089 -806793.29430
[28,] 493285.69665 -91537.29089
[29,] 102851.24420 493285.69665
[30,] 321503.16649 102851.24420
[31,] -104258.33967 321503.16649
[32,] 255545.92790 -104258.33967
[33,] 160686.08391 255545.92790
[34,] -132980.09112 160686.08391
[35,] -66128.03087 -132980.09112
[36,] 238335.88236 -66128.03087
[37,] 142898.18103 238335.88236
[38,] 44953.18110 142898.18103
[39,] -135648.29772 44953.18110
[40,] -143145.37863 -135648.29772
[41,] -91583.43595 -143145.37863
[42,] 121975.57705 -91583.43595
[43,] -56834.07852 121975.57705
[44,] 44665.36706 -56834.07852
[45,] 2854.95073 44665.36706
[46,] -109747.17134 2854.95073
[47,] -353406.30019 -109747.17134
[48,] 28280.34861 -353406.30019
[49,] -149859.43355 28280.34861
[50,] 1830.56568 -149859.43355
[51,] -154893.59582 1830.56568
[52,] -163816.82868 -154893.59582
[53,] 108069.30970 -163816.82868
[54,] -398627.52899 108069.30970
[55,] 96414.42623 -398627.52899
[56,] 344119.88167 96414.42623
[57,] 65897.34530 344119.88167
[58,] -39937.02809 65897.34530
[59,] -268342.71677 -39937.02809
[60,] 30684.17243 -268342.71677
[61,] -165191.81686 30684.17243
[62,] 54293.11167 -165191.81686
[63,] -229559.34975 54293.11167
[64,] 44287.34309 -229559.34975
[65,] -30222.19008 44287.34309
[66,] -89225.87105 -30222.19008
[67,] 12463.09493 -89225.87105
[68,] -62985.93646 12463.09493
[69,] -110309.30452 -62985.93646
[70,] -206946.64271 -110309.30452
[71,] 90240.21711 -206946.64271
[72,] -88273.72819 90240.21711
[73,] -21316.72675 -88273.72819
[74,] -30243.60472 -21316.72675
[75,] -66184.57835 -30243.60472
[76,] -44934.51284 -66184.57835
[77,] -101207.15920 -44934.51284
[78,] -11678.29645 -101207.15920
[79,] 63643.05378 -11678.29645
[80,] 34641.61231 63643.05378
[81,] -73790.49159 34641.61231
[82,] -104685.81508 -73790.49159
[83,] 29937.67397 -104685.81508
[84,] -12725.14358 29937.67397
[85,] -167364.74107 -12725.14358
[86,] 91819.60560 -167364.74107
[87,] -38862.44579 91819.60560
[88,] -91394.42511 -38862.44579
[89,] -5660.57195 -91394.42511
[90,] 13129.61159 -5660.57195
[91,] -14033.50196 13129.61159
[92,] 137631.86438 -14033.50196
[93,] -89777.21321 137631.86438
[94,] 90481.59227 -89777.21321
[95,] 10678.39226 90481.59227
[96,] 19473.69877 10678.39226
[97,] 21519.12241 19473.69877
[98,] -48344.28494 21519.12241
[99,] 42245.68748 -48344.28494
[100,] -21816.82929 42245.68748
[101,] -133542.53336 -21816.82929
[102,] -108868.71359 -133542.53336
[103,] -98741.69172 -108868.71359
[104,] 39063.39389 -98741.69172
[105,] -59967.40542 39063.39389
[106,] -289570.08055 -59967.40542
[107,] 198107.49533 -289570.08055
[108,] -86218.57294 198107.49533
[109,] -66492.21629 -86218.57294
[110,] 83864.87293 -66492.21629
[111,] -222070.80804 83864.87293
[112,] -84451.81816 -222070.80804
[113,] -19997.95631 -84451.81816
[114,] 112369.79821 -19997.95631
[115,] 39988.90315 112369.79821
[116,] -24563.04824 39988.90315
[117,] 71140.54934 -24563.04824
[118,] -86911.82797 71140.54934
[119,] 3663.56051 -86911.82797
[120,] -65922.30668 3663.56051
[121,] 9489.98057 -65922.30668
[122,] -3418.89374 9489.98057
[123,] 152158.54615 -3418.89374
[124,] 29286.34957 152158.54615
[125,] -138918.74718 29286.34957
[126,] 4695.80258 -138918.74718
[127,] 5100.36378 4695.80258
[128,] 14817.85613 5100.36378
[129,] 4695.80258 14817.85613
[130,] 4973.42744 4695.80258
[131,] 65208.58145 4973.42744
[132,] 4695.80258 65208.58145
[133,] 5357.41047 4695.80258
[134,] -1120.95879 5357.41047
[135,] 3882.97829 -1120.95879
[136,] 9075.68325 3882.97829
[137,] 4364.04313 9075.68325
[138,] 444.25026 4364.04313
[139,] 4695.80258 444.25026
[140,] 17633.90268 4695.80258
[141,] 4695.80258 17633.90268
[142,] -7587.46063 4695.80258
[143,] 5383.56487 -7587.46063
[144,] 14587.67173 5383.56487
[145,] 10257.96358 14587.67173
[146,] 4695.80258 10257.96358
[147,] 5344.44639 4695.80258
[148,] 13569.03592 5344.44639
[149,] 4560.05297 13569.03592
[150,] -16424.65558 4560.05297
[151,] -4530.63985 -16424.65558
[152,] 5628.27428 -4530.63985
[153,] -9949.03909 5628.27428
[154,] 4418.84197 -9949.03909
[155,] 4695.80258 4418.84197
[156,] 4695.80258 4695.80258
[157,] 2967.75607 4695.80258
[158,] 24390.36985 2967.75607
[159,] 4695.80258 24390.36985
[160,] 4695.80258 4695.80258
[161,] 4695.80258 4695.80258
[162,] 4695.80258 4695.80258
[163,] -5759.11887 4695.80258
[164,] -11139.17833 -5759.11887
[165,] 4146.91131 -11139.17833
[166,] 4695.80258 4146.91131
[167,] 100101.84568 4695.80258
[168,] 17631.71406 100101.84568
[169,] 4695.80258 17631.71406
[170,] 4695.80258 4695.80258
[171,] 502.88764 4695.80258
[172,] 4695.80258 502.88764
[173,] -2747.76674 4695.80258
[174,] 4695.80258 -2747.76674
[175,] 5279.74355 4695.80258
[176,] 21636.35992 5279.74355
[177,] 95888.64101 21636.35992
[178,] 4695.80258 95888.64101
[179,] 4695.80258 4695.80258
[180,] 26546.37872 4695.80258
[181,] 69353.44777 26546.37872
[182,] 5194.43625 69353.44777
[183,] 27469.56632 5194.43625
[184,] 82604.48615 27469.56632
[185,] 71777.73429 82604.48615
[186,] 20793.09613 71777.73429
[187,] 505.91678 20793.09613
[188,] -28647.56861 505.91678
[189,] 4695.80258 -28647.56861
[190,] 10293.81670 4695.80258
[191,] 4695.80258 10293.81670
[192,] 3547.21593 4695.80258
[193,] 4695.80258 3547.21593
[194,] 5099.07171 4695.80258
[195,] 96745.06614 5099.07171
[196,] -159756.58648 96745.06614
[197,] 96745.06614 -159756.58648
[198,] 96745.06614 96745.06614
[199,] -11466.83054 96745.06614
[200,] 96745.06614 -11466.83054
[201,] 96745.06614 96745.06614
[202,] 96745.06614 96745.06614
[203,] 48580.16067 96745.06614
[204,] 60559.55258 48580.16067
[205,] 84396.23774 60559.55258
[206,] 92599.63230 84396.23774
[207,] 96745.06614 92599.63230
[208,] 34240.26905 96745.06614
[209,] 74513.23008 34240.26905
[210,] 4695.80258 74513.23008
[211,] 4695.80258 4695.80258
[212,] 25925.93172 4695.80258
[213,] -24615.08018 25925.93172
[214,] -11571.51316 -24615.08018
[215,] 5068.79126 -11571.51316
[216,] 5836.24263 5068.79126
[217,] 6492.62189 5836.24263
[218,] 5068.79126 6492.62189
[219,] 15978.94819 5068.79126
[220,] -14524.58874 15978.94819
[221,] 4695.80258 -14524.58874
[222,] 4882.91782 4695.80258
[223,] 2842.93476 4882.91782
[224,] 4815.88543 2842.93476
[225,] 5637.80160 4815.88543
[226,] 5237.15542 5637.80160
[227,] 4695.80258 5237.15542
[228,] -7592.79784 4695.80258
[229,] 4695.80258 -7592.79784
[230,] -18346.43070 4695.80258
[231,] 8642.79244 -18346.43070
[232,] 4695.80258 8642.79244
[233,] -8722.67173 4695.80258
[234,] 4695.80258 -8722.67173
[235,] 4749.02911 4695.80258
[236,] 4695.80258 4749.02911
[237,] 6664.52600 4695.80258
[238,] 4695.80258 6664.52600
[239,] 4695.80258 4695.80258
[240,] 4695.80258 4695.80258
[241,] 1072.29490 4695.80258
[242,] 4695.80258 1072.29490
[243,] 7849.74472 4695.80258
[244,] -20413.98966 7849.74472
[245,] 96745.06614 -20413.98966
[246,] 96745.06614 96745.06614
[247,] -1625.30469 96745.06614
[248,] 96745.06614 -1625.30469
[249,] 96745.06614 96745.06614
[250,] 1105.51459 96745.06614
[251,] 96745.06614 1105.51459
[252,] -10712.38062 96745.06614
[253,] -25693.29890 -10712.38062
[254,] -100020.13876 -25693.29890
[255,] 4695.80258 -100020.13876
[256,] 9540.54090 4695.80258
[257,] 15504.80685 9540.54090
[258,] 96745.06614 15504.80685
[259,] 4695.80258 96745.06614
[260,] 96745.06614 4695.80258
[261,] 96745.06614 96745.06614
[262,] 11524.03289 96745.06614
[263,] 25124.49007 11524.03289
[264,] 96745.06614 25124.49007
[265,] 4695.80258 96745.06614
[266,] 96745.06614 4695.80258
[267,] 39068.86267 96745.06614
[268,] 96745.06614 39068.86267
[269,] 7614.44532 96745.06614
[270,] 96745.06614 7614.44532
[271,] 96745.06614 96745.06614
[272,] 96745.06614 96745.06614
[273,] 2280.74545 96745.06614
[274,] 4695.80258 2280.74545
[275,] 6856.73328 4695.80258
[276,] 5750.88219 6856.73328
[277,] 4695.80258 5750.88219
[278,] -32509.02284 4695.80258
[279,] -13492.53170 -32509.02284
[280,] 96745.06614 -13492.53170
[281,] 4695.80258 96745.06614
[282,] 13981.74577 4695.80258
[283,] 4695.80258 13981.74577
[284,] 1474.75604 4695.80258
[285,] 4695.80258 1474.75604
[286,] 73875.87396 4695.80258
[287,] 5305.87880 73875.87396
[288,] 96745.06614 5305.87880
[289,] -4616.91129 96745.06614
[290,] 96745.06614 -4616.91129
[291,] 11869.04635 96745.06614
[292,] 17785.58720 11869.04635
[293,] 4695.80258 17785.58720
[294,] -4372.82474 4695.80258
[295,] 90735.14542 -4372.82474
[296,] 4695.80258 90735.14542
[297,] 96745.06614 4695.80258
[298,] 4695.80258 96745.06614
[299,] 6374.25164 4695.80258
[300,] 83624.34209 6374.25164
[301,] -7133.77660 83624.34209
[302,] 4882.29692 -7133.77660
[303,] 96745.06614 4882.29692
[304,] 32437.27291 96745.06614
[305,] 6909.15144 32437.27291
[306,] 96745.06614 6909.15144
[307,] 44212.05015 96745.06614
[308,] 11845.65282 44212.05015
[309,] 96745.06614 11845.65282
[310,] 96745.06614 96745.06614
[311,] -10213.42054 96745.06614
[312,] 96745.06614 -10213.42054
[313,] 66041.66334 96745.06614
[314,] 6370.75503 66041.66334
[315,] 18828.45632 6370.75503
[316,] 5508.08560 18828.45632
[317,] -12520.74431 5508.08560
[318,] 5617.12230 -12520.74431
[319,] 4695.80258 5617.12230
[320,] 4695.80258 4695.80258
[321,] 4695.80258 4695.80258
[322,] 4695.80258 4695.80258
[323,] 96745.06614 4695.80258
[324,] 4695.80258 96745.06614
[325,] 1889.38269 4695.80258
[326,] 96745.06614 1889.38269
[327,] 96745.06614 96745.06614
[328,] 4695.80258 96745.06614
[329,] 4695.80258 4695.80258
[330,] 4695.80258 4695.80258
[331,] 4695.80258 4695.80258
[332,] 4695.80258 4695.80258
[333,] 3953.39132 4695.80258
[334,] 4871.95470 3953.39132
[335,] 3846.84201 4871.95470
[336,] -15467.40013 3846.84201
[337,] 4695.80258 -15467.40013
[338,] 2916.79186 4695.80258
[339,] 10119.59876 2916.79186
[340,] -18672.05287 10119.59876
[341,] 4695.80258 -18672.05287
[342,] 4695.80258 4695.80258
[343,] 4695.80258 4695.80258
[344,] -3140.29155 4695.80258
[345,] 4695.80258 -3140.29155
[346,] -4006.90833 4695.80258
[347,] 4695.80258 -4006.90833
[348,] 4695.80258 4695.80258
[349,] 2493.78740 4695.80258
[350,] 11018.66125 2493.78740
[351,] 3600.11020 11018.66125
[352,] 16051.07917 3600.11020
[353,] 476.58495 16051.07917
[354,] 4695.80258 476.58495
[355,] 4695.80258 4695.80258
[356,] -22396.84628 4695.80258
[357,] 4695.80258 -22396.84628
[358,] 4695.80258 4695.80258
[359,] 4695.80258 4695.80258
[360,] 5450.93162 4695.80258
[361,] 4695.80258 5450.93162
[362,] 4695.80258 4695.80258
[363,] -22742.53307 4695.80258
[364,] 8531.29201 -22742.53307
[365,] 4695.80258 8531.29201
[366,] 4695.80258 4695.80258
[367,] 4695.80258 4695.80258
[368,] 4695.80258 4695.80258
[369,] 5073.19778 4695.80258
[370,] 4695.80258 5073.19778
[371,] 13502.47832 4695.80258
[372,] -4122.43855 13502.47832
[373,] -11364.54599 -4122.43855
[374,] 4695.80258 -11364.54599
[375,] 4695.80258 4695.80258
[376,] -5249.64067 4695.80258
[377,] 10434.70666 -5249.64067
[378,] -127.53011 10434.70666
[379,] 4695.80258 -127.53011
[380,] 4695.80258 4695.80258
[381,] -16852.39645 4695.80258
[382,] 5978.97241 -16852.39645
[383,] 4695.80258 5978.97241
[384,] 4730.19031 4695.80258
[385,] -93173.15808 4730.19031
[386,] 18295.44481 -93173.15808
[387,] -1201.64187 18295.44481
[388,] 4333.03834 -1201.64187
[389,] -87.88601 4333.03834
[390,] 7724.02211 -87.88601
[391,] -17053.27366 7724.02211
[392,] -49494.53350 -17053.27366
[393,] 40647.49443 -49494.53350
[394,] -40778.26494 40647.49443
[395,] -24144.92966 -40778.26494
[396,] -72002.46594 -24144.92966
[397,] -432306.25565 -72002.46594
[398,] -21770.74622 -432306.25565
[399,] -135659.00978 -21770.74622
[400,] -29640.24458 -135659.00978
[401,] 13728.65001 -29640.24458
[402,] -12827.11030 13728.65001
[403,] -29957.78483 -12827.11030
[404,] -42211.45613 -29957.78483
[405,] -6534.19128 -42211.45613
[406,] -8555.89563 -6534.19128
[407,] -5487.64777 -8555.89563
[408,] -8564.90722 -5487.64777
[409,] -99513.27125 -8564.90722
[410,] -14573.70679 -99513.27125
[411,] -157855.32051 -14573.70679
[412,] -13610.02368 -157855.32051
[413,] 68913.23945 -13610.02368
[414,] -188433.40848 68913.23945
[415,] -135107.61599 -188433.40848
[416,] 132680.84269 -135107.61599
[417,] -5736.85380 132680.84269
[418,] -53424.84627 -5736.85380
[419,] -77168.91215 -53424.84627
[420,] -79330.08351 -77168.91215
[421,] -93066.08147 -79330.08351
[422,] -731095.45915 -93066.08147
[423,] -372773.51927 -731095.45915
[424,] -35145.58987 -372773.51927
[425,] -331902.94857 -35145.58987
[426,] -266234.14336 -331902.94857
[427,] -322175.21212 -266234.14336
[428,] -133587.07324 -322175.21212
[429,] -30715.63764 -133587.07324
[430,] -332780.81006 -30715.63764
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 2275662.88080 208308.99237
2 795828.19898 2275662.88080
3 -270987.19826 795828.19898
4 -809715.40781 -270987.19826
5 -22266.82622 -809715.40781
6 800116.23761 -22266.82622
7 -49155.83317 800116.23761
8 -691527.85593 -49155.83317
9 520274.61695 -691527.85593
10 300395.15083 520274.61695
11 657029.06624 300395.15083
12 -223632.43021 657029.06624
13 -214806.51231 -223632.43021
14 -584732.95539 -214806.51231
15 567057.94638 -584732.95539
16 -120685.92990 567057.94638
17 -555208.39746 -120685.92990
18 -692034.64645 -555208.39746
19 -22517.65388 -692034.64645
20 147556.65102 -22517.65388
21 -198905.77113 147556.65102
22 144552.96763 -198905.77113
23 -85472.26929 144552.96763
24 -172515.67738 -85472.26929
25 199630.89796 -172515.67738
26 -806793.29430 199630.89796
27 -91537.29089 -806793.29430
28 493285.69665 -91537.29089
29 102851.24420 493285.69665
30 321503.16649 102851.24420
31 -104258.33967 321503.16649
32 255545.92790 -104258.33967
33 160686.08391 255545.92790
34 -132980.09112 160686.08391
35 -66128.03087 -132980.09112
36 238335.88236 -66128.03087
37 142898.18103 238335.88236
38 44953.18110 142898.18103
39 -135648.29772 44953.18110
40 -143145.37863 -135648.29772
41 -91583.43595 -143145.37863
42 121975.57705 -91583.43595
43 -56834.07852 121975.57705
44 44665.36706 -56834.07852
45 2854.95073 44665.36706
46 -109747.17134 2854.95073
47 -353406.30019 -109747.17134
48 28280.34861 -353406.30019
49 -149859.43355 28280.34861
50 1830.56568 -149859.43355
51 -154893.59582 1830.56568
52 -163816.82868 -154893.59582
53 108069.30970 -163816.82868
54 -398627.52899 108069.30970
55 96414.42623 -398627.52899
56 344119.88167 96414.42623
57 65897.34530 344119.88167
58 -39937.02809 65897.34530
59 -268342.71677 -39937.02809
60 30684.17243 -268342.71677
61 -165191.81686 30684.17243
62 54293.11167 -165191.81686
63 -229559.34975 54293.11167
64 44287.34309 -229559.34975
65 -30222.19008 44287.34309
66 -89225.87105 -30222.19008
67 12463.09493 -89225.87105
68 -62985.93646 12463.09493
69 -110309.30452 -62985.93646
70 -206946.64271 -110309.30452
71 90240.21711 -206946.64271
72 -88273.72819 90240.21711
73 -21316.72675 -88273.72819
74 -30243.60472 -21316.72675
75 -66184.57835 -30243.60472
76 -44934.51284 -66184.57835
77 -101207.15920 -44934.51284
78 -11678.29645 -101207.15920
79 63643.05378 -11678.29645
80 34641.61231 63643.05378
81 -73790.49159 34641.61231
82 -104685.81508 -73790.49159
83 29937.67397 -104685.81508
84 -12725.14358 29937.67397
85 -167364.74107 -12725.14358
86 91819.60560 -167364.74107
87 -38862.44579 91819.60560
88 -91394.42511 -38862.44579
89 -5660.57195 -91394.42511
90 13129.61159 -5660.57195
91 -14033.50196 13129.61159
92 137631.86438 -14033.50196
93 -89777.21321 137631.86438
94 90481.59227 -89777.21321
95 10678.39226 90481.59227
96 19473.69877 10678.39226
97 21519.12241 19473.69877
98 -48344.28494 21519.12241
99 42245.68748 -48344.28494
100 -21816.82929 42245.68748
101 -133542.53336 -21816.82929
102 -108868.71359 -133542.53336
103 -98741.69172 -108868.71359
104 39063.39389 -98741.69172
105 -59967.40542 39063.39389
106 -289570.08055 -59967.40542
107 198107.49533 -289570.08055
108 -86218.57294 198107.49533
109 -66492.21629 -86218.57294
110 83864.87293 -66492.21629
111 -222070.80804 83864.87293
112 -84451.81816 -222070.80804
113 -19997.95631 -84451.81816
114 112369.79821 -19997.95631
115 39988.90315 112369.79821
116 -24563.04824 39988.90315
117 71140.54934 -24563.04824
118 -86911.82797 71140.54934
119 3663.56051 -86911.82797
120 -65922.30668 3663.56051
121 9489.98057 -65922.30668
122 -3418.89374 9489.98057
123 152158.54615 -3418.89374
124 29286.34957 152158.54615
125 -138918.74718 29286.34957
126 4695.80258 -138918.74718
127 5100.36378 4695.80258
128 14817.85613 5100.36378
129 4695.80258 14817.85613
130 4973.42744 4695.80258
131 65208.58145 4973.42744
132 4695.80258 65208.58145
133 5357.41047 4695.80258
134 -1120.95879 5357.41047
135 3882.97829 -1120.95879
136 9075.68325 3882.97829
137 4364.04313 9075.68325
138 444.25026 4364.04313
139 4695.80258 444.25026
140 17633.90268 4695.80258
141 4695.80258 17633.90268
142 -7587.46063 4695.80258
143 5383.56487 -7587.46063
144 14587.67173 5383.56487
145 10257.96358 14587.67173
146 4695.80258 10257.96358
147 5344.44639 4695.80258
148 13569.03592 5344.44639
149 4560.05297 13569.03592
150 -16424.65558 4560.05297
151 -4530.63985 -16424.65558
152 5628.27428 -4530.63985
153 -9949.03909 5628.27428
154 4418.84197 -9949.03909
155 4695.80258 4418.84197
156 4695.80258 4695.80258
157 2967.75607 4695.80258
158 24390.36985 2967.75607
159 4695.80258 24390.36985
160 4695.80258 4695.80258
161 4695.80258 4695.80258
162 4695.80258 4695.80258
163 -5759.11887 4695.80258
164 -11139.17833 -5759.11887
165 4146.91131 -11139.17833
166 4695.80258 4146.91131
167 100101.84568 4695.80258
168 17631.71406 100101.84568
169 4695.80258 17631.71406
170 4695.80258 4695.80258
171 502.88764 4695.80258
172 4695.80258 502.88764
173 -2747.76674 4695.80258
174 4695.80258 -2747.76674
175 5279.74355 4695.80258
176 21636.35992 5279.74355
177 95888.64101 21636.35992
178 4695.80258 95888.64101
179 4695.80258 4695.80258
180 26546.37872 4695.80258
181 69353.44777 26546.37872
182 5194.43625 69353.44777
183 27469.56632 5194.43625
184 82604.48615 27469.56632
185 71777.73429 82604.48615
186 20793.09613 71777.73429
187 505.91678 20793.09613
188 -28647.56861 505.91678
189 4695.80258 -28647.56861
190 10293.81670 4695.80258
191 4695.80258 10293.81670
192 3547.21593 4695.80258
193 4695.80258 3547.21593
194 5099.07171 4695.80258
195 96745.06614 5099.07171
196 -159756.58648 96745.06614
197 96745.06614 -159756.58648
198 96745.06614 96745.06614
199 -11466.83054 96745.06614
200 96745.06614 -11466.83054
201 96745.06614 96745.06614
202 96745.06614 96745.06614
203 48580.16067 96745.06614
204 60559.55258 48580.16067
205 84396.23774 60559.55258
206 92599.63230 84396.23774
207 96745.06614 92599.63230
208 34240.26905 96745.06614
209 74513.23008 34240.26905
210 4695.80258 74513.23008
211 4695.80258 4695.80258
212 25925.93172 4695.80258
213 -24615.08018 25925.93172
214 -11571.51316 -24615.08018
215 5068.79126 -11571.51316
216 5836.24263 5068.79126
217 6492.62189 5836.24263
218 5068.79126 6492.62189
219 15978.94819 5068.79126
220 -14524.58874 15978.94819
221 4695.80258 -14524.58874
222 4882.91782 4695.80258
223 2842.93476 4882.91782
224 4815.88543 2842.93476
225 5637.80160 4815.88543
226 5237.15542 5637.80160
227 4695.80258 5237.15542
228 -7592.79784 4695.80258
229 4695.80258 -7592.79784
230 -18346.43070 4695.80258
231 8642.79244 -18346.43070
232 4695.80258 8642.79244
233 -8722.67173 4695.80258
234 4695.80258 -8722.67173
235 4749.02911 4695.80258
236 4695.80258 4749.02911
237 6664.52600 4695.80258
238 4695.80258 6664.52600
239 4695.80258 4695.80258
240 4695.80258 4695.80258
241 1072.29490 4695.80258
242 4695.80258 1072.29490
243 7849.74472 4695.80258
244 -20413.98966 7849.74472
245 96745.06614 -20413.98966
246 96745.06614 96745.06614
247 -1625.30469 96745.06614
248 96745.06614 -1625.30469
249 96745.06614 96745.06614
250 1105.51459 96745.06614
251 96745.06614 1105.51459
252 -10712.38062 96745.06614
253 -25693.29890 -10712.38062
254 -100020.13876 -25693.29890
255 4695.80258 -100020.13876
256 9540.54090 4695.80258
257 15504.80685 9540.54090
258 96745.06614 15504.80685
259 4695.80258 96745.06614
260 96745.06614 4695.80258
261 96745.06614 96745.06614
262 11524.03289 96745.06614
263 25124.49007 11524.03289
264 96745.06614 25124.49007
265 4695.80258 96745.06614
266 96745.06614 4695.80258
267 39068.86267 96745.06614
268 96745.06614 39068.86267
269 7614.44532 96745.06614
270 96745.06614 7614.44532
271 96745.06614 96745.06614
272 96745.06614 96745.06614
273 2280.74545 96745.06614
274 4695.80258 2280.74545
275 6856.73328 4695.80258
276 5750.88219 6856.73328
277 4695.80258 5750.88219
278 -32509.02284 4695.80258
279 -13492.53170 -32509.02284
280 96745.06614 -13492.53170
281 4695.80258 96745.06614
282 13981.74577 4695.80258
283 4695.80258 13981.74577
284 1474.75604 4695.80258
285 4695.80258 1474.75604
286 73875.87396 4695.80258
287 5305.87880 73875.87396
288 96745.06614 5305.87880
289 -4616.91129 96745.06614
290 96745.06614 -4616.91129
291 11869.04635 96745.06614
292 17785.58720 11869.04635
293 4695.80258 17785.58720
294 -4372.82474 4695.80258
295 90735.14542 -4372.82474
296 4695.80258 90735.14542
297 96745.06614 4695.80258
298 4695.80258 96745.06614
299 6374.25164 4695.80258
300 83624.34209 6374.25164
301 -7133.77660 83624.34209
302 4882.29692 -7133.77660
303 96745.06614 4882.29692
304 32437.27291 96745.06614
305 6909.15144 32437.27291
306 96745.06614 6909.15144
307 44212.05015 96745.06614
308 11845.65282 44212.05015
309 96745.06614 11845.65282
310 96745.06614 96745.06614
311 -10213.42054 96745.06614
312 96745.06614 -10213.42054
313 66041.66334 96745.06614
314 6370.75503 66041.66334
315 18828.45632 6370.75503
316 5508.08560 18828.45632
317 -12520.74431 5508.08560
318 5617.12230 -12520.74431
319 4695.80258 5617.12230
320 4695.80258 4695.80258
321 4695.80258 4695.80258
322 4695.80258 4695.80258
323 96745.06614 4695.80258
324 4695.80258 96745.06614
325 1889.38269 4695.80258
326 96745.06614 1889.38269
327 96745.06614 96745.06614
328 4695.80258 96745.06614
329 4695.80258 4695.80258
330 4695.80258 4695.80258
331 4695.80258 4695.80258
332 4695.80258 4695.80258
333 3953.39132 4695.80258
334 4871.95470 3953.39132
335 3846.84201 4871.95470
336 -15467.40013 3846.84201
337 4695.80258 -15467.40013
338 2916.79186 4695.80258
339 10119.59876 2916.79186
340 -18672.05287 10119.59876
341 4695.80258 -18672.05287
342 4695.80258 4695.80258
343 4695.80258 4695.80258
344 -3140.29155 4695.80258
345 4695.80258 -3140.29155
346 -4006.90833 4695.80258
347 4695.80258 -4006.90833
348 4695.80258 4695.80258
349 2493.78740 4695.80258
350 11018.66125 2493.78740
351 3600.11020 11018.66125
352 16051.07917 3600.11020
353 476.58495 16051.07917
354 4695.80258 476.58495
355 4695.80258 4695.80258
356 -22396.84628 4695.80258
357 4695.80258 -22396.84628
358 4695.80258 4695.80258
359 4695.80258 4695.80258
360 5450.93162 4695.80258
361 4695.80258 5450.93162
362 4695.80258 4695.80258
363 -22742.53307 4695.80258
364 8531.29201 -22742.53307
365 4695.80258 8531.29201
366 4695.80258 4695.80258
367 4695.80258 4695.80258
368 4695.80258 4695.80258
369 5073.19778 4695.80258
370 4695.80258 5073.19778
371 13502.47832 4695.80258
372 -4122.43855 13502.47832
373 -11364.54599 -4122.43855
374 4695.80258 -11364.54599
375 4695.80258 4695.80258
376 -5249.64067 4695.80258
377 10434.70666 -5249.64067
378 -127.53011 10434.70666
379 4695.80258 -127.53011
380 4695.80258 4695.80258
381 -16852.39645 4695.80258
382 5978.97241 -16852.39645
383 4695.80258 5978.97241
384 4730.19031 4695.80258
385 -93173.15808 4730.19031
386 18295.44481 -93173.15808
387 -1201.64187 18295.44481
388 4333.03834 -1201.64187
389 -87.88601 4333.03834
390 7724.02211 -87.88601
391 -17053.27366 7724.02211
392 -49494.53350 -17053.27366
393 40647.49443 -49494.53350
394 -40778.26494 40647.49443
395 -24144.92966 -40778.26494
396 -72002.46594 -24144.92966
397 -432306.25565 -72002.46594
398 -21770.74622 -432306.25565
399 -135659.00978 -21770.74622
400 -29640.24458 -135659.00978
401 13728.65001 -29640.24458
402 -12827.11030 13728.65001
403 -29957.78483 -12827.11030
404 -42211.45613 -29957.78483
405 -6534.19128 -42211.45613
406 -8555.89563 -6534.19128
407 -5487.64777 -8555.89563
408 -8564.90722 -5487.64777
409 -99513.27125 -8564.90722
410 -14573.70679 -99513.27125
411 -157855.32051 -14573.70679
412 -13610.02368 -157855.32051
413 68913.23945 -13610.02368
414 -188433.40848 68913.23945
415 -135107.61599 -188433.40848
416 132680.84269 -135107.61599
417 -5736.85380 132680.84269
418 -53424.84627 -5736.85380
419 -77168.91215 -53424.84627
420 -79330.08351 -77168.91215
421 -93066.08147 -79330.08351
422 -731095.45915 -93066.08147
423 -372773.51927 -731095.45915
424 -35145.58987 -372773.51927
425 -331902.94857 -35145.58987
426 -266234.14336 -331902.94857
427 -322175.21212 -266234.14336
428 -133587.07324 -322175.21212
429 -30715.63764 -133587.07324
430 -332780.81006 -30715.63764
> 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/73qfr1291225443.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/83qfr1291225443.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/93qfr1291225443.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/10w0fd1291225443.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/11zivi1291225443.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/12k1u61291225443.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/13yssx1291225443.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/142t831291225443.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/155tp91291225443.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/169c5f1291225443.tab")
+ }
>
> try(system("convert tmp/1ph0j1291225443.ps tmp/1ph0j1291225443.png",intern=TRUE))
character(0)
> try(system("convert tmp/2ph0j1291225443.ps tmp/2ph0j1291225443.png",intern=TRUE))
character(0)
> try(system("convert tmp/3ph0j1291225443.ps tmp/3ph0j1291225443.png",intern=TRUE))
character(0)
> try(system("convert tmp/40qhm1291225443.ps tmp/40qhm1291225443.png",intern=TRUE))
character(0)
> try(system("convert tmp/50qhm1291225443.ps tmp/50qhm1291225443.png",intern=TRUE))
character(0)
> try(system("convert tmp/6shy71291225443.ps tmp/6shy71291225443.png",intern=TRUE))
character(0)
> try(system("convert tmp/73qfr1291225443.ps tmp/73qfr1291225443.png",intern=TRUE))
character(0)
> try(system("convert tmp/83qfr1291225443.ps tmp/83qfr1291225443.png",intern=TRUE))
character(0)
> try(system("convert tmp/93qfr1291225443.ps tmp/93qfr1291225443.png",intern=TRUE))
character(0)
> try(system("convert tmp/10w0fd1291225443.ps tmp/10w0fd1291225443.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
14.328 2.371 61.470