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(0
+ ,162556
+ ,0
+ ,807
+ ,0
+ ,213118
+ ,0
+ ,6282154
+ ,0
+ ,29790
+ ,0
+ ,444
+ ,0
+ ,81767
+ ,0
+ ,4321023
+ ,0
+ ,87550
+ ,0
+ ,412
+ ,0
+ ,153198
+ ,0
+ ,4111912
+ ,1
+ ,84738
+ ,254214
+ ,428
+ ,2140
+ ,-26007
+ ,-52014
+ ,223193
+ ,0
+ ,54660
+ ,0
+ ,315
+ ,0
+ ,126942
+ ,0
+ ,1491348
+ ,0
+ ,42634
+ ,0
+ ,168
+ ,0
+ ,157214
+ ,0
+ ,1629616
+ ,1
+ ,40949
+ ,122847
+ ,263
+ ,1315
+ ,129352
+ ,258704
+ ,1398893
+ ,0
+ ,45187
+ ,0
+ ,267
+ ,0
+ ,234817
+ ,0
+ ,1926517
+ ,0
+ ,37704
+ ,0
+ ,228
+ ,0
+ ,60448
+ ,0
+ ,983660
+ ,0
+ ,16275
+ ,0
+ ,129
+ ,0
+ ,47818
+ ,0
+ ,1443586
+ ,1
+ ,25830
+ ,77490
+ ,104
+ ,520
+ ,245546
+ ,491092
+ ,1073089
+ ,1
+ ,12679
+ ,38037
+ ,122
+ ,610
+ ,48020
+ ,96040
+ ,984885
+ ,0
+ ,18014
+ ,0
+ ,393
+ ,0
+ ,-1710
+ ,0
+ ,1405225
+ ,1
+ ,43556
+ ,130668
+ ,190
+ ,950
+ ,32648
+ ,65296
+ ,227132
+ ,0
+ ,24811
+ ,0
+ ,280
+ ,0
+ ,95350
+ ,0
+ ,929118
+ ,1
+ ,6575
+ ,19725
+ ,63
+ ,315
+ ,151352
+ ,302704
+ ,1071292
+ ,1
+ ,7123
+ ,21369
+ ,102
+ ,510
+ ,288170
+ ,576340
+ ,638830
+ ,0
+ ,21950
+ ,0
+ ,265
+ ,0
+ ,114337
+ ,0
+ ,856956
+ ,0
+ ,37597
+ ,0
+ ,234
+ ,0
+ ,37884
+ ,0
+ ,992426
+ ,1
+ ,17821
+ ,53463
+ ,277
+ ,1385
+ ,122844
+ ,245688
+ ,444477
+ ,0
+ ,12988
+ ,0
+ ,73
+ ,0
+ ,82340
+ ,0
+ ,857217
+ ,0
+ ,22330
+ ,0
+ ,67
+ ,0
+ ,79801
+ ,0
+ ,711969
+ ,1
+ ,13326
+ ,39978
+ ,103
+ ,515
+ ,165548
+ ,331096
+ ,702380
+ ,1
+ ,16189
+ ,48567
+ ,290
+ ,1450
+ ,116384
+ ,232768
+ ,358589
+ ,1
+ ,7146
+ ,21438
+ ,83
+ ,415
+ ,134028
+ ,268056
+ ,297978
+ ,1
+ ,15824
+ ,47472
+ ,56
+ ,280
+ ,63838
+ ,127676
+ ,585715
+ ,0
+ ,27664
+ ,0
+ ,236
+ ,0
+ ,74996
+ ,0
+ ,657954
+ ,1
+ ,11920
+ ,35760
+ ,73
+ ,365
+ ,31080
+ ,62160
+ ,209458
+ ,1
+ ,8568
+ ,25704
+ ,34
+ ,170
+ ,32168
+ ,64336
+ ,786690
+ ,1
+ ,14416
+ ,43248
+ ,139
+ ,695
+ ,49857
+ ,99714
+ ,439798
+ ,0
+ ,3369
+ ,0
+ ,26
+ ,0
+ ,87161
+ ,0
+ ,688779
+ ,0
+ ,11819
+ ,0
+ ,70
+ ,0
+ ,106113
+ ,0
+ ,574339
+ ,0
+ ,6984
+ ,0
+ ,40
+ ,0
+ ,80570
+ ,0
+ ,741409
+ ,0
+ ,4519
+ ,0
+ ,42
+ ,0
+ ,102129
+ ,0
+ ,597793
+ ,1
+ ,2220
+ ,6660
+ ,12
+ ,60
+ ,301670
+ ,603340
+ ,644190
+ ,1
+ ,18562
+ ,55686
+ ,211
+ ,1055
+ ,102313
+ ,204626
+ ,377934
+ ,1
+ ,10327
+ ,30981
+ ,74
+ ,370
+ ,88577
+ ,177154
+ ,640273
+ ,0
+ ,5336
+ ,0
+ ,80
+ ,0
+ ,112477
+ ,0
+ ,697458
+ ,0
+ ,2365
+ ,0
+ ,83
+ ,0
+ ,191778
+ ,0
+ ,550608
+ ,1
+ ,4069
+ ,12207
+ ,131
+ ,655
+ ,79804
+ ,159608
+ ,207393
+ ,1
+ ,8636
+ ,25908
+ ,203
+ ,1015
+ ,128294
+ ,256588
+ ,301607
+ ,1
+ ,13718
+ ,41154
+ ,56
+ ,280
+ ,96448
+ ,192896
+ ,345783
+ ,1
+ ,4525
+ ,13575
+ ,89
+ ,445
+ ,93811
+ ,187622
+ ,501749
+ ,1
+ ,6869
+ ,20607
+ ,88
+ ,440
+ ,117520
+ ,235040
+ ,379983
+ ,1
+ ,4628
+ ,13884
+ ,39
+ ,195
+ ,69159
+ ,138318
+ ,387475
+ ,0
+ ,3689
+ ,0
+ ,25
+ ,0
+ ,101792
+ ,0
+ ,377305
+ ,0
+ ,4891
+ ,0
+ ,49
+ ,0
+ ,210568
+ ,0
+ ,370837
+ ,0
+ ,7489
+ ,0
+ ,149
+ ,0
+ ,136996
+ ,0
+ ,430866
+ ,1
+ ,4901
+ ,14703
+ ,58
+ ,290
+ ,121920
+ ,243840
+ ,469107
+ ,1
+ ,2284
+ ,6852
+ ,41
+ ,205
+ ,76403
+ ,152806
+ ,194493
+ ,0
+ ,3160
+ ,0
+ ,90
+ ,0
+ ,108094
+ ,0
+ ,530670
+ ,0
+ ,4150
+ ,0
+ ,136
+ ,0
+ ,134759
+ ,0
+ ,518365
+ ,0
+ ,7285
+ ,0
+ ,97
+ ,0
+ ,188873
+ ,0
+ ,491303
+ ,0
+ ,1134
+ ,0
+ ,63
+ ,0
+ ,146216
+ ,0
+ ,527021
+ ,0
+ ,4658
+ ,0
+ ,114
+ ,0
+ ,156608
+ ,0
+ ,233773
+ ,1
+ ,2384
+ ,7152
+ ,77
+ ,385
+ ,61348
+ ,122696
+ ,405972
+ ,1
+ ,3748
+ ,11244
+ ,6
+ ,30
+ ,50350
+ ,100700
+ ,652925
+ ,1
+ ,5371
+ ,16113
+ ,47
+ ,235
+ ,87720
+ ,175440
+ ,446211
+ ,1
+ ,1285
+ ,3855
+ ,51
+ ,255
+ ,99489
+ ,198978
+ ,341340
+ ,0
+ ,9327
+ ,0
+ ,85
+ ,0
+ ,87419
+ ,0
+ ,387699
+ ,0
+ ,5565
+ ,0
+ ,43
+ ,0
+ ,94355
+ ,0
+ ,493408
+ ,1
+ ,1528
+ ,4584
+ ,32
+ ,160
+ ,60326
+ ,120652
+ ,146494
+ ,0
+ ,3122
+ ,0
+ ,25
+ ,0
+ ,94670
+ ,0
+ ,414462
+ ,0
+ ,7561
+ ,0
+ ,77
+ ,0
+ ,82425
+ ,0
+ ,364304
+ ,1
+ ,2675
+ ,8025
+ ,54
+ ,270
+ ,59017
+ ,118034
+ ,355178
+ ,1
+ ,13253
+ ,39759
+ ,251
+ ,1255
+ ,90829
+ ,181658
+ ,357760
+ ,1
+ ,880
+ ,2640
+ ,15
+ ,75
+ ,80791
+ ,161582
+ ,261216
+ ,0
+ ,2053
+ ,0
+ ,44
+ ,0
+ ,100423
+ ,0
+ ,397144
+ ,1
+ ,1424
+ ,4272
+ ,73
+ ,365
+ ,131116
+ ,262232
+ ,374943
+ ,0
+ ,4036
+ ,0
+ ,85
+ ,0
+ ,100269
+ ,0
+ ,424898
+ ,0
+ ,3045
+ ,0
+ ,49
+ ,0
+ ,27330
+ ,0
+ ,202055
+ ,1
+ ,5119
+ ,15357
+ ,38
+ ,190
+ ,39039
+ ,78078
+ ,378525
+ ,1
+ ,1431
+ ,4293
+ ,35
+ ,175
+ ,106885
+ ,213770
+ ,310768
+ ,1
+ ,554
+ ,1662
+ ,9
+ ,45
+ ,79285
+ ,158570
+ ,325738
+ ,1
+ ,1975
+ ,5925
+ ,34
+ ,170
+ ,118881
+ ,237762
+ ,394510
+ ,0
+ ,1765
+ ,0
+ ,20
+ ,0
+ ,77623
+ ,0
+ ,247060
+ ,1
+ ,1012
+ ,3036
+ ,29
+ ,145
+ ,114768
+ ,229536
+ ,368078
+ ,1
+ ,810
+ ,2430
+ ,11
+ ,55
+ ,74015
+ ,148030
+ ,236761
+ ,1
+ ,1280
+ ,3840
+ ,52
+ ,260
+ ,69465
+ ,138930
+ ,312378
+ ,0
+ ,666
+ ,0
+ ,13
+ ,0
+ ,117869
+ ,0
+ ,339836
+ ,1
+ ,1380
+ ,4140
+ ,29
+ ,145
+ ,60982
+ ,121964
+ ,347385
+ ,0
+ ,4677
+ ,0
+ ,66
+ ,0
+ ,90131
+ ,0
+ ,426280
+ ,1
+ ,876
+ ,2628
+ ,33
+ ,165
+ ,138971
+ ,277942
+ ,352850
+ ,1
+ ,814
+ ,2442
+ ,15
+ ,75
+ ,39625
+ ,79250
+ ,301881
+ ,1
+ ,514
+ ,1542
+ ,15
+ ,75
+ ,102725
+ ,205450
+ ,377516
+ ,0
+ ,5692
+ ,0
+ ,68
+ ,0
+ ,64239
+ ,0
+ ,357312
+ ,1
+ ,3642
+ ,10926
+ ,100
+ ,500
+ ,90262
+ ,180524
+ ,458343
+ ,1
+ ,540
+ ,1620
+ ,13
+ ,65
+ ,103960
+ ,207920
+ ,354228
+ ,1
+ ,2099
+ ,6297
+ ,45
+ ,225
+ ,106611
+ ,213222
+ ,308636
+ ,1
+ ,567
+ ,1701
+ ,14
+ ,70
+ ,103345
+ ,206690
+ ,386212
+ ,1
+ ,2001
+ ,6003
+ ,36
+ ,180
+ ,95551
+ ,191102
+ ,393343
+ ,0
+ ,2949
+ ,0
+ ,40
+ ,0
+ ,82903
+ ,0
+ ,378509
+ ,1
+ ,2253
+ ,6759
+ ,68
+ ,340
+ ,63593
+ ,127186
+ ,452469
+ ,0
+ ,6533
+ ,0
+ ,29
+ ,0
+ ,126910
+ ,0
+ ,364839
+ ,1
+ ,1889
+ ,5667
+ ,43
+ ,215
+ ,37527
+ ,75054
+ ,358649
+ ,0
+ ,3055
+ ,0
+ ,30
+ ,0
+ ,60247
+ ,0
+ ,376641
+ ,1
+ ,272
+ ,816
+ ,9
+ ,45
+ ,112995
+ ,225990
+ ,429112
+ ,0
+ ,1414
+ ,0
+ ,22
+ ,0
+ ,70184
+ ,0
+ ,330546
+ ,1
+ ,2564
+ ,7692
+ ,19
+ ,95
+ ,130140
+ ,260280
+ ,403560
+ ,0
+ ,1383
+ ,0
+ ,9
+ ,0
+ ,73221
+ ,0
+ ,317892
+ ,0
+ ,1261
+ ,0
+ ,31
+ ,0
+ ,76114
+ ,0
+ ,307528
+ ,1
+ ,975
+ ,2925
+ ,19
+ ,95
+ ,90534
+ ,181068
+ ,235133
+ ,1
+ ,3366
+ ,10098
+ ,55
+ ,275
+ ,108479
+ ,216958
+ ,299243
+ ,1
+ ,576
+ ,1728
+ ,8
+ ,40
+ ,113761
+ ,227522
+ ,314073
+ ,1
+ ,1686
+ ,5058
+ ,28
+ ,140
+ ,68696
+ ,137392
+ ,368186
+ ,1
+ ,746
+ ,2238
+ ,29
+ ,145
+ ,71561
+ ,143122
+ ,269661
+ ,0
+ ,3192
+ ,0
+ ,48
+ ,0
+ ,59831
+ ,0
+ ,125390
+ ,0
+ ,2045
+ ,0
+ ,16
+ ,0
+ ,97890
+ ,0
+ ,510834
+ ,1
+ ,5702
+ ,17106
+ ,47
+ ,235
+ ,101481
+ ,202962
+ ,321896
+ ,0
+ ,1932
+ ,0
+ ,20
+ ,0
+ ,72954
+ ,0
+ ,249898
+ ,1
+ ,936
+ ,2808
+ ,22
+ ,110
+ ,67939
+ ,135878
+ ,408881
+ ,0
+ ,3437
+ ,0
+ ,33
+ ,0
+ ,48022
+ ,0
+ ,158492
+ ,1
+ ,5131
+ ,15393
+ ,44
+ ,220
+ ,86111
+ ,172222
+ ,292154
+ ,0
+ ,2397
+ ,0
+ ,13
+ ,0
+ ,74020
+ ,0
+ ,289513
+ ,0
+ ,1389
+ ,0
+ ,6
+ ,0
+ ,57530
+ ,0
+ ,378049
+ ,1
+ ,1503
+ ,4509
+ ,35
+ ,175
+ ,56364
+ ,112728
+ ,343466
+ ,1
+ ,402
+ ,1206
+ ,8
+ ,40
+ ,84990
+ ,169980
+ ,332743
+ ,1
+ ,2239
+ ,6717
+ ,17
+ ,85
+ ,88590
+ ,177180
+ ,442882
+ ,0
+ ,2234
+ ,0
+ ,11
+ ,0
+ ,77200
+ ,0
+ ,214215
+ ,1
+ ,837
+ ,2511
+ ,21
+ ,105
+ ,61262
+ ,122524
+ ,315688
+ ,1
+ ,10579
+ ,31737
+ ,92
+ ,460
+ ,110309
+ ,220618
+ ,375195
+ ,1
+ ,875
+ ,2625
+ ,12
+ ,60
+ ,67000
+ ,134000
+ ,334280
+ ,1
+ ,1585
+ ,4755
+ ,112
+ ,560
+ ,93099
+ ,186198
+ ,355864
+ ,0
+ ,1659
+ ,0
+ ,25
+ ,0
+ ,107577
+ ,0
+ ,480382
+ ,0
+ ,2647
+ ,0
+ ,17
+ ,0
+ ,62920
+ ,0
+ ,353058
+ ,0
+ ,3294
+ ,0
+ ,23
+ ,0
+ ,75832
+ ,0
+ ,217193
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,94
+ ,282
+ ,10
+ ,50
+ ,60793
+ ,121586
+ ,314533
+ ,1
+ ,422
+ ,1266
+ ,23
+ ,115
+ ,57935
+ ,115870
+ ,318056
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,34
+ ,102
+ ,7
+ ,35
+ ,60630
+ ,121260
+ ,314353
+ ,1
+ ,1558
+ ,4674
+ ,25
+ ,125
+ ,55637
+ ,111274
+ ,369448
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,43
+ ,129
+ ,20
+ ,100
+ ,60887
+ ,121774
+ ,312846
+ ,1
+ ,645
+ ,1935
+ ,4
+ ,20
+ ,60720
+ ,121440
+ ,312075
+ ,1
+ ,316
+ ,948
+ ,4
+ ,20
+ ,60505
+ ,121010
+ ,315009
+ ,1
+ ,115
+ ,345
+ ,10
+ ,50
+ ,60945
+ ,121890
+ ,318903
+ ,1
+ ,5
+ ,15
+ ,1
+ ,5
+ ,60720
+ ,121440
+ ,314887
+ ,1
+ ,897
+ ,2691
+ ,4
+ ,20
+ ,60720
+ ,121440
+ ,314913
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,389
+ ,1167
+ ,8
+ ,40
+ ,58990
+ ,117980
+ ,325506
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,1002
+ ,3006
+ ,11
+ ,55
+ ,56750
+ ,113500
+ ,298568
+ ,1
+ ,36
+ ,108
+ ,4
+ ,20
+ ,60894
+ ,121788
+ ,315834
+ ,1
+ ,460
+ ,1380
+ ,15
+ ,75
+ ,63346
+ ,126692
+ ,329784
+ ,1
+ ,309
+ ,927
+ ,9
+ ,45
+ ,56535
+ ,113070
+ ,312878
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,9
+ ,27
+ ,7
+ ,35
+ ,60835
+ ,121670
+ ,314987
+ ,1
+ ,271
+ ,813
+ ,2
+ ,10
+ ,60720
+ ,121440
+ ,325249
+ ,1
+ ,14
+ ,42
+ ,0
+ ,0
+ ,61016
+ ,122032
+ ,315877
+ ,1
+ ,520
+ ,1560
+ ,7
+ ,35
+ ,58650
+ ,117300
+ ,291650
+ ,1
+ ,1766
+ ,5298
+ ,46
+ ,230
+ ,60438
+ ,120876
+ ,305959
+ ,1
+ ,0
+ ,0
+ ,5
+ ,25
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,458
+ ,1374
+ ,7
+ ,35
+ ,58625
+ ,117250
+ ,297765
+ ,1
+ ,20
+ ,60
+ ,2
+ ,10
+ ,60938
+ ,121876
+ ,315245
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,98
+ ,294
+ ,2
+ ,10
+ ,61490
+ ,122980
+ ,315236
+ ,1
+ ,405
+ ,1215
+ ,5
+ ,25
+ ,60845
+ ,121690
+ ,336425
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,483
+ ,1449
+ ,7
+ ,35
+ ,60830
+ ,121660
+ ,306268
+ ,1
+ ,454
+ ,1362
+ ,24
+ ,120
+ ,63261
+ ,126522
+ ,302187
+ ,1
+ ,47
+ ,141
+ ,1
+ ,5
+ ,60720
+ ,121440
+ ,314882
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,757
+ ,2271
+ ,18
+ ,90
+ ,45689
+ ,91378
+ ,382712
+ ,1
+ ,4655
+ ,13965
+ ,55
+ ,275
+ ,60720
+ ,121440
+ ,341570
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,36
+ ,108
+ ,3
+ ,15
+ ,61564
+ ,123128
+ ,312412
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,203
+ ,609
+ ,9
+ ,45
+ ,61938
+ ,123876
+ ,309596
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,126
+ ,378
+ ,8
+ ,40
+ ,60951
+ ,121902
+ ,315547
+ ,1
+ ,400
+ ,1200
+ ,113
+ ,565
+ ,60720
+ ,121440
+ ,313267
+ ,0
+ ,71
+ ,0
+ ,0
+ ,0
+ ,60745
+ ,0
+ ,316176
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,972
+ ,2916
+ ,19
+ ,95
+ ,71642
+ ,143284
+ ,359335
+ ,0
+ ,531
+ ,0
+ ,11
+ ,0
+ ,71641
+ ,0
+ ,330068
+ ,1
+ ,2461
+ ,7383
+ ,25
+ ,125
+ ,55792
+ ,111584
+ ,314289
+ ,0
+ ,378
+ ,0
+ ,16
+ ,0
+ ,71873
+ ,0
+ ,297413
+ ,0
+ ,23
+ ,0
+ ,5
+ ,0
+ ,62555
+ ,0
+ ,314806
+ ,0
+ ,638
+ ,0
+ ,11
+ ,0
+ ,60370
+ ,0
+ ,333210
+ ,0
+ ,2300
+ ,0
+ ,23
+ ,0
+ ,64873
+ ,0
+ ,352108
+ ,1
+ ,149
+ ,447
+ ,6
+ ,30
+ ,62041
+ ,124082
+ ,313332
+ ,1
+ ,226
+ ,678
+ ,5
+ ,25
+ ,65745
+ ,131490
+ ,291787
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,275
+ ,825
+ ,7
+ ,35
+ ,59500
+ ,119000
+ ,318745
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,141
+ ,423
+ ,7
+ ,35
+ ,61630
+ ,123260
+ ,315366
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,28
+ ,84
+ ,3
+ ,15
+ ,60890
+ ,121780
+ ,315688
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,4980
+ ,0
+ ,89
+ ,0
+ ,113521
+ ,0
+ ,409642
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,472
+ ,0
+ ,19
+ ,0
+ ,80045
+ ,0
+ ,269587
+ ,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
+ ,203
+ ,0
+ ,12
+ ,0
+ ,50804
+ ,0
+ ,300962
+ ,0
+ ,496
+ ,0
+ ,12
+ ,0
+ ,87390
+ ,0
+ ,325479
+ ,0
+ ,10
+ ,0
+ ,5
+ ,0
+ ,61656
+ ,0
+ ,316155
+ ,0
+ ,63
+ ,0
+ ,2
+ ,0
+ ,65688
+ ,0
+ ,318574
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,1136
+ ,0
+ ,26
+ ,0
+ ,48522
+ ,0
+ ,343613
+ ,0
+ ,265
+ ,0
+ ,3
+ ,0
+ ,60720
+ ,0
+ ,306948
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,267
+ ,801
+ ,11
+ ,55
+ ,57640
+ ,115280
+ ,330059
+ ,1
+ ,474
+ ,1422
+ ,10
+ ,50
+ ,61977
+ ,123954
+ ,288985
+ ,1
+ ,534
+ ,1602
+ ,5
+ ,25
+ ,62620
+ ,125240
+ ,304485
+ ,1
+ ,0
+ ,0
+ ,2
+ ,10
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,15
+ ,45
+ ,6
+ ,30
+ ,60831
+ ,121662
+ ,315688
+ ,1
+ ,397
+ ,1191
+ ,7
+ ,35
+ ,60646
+ ,121292
+ ,317736
+ ,1
+ ,0
+ ,0
+ ,2
+ ,10
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,1866
+ ,5598
+ ,28
+ ,140
+ ,56225
+ ,112450
+ ,322331
+ ,1
+ ,288
+ ,864
+ ,3
+ ,15
+ ,60510
+ ,121020
+ ,296656
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,3
+ ,9
+ ,1
+ ,5
+ ,60698
+ ,121396
+ ,315354
+ ,1
+ ,468
+ ,1404
+ ,20
+ ,100
+ ,60720
+ ,121440
+ ,312161
+ ,1
+ ,20
+ ,60
+ ,1
+ ,5
+ ,60805
+ ,121610
+ ,315576
+ ,1
+ ,278
+ ,834
+ ,22
+ ,110
+ ,61404
+ ,122808
+ ,314922
+ ,1
+ ,61
+ ,183
+ ,9
+ ,45
+ ,60720
+ ,121440
+ ,314551
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,192
+ ,576
+ ,2
+ ,10
+ ,65276
+ ,130552
+ ,312339
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,317
+ ,951
+ ,7
+ ,35
+ ,63915
+ ,127830
+ ,298700
+ ,1
+ ,738
+ ,2214
+ ,9
+ ,45
+ ,60720
+ ,121440
+ ,321376
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,368
+ ,1104
+ ,13
+ ,65
+ ,61686
+ ,123372
+ ,303230
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,2
+ ,6
+ ,0
+ ,0
+ ,60743
+ ,121486
+ ,315487
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,53
+ ,159
+ ,6
+ ,30
+ ,60349
+ ,120698
+ ,315793
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,94
+ ,282
+ ,3
+ ,15
+ ,61360
+ ,122720
+ ,312887
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,24
+ ,72
+ ,7
+ ,35
+ ,59818
+ ,119636
+ ,315637
+ ,1
+ ,2332
+ ,6996
+ ,2
+ ,10
+ ,72680
+ ,145360
+ ,324385
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,131
+ ,393
+ ,15
+ ,75
+ ,61808
+ ,123616
+ ,308989
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,206
+ ,618
+ ,9
+ ,45
+ ,53110
+ ,106220
+ ,296702
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,167
+ ,501
+ ,1
+ ,5
+ ,64245
+ ,128490
+ ,307322
+ ,1
+ ,622
+ ,1866
+ ,38
+ ,190
+ ,73007
+ ,146014
+ ,304376
+ ,1
+ ,2328
+ ,6984
+ ,57
+ ,285
+ ,82732
+ ,165464
+ ,253588
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,365
+ ,1095
+ ,7
+ ,35
+ ,54820
+ ,109640
+ ,309560
+ ,1
+ ,364
+ ,1092
+ ,26
+ ,130
+ ,47705
+ ,95410
+ ,298466
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,226
+ ,678
+ ,13
+ ,65
+ ,72835
+ ,145670
+ ,343929
+ ,1
+ ,307
+ ,921
+ ,10
+ ,50
+ ,58856
+ ,117712
+ ,331955
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,188
+ ,564
+ ,9
+ ,45
+ ,77655
+ ,155310
+ ,381180
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,138
+ ,414
+ ,26
+ ,130
+ ,69817
+ ,139634
+ ,331420
+ ,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
+ ,125
+ ,375
+ ,19
+ ,95
+ ,60798
+ ,121596
+ ,310201
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,282
+ ,846
+ ,12
+ ,60
+ ,62452
+ ,124904
+ ,320016
+ ,1
+ ,335
+ ,1005
+ ,23
+ ,115
+ ,64175
+ ,128350
+ ,320398
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,0
+ ,1324
+ ,0
+ ,29
+ ,0
+ ,67440
+ ,0
+ ,291841
+ ,1
+ ,176
+ ,528
+ ,8
+ ,40
+ ,68136
+ ,136272
+ ,310670
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,249
+ ,747
+ ,26
+ ,130
+ ,56726
+ ,113452
+ ,313491
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,333
+ ,999
+ ,9
+ ,45
+ ,70811
+ ,141622
+ ,331323
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,0
+ ,601
+ ,0
+ ,5
+ ,0
+ ,60720
+ ,0
+ ,319210
+ ,1
+ ,30
+ ,90
+ ,3
+ ,15
+ ,62045
+ ,124090
+ ,318098
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,249
+ ,747
+ ,13
+ ,65
+ ,54323
+ ,108646
+ ,292754
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,165
+ ,495
+ ,12
+ ,60
+ ,62841
+ ,125682
+ ,325176
+ ,1
+ ,453
+ ,1359
+ ,19
+ ,95
+ ,81125
+ ,162250
+ ,365959
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,53
+ ,159
+ ,10
+ ,50
+ ,59506
+ ,119012
+ ,302409
+ ,0
+ ,382
+ ,0
+ ,9
+ ,0
+ ,59365
+ ,0
+ ,340968
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,9
+ ,45
+ ,60720
+ ,121440
+ ,315380
+ ,0
+ ,30
+ ,0
+ ,4
+ ,0
+ ,60798
+ ,0
+ ,313164
+ ,1
+ ,290
+ ,870
+ ,1
+ ,5
+ ,58790
+ ,117580
+ ,301164
+ ,1
+ ,0
+ ,0
+ ,1
+ ,5
+ ,60720
+ ,121440
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,366
+ ,1098
+ ,14
+ ,70
+ ,61808
+ ,123616
+ ,344425
+ ,1
+ ,2
+ ,6
+ ,12
+ ,60
+ ,60735
+ ,121470
+ ,315394
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,209
+ ,0
+ ,19
+ ,0
+ ,64016
+ ,0
+ ,316647
+ ,1
+ ,384
+ ,1152
+ ,17
+ ,85
+ ,54683
+ ,109366
+ ,309836
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,365
+ ,1095
+ ,32
+ ,160
+ ,87192
+ ,174384
+ ,346611
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,49
+ ,0
+ ,14
+ ,0
+ ,64107
+ ,0
+ ,322031
+ ,1
+ ,3
+ ,9
+ ,8
+ ,40
+ ,60761
+ ,121522
+ ,315656
+ ,1
+ ,133
+ ,399
+ ,4
+ ,20
+ ,65990
+ ,131980
+ ,339445
+ ,1
+ ,32
+ ,96
+ ,0
+ ,0
+ ,59988
+ ,119976
+ ,314964
+ ,1
+ ,368
+ ,1104
+ ,20
+ ,100
+ ,61167
+ ,122334
+ ,297141
+ ,1
+ ,1
+ ,3
+ ,5
+ ,25
+ ,60719
+ ,121438
+ ,315372
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,22
+ ,66
+ ,1
+ ,5
+ ,60722
+ ,121444
+ ,312502
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,96
+ ,288
+ ,4
+ ,20
+ ,60379
+ ,120758
+ ,313729
+ ,1
+ ,1
+ ,3
+ ,1
+ ,5
+ ,60727
+ ,121454
+ ,315388
+ ,1
+ ,314
+ ,942
+ ,4
+ ,20
+ ,60720
+ ,121440
+ ,315371
+ ,1
+ ,844
+ ,2532
+ ,20
+ ,100
+ ,60925
+ ,121850
+ ,296139
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,26
+ ,78
+ ,1
+ ,5
+ ,60896
+ ,121792
+ ,313880
+ ,1
+ ,125
+ ,375
+ ,10
+ ,50
+ ,59734
+ ,119468
+ ,317698
+ ,1
+ ,304
+ ,912
+ ,12
+ ,60
+ ,62969
+ ,125938
+ ,295580
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,621
+ ,1863
+ ,13
+ ,65
+ ,60720
+ ,121440
+ ,308256
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,119
+ ,357
+ ,3
+ ,15
+ ,59118
+ ,118236
+ ,303677
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,1595
+ ,4785
+ ,10
+ ,50
+ ,60720
+ ,121440
+ ,319369
+ ,1
+ ,312
+ ,936
+ ,3
+ ,15
+ ,58598
+ ,117196
+ ,318690
+ ,1
+ ,60
+ ,180
+ ,7
+ ,35
+ ,61124
+ ,122248
+ ,314049
+ ,1
+ ,587
+ ,1761
+ ,10
+ ,50
+ ,59595
+ ,119190
+ ,325699
+ ,1
+ ,135
+ ,405
+ ,1
+ ,5
+ ,62065
+ ,124130
+ ,314210
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,514
+ ,1542
+ ,15
+ ,75
+ ,78780
+ ,157560
+ ,322378
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,1
+ ,3
+ ,4
+ ,20
+ ,60722
+ ,121444
+ ,315398
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,0
+ ,1763
+ ,0
+ ,28
+ ,0
+ ,61600
+ ,0
+ ,308336
+ ,1
+ ,180
+ ,540
+ ,9
+ ,45
+ ,59635
+ ,119270
+ ,316386
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,218
+ ,654
+ ,7
+ ,35
+ ,60720
+ ,121440
+ ,315553
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,448
+ ,1344
+ ,7
+ ,35
+ ,59781
+ ,119562
+ ,323361
+ ,1
+ ,227
+ ,681
+ ,7
+ ,35
+ ,76644
+ ,153288
+ ,336639
+ ,1
+ ,174
+ ,522
+ ,3
+ ,15
+ ,64820
+ ,129640
+ ,307424
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,121
+ ,363
+ ,11
+ ,55
+ ,56178
+ ,112356
+ ,295370
+ ,1
+ ,607
+ ,1821
+ ,7
+ ,35
+ ,60436
+ ,120872
+ ,322340
+ ,1
+ ,2212
+ ,6636
+ ,10
+ ,50
+ ,60720
+ ,121440
+ ,319864
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,530
+ ,1590
+ ,18
+ ,90
+ ,73433
+ ,146866
+ ,317291
+ ,1
+ ,571
+ ,1713
+ ,14
+ ,70
+ ,41477
+ ,82954
+ ,280398
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,121440
+ ,315380
+ ,1
+ ,78
+ ,234
+ ,12
+ ,60
+ ,62700
+ ,125400
+ ,317330
+ ,1
+ ,2489
+ ,7467
+ ,29
+ ,145
+ ,67804
+ ,135608
+ ,238125
+ ,1
+ ,131
+ ,393
+ ,3
+ ,15
+ ,59661
+ ,119322
+ ,327071
+ ,1
+ ,923
+ ,2769
+ ,6
+ ,30
+ ,58620
+ ,117240
+ ,309038
+ ,1
+ ,72
+ ,216
+ ,3
+ ,15
+ ,60398
+ ,120796
+ ,314210
+ ,1
+ ,572
+ ,1716
+ ,8
+ ,40
+ ,58580
+ ,117160
+ ,307930
+ ,1
+ ,397
+ ,1191
+ ,10
+ ,50
+ ,62710
+ ,125420
+ ,322327
+ ,1
+ ,450
+ ,1350
+ ,6
+ ,30
+ ,59325
+ ,118650
+ ,292136
+ ,1
+ ,622
+ ,1866
+ ,8
+ ,40
+ ,60950
+ ,121900
+ ,263276
+ ,1
+ ,694
+ ,2082
+ ,6
+ ,30
+ ,68060
+ ,136120
+ ,367655
+ ,0
+ ,3425
+ ,0
+ ,9
+ ,0
+ ,83620
+ ,0
+ ,283910
+ ,1
+ ,562
+ ,1686
+ ,8
+ ,40
+ ,58456
+ ,116912
+ ,283587
+ ,1
+ ,4917
+ ,14751
+ ,26
+ ,130
+ ,52811
+ ,105622
+ ,243650
+ ,0
+ ,1442
+ ,0
+ ,239
+ ,0
+ ,121173
+ ,0
+ ,438493
+ ,1
+ ,529
+ ,1587
+ ,7
+ ,35
+ ,63870
+ ,127740
+ ,296261
+ ,0
+ ,2126
+ ,0
+ ,41
+ ,0
+ ,21001
+ ,0
+ ,230621
+ ,1
+ ,1061
+ ,3183
+ ,3
+ ,15
+ ,70415
+ ,140830
+ ,304252
+ ,1
+ ,776
+ ,2328
+ ,8
+ ,40
+ ,64230
+ ,128460
+ ,333505
+ ,1
+ ,611
+ ,1833
+ ,6
+ ,30
+ ,59190
+ ,118380
+ ,296919
+ ,0
+ ,1526
+ ,0
+ ,21
+ ,0
+ ,69351
+ ,0
+ ,278990
+ ,1
+ ,592
+ ,1776
+ ,7
+ ,35
+ ,64270
+ ,128540
+ ,276898
+ ,1
+ ,1182
+ ,3546
+ ,11
+ ,55
+ ,70694
+ ,141388
+ ,327007
+ ,1
+ ,621
+ ,1863
+ ,11
+ ,55
+ ,68005
+ ,136010
+ ,317046
+ ,1
+ ,989
+ ,2967
+ ,12
+ ,60
+ ,58930
+ ,117860
+ ,304555
+ ,1
+ ,438
+ ,1314
+ ,9
+ ,45
+ ,58320
+ ,116640
+ ,298096
+ ,1
+ ,726
+ ,2178
+ ,3
+ ,15
+ ,69980
+ ,139960
+ ,231861
+ ,1
+ ,1303
+ ,3909
+ ,57
+ ,285
+ ,69863
+ ,139726
+ ,309422
+ ,0
+ ,7419
+ ,0
+ ,21
+ ,0
+ ,63255
+ ,0
+ ,286963
+ ,0
+ ,1164
+ ,0
+ ,15
+ ,0
+ ,57320
+ ,0
+ ,269753
+ ,0
+ ,3310
+ ,0
+ ,32
+ ,0
+ ,75230
+ ,0
+ ,448243
+ ,1
+ ,1920
+ ,5760
+ ,11
+ ,55
+ ,79420
+ ,158840
+ ,165404
+ ,1
+ ,965
+ ,2895
+ ,2
+ ,10
+ ,73490
+ ,146980
+ ,204325
+ ,1
+ ,3256
+ ,9768
+ ,23
+ ,115
+ ,35250
+ ,70500
+ ,407159
+ ,0
+ ,1135
+ ,0
+ ,20
+ ,0
+ ,62285
+ ,0
+ ,290476
+ ,1
+ ,1270
+ ,3810
+ ,24
+ ,120
+ ,69206
+ ,138412
+ ,275311
+ ,1
+ ,661
+ ,1983
+ ,1
+ ,5
+ ,65920
+ ,131840
+ ,246541
+ ,1
+ ,1013
+ ,3039
+ ,1
+ ,5
+ ,69770
+ ,139540
+ ,253468
+ ,1
+ ,2844
+ ,8532
+ ,74
+ ,370
+ ,72683
+ ,145366
+ ,240897
+ ,0
+ ,11528
+ ,0
+ ,68
+ ,0
+ ,-14545
+ ,0
+ ,-83265
+ ,1
+ ,6526
+ ,19578
+ ,20
+ ,100
+ ,55830
+ ,111660
+ ,-42143
+ ,1
+ ,2264
+ ,6792
+ ,20
+ ,100
+ ,55174
+ ,110348
+ ,272713
+ ,0
+ ,5109
+ ,0
+ ,82
+ ,0
+ ,67038
+ ,0
+ ,215362
+ ,1
+ ,3999
+ ,11997
+ ,21
+ ,105
+ ,51252
+ ,102504
+ ,42754
+ ,1
+ ,35624
+ ,106872
+ ,244
+ ,1220
+ ,157278
+ ,314556
+ ,306275
+ ,1
+ ,9252
+ ,27756
+ ,32
+ ,160
+ ,79510
+ ,159020
+ ,253537
+ ,1
+ ,15236
+ ,45708
+ ,86
+ ,430
+ ,77440
+ ,154880
+ ,372631
+ ,1
+ ,18073
+ ,54219
+ ,69
+ ,345
+ ,27284
+ ,54568
+ ,-7170)
+ ,dim=c(8
+ ,431)
+ ,dimnames=list(c('Group'
+ ,'Costs'
+ ,'Costs_p'
+ ,'Orders'
+ ,'Orders_p'
+ ,'Dividends'
+ ,'Dividends_p'
+ ,'Wealth
')
+ ,1:431))
> y <- array(NA,dim=c(8,431),dimnames=list(c('Group','Costs','Costs_p','Orders','Orders_p','Dividends','Dividends_p','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\r Group Costs Costs_p Orders Orders_p Dividends Dividends_p
1 6282154 0 162556 0 807 0 213118 0
2 4321023 0 29790 0 444 0 81767 0
3 4111912 0 87550 0 412 0 153198 0
4 223193 1 84738 254214 428 2140 -26007 -52014
5 1491348 0 54660 0 315 0 126942 0
6 1629616 0 42634 0 168 0 157214 0
7 1398893 1 40949 122847 263 1315 129352 258704
8 1926517 0 45187 0 267 0 234817 0
9 983660 0 37704 0 228 0 60448 0
10 1443586 0 16275 0 129 0 47818 0
11 1073089 1 25830 77490 104 520 245546 491092
12 984885 1 12679 38037 122 610 48020 96040
13 1405225 0 18014 0 393 0 -1710 0
14 227132 1 43556 130668 190 950 32648 65296
15 929118 0 24811 0 280 0 95350 0
16 1071292 1 6575 19725 63 315 151352 302704
17 638830 1 7123 21369 102 510 288170 576340
18 856956 0 21950 0 265 0 114337 0
19 992426 0 37597 0 234 0 37884 0
20 444477 1 17821 53463 277 1385 122844 245688
21 857217 0 12988 0 73 0 82340 0
22 711969 0 22330 0 67 0 79801 0
23 702380 1 13326 39978 103 515 165548 331096
24 358589 1 16189 48567 290 1450 116384 232768
25 297978 1 7146 21438 83 415 134028 268056
26 585715 1 15824 47472 56 280 63838 127676
27 657954 0 27664 0 236 0 74996 0
28 209458 1 11920 35760 73 365 31080 62160
29 786690 1 8568 25704 34 170 32168 64336
30 439798 1 14416 43248 139 695 49857 99714
31 688779 0 3369 0 26 0 87161 0
32 574339 0 11819 0 70 0 106113 0
33 741409 0 6984 0 40 0 80570 0
34 597793 0 4519 0 42 0 102129 0
35 644190 1 2220 6660 12 60 301670 603340
36 377934 1 18562 55686 211 1055 102313 204626
37 640273 1 10327 30981 74 370 88577 177154
38 697458 0 5336 0 80 0 112477 0
39 550608 0 2365 0 83 0 191778 0
40 207393 1 4069 12207 131 655 79804 159608
41 301607 1 8636 25908 203 1015 128294 256588
42 345783 1 13718 41154 56 280 96448 192896
43 501749 1 4525 13575 89 445 93811 187622
44 379983 1 6869 20607 88 440 117520 235040
45 387475 1 4628 13884 39 195 69159 138318
46 377305 0 3689 0 25 0 101792 0
47 370837 0 4891 0 49 0 210568 0
48 430866 0 7489 0 149 0 136996 0
49 469107 1 4901 14703 58 290 121920 243840
50 194493 1 2284 6852 41 205 76403 152806
51 530670 0 3160 0 90 0 108094 0
52 518365 0 4150 0 136 0 134759 0
53 491303 0 7285 0 97 0 188873 0
54 527021 0 1134 0 63 0 146216 0
55 233773 0 4658 0 114 0 156608 0
56 405972 1 2384 7152 77 385 61348 122696
57 652925 1 3748 11244 6 30 50350 100700
58 446211 1 5371 16113 47 235 87720 175440
59 341340 1 1285 3855 51 255 99489 198978
60 387699 0 9327 0 85 0 87419 0
61 493408 0 5565 0 43 0 94355 0
62 146494 1 1528 4584 32 160 60326 120652
63 414462 0 3122 0 25 0 94670 0
64 364304 0 7561 0 77 0 82425 0
65 355178 1 2675 8025 54 270 59017 118034
66 357760 1 13253 39759 251 1255 90829 181658
67 261216 1 880 2640 15 75 80791 161582
68 397144 0 2053 0 44 0 100423 0
69 374943 1 1424 4272 73 365 131116 262232
70 424898 0 4036 0 85 0 100269 0
71 202055 0 3045 0 49 0 27330 0
72 378525 1 5119 15357 38 190 39039 78078
73 310768 1 1431 4293 35 175 106885 213770
74 325738 1 554 1662 9 45 79285 158570
75 394510 1 1975 5925 34 170 118881 237762
76 247060 0 1765 0 20 0 77623 0
77 368078 1 1012 3036 29 145 114768 229536
78 236761 1 810 2430 11 55 74015 148030
79 312378 1 1280 3840 52 260 69465 138930
80 339836 0 666 0 13 0 117869 0
81 347385 1 1380 4140 29 145 60982 121964
82 426280 0 4677 0 66 0 90131 0
83 352850 1 876 2628 33 165 138971 277942
84 301881 1 814 2442 15 75 39625 79250
85 377516 1 514 1542 15 75 102725 205450
86 357312 0 5692 0 68 0 64239 0
87 458343 1 3642 10926 100 500 90262 180524
88 354228 1 540 1620 13 65 103960 207920
89 308636 1 2099 6297 45 225 106611 213222
90 386212 1 567 1701 14 70 103345 206690
91 393343 1 2001 6003 36 180 95551 191102
92 378509 0 2949 0 40 0 82903 0
93 452469 1 2253 6759 68 340 63593 127186
94 364839 0 6533 0 29 0 126910 0
95 358649 1 1889 5667 43 215 37527 75054
96 376641 0 3055 0 30 0 60247 0
97 429112 1 272 816 9 45 112995 225990
98 330546 0 1414 0 22 0 70184 0
99 403560 1 2564 7692 19 95 130140 260280
100 317892 0 1383 0 9 0 73221 0
101 307528 0 1261 0 31 0 76114 0
102 235133 1 975 2925 19 95 90534 181068
103 299243 1 3366 10098 55 275 108479 216958
104 314073 1 576 1728 8 40 113761 227522
105 368186 1 1686 5058 28 140 68696 137392
106 269661 1 746 2238 29 145 71561 143122
107 125390 0 3192 0 48 0 59831 0
108 510834 0 2045 0 16 0 97890 0
109 321896 1 5702 17106 47 235 101481 202962
110 249898 0 1932 0 20 0 72954 0
111 408881 1 936 2808 22 110 67939 135878
112 158492 0 3437 0 33 0 48022 0
113 292154 1 5131 15393 44 220 86111 172222
114 289513 0 2397 0 13 0 74020 0
115 378049 0 1389 0 6 0 57530 0
116 343466 1 1503 4509 35 175 56364 112728
117 332743 1 402 1206 8 40 84990 169980
118 442882 1 2239 6717 17 85 88590 177180
119 214215 0 2234 0 11 0 77200 0
120 315688 1 837 2511 21 105 61262 122524
121 375195 1 10579 31737 92 460 110309 220618
122 334280 1 875 2625 12 60 67000 134000
123 355864 1 1585 4755 112 560 93099 186198
124 480382 0 1659 0 25 0 107577 0
125 353058 0 2647 0 17 0 62920 0
126 217193 0 3294 0 23 0 75832 0
127 315380 1 0 0 0 0 60720 121440
128 314533 1 94 282 10 50 60793 121586
129 318056 1 422 1266 23 115 57935 115870
130 315380 1 0 0 0 0 60720 121440
131 314353 1 34 102 7 35 60630 121260
132 369448 1 1558 4674 25 125 55637 111274
133 315380 1 0 0 0 0 60720 121440
134 312846 1 43 129 20 100 60887 121774
135 312075 1 645 1935 4 20 60720 121440
136 315009 1 316 948 4 20 60505 121010
137 318903 1 115 345 10 50 60945 121890
138 314887 1 5 15 1 5 60720 121440
139 314913 1 897 2691 4 20 60720 121440
140 315380 1 0 0 0 0 60720 121440
141 325506 1 389 1167 8 40 58990 117980
142 315380 1 0 0 0 0 60720 121440
143 298568 1 1002 3006 11 55 56750 113500
144 315834 1 36 108 4 20 60894 121788
145 329784 1 460 1380 15 75 63346 126692
146 312878 1 309 927 9 45 56535 113070
147 315380 1 0 0 0 0 60720 121440
148 314987 1 9 27 7 35 60835 121670
149 325249 1 271 813 2 10 60720 121440
150 315877 1 14 42 0 0 61016 122032
151 291650 1 520 1560 7 35 58650 117300
152 305959 1 1766 5298 46 230 60438 120876
153 315380 1 0 0 5 25 60720 121440
154 297765 1 458 1374 7 35 58625 117250
155 315245 1 20 60 2 10 60938 121876
156 315380 1 0 0 0 0 60720 121440
157 315380 1 0 0 0 0 60720 121440
158 315236 1 98 294 2 10 61490 122980
159 336425 1 405 1215 5 25 60845 121690
160 315380 1 0 0 0 0 60720 121440
161 315380 1 0 0 0 0 60720 121440
162 315380 1 0 0 0 0 60720 121440
163 315380 1 0 0 0 0 60720 121440
164 306268 1 483 1449 7 35 60830 121660
165 302187 1 454 1362 24 120 63261 126522
166 314882 1 47 141 1 5 60720 121440
167 315380 1 0 0 0 0 60720 121440
168 382712 1 757 2271 18 90 45689 91378
169 341570 1 4655 13965 55 275 60720 121440
170 315380 1 0 0 0 0 60720 121440
171 315380 1 0 0 0 0 60720 121440
172 312412 1 36 108 3 15 61564 123128
173 315380 1 0 0 0 0 60720 121440
174 309596 1 203 609 9 45 61938 123876
175 315380 1 0 0 0 0 60720 121440
176 315547 1 126 378 8 40 60951 121902
177 313267 1 400 1200 113 565 60720 121440
178 316176 0 71 0 0 0 60745 0
179 315380 1 0 0 0 0 60720 121440
180 315380 1 0 0 0 0 60720 121440
181 359335 1 972 2916 19 95 71642 143284
182 330068 0 531 0 11 0 71641 0
183 314289 1 2461 7383 25 125 55792 111584
184 297413 0 378 0 16 0 71873 0
185 314806 0 23 0 5 0 62555 0
186 333210 0 638 0 11 0 60370 0
187 352108 0 2300 0 23 0 64873 0
188 313332 1 149 447 6 30 62041 124082
189 291787 1 226 678 5 25 65745 131490
190 315380 1 0 0 0 0 60720 121440
191 318745 1 275 825 7 35 59500 119000
192 315380 1 0 0 0 0 60720 121440
193 315366 1 141 423 7 35 61630 123260
194 315380 1 0 0 0 0 60720 121440
195 315688 1 28 84 3 15 60890 121780
196 315380 0 0 0 0 0 60720 0
197 409642 0 4980 0 89 0 113521 0
198 315380 0 0 0 0 0 60720 0
199 315380 0 0 0 0 0 60720 0
200 269587 0 472 0 19 0 80045 0
201 315380 0 0 0 0 0 60720 0
202 315380 0 0 0 0 0 60720 0
203 315380 0 0 0 0 0 60720 0
204 300962 0 203 0 12 0 50804 0
205 325479 0 496 0 12 0 87390 0
206 316155 0 10 0 5 0 61656 0
207 318574 0 63 0 2 0 65688 0
208 315380 0 0 0 0 0 60720 0
209 343613 0 1136 0 26 0 48522 0
210 306948 0 265 0 3 0 60720 0
211 315380 1 0 0 0 0 60720 121440
212 315380 1 0 0 0 0 60720 121440
213 330059 1 267 801 11 55 57640 115280
214 288985 1 474 1422 10 50 61977 123954
215 304485 1 534 1602 5 25 62620 125240
216 315380 1 0 0 2 10 60720 121440
217 315688 1 15 45 6 30 60831 121662
218 317736 1 397 1191 7 35 60646 121292
219 315380 1 0 0 2 10 60720 121440
220 322331 1 1866 5598 28 140 56225 112450
221 296656 1 288 864 3 15 60510 121020
222 315380 1 0 0 0 0 60720 121440
223 315354 1 3 9 1 5 60698 121396
224 312161 1 468 1404 20 100 60720 121440
225 315576 1 20 60 1 5 60805 121610
226 314922 1 278 834 22 110 61404 122808
227 314551 1 61 183 9 45 60720 121440
228 315380 1 0 0 0 0 60720 121440
229 312339 1 192 576 2 10 65276 130552
230 315380 1 0 0 0 0 60720 121440
231 298700 1 317 951 7 35 63915 127830
232 321376 1 738 2214 9 45 60720 121440
233 315380 1 0 0 0 0 60720 121440
234 303230 1 368 1104 13 65 61686 123372
235 315380 1 0 0 0 0 60720 121440
236 315487 1 2 6 0 0 60743 121486
237 315380 1 0 0 0 0 60720 121440
238 315793 1 53 159 6 30 60349 120698
239 315380 1 0 0 0 0 60720 121440
240 315380 1 0 0 0 0 60720 121440
241 315380 1 0 0 0 0 60720 121440
242 312887 1 94 282 3 15 61360 122720
243 315380 1 0 0 0 0 60720 121440
244 315637 1 24 72 7 35 59818 119636
245 324385 1 2332 6996 2 10 72680 145360
246 315380 0 0 0 0 0 60720 0
247 315380 0 0 0 0 0 60720 0
248 308989 1 131 393 15 75 61808 123616
249 315380 0 0 0 0 0 60720 0
250 315380 0 0 0 0 0 60720 0
251 296702 1 206 618 9 45 53110 106220
252 315380 0 0 0 0 0 60720 0
253 307322 1 167 501 1 5 64245 128490
254 304376 1 622 1866 38 190 73007 146014
255 253588 1 2328 6984 57 285 82732 165464
256 315380 1 0 0 0 0 60720 121440
257 309560 1 365 1095 7 35 54820 109640
258 298466 1 364 1092 26 130 47705 95410
259 315380 0 0 0 0 0 60720 0
260 315380 1 0 0 0 0 60720 121440
261 315380 0 0 0 0 0 60720 0
262 315380 0 0 0 0 0 60720 0
263 343929 1 226 678 13 65 72835 145670
264 331955 1 307 921 10 50 58856 117712
265 315380 0 0 0 0 0 60720 0
266 315380 1 0 0 0 0 60720 121440
267 315380 0 0 0 0 0 60720 0
268 381180 1 188 564 9 45 77655 155310
269 315380 0 0 0 0 0 60720 0
270 331420 1 138 414 26 130 69817 139634
271 315380 0 0 0 0 0 60720 0
272 315380 0 0 0 0 0 60720 0
273 315380 0 0 0 0 0 60720 0
274 310201 1 125 375 19 95 60798 121596
275 315380 1 0 0 0 0 60720 121440
276 320016 1 282 846 12 60 62452 124904
277 320398 1 335 1005 23 115 64175 128350
278 315380 1 0 0 0 0 60720 121440
279 291841 0 1324 0 29 0 67440 0
280 310670 1 176 528 8 40 68136 136272
281 315380 0 0 0 0 0 60720 0
282 315380 1 0 0 0 0 60720 121440
283 313491 1 249 747 26 130 56726 113452
284 315380 1 0 0 0 0 60720 121440
285 331323 1 333 999 9 45 70811 141622
286 315380 1 0 0 0 0 60720 121440
287 319210 0 601 0 5 0 60720 0
288 318098 1 30 90 3 15 62045 124090
289 315380 0 0 0 0 0 60720 0
290 292754 1 249 747 13 65 54323 108646
291 315380 0 0 0 0 0 60720 0
292 325176 1 165 495 12 60 62841 125682
293 365959 1 453 1359 19 95 81125 162250
294 315380 1 0 0 0 0 60720 121440
295 302409 1 53 159 10 50 59506 119012
296 340968 0 382 0 9 0 59365 0
297 315380 1 0 0 0 0 60720 121440
298 315380 0 0 0 0 0 60720 0
299 315380 1 0 0 0 0 60720 121440
300 315380 1 0 0 9 45 60720 121440
301 313164 0 30 0 4 0 60798 0
302 301164 1 290 870 1 5 58790 117580
303 315380 1 0 0 1 5 60720 121440
304 315380 0 0 0 0 0 60720 0
305 344425 1 366 1098 14 70 61808 123616
306 315394 1 2 6 12 60 60735 121470
307 315380 0 0 0 0 0 60720 0
308 316647 0 209 0 19 0 64016 0
309 309836 1 384 1152 17 85 54683 109366
310 315380 0 0 0 0 0 60720 0
311 315380 0 0 0 0 0 60720 0
312 346611 1 365 1095 32 160 87192 174384
313 315380 0 0 0 0 0 60720 0
314 322031 0 49 0 14 0 64107 0
315 315656 1 3 9 8 40 60761 121522
316 339445 1 133 399 4 20 65990 131980
317 314964 1 32 96 0 0 59988 119976
318 297141 1 368 1104 20 100 61167 122334
319 315372 1 1 3 5 25 60719 121438
320 315380 1 0 0 0 0 60720 121440
321 315380 1 0 0 0 0 60720 121440
322 315380 1 0 0 0 0 60720 121440
323 315380 1 0 0 0 0 60720 121440
324 315380 0 0 0 0 0 60720 0
325 315380 1 0 0 0 0 60720 121440
326 312502 1 22 66 1 5 60722 121444
327 315380 0 0 0 0 0 60720 0
328 315380 0 0 0 0 0 60720 0
329 315380 1 0 0 0 0 60720 121440
330 315380 1 0 0 0 0 60720 121440
331 315380 1 0 0 0 0 60720 121440
332 315380 1 0 0 0 0 60720 121440
333 315380 1 0 0 0 0 60720 121440
334 313729 1 96 288 4 20 60379 120758
335 315388 1 1 3 1 5 60727 121454
336 315371 1 314 942 4 20 60720 121440
337 296139 1 844 2532 20 100 60925 121850
338 315380 1 0 0 0 0 60720 121440
339 313880 1 26 78 1 5 60896 121792
340 317698 1 125 375 10 50 59734 119468
341 295580 1 304 912 12 60 62969 125938
342 315380 1 0 0 0 0 60720 121440
343 315380 1 0 0 0 0 60720 121440
344 315380 1 0 0 0 0 60720 121440
345 308256 1 621 1863 13 65 60720 121440
346 315380 1 0 0 0 0 60720 121440
347 303677 1 119 357 3 15 59118 118236
348 315380 1 0 0 0 0 60720 121440
349 315380 1 0 0 0 0 60720 121440
350 319369 1 1595 4785 10 50 60720 121440
351 318690 1 312 936 3 15 58598 117196
352 314049 1 60 180 7 35 61124 122248
353 325699 1 587 1761 10 50 59595 119190
354 314210 1 135 405 1 5 62065 124130
355 315380 1 0 0 0 0 60720 121440
356 315380 1 0 0 0 0 60720 121440
357 322378 1 514 1542 15 75 78780 157560
358 315380 1 0 0 0 0 60720 121440
359 315380 1 0 0 0 0 60720 121440
360 315380 1 0 0 0 0 60720 121440
361 315398 1 1 3 4 20 60722 121444
362 315380 1 0 0 0 0 60720 121440
363 315380 1 0 0 0 0 60720 121440
364 308336 0 1763 0 28 0 61600 0
365 316386 1 180 540 9 45 59635 119270
366 315380 1 0 0 0 0 60720 121440
367 315380 1 0 0 0 0 60720 121440
368 315380 1 0 0 0 0 60720 121440
369 315380 1 0 0 0 0 60720 121440
370 315553 1 218 654 7 35 60720 121440
371 315380 1 0 0 0 0 60720 121440
372 323361 1 448 1344 7 35 59781 119562
373 336639 1 227 681 7 35 76644 153288
374 307424 1 174 522 3 15 64820 129640
375 315380 1 0 0 0 0 60720 121440
376 315380 1 0 0 0 0 60720 121440
377 295370 1 121 363 11 55 56178 112356
378 322340 1 607 1821 7 35 60436 120872
379 319864 1 2212 6636 10 50 60720 121440
380 315380 1 0 0 0 0 60720 121440
381 315380 1 0 0 0 0 60720 121440
382 317291 1 530 1590 18 90 73433 146866
383 280398 1 571 1713 14 70 41477 82954
384 315380 1 0 0 0 0 60720 121440
385 317330 1 78 234 12 60 62700 125400
386 238125 1 2489 7467 29 145 67804 135608
387 327071 1 131 393 3 15 59661 119322
388 309038 1 923 2769 6 30 58620 117240
389 314210 1 72 216 3 15 60398 120796
390 307930 1 572 1716 8 40 58580 117160
391 322327 1 397 1191 10 50 62710 125420
392 292136 1 450 1350 6 30 59325 118650
393 263276 1 622 1866 8 40 60950 121900
394 367655 1 694 2082 6 30 68060 136120
395 283910 0 3425 0 9 0 83620 0
396 283587 1 562 1686 8 40 58456 116912
397 243650 1 4917 14751 26 130 52811 105622
398 438493 0 1442 0 239 0 121173 0
399 296261 1 529 1587 7 35 63870 127740
400 230621 0 2126 0 41 0 21001 0
401 304252 1 1061 3183 3 15 70415 140830
402 333505 1 776 2328 8 40 64230 128460
403 296919 1 611 1833 6 30 59190 118380
404 278990 0 1526 0 21 0 69351 0
405 276898 1 592 1776 7 35 64270 128540
406 327007 1 1182 3546 11 55 70694 141388
407 317046 1 621 1863 11 55 68005 136010
408 304555 1 989 2967 12 60 58930 117860
409 298096 1 438 1314 9 45 58320 116640
410 231861 1 726 2178 3 15 69980 139960
411 309422 1 1303 3909 57 285 69863 139726
412 286963 0 7419 0 21 0 63255 0
413 269753 0 1164 0 15 0 57320 0
414 448243 0 3310 0 32 0 75230 0
415 165404 1 1920 5760 11 55 79420 158840
416 204325 1 965 2895 2 10 73490 146980
417 407159 1 3256 9768 23 115 35250 70500
418 290476 0 1135 0 20 0 62285 0
419 275311 1 1270 3810 24 120 69206 138412
420 246541 1 661 1983 1 5 65920 131840
421 253468 1 1013 3039 1 5 69770 139540
422 240897 1 2844 8532 74 370 72683 145366
423 -83265 0 11528 0 68 0 -14545 0
424 -42143 1 6526 19578 20 100 55830 111660
425 272713 1 2264 6792 20 100 55174 110348
426 215362 0 5109 0 82 0 67038 0
427 42754 1 3999 11997 21 105 51252 102504
428 306275 1 35624 106872 244 1220 157278 314556
429 253537 1 9252 27756 32 160 79510 159020
430 372631 1 15236 45708 86 430 77440 154880
431 -7170 1 18073 54219 69 345 27284 54568
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) Group Costs Costs_p Orders Orders_p
2.091e+05 -1.373e+04 2.322e+01 -6.056e+00 2.549e+03 -5.471e+02
Dividends Dividends_p
1.567e-01 8.710e-01
> (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) 2.091e+05 4.111e+04 5.087 5.49e-07 ***
Group -1.373e+04 4.962e+04 -0.277 0.78218
Costs 2.322e+01 1.992e+00 11.654 < 2e-16 ***
Costs_p -6.056e+00 1.160e+00 -5.219 2.83e-07 ***
Orders 2.549e+03 3.336e+02 7.641 1.46e-13 ***
Orders_p -5.471e+02 1.065e+02 -5.137 4.26e-07 ***
Dividends 1.567e-01 5.031e-01 0.312 0.75555
Dividends_p 8.710e-01 3.181e-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/1p6jf1291284860.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/2p6jf1291284860.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/3p6jf1291284860.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/40x0h1291284860.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/50x0h1291284860.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/6bphk1291284860.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/73gg51291284860.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/83gg51291284860.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/93gg51291284860.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/10wpx81291284860.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/11z7ew1291284860.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/12d0fx1291284861.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/13rsdo1291284861.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/14cabt1291284861.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/15gbah1291284861.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/16ul881291284861.tab")
+ }
>
> try(system("convert tmp/1p6jf1291284860.ps tmp/1p6jf1291284860.png",intern=TRUE))
character(0)
> try(system("convert tmp/2p6jf1291284860.ps tmp/2p6jf1291284860.png",intern=TRUE))
character(0)
> try(system("convert tmp/3p6jf1291284860.ps tmp/3p6jf1291284860.png",intern=TRUE))
character(0)
> try(system("convert tmp/40x0h1291284860.ps tmp/40x0h1291284860.png",intern=TRUE))
character(0)
> try(system("convert tmp/50x0h1291284860.ps tmp/50x0h1291284860.png",intern=TRUE))
character(0)
> try(system("convert tmp/6bphk1291284860.ps tmp/6bphk1291284860.png",intern=TRUE))
character(0)
> try(system("convert tmp/73gg51291284860.ps tmp/73gg51291284860.png",intern=TRUE))
character(0)
> try(system("convert tmp/83gg51291284860.ps tmp/83gg51291284860.png",intern=TRUE))
character(0)
> try(system("convert tmp/93gg51291284860.ps tmp/93gg51291284860.png",intern=TRUE))
character(0)
> try(system("convert tmp/10wpx81291284860.ps tmp/10wpx81291284860.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
14.077 2.290 35.644