R version 2.8.0 (2008-10-20)
Copyright (C) 2008 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
Natural language support but running in an English locale
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> x <- array(list(1
+ ,162556
+ ,162556
+ ,1081
+ ,1081
+ ,213118
+ ,213118
+ ,230380558
+ ,6282929
+ ,1
+ ,29790
+ ,29790
+ ,309
+ ,309
+ ,81767
+ ,81767
+ ,25266003
+ ,4324047
+ ,1
+ ,87550
+ ,87550
+ ,458
+ ,458
+ ,153198
+ ,153198
+ ,70164684
+ ,4108272
+ ,0
+ ,84738
+ ,0
+ ,588
+ ,0
+ ,-26007
+ ,0
+ ,-15292116
+ ,-1212617
+ ,1
+ ,54660
+ ,54660
+ ,299
+ ,299
+ ,126942
+ ,126942
+ ,37955658
+ ,1485329
+ ,1
+ ,42634
+ ,42634
+ ,156
+ ,156
+ ,157214
+ ,157214
+ ,24525384
+ ,1779876
+ ,0
+ ,40949
+ ,0
+ ,481
+ ,0
+ ,129352
+ ,0
+ ,62218312
+ ,1367203
+ ,1
+ ,42312
+ ,42312
+ ,323
+ ,323
+ ,234817
+ ,234817
+ ,75845891
+ ,2519076
+ ,1
+ ,37704
+ ,37704
+ ,452
+ ,452
+ ,60448
+ ,60448
+ ,27322496
+ ,912684
+ ,1
+ ,16275
+ ,16275
+ ,109
+ ,109
+ ,47818
+ ,47818
+ ,5212162
+ ,1443586
+ ,0
+ ,25830
+ ,0
+ ,115
+ ,0
+ ,245546
+ ,0
+ ,28237790
+ ,1220017
+ ,0
+ ,12679
+ ,0
+ ,110
+ ,0
+ ,48020
+ ,0
+ ,5282200
+ ,984885
+ ,1
+ ,18014
+ ,18014
+ ,239
+ ,239
+ ,-1710
+ ,-1710
+ ,-408690
+ ,1457425
+ ,0
+ ,43556
+ ,0
+ ,247
+ ,0
+ ,32648
+ ,0
+ ,8064056
+ ,-572920
+ ,1
+ ,24524
+ ,24524
+ ,497
+ ,497
+ ,95350
+ ,95350
+ ,47388950
+ ,929144
+ ,0
+ ,6532
+ ,0
+ ,103
+ ,0
+ ,151352
+ ,0
+ ,15589256
+ ,1151176
+ ,0
+ ,7123
+ ,0
+ ,109
+ ,0
+ ,288170
+ ,0
+ ,31410530
+ ,790090
+ ,1
+ ,20813
+ ,20813
+ ,502
+ ,502
+ ,114337
+ ,114337
+ ,57397174
+ ,774497
+ ,1
+ ,37597
+ ,37597
+ ,248
+ ,248
+ ,37884
+ ,37884
+ ,9395232
+ ,990576
+ ,0
+ ,17821
+ ,0
+ ,373
+ ,0
+ ,122844
+ ,0
+ ,45820812
+ ,454195
+ ,1
+ ,12988
+ ,12988
+ ,119
+ ,119
+ ,82340
+ ,82340
+ ,9798460
+ ,876607
+ ,1
+ ,22330
+ ,22330
+ ,84
+ ,84
+ ,79801
+ ,79801
+ ,6703284
+ ,711969
+ ,0
+ ,13326
+ ,0
+ ,102
+ ,0
+ ,165548
+ ,0
+ ,16885896
+ ,702380
+ ,0
+ ,16189
+ ,0
+ ,295
+ ,0
+ ,116384
+ ,0
+ ,34333280
+ ,264449
+ ,0
+ ,7146
+ ,0
+ ,105
+ ,0
+ ,134028
+ ,0
+ ,14072940
+ ,450033
+ ,0
+ ,15824
+ ,0
+ ,64
+ ,0
+ ,63838
+ ,0
+ ,4085632
+ ,541063
+ ,1
+ ,26088
+ ,26088
+ ,267
+ ,267
+ ,74996
+ ,74996
+ ,20023932
+ ,588864
+ ,0
+ ,11326
+ ,0
+ ,129
+ ,0
+ ,31080
+ ,0
+ ,4009320
+ ,-37216
+ ,0
+ ,8568
+ ,0
+ ,37
+ ,0
+ ,32168
+ ,0
+ ,1190216
+ ,783310
+ ,0
+ ,14416
+ ,0
+ ,361
+ ,0
+ ,49857
+ ,0
+ ,17998377
+ ,467359
+ ,1
+ ,3369
+ ,3369
+ ,28
+ ,28
+ ,87161
+ ,87161
+ ,2440508
+ ,688779
+ ,1
+ ,11819
+ ,11819
+ ,85
+ ,85
+ ,106113
+ ,106113
+ ,9019605
+ ,608419
+ ,1
+ ,6620
+ ,6620
+ ,44
+ ,44
+ ,80570
+ ,80570
+ ,3545080
+ ,696348
+ ,1
+ ,4519
+ ,4519
+ ,49
+ ,49
+ ,102129
+ ,102129
+ ,5004321
+ ,597793
+ ,0
+ ,2220
+ ,0
+ ,22
+ ,0
+ ,301670
+ ,0
+ ,6636740
+ ,821730
+ ,0
+ ,18562
+ ,0
+ ,155
+ ,0
+ ,102313
+ ,0
+ ,15858515
+ ,377934
+ ,0
+ ,10327
+ ,0
+ ,91
+ ,0
+ ,88577
+ ,0
+ ,8060507
+ ,651939
+ ,1
+ ,5336
+ ,5336
+ ,81
+ ,81
+ ,112477
+ ,112477
+ ,9110637
+ ,697458
+ ,1
+ ,2365
+ ,2365
+ ,79
+ ,79
+ ,191778
+ ,191778
+ ,15150462
+ ,700368
+ ,0
+ ,4069
+ ,0
+ ,145
+ ,0
+ ,79804
+ ,0
+ ,11571580
+ ,225986
+ ,0
+ ,7710
+ ,0
+ ,816
+ ,0
+ ,128294
+ ,0
+ ,104687904
+ ,348695
+ ,0
+ ,13718
+ ,0
+ ,61
+ ,0
+ ,96448
+ ,0
+ ,5883328
+ ,373683
+ ,0
+ ,4525
+ ,0
+ ,226
+ ,0
+ ,93811
+ ,0
+ ,21201286
+ ,501709
+ ,0
+ ,6869
+ ,0
+ ,105
+ ,0
+ ,117520
+ ,0
+ ,12339600
+ ,413743
+ ,0
+ ,4628
+ ,0
+ ,62
+ ,0
+ ,69159
+ ,0
+ ,4287858
+ ,379825
+ ,1
+ ,3653
+ ,3653
+ ,24
+ ,24
+ ,101792
+ ,101792
+ ,2443008
+ ,336260
+ ,1
+ ,1265
+ ,1265
+ ,26
+ ,26
+ ,210568
+ ,210568
+ ,5474768
+ ,636765
+ ,1
+ ,7489
+ ,7489
+ ,322
+ ,322
+ ,136996
+ ,136996
+ ,44112712
+ ,481231
+ ,0
+ ,4901
+ ,0
+ ,84
+ ,0
+ ,121920
+ ,0
+ ,10241280
+ ,469107
+ ,0
+ ,2284
+ ,0
+ ,33
+ ,0
+ ,76403
+ ,0
+ ,2521299
+ ,211928
+ ,1
+ ,3160
+ ,3160
+ ,108
+ ,108
+ ,108094
+ ,108094
+ ,11674152
+ ,563925
+ ,1
+ ,4150
+ ,4150
+ ,150
+ ,150
+ ,134759
+ ,134759
+ ,20213850
+ ,511939
+ ,1
+ ,7285
+ ,7285
+ ,115
+ ,115
+ ,188873
+ ,188873
+ ,21720395
+ ,521016
+ ,1
+ ,1134
+ ,1134
+ ,162
+ ,162
+ ,146216
+ ,146216
+ ,23686992
+ ,543856
+ ,1
+ ,4658
+ ,4658
+ ,158
+ ,158
+ ,156608
+ ,156608
+ ,24744064
+ ,329304
+ ,0
+ ,2384
+ ,0
+ ,97
+ ,0
+ ,61348
+ ,0
+ ,5950756
+ ,423262
+ ,0
+ ,3748
+ ,0
+ ,9
+ ,0
+ ,50350
+ ,0
+ ,453150
+ ,509665
+ ,0
+ ,5371
+ ,0
+ ,66
+ ,0
+ ,87720
+ ,0
+ ,5789520
+ ,455881
+ ,0
+ ,1285
+ ,0
+ ,107
+ ,0
+ ,99489
+ ,0
+ ,10645323
+ ,367772
+ ,1
+ ,9327
+ ,9327
+ ,101
+ ,101
+ ,87419
+ ,87419
+ ,8829319
+ ,406339
+ ,1
+ ,5565
+ ,5565
+ ,47
+ ,47
+ ,94355
+ ,94355
+ ,4434685
+ ,493408
+ ,0
+ ,1528
+ ,0
+ ,38
+ ,0
+ ,60326
+ ,0
+ ,2292388
+ ,232942
+ ,1
+ ,3122
+ ,3122
+ ,34
+ ,34
+ ,94670
+ ,94670
+ ,3218780
+ ,416002
+ ,1
+ ,7317
+ ,7317
+ ,84
+ ,84
+ ,82425
+ ,82425
+ ,6923700
+ ,337430
+ ,0
+ ,2675
+ ,0
+ ,79
+ ,0
+ ,59017
+ ,0
+ ,4662343
+ ,361517
+ ,0
+ ,13253
+ ,0
+ ,947
+ ,0
+ ,90829
+ ,0
+ ,86015063
+ ,360962
+ ,0
+ ,880
+ ,0
+ ,74
+ ,0
+ ,80791
+ ,0
+ ,5978534
+ ,235561
+ ,1
+ ,2053
+ ,2053
+ ,53
+ ,53
+ ,100423
+ ,100423
+ ,5322419
+ ,408247
+ ,0
+ ,1424
+ ,0
+ ,94
+ ,0
+ ,131116
+ ,0
+ ,12324904
+ ,450296
+ ,1
+ ,4036
+ ,4036
+ ,63
+ ,63
+ ,100269
+ ,100269
+ ,6316947
+ ,418799
+ ,1
+ ,3045
+ ,3045
+ ,58
+ ,58
+ ,27330
+ ,27330
+ ,1585140
+ ,247405
+ ,0
+ ,5119
+ ,0
+ ,49
+ ,0
+ ,39039
+ ,0
+ ,1912911
+ ,378519
+ ,0
+ ,1431
+ ,0
+ ,34
+ ,0
+ ,106885
+ ,0
+ ,3634090
+ ,326638
+ ,0
+ ,554
+ ,0
+ ,11
+ ,0
+ ,79285
+ ,0
+ ,872135
+ ,328233
+ ,0
+ ,1975
+ ,0
+ ,35
+ ,0
+ ,118881
+ ,0
+ ,4160835
+ ,386225
+ ,1
+ ,1286
+ ,1286
+ ,17
+ ,17
+ ,77623
+ ,77623
+ ,1319591
+ ,283662
+ ,0
+ ,1012
+ ,0
+ ,47
+ ,0
+ ,114768
+ ,0
+ ,5394096
+ ,370225
+ ,0
+ ,810
+ ,0
+ ,43
+ ,0
+ ,74015
+ ,0
+ ,3182645
+ ,269236
+ ,0
+ ,1280
+ ,0
+ ,117
+ ,0
+ ,69465
+ ,0
+ ,8127405
+ ,365732
+ ,1
+ ,666
+ ,666
+ ,171
+ ,171
+ ,117869
+ ,117869
+ ,20155599
+ ,420383
+ ,0
+ ,1380
+ ,0
+ ,26
+ ,0
+ ,60982
+ ,0
+ ,1585532
+ ,345811
+ ,1
+ ,4608
+ ,4608
+ ,73
+ ,73
+ ,90131
+ ,90131
+ ,6579563
+ ,431809
+ ,0
+ ,876
+ ,0
+ ,59
+ ,0
+ ,138971
+ ,0
+ ,8199289
+ ,418876
+ ,0
+ ,814
+ ,0
+ ,18
+ ,0
+ ,39625
+ ,0
+ ,713250
+ ,297476
+ ,0
+ ,514
+ ,0
+ ,15
+ ,0
+ ,102725
+ ,0
+ ,1540875
+ ,416776
+ ,1
+ ,5692
+ ,5692
+ ,72
+ ,72
+ ,64239
+ ,64239
+ ,4625208
+ ,357257
+ ,0
+ ,3642
+ ,0
+ ,86
+ ,0
+ ,90262
+ ,0
+ ,7762532
+ ,458343
+ ,0
+ ,540
+ ,0
+ ,14
+ ,0
+ ,103960
+ ,0
+ ,1455440
+ ,388386
+ ,0
+ ,2099
+ ,0
+ ,64
+ ,0
+ ,106611
+ ,0
+ ,6823104
+ ,358934
+ ,0
+ ,567
+ ,0
+ ,11
+ ,0
+ ,103345
+ ,0
+ ,1136795
+ ,407560
+ ,0
+ ,2001
+ ,0
+ ,52
+ ,0
+ ,95551
+ ,0
+ ,4968652
+ ,392558
+ ,1
+ ,2949
+ ,2949
+ ,41
+ ,41
+ ,82903
+ ,82903
+ ,3399023
+ ,373177
+ ,0
+ ,2253
+ ,0
+ ,99
+ ,0
+ ,63593
+ ,0
+ ,6295707
+ ,428370
+ ,1
+ ,6533
+ ,6533
+ ,75
+ ,75
+ ,126910
+ ,126910
+ ,9518250
+ ,369419
+ ,0
+ ,1889
+ ,0
+ ,45
+ ,0
+ ,37527
+ ,0
+ ,1688715
+ ,358649
+ ,1
+ ,3055
+ ,3055
+ ,43
+ ,43
+ ,60247
+ ,60247
+ ,2590621
+ ,376641
+ ,0
+ ,272
+ ,0
+ ,8
+ ,0
+ ,112995
+ ,0
+ ,903960
+ ,467427
+ ,1
+ ,1414
+ ,1414
+ ,198
+ ,198
+ ,70184
+ ,70184
+ ,13896432
+ ,364885
+ ,0
+ ,2564
+ ,0
+ ,22
+ ,0
+ ,130140
+ ,0
+ ,2863080
+ ,436230
+ ,1
+ ,1383
+ ,1383
+ ,11
+ ,11
+ ,73221
+ ,73221
+ ,805431
+ ,329118
+ ,1
+ ,1261
+ ,1261
+ ,33
+ ,33
+ ,76114
+ ,76114
+ ,2511762
+ ,317365
+ ,0
+ ,975
+ ,0
+ ,23
+ ,0
+ ,90534
+ ,0
+ ,2082282
+ ,286849
+ ,0
+ ,3366
+ ,0
+ ,80
+ ,0
+ ,108479
+ ,0
+ ,8678320
+ ,376685
+ ,0
+ ,576
+ ,0
+ ,18
+ ,0
+ ,113761
+ ,0
+ ,2047698
+ ,407198
+ ,0
+ ,1306
+ ,0
+ ,28
+ ,0
+ ,68696
+ ,0
+ ,1923488
+ ,377772
+ ,0
+ ,746
+ ,0
+ ,23
+ ,0
+ ,71561
+ ,0
+ ,1645903
+ ,271483
+ ,1
+ ,3192
+ ,3192
+ ,60
+ ,60
+ ,59831
+ ,59831
+ ,3589860
+ ,153661
+ ,1
+ ,2045
+ ,2045
+ ,20
+ ,20
+ ,97890
+ ,97890
+ ,1957800
+ ,513294
+ ,0
+ ,5477
+ ,0
+ ,59
+ ,0
+ ,101481
+ ,0
+ ,5987379
+ ,324881
+ ,1
+ ,1932
+ ,1932
+ ,36
+ ,36
+ ,72954
+ ,72954
+ ,2626344
+ ,264512
+ ,0
+ ,936
+ ,0
+ ,30
+ ,0
+ ,67939
+ ,0
+ ,2038170
+ ,420968
+ ,1
+ ,3437
+ ,3437
+ ,47
+ ,47
+ ,48022
+ ,48022
+ ,2257034
+ ,129302
+ ,0
+ ,5131
+ ,0
+ ,71
+ ,0
+ ,86111
+ ,0
+ ,6113881
+ ,191521
+ ,1
+ ,2397
+ ,2397
+ ,14
+ ,14
+ ,74020
+ ,74020
+ ,1036280
+ ,268673
+ ,1
+ ,1389
+ ,1389
+ ,9
+ ,9
+ ,57530
+ ,57530
+ ,517770
+ ,353179
+ ,0
+ ,1503
+ ,0
+ ,39
+ ,0
+ ,56364
+ ,0
+ ,2198196
+ ,354624
+ ,0
+ ,402
+ ,0
+ ,26
+ ,0
+ ,84990
+ ,0
+ ,2209740
+ ,363713
+ ,0
+ ,2239
+ ,0
+ ,21
+ ,0
+ ,88590
+ ,0
+ ,1860390
+ ,456657
+ ,1
+ ,2234
+ ,2234
+ ,16
+ ,16
+ ,77200
+ ,77200
+ ,1235200
+ ,211742
+ ,0
+ ,837
+ ,0
+ ,69
+ ,0
+ ,61262
+ ,0
+ ,4227078
+ ,338381
+ ,0
+ ,10579
+ ,0
+ ,92
+ ,0
+ ,110309
+ ,0
+ ,10148428
+ ,418530
+ ,0
+ ,875
+ ,0
+ ,14
+ ,0
+ ,67000
+ ,0
+ ,938000
+ ,351483
+ ,0
+ ,1395
+ ,0
+ ,103
+ ,0
+ ,93099
+ ,0
+ ,9589197
+ ,372928
+ ,1
+ ,1659
+ ,1659
+ ,29
+ ,29
+ ,107577
+ ,107577
+ ,3119733
+ ,485538
+ ,1
+ ,2647
+ ,2647
+ ,37
+ ,37
+ ,62920
+ ,62920
+ ,2328040
+ ,279268
+ ,1
+ ,3294
+ ,3294
+ ,23
+ ,23
+ ,75832
+ ,75832
+ ,1744136
+ ,219060
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,94
+ ,0
+ ,7
+ ,0
+ ,60793
+ ,0
+ ,425551
+ ,325314
+ ,0
+ ,422
+ ,0
+ ,28
+ ,0
+ ,57935
+ ,0
+ ,1622180
+ ,322046
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,34
+ ,0
+ ,8
+ ,0
+ ,60630
+ ,0
+ ,485040
+ ,325599
+ ,0
+ ,1558
+ ,0
+ ,63
+ ,0
+ ,55637
+ ,0
+ ,3505131
+ ,377028
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,43
+ ,0
+ ,3
+ ,0
+ ,60887
+ ,0
+ ,182661
+ ,323850
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,316
+ ,0
+ ,9
+ ,0
+ ,60505
+ ,0
+ ,544545
+ ,331514
+ ,0
+ ,115
+ ,0
+ ,13
+ ,0
+ ,60945
+ ,0
+ ,792285
+ ,325632
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,389
+ ,0
+ ,14
+ ,0
+ ,58990
+ ,0
+ ,825860
+ ,322265
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,1002
+ ,0
+ ,15
+ ,0
+ ,56750
+ ,0
+ ,851250
+ ,325906
+ ,0
+ ,36
+ ,0
+ ,3
+ ,0
+ ,60894
+ ,0
+ ,182682
+ ,325985
+ ,0
+ ,460
+ ,0
+ ,15
+ ,0
+ ,63346
+ ,0
+ ,950190
+ ,346145
+ ,0
+ ,309
+ ,0
+ ,11
+ ,0
+ ,56535
+ ,0
+ ,621885
+ ,325898
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,9
+ ,0
+ ,6
+ ,0
+ ,60835
+ ,0
+ ,365010
+ ,325356
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,14
+ ,0
+ ,1
+ ,0
+ ,61016
+ ,0
+ ,61016
+ ,325930
+ ,0
+ ,520
+ ,0
+ ,10
+ ,0
+ ,58650
+ ,0
+ ,586500
+ ,318020
+ ,0
+ ,1766
+ ,0
+ ,73
+ ,0
+ ,60438
+ ,0
+ ,4411974
+ ,326389
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,458
+ ,0
+ ,11
+ ,0
+ ,58625
+ ,0
+ ,644875
+ ,302925
+ ,0
+ ,20
+ ,0
+ ,3
+ ,0
+ ,60938
+ ,0
+ ,182814
+ ,325540
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,98
+ ,0
+ ,2
+ ,0
+ ,61490
+ ,0
+ ,122980
+ ,326736
+ ,0
+ ,405
+ ,0
+ ,7
+ ,0
+ ,60845
+ ,0
+ ,425915
+ ,340580
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,483
+ ,0
+ ,27
+ ,0
+ ,60830
+ ,0
+ ,1642410
+ ,331828
+ ,0
+ ,454
+ ,0
+ ,51
+ ,0
+ ,63261
+ ,0
+ ,3226311
+ ,323299
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,757
+ ,0
+ ,19
+ ,0
+ ,45689
+ ,0
+ ,868091
+ ,387722
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,36
+ ,0
+ ,4
+ ,0
+ ,61564
+ ,0
+ ,246256
+ ,324598
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,203
+ ,0
+ ,9
+ ,0
+ ,61938
+ ,0
+ ,557442
+ ,328726
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,90
+ ,0
+ ,8
+ ,0
+ ,60951
+ ,0
+ ,487608
+ ,325043
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,1
+ ,71
+ ,71
+ ,1
+ ,1
+ ,60745
+ ,60745
+ ,60745
+ ,325806
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,972
+ ,0
+ ,34
+ ,0
+ ,71642
+ ,0
+ ,2435828
+ ,387732
+ ,1
+ ,531
+ ,531
+ ,10
+ ,10
+ ,71641
+ ,71641
+ ,716410
+ ,349729
+ ,0
+ ,604
+ ,0
+ ,38
+ ,0
+ ,55792
+ ,0
+ ,2120096
+ ,332202
+ ,1
+ ,283
+ ,283
+ ,10
+ ,10
+ ,71873
+ ,71873
+ ,718730
+ ,305442
+ ,1
+ ,23
+ ,23
+ ,5
+ ,5
+ ,62555
+ ,62555
+ ,312775
+ ,329537
+ ,1
+ ,638
+ ,638
+ ,14
+ ,14
+ ,60370
+ ,60370
+ ,845180
+ ,327055
+ ,1
+ ,699
+ ,699
+ ,16
+ ,16
+ ,64873
+ ,64873
+ ,1037968
+ ,356245
+ ,0
+ ,149
+ ,0
+ ,5
+ ,0
+ ,62041
+ ,0
+ ,310205
+ ,328451
+ ,0
+ ,226
+ ,0
+ ,5
+ ,0
+ ,65745
+ ,0
+ ,328725
+ ,307062
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,275
+ ,0
+ ,4
+ ,0
+ ,59500
+ ,0
+ ,238000
+ ,331345
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,141
+ ,0
+ ,6
+ ,0
+ ,61630
+ ,0
+ ,369780
+ ,331824
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,28
+ ,0
+ ,2
+ ,0
+ ,60890
+ ,0
+ ,121780
+ ,325685
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,2566
+ ,2566
+ ,80
+ ,80
+ ,113521
+ ,113521
+ ,9081680
+ ,404480
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,472
+ ,472
+ ,20
+ ,20
+ ,80045
+ ,80045
+ ,1600900
+ ,318314
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,203
+ ,203
+ ,27
+ ,27
+ ,50804
+ ,50804
+ ,1371708
+ ,311807
+ ,1
+ ,496
+ ,496
+ ,17
+ ,17
+ ,87390
+ ,87390
+ ,1485630
+ ,337724
+ ,1
+ ,10
+ ,10
+ ,2
+ ,2
+ ,61656
+ ,61656
+ ,123312
+ ,326431
+ ,1
+ ,63
+ ,63
+ ,4
+ ,4
+ ,65688
+ ,65688
+ ,262752
+ ,327556
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,1136
+ ,1136
+ ,32
+ ,32
+ ,48522
+ ,48522
+ ,1552704
+ ,356850
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,267
+ ,0
+ ,32
+ ,0
+ ,57640
+ ,0
+ ,1844480
+ ,322741
+ ,0
+ ,474
+ ,0
+ ,20
+ ,0
+ ,61977
+ ,0
+ ,1239540
+ ,310902
+ ,0
+ ,534
+ ,0
+ ,7
+ ,0
+ ,62620
+ ,0
+ ,438340
+ ,324295
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,15
+ ,0
+ ,8
+ ,0
+ ,60831
+ ,0
+ ,486648
+ ,326156
+ ,0
+ ,397
+ ,0
+ ,28
+ ,0
+ ,60646
+ ,0
+ ,1698088
+ ,326960
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,1061
+ ,0
+ ,20
+ ,0
+ ,56225
+ ,0
+ ,1124500
+ ,333411
+ ,0
+ ,288
+ ,0
+ ,4
+ ,0
+ ,60510
+ ,0
+ ,242040
+ ,297761
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,3
+ ,0
+ ,2
+ ,0
+ ,60698
+ ,0
+ ,121396
+ ,325536
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,20
+ ,0
+ ,2
+ ,0
+ ,60805
+ ,0
+ ,121610
+ ,325762
+ ,0
+ ,278
+ ,0
+ ,26
+ ,0
+ ,61404
+ ,0
+ ,1596504
+ ,327957
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,192
+ ,0
+ ,4
+ ,0
+ ,65276
+ ,0
+ ,261104
+ ,318521
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,317
+ ,0
+ ,9
+ ,0
+ ,63915
+ ,0
+ ,575235
+ ,319775
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,368
+ ,0
+ ,17
+ ,0
+ ,61686
+ ,0
+ ,1048662
+ ,332128
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,2
+ ,0
+ ,1
+ ,0
+ ,60743
+ ,0
+ ,60743
+ ,325486
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,53
+ ,0
+ ,6
+ ,0
+ ,60349
+ ,0
+ ,362094
+ ,325838
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,94
+ ,0
+ ,3
+ ,0
+ ,61360
+ ,0
+ ,184080
+ ,331767
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,24
+ ,0
+ ,8
+ ,0
+ ,59818
+ ,0
+ ,478544
+ ,324523
+ ,0
+ ,2332
+ ,0
+ ,4
+ ,0
+ ,72680
+ ,0
+ ,290720
+ ,339995
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,131
+ ,0
+ ,11
+ ,0
+ ,61808
+ ,0
+ ,679888
+ ,319582
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,206
+ ,0
+ ,9
+ ,0
+ ,53110
+ ,0
+ ,477990
+ ,307245
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,167
+ ,0
+ ,2
+ ,0
+ ,64245
+ ,0
+ ,128490
+ ,317967
+ ,0
+ ,622
+ ,0
+ ,73
+ ,0
+ ,73007
+ ,0
+ ,5329511
+ ,331488
+ ,0
+ ,885
+ ,0
+ ,85
+ ,0
+ ,82732
+ ,0
+ ,7032220
+ ,335452
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,365
+ ,0
+ ,8
+ ,0
+ ,54820
+ ,0
+ ,438560
+ ,334184
+ ,0
+ ,364
+ ,0
+ ,35
+ ,0
+ ,47705
+ ,0
+ ,1669675
+ ,313213
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,226
+ ,0
+ ,12
+ ,0
+ ,72835
+ ,0
+ ,874020
+ ,348678
+ ,0
+ ,307
+ ,0
+ ,15
+ ,0
+ ,58856
+ ,0
+ ,882840
+ ,328727
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,188
+ ,0
+ ,11
+ ,0
+ ,77655
+ ,0
+ ,854205
+ ,387978
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,138
+ ,0
+ ,6
+ ,0
+ ,69817
+ ,0
+ ,418902
+ ,336704
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,125
+ ,0
+ ,12
+ ,0
+ ,60798
+ ,0
+ ,729576
+ ,322076
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,282
+ ,0
+ ,30
+ ,0
+ ,62452
+ ,0
+ ,1873560
+ ,334272
+ ,0
+ ,335
+ ,0
+ ,33
+ ,0
+ ,64175
+ ,0
+ ,2117775
+ ,338197
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,1
+ ,813
+ ,813
+ ,82
+ ,82
+ ,67440
+ ,67440
+ ,5530080
+ ,321024
+ ,0
+ ,176
+ ,0
+ ,28
+ ,0
+ ,68136
+ ,0
+ ,1907808
+ ,322145
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,249
+ ,0
+ ,72
+ ,0
+ ,56726
+ ,0
+ ,4084272
+ ,323351
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,333
+ ,0
+ ,13
+ ,0
+ ,70811
+ ,0
+ ,920543
+ ,327748
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,30
+ ,0
+ ,4
+ ,0
+ ,62045
+ ,0
+ ,248180
+ ,328157
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,249
+ ,0
+ ,62
+ ,0
+ ,54323
+ ,0
+ ,3368026
+ ,311594
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,165
+ ,0
+ ,24
+ ,0
+ ,62841
+ ,0
+ ,1508184
+ ,335962
+ ,0
+ ,453
+ ,0
+ ,21
+ ,0
+ ,81125
+ ,0
+ ,1703625
+ ,372426
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,53
+ ,0
+ ,14
+ ,0
+ ,59506
+ ,0
+ ,833084
+ ,319844
+ ,1
+ ,382
+ ,382
+ ,21
+ ,21
+ ,59365
+ ,59365
+ ,1246665
+ ,355822
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,1
+ ,30
+ ,30
+ ,4
+ ,4
+ ,60798
+ ,60798
+ ,243192
+ ,324047
+ ,0
+ ,290
+ ,0
+ ,2
+ ,0
+ ,58790
+ ,0
+ ,117580
+ ,311464
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,366
+ ,0
+ ,53
+ ,0
+ ,61808
+ ,0
+ ,3275824
+ ,353417
+ ,0
+ ,2
+ ,0
+ ,9
+ ,0
+ ,60735
+ ,0
+ ,546615
+ ,325590
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,209
+ ,209
+ ,13
+ ,13
+ ,64016
+ ,64016
+ ,832208
+ ,328576
+ ,0
+ ,384
+ ,0
+ ,22
+ ,0
+ ,54683
+ ,0
+ ,1203026
+ ,326126
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,365
+ ,0
+ ,83
+ ,0
+ ,87192
+ ,0
+ ,7236936
+ ,369376
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,49
+ ,49
+ ,8
+ ,8
+ ,64107
+ ,64107
+ ,512856
+ ,332013
+ ,0
+ ,3
+ ,0
+ ,4
+ ,0
+ ,60761
+ ,0
+ ,243044
+ ,325871
+ ,0
+ ,133
+ ,0
+ ,14
+ ,0
+ ,65990
+ ,0
+ ,923860
+ ,342165
+ ,0
+ ,32
+ ,0
+ ,1
+ ,0
+ ,59988
+ ,0
+ ,59988
+ ,324967
+ ,0
+ ,368
+ ,0
+ ,17
+ ,0
+ ,61167
+ ,0
+ ,1039839
+ ,314832
+ ,0
+ ,1
+ ,0
+ ,6
+ ,0
+ ,60719
+ ,0
+ ,364314
+ ,325557
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,22
+ ,0
+ ,2
+ ,0
+ ,60722
+ ,0
+ ,121444
+ ,322649
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,96
+ ,0
+ ,5
+ ,0
+ ,60379
+ ,0
+ ,301895
+ ,324598
+ ,0
+ ,1
+ ,0
+ ,2
+ ,0
+ ,60727
+ ,0
+ ,121454
+ ,325567
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,81
+ ,0
+ ,7
+ ,0
+ ,60925
+ ,0
+ ,426475
+ ,324005
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,26
+ ,0
+ ,1
+ ,0
+ ,60896
+ ,0
+ ,60896
+ ,325748
+ ,0
+ ,125
+ ,0
+ ,13
+ ,0
+ ,59734
+ ,0
+ ,776542
+ ,323385
+ ,0
+ ,304
+ ,0
+ ,15
+ ,0
+ ,62969
+ ,0
+ ,944535
+ ,315409
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,119
+ ,0
+ ,6
+ ,0
+ ,59118
+ ,0
+ ,354708
+ ,312275
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,312
+ ,0
+ ,14
+ ,0
+ ,58598
+ ,0
+ ,820372
+ ,320576
+ ,0
+ ,60
+ ,0
+ ,10
+ ,0
+ ,61124
+ ,0
+ ,611240
+ ,325246
+ ,0
+ ,587
+ ,0
+ ,12
+ ,0
+ ,59595
+ ,0
+ ,715140
+ ,332961
+ ,0
+ ,135
+ ,0
+ ,2
+ ,0
+ ,62065
+ ,0
+ ,124130
+ ,323010
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,514
+ ,0
+ ,52
+ ,0
+ ,78780
+ ,0
+ ,4096560
+ ,345253
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,1
+ ,0
+ ,4
+ ,0
+ ,60722
+ ,0
+ ,242888
+ ,325559
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,1
+ ,58
+ ,58
+ ,3
+ ,3
+ ,61600
+ ,61600
+ ,184800
+ ,319634
+ ,0
+ ,180
+ ,0
+ ,11
+ ,0
+ ,59635
+ ,0
+ ,655985
+ ,319951
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,448
+ ,0
+ ,40
+ ,0
+ ,59781
+ ,0
+ ,2391240
+ ,318519
+ ,0
+ ,227
+ ,0
+ ,9
+ ,0
+ ,76644
+ ,0
+ ,689796
+ ,343222
+ ,0
+ ,174
+ ,0
+ ,1
+ ,0
+ ,64820
+ ,0
+ ,64820
+ ,317234
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,121
+ ,0
+ ,24
+ ,0
+ ,56178
+ ,0
+ ,1348272
+ ,314025
+ ,0
+ ,607
+ ,0
+ ,11
+ ,0
+ ,60436
+ ,0
+ ,664796
+ ,320249
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,530
+ ,0
+ ,60
+ ,0
+ ,73433
+ ,0
+ ,4405980
+ ,349365
+ ,0
+ ,571
+ ,0
+ ,80
+ ,0
+ ,41477
+ ,0
+ ,3318160
+ ,289197
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,78
+ ,0
+ ,16
+ ,0
+ ,62700
+ ,0
+ ,1003200
+ ,329245
+ ,0
+ ,2489
+ ,0
+ ,40
+ ,0
+ ,67804
+ ,0
+ ,2712160
+ ,240869
+ ,0
+ ,131
+ ,0
+ ,6
+ ,0
+ ,59661
+ ,0
+ ,357966
+ ,327182
+ ,0
+ ,923
+ ,0
+ ,8
+ ,0
+ ,58620
+ ,0
+ ,468960
+ ,322876
+ ,0
+ ,72
+ ,0
+ ,3
+ ,0
+ ,60398
+ ,0
+ ,181194
+ ,323117
+ ,0
+ ,572
+ ,0
+ ,16
+ ,0
+ ,58580
+ ,0
+ ,937280
+ ,306351
+ ,0
+ ,397
+ ,0
+ ,10
+ ,0
+ ,62710
+ ,0
+ ,627100
+ ,335137
+ ,0
+ ,450
+ ,0
+ ,8
+ ,0
+ ,59325
+ ,0
+ ,474600
+ ,308271
+ ,0
+ ,622
+ ,0
+ ,7
+ ,0
+ ,60950
+ ,0
+ ,426650
+ ,301731
+ ,0
+ ,694
+ ,0
+ ,8
+ ,0
+ ,68060
+ ,0
+ ,544480
+ ,382409
+ ,1
+ ,3425
+ ,3425
+ ,12
+ ,12
+ ,83620
+ ,83620
+ ,1003440
+ ,279230
+ ,0
+ ,562
+ ,0
+ ,13
+ ,0
+ ,58456
+ ,0
+ ,759928
+ ,298731
+ ,0
+ ,4917
+ ,0
+ ,42
+ ,0
+ ,52811
+ ,0
+ ,2218062
+ ,243650
+ ,1
+ ,1442
+ ,1442
+ ,118
+ ,118
+ ,121173
+ ,121173
+ ,14298414
+ ,532682
+ ,0
+ ,529
+ ,0
+ ,9
+ ,0
+ ,63870
+ ,0
+ ,574830
+ ,319771
+ ,1
+ ,2126
+ ,2126
+ ,138
+ ,138
+ ,21001
+ ,21001
+ ,2898138
+ ,171493
+ ,0
+ ,1061
+ ,0
+ ,5
+ ,0
+ ,70415
+ ,0
+ ,352075
+ ,347262
+ ,0
+ ,776
+ ,0
+ ,9
+ ,0
+ ,64230
+ ,0
+ ,578070
+ ,343945
+ ,0
+ ,611
+ ,0
+ ,8
+ ,0
+ ,59190
+ ,0
+ ,473520
+ ,311874
+ ,1
+ ,1526
+ ,1526
+ ,25
+ ,25
+ ,69351
+ ,69351
+ ,1733775
+ ,302211
+ ,0
+ ,592
+ ,0
+ ,7
+ ,0
+ ,64270
+ ,0
+ ,449890
+ ,316708
+ ,0
+ ,1182
+ ,0
+ ,13
+ ,0
+ ,70694
+ ,0
+ ,919022
+ ,333463
+ ,0
+ ,621
+ ,0
+ ,16
+ ,0
+ ,68005
+ ,0
+ ,1088080
+ ,344282
+ ,0
+ ,989
+ ,0
+ ,11
+ ,0
+ ,58930
+ ,0
+ ,648230
+ ,319635
+ ,0
+ ,438
+ ,0
+ ,11
+ ,0
+ ,58320
+ ,0
+ ,641520
+ ,301186
+ ,0
+ ,726
+ ,0
+ ,3
+ ,0
+ ,69980
+ ,0
+ ,209940
+ ,300381
+ ,0
+ ,1303
+ ,0
+ ,61
+ ,0
+ ,69863
+ ,0
+ ,4261643
+ ,318765
+ ,1
+ ,6341
+ ,6341
+ ,24
+ ,24
+ ,63255
+ ,63255
+ ,1518120
+ ,286146
+ ,1
+ ,1164
+ ,1164
+ ,17
+ ,17
+ ,57320
+ ,57320
+ ,974440
+ ,306844
+ ,1
+ ,3310
+ ,3310
+ ,33
+ ,33
+ ,75230
+ ,75230
+ ,2482590
+ ,307705
+ ,0
+ ,1366
+ ,0
+ ,7
+ ,0
+ ,79420
+ ,0
+ ,555940
+ ,312448
+ ,0
+ ,965
+ ,0
+ ,3
+ ,0
+ ,73490
+ ,0
+ ,220470
+ ,299715
+ ,0
+ ,3256
+ ,0
+ ,66
+ ,0
+ ,35250
+ ,0
+ ,2326500
+ ,373399
+ ,1
+ ,1135
+ ,1135
+ ,17
+ ,17
+ ,62285
+ ,62285
+ ,1058845
+ ,299446
+ ,0
+ ,1270
+ ,0
+ ,26
+ ,0
+ ,69206
+ ,0
+ ,1799356
+ ,325586
+ ,0
+ ,661
+ ,0
+ ,3
+ ,0
+ ,65920
+ ,0
+ ,197760
+ ,291221
+ ,0
+ ,1013
+ ,0
+ ,2
+ ,0
+ ,69770
+ ,0
+ ,139540
+ ,261173
+ ,0
+ ,2844
+ ,0
+ ,67
+ ,0
+ ,72683
+ ,0
+ ,4869761
+ ,255027
+ ,1
+ ,11528
+ ,11528
+ ,70
+ ,70
+ ,-14545
+ ,-14545
+ ,-1018150
+ ,-78375
+ ,0
+ ,6526
+ ,0
+ ,26
+ ,0
+ ,55830
+ ,0
+ ,1451580
+ ,-58143
+ ,0
+ ,2264
+ ,0
+ ,24
+ ,0
+ ,55174
+ ,0
+ ,1324176
+ ,227033
+ ,1
+ ,4461
+ ,4461
+ ,94
+ ,94
+ ,67038
+ ,67038
+ ,6301572
+ ,235098
+ ,0
+ ,3999
+ ,0
+ ,30
+ ,0
+ ,51252
+ ,0
+ ,1537560
+ ,21267
+ ,0
+ ,35624
+ ,0
+ ,223
+ ,0
+ ,157278
+ ,0
+ ,35072994
+ ,238675
+ ,0
+ ,9252
+ ,0
+ ,48
+ ,0
+ ,79510
+ ,0
+ ,3816480
+ ,197687
+ ,0
+ ,15236
+ ,0
+ ,90
+ ,0
+ ,77440
+ ,0
+ ,6969600
+ ,418341
+ ,0
+ ,18073
+ ,0
+ ,180
+ ,0
+ ,27284
+ ,0
+ ,4911120
+ ,-297706)
+ ,dim=c(9
+ ,431)
+ ,dimnames=list(c('Group'
+ ,'Costs'
+ ,'GrCosts'
+ ,'Trades'
+ ,'GrTrades'
+ ,'Dividends'
+ ,'GrDiv'
+ ,'TrDiv'
+ ,'Wealth
')
+ ,1:431))
> y <- array(NA,dim=c(9,431),dimnames=list(c('Group','Costs','GrCosts','Trades','GrTrades','Dividends','GrDiv','TrDiv','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 = '9'
> #'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 GrCosts Trades GrTrades Dividends GrDiv TrDiv
1 6282929 1 162556 162556 1081 1081 213118 213118 230380558
2 4324047 1 29790 29790 309 309 81767 81767 25266003
3 4108272 1 87550 87550 458 458 153198 153198 70164684
4 -1212617 0 84738 0 588 0 -26007 0 -15292116
5 1485329 1 54660 54660 299 299 126942 126942 37955658
6 1779876 1 42634 42634 156 156 157214 157214 24525384
7 1367203 0 40949 0 481 0 129352 0 62218312
8 2519076 1 42312 42312 323 323 234817 234817 75845891
9 912684 1 37704 37704 452 452 60448 60448 27322496
10 1443586 1 16275 16275 109 109 47818 47818 5212162
11 1220017 0 25830 0 115 0 245546 0 28237790
12 984885 0 12679 0 110 0 48020 0 5282200
13 1457425 1 18014 18014 239 239 -1710 -1710 -408690
14 -572920 0 43556 0 247 0 32648 0 8064056
15 929144 1 24524 24524 497 497 95350 95350 47388950
16 1151176 0 6532 0 103 0 151352 0 15589256
17 790090 0 7123 0 109 0 288170 0 31410530
18 774497 1 20813 20813 502 502 114337 114337 57397174
19 990576 1 37597 37597 248 248 37884 37884 9395232
20 454195 0 17821 0 373 0 122844 0 45820812
21 876607 1 12988 12988 119 119 82340 82340 9798460
22 711969 1 22330 22330 84 84 79801 79801 6703284
23 702380 0 13326 0 102 0 165548 0 16885896
24 264449 0 16189 0 295 0 116384 0 34333280
25 450033 0 7146 0 105 0 134028 0 14072940
26 541063 0 15824 0 64 0 63838 0 4085632
27 588864 1 26088 26088 267 267 74996 74996 20023932
28 -37216 0 11326 0 129 0 31080 0 4009320
29 783310 0 8568 0 37 0 32168 0 1190216
30 467359 0 14416 0 361 0 49857 0 17998377
31 688779 1 3369 3369 28 28 87161 87161 2440508
32 608419 1 11819 11819 85 85 106113 106113 9019605
33 696348 1 6620 6620 44 44 80570 80570 3545080
34 597793 1 4519 4519 49 49 102129 102129 5004321
35 821730 0 2220 0 22 0 301670 0 6636740
36 377934 0 18562 0 155 0 102313 0 15858515
37 651939 0 10327 0 91 0 88577 0 8060507
38 697458 1 5336 5336 81 81 112477 112477 9110637
39 700368 1 2365 2365 79 79 191778 191778 15150462
40 225986 0 4069 0 145 0 79804 0 11571580
41 348695 0 7710 0 816 0 128294 0 104687904
42 373683 0 13718 0 61 0 96448 0 5883328
43 501709 0 4525 0 226 0 93811 0 21201286
44 413743 0 6869 0 105 0 117520 0 12339600
45 379825 0 4628 0 62 0 69159 0 4287858
46 336260 1 3653 3653 24 24 101792 101792 2443008
47 636765 1 1265 1265 26 26 210568 210568 5474768
48 481231 1 7489 7489 322 322 136996 136996 44112712
49 469107 0 4901 0 84 0 121920 0 10241280
50 211928 0 2284 0 33 0 76403 0 2521299
51 563925 1 3160 3160 108 108 108094 108094 11674152
52 511939 1 4150 4150 150 150 134759 134759 20213850
53 521016 1 7285 7285 115 115 188873 188873 21720395
54 543856 1 1134 1134 162 162 146216 146216 23686992
55 329304 1 4658 4658 158 158 156608 156608 24744064
56 423262 0 2384 0 97 0 61348 0 5950756
57 509665 0 3748 0 9 0 50350 0 453150
58 455881 0 5371 0 66 0 87720 0 5789520
59 367772 0 1285 0 107 0 99489 0 10645323
60 406339 1 9327 9327 101 101 87419 87419 8829319
61 493408 1 5565 5565 47 47 94355 94355 4434685
62 232942 0 1528 0 38 0 60326 0 2292388
63 416002 1 3122 3122 34 34 94670 94670 3218780
64 337430 1 7317 7317 84 84 82425 82425 6923700
65 361517 0 2675 0 79 0 59017 0 4662343
66 360962 0 13253 0 947 0 90829 0 86015063
67 235561 0 880 0 74 0 80791 0 5978534
68 408247 1 2053 2053 53 53 100423 100423 5322419
69 450296 0 1424 0 94 0 131116 0 12324904
70 418799 1 4036 4036 63 63 100269 100269 6316947
71 247405 1 3045 3045 58 58 27330 27330 1585140
72 378519 0 5119 0 49 0 39039 0 1912911
73 326638 0 1431 0 34 0 106885 0 3634090
74 328233 0 554 0 11 0 79285 0 872135
75 386225 0 1975 0 35 0 118881 0 4160835
76 283662 1 1286 1286 17 17 77623 77623 1319591
77 370225 0 1012 0 47 0 114768 0 5394096
78 269236 0 810 0 43 0 74015 0 3182645
79 365732 0 1280 0 117 0 69465 0 8127405
80 420383 1 666 666 171 171 117869 117869 20155599
81 345811 0 1380 0 26 0 60982 0 1585532
82 431809 1 4608 4608 73 73 90131 90131 6579563
83 418876 0 876 0 59 0 138971 0 8199289
84 297476 0 814 0 18 0 39625 0 713250
85 416776 0 514 0 15 0 102725 0 1540875
86 357257 1 5692 5692 72 72 64239 64239 4625208
87 458343 0 3642 0 86 0 90262 0 7762532
88 388386 0 540 0 14 0 103960 0 1455440
89 358934 0 2099 0 64 0 106611 0 6823104
90 407560 0 567 0 11 0 103345 0 1136795
91 392558 0 2001 0 52 0 95551 0 4968652
92 373177 1 2949 2949 41 41 82903 82903 3399023
93 428370 0 2253 0 99 0 63593 0 6295707
94 369419 1 6533 6533 75 75 126910 126910 9518250
95 358649 0 1889 0 45 0 37527 0 1688715
96 376641 1 3055 3055 43 43 60247 60247 2590621
97 467427 0 272 0 8 0 112995 0 903960
98 364885 1 1414 1414 198 198 70184 70184 13896432
99 436230 0 2564 0 22 0 130140 0 2863080
100 329118 1 1383 1383 11 11 73221 73221 805431
101 317365 1 1261 1261 33 33 76114 76114 2511762
102 286849 0 975 0 23 0 90534 0 2082282
103 376685 0 3366 0 80 0 108479 0 8678320
104 407198 0 576 0 18 0 113761 0 2047698
105 377772 0 1306 0 28 0 68696 0 1923488
106 271483 0 746 0 23 0 71561 0 1645903
107 153661 1 3192 3192 60 60 59831 59831 3589860
108 513294 1 2045 2045 20 20 97890 97890 1957800
109 324881 0 5477 0 59 0 101481 0 5987379
110 264512 1 1932 1932 36 36 72954 72954 2626344
111 420968 0 936 0 30 0 67939 0 2038170
112 129302 1 3437 3437 47 47 48022 48022 2257034
113 191521 0 5131 0 71 0 86111 0 6113881
114 268673 1 2397 2397 14 14 74020 74020 1036280
115 353179 1 1389 1389 9 9 57530 57530 517770
116 354624 0 1503 0 39 0 56364 0 2198196
117 363713 0 402 0 26 0 84990 0 2209740
118 456657 0 2239 0 21 0 88590 0 1860390
119 211742 1 2234 2234 16 16 77200 77200 1235200
120 338381 0 837 0 69 0 61262 0 4227078
121 418530 0 10579 0 92 0 110309 0 10148428
122 351483 0 875 0 14 0 67000 0 938000
123 372928 0 1395 0 103 0 93099 0 9589197
124 485538 1 1659 1659 29 29 107577 107577 3119733
125 279268 1 2647 2647 37 37 62920 62920 2328040
126 219060 1 3294 3294 23 23 75832 75832 1744136
127 325560 0 0 0 0 0 60720 0 0
128 325314 0 94 0 7 0 60793 0 425551
129 322046 0 422 0 28 0 57935 0 1622180
130 325560 0 0 0 0 0 60720 0 0
131 325599 0 34 0 8 0 60630 0 485040
132 377028 0 1558 0 63 0 55637 0 3505131
133 325560 0 0 0 0 0 60720 0 0
134 323850 0 43 0 3 0 60887 0 182661
135 325560 0 0 0 0 0 60720 0 0
136 331514 0 316 0 9 0 60505 0 544545
137 325632 0 115 0 13 0 60945 0 792285
138 325560 0 0 0 0 0 60720 0 0
139 325560 0 0 0 0 0 60720 0 0
140 325560 0 0 0 0 0 60720 0 0
141 322265 0 389 0 14 0 58990 0 825860
142 325560 0 0 0 0 0 60720 0 0
143 325906 0 1002 0 15 0 56750 0 851250
144 325985 0 36 0 3 0 60894 0 182682
145 346145 0 460 0 15 0 63346 0 950190
146 325898 0 309 0 11 0 56535 0 621885
147 325560 0 0 0 0 0 60720 0 0
148 325356 0 9 0 6 0 60835 0 365010
149 325560 0 0 0 0 0 60720 0 0
150 325930 0 14 0 1 0 61016 0 61016
151 318020 0 520 0 10 0 58650 0 586500
152 326389 0 1766 0 73 0 60438 0 4411974
153 325560 0 0 0 0 0 60720 0 0
154 302925 0 458 0 11 0 58625 0 644875
155 325540 0 20 0 3 0 60938 0 182814
156 325560 0 0 0 0 0 60720 0 0
157 325560 0 0 0 0 0 60720 0 0
158 326736 0 98 0 2 0 61490 0 122980
159 340580 0 405 0 7 0 60845 0 425915
160 325560 0 0 0 0 0 60720 0 0
161 325560 0 0 0 0 0 60720 0 0
162 325560 0 0 0 0 0 60720 0 0
163 325560 0 0 0 0 0 60720 0 0
164 331828 0 483 0 27 0 60830 0 1642410
165 323299 0 454 0 51 0 63261 0 3226311
166 325560 0 0 0 0 0 60720 0 0
167 325560 0 0 0 0 0 60720 0 0
168 387722 0 757 0 19 0 45689 0 868091
169 325560 0 0 0 0 0 60720 0 0
170 325560 0 0 0 0 0 60720 0 0
171 325560 0 0 0 0 0 60720 0 0
172 324598 0 36 0 4 0 61564 0 246256
173 325560 0 0 0 0 0 60720 0 0
174 328726 0 203 0 9 0 61938 0 557442
175 325560 0 0 0 0 0 60720 0 0
176 325043 0 90 0 8 0 60951 0 487608
177 325560 0 0 0 0 0 60720 0 0
178 325806 1 71 71 1 1 60745 60745 60745
179 325560 0 0 0 0 0 60720 0 0
180 325560 0 0 0 0 0 60720 0 0
181 387732 0 972 0 34 0 71642 0 2435828
182 349729 1 531 531 10 10 71641 71641 716410
183 332202 0 604 0 38 0 55792 0 2120096
184 305442 1 283 283 10 10 71873 71873 718730
185 329537 1 23 23 5 5 62555 62555 312775
186 327055 1 638 638 14 14 60370 60370 845180
187 356245 1 699 699 16 16 64873 64873 1037968
188 328451 0 149 0 5 0 62041 0 310205
189 307062 0 226 0 5 0 65745 0 328725
190 325560 0 0 0 0 0 60720 0 0
191 331345 0 275 0 4 0 59500 0 238000
192 325560 0 0 0 0 0 60720 0 0
193 331824 0 141 0 6 0 61630 0 369780
194 325560 0 0 0 0 0 60720 0 0
195 325685 0 28 0 2 0 60890 0 121780
196 325560 1 0 0 0 0 60720 60720 0
197 404480 1 2566 2566 80 80 113521 113521 9081680
198 325560 1 0 0 0 0 60720 60720 0
199 325560 1 0 0 0 0 60720 60720 0
200 318314 1 472 472 20 20 80045 80045 1600900
201 325560 1 0 0 0 0 60720 60720 0
202 325560 1 0 0 0 0 60720 60720 0
203 325560 1 0 0 0 0 60720 60720 0
204 311807 1 203 203 27 27 50804 50804 1371708
205 337724 1 496 496 17 17 87390 87390 1485630
206 326431 1 10 10 2 2 61656 61656 123312
207 327556 1 63 63 4 4 65688 65688 262752
208 325560 1 0 0 0 0 60720 60720 0
209 356850 1 1136 1136 32 32 48522 48522 1552704
210 325560 1 0 0 0 0 60720 60720 0
211 325560 0 0 0 0 0 60720 0 0
212 325560 0 0 0 0 0 60720 0 0
213 322741 0 267 0 32 0 57640 0 1844480
214 310902 0 474 0 20 0 61977 0 1239540
215 324295 0 534 0 7 0 62620 0 438340
216 325560 0 0 0 0 0 60720 0 0
217 326156 0 15 0 8 0 60831 0 486648
218 326960 0 397 0 28 0 60646 0 1698088
219 325560 0 0 0 0 0 60720 0 0
220 333411 0 1061 0 20 0 56225 0 1124500
221 297761 0 288 0 4 0 60510 0 242040
222 325560 0 0 0 0 0 60720 0 0
223 325536 0 3 0 2 0 60698 0 121396
224 325560 0 0 0 0 0 60720 0 0
225 325762 0 20 0 2 0 60805 0 121610
226 327957 0 278 0 26 0 61404 0 1596504
227 325560 0 0 0 0 0 60720 0 0
228 325560 0 0 0 0 0 60720 0 0
229 318521 0 192 0 4 0 65276 0 261104
230 325560 0 0 0 0 0 60720 0 0
231 319775 0 317 0 9 0 63915 0 575235
232 325560 0 0 0 0 0 60720 0 0
233 325560 0 0 0 0 0 60720 0 0
234 332128 0 368 0 17 0 61686 0 1048662
235 325560 0 0 0 0 0 60720 0 0
236 325486 0 2 0 1 0 60743 0 60743
237 325560 0 0 0 0 0 60720 0 0
238 325838 0 53 0 6 0 60349 0 362094
239 325560 0 0 0 0 0 60720 0 0
240 325560 0 0 0 0 0 60720 0 0
241 325560 0 0 0 0 0 60720 0 0
242 331767 0 94 0 3 0 61360 0 184080
243 325560 0 0 0 0 0 60720 0 0
244 324523 0 24 0 8 0 59818 0 478544
245 339995 0 2332 0 4 0 72680 0 290720
246 325560 1 0 0 0 0 60720 60720 0
247 325560 1 0 0 0 0 60720 60720 0
248 319582 0 131 0 11 0 61808 0 679888
249 325560 1 0 0 0 0 60720 60720 0
250 325560 1 0 0 0 0 60720 60720 0
251 307245 0 206 0 9 0 53110 0 477990
252 325560 1 0 0 0 0 60720 60720 0
253 317967 0 167 0 2 0 64245 0 128490
254 331488 0 622 0 73 0 73007 0 5329511
255 335452 0 885 0 85 0 82732 0 7032220
256 325560 0 0 0 0 0 60720 0 0
257 334184 0 365 0 8 0 54820 0 438560
258 313213 0 364 0 35 0 47705 0 1669675
259 325560 1 0 0 0 0 60720 60720 0
260 325560 0 0 0 0 0 60720 0 0
261 325560 1 0 0 0 0 60720 60720 0
262 325560 1 0 0 0 0 60720 60720 0
263 348678 0 226 0 12 0 72835 0 874020
264 328727 0 307 0 15 0 58856 0 882840
265 325560 1 0 0 0 0 60720 60720 0
266 325560 0 0 0 0 0 60720 0 0
267 325560 1 0 0 0 0 60720 60720 0
268 387978 0 188 0 11 0 77655 0 854205
269 325560 1 0 0 0 0 60720 60720 0
270 336704 0 138 0 6 0 69817 0 418902
271 325560 1 0 0 0 0 60720 60720 0
272 325560 1 0 0 0 0 60720 60720 0
273 325560 1 0 0 0 0 60720 60720 0
274 322076 0 125 0 12 0 60798 0 729576
275 325560 0 0 0 0 0 60720 0 0
276 334272 0 282 0 30 0 62452 0 1873560
277 338197 0 335 0 33 0 64175 0 2117775
278 325560 0 0 0 0 0 60720 0 0
279 321024 1 813 813 82 82 67440 67440 5530080
280 322145 0 176 0 28 0 68136 0 1907808
281 325560 1 0 0 0 0 60720 60720 0
282 325560 0 0 0 0 0 60720 0 0
283 323351 0 249 0 72 0 56726 0 4084272
284 325560 0 0 0 0 0 60720 0 0
285 327748 0 333 0 13 0 70811 0 920543
286 325560 0 0 0 0 0 60720 0 0
287 325560 1 0 0 0 0 60720 60720 0
288 328157 0 30 0 4 0 62045 0 248180
289 325560 1 0 0 0 0 60720 60720 0
290 311594 0 249 0 62 0 54323 0 3368026
291 325560 1 0 0 0 0 60720 60720 0
292 335962 0 165 0 24 0 62841 0 1508184
293 372426 0 453 0 21 0 81125 0 1703625
294 325560 0 0 0 0 0 60720 0 0
295 319844 0 53 0 14 0 59506 0 833084
296 355822 1 382 382 21 21 59365 59365 1246665
297 325560 0 0 0 0 0 60720 0 0
298 325560 1 0 0 0 0 60720 60720 0
299 325560 0 0 0 0 0 60720 0 0
300 325560 0 0 0 0 0 60720 0 0
301 324047 1 30 30 4 4 60798 60798 243192
302 311464 0 290 0 2 0 58790 0 117580
303 325560 0 0 0 0 0 60720 0 0
304 325560 1 0 0 0 0 60720 60720 0
305 353417 0 366 0 53 0 61808 0 3275824
306 325590 0 2 0 9 0 60735 0 546615
307 325560 1 0 0 0 0 60720 60720 0
308 328576 1 209 209 13 13 64016 64016 832208
309 326126 0 384 0 22 0 54683 0 1203026
310 325560 1 0 0 0 0 60720 60720 0
311 325560 1 0 0 0 0 60720 60720 0
312 369376 0 365 0 83 0 87192 0 7236936
313 325560 1 0 0 0 0 60720 60720 0
314 332013 1 49 49 8 8 64107 64107 512856
315 325871 0 3 0 4 0 60761 0 243044
316 342165 0 133 0 14 0 65990 0 923860
317 324967 0 32 0 1 0 59988 0 59988
318 314832 0 368 0 17 0 61167 0 1039839
319 325557 0 1 0 6 0 60719 0 364314
320 325560 0 0 0 0 0 60720 0 0
321 325560 0 0 0 0 0 60720 0 0
322 325560 0 0 0 0 0 60720 0 0
323 325560 0 0 0 0 0 60720 0 0
324 325560 1 0 0 0 0 60720 60720 0
325 325560 0 0 0 0 0 60720 0 0
326 322649 0 22 0 2 0 60722 0 121444
327 325560 1 0 0 0 0 60720 60720 0
328 325560 1 0 0 0 0 60720 60720 0
329 325560 0 0 0 0 0 60720 0 0
330 325560 0 0 0 0 0 60720 0 0
331 325560 0 0 0 0 0 60720 0 0
332 325560 0 0 0 0 0 60720 0 0
333 325560 0 0 0 0 0 60720 0 0
334 324598 0 96 0 5 0 60379 0 301895
335 325567 0 1 0 2 0 60727 0 121454
336 325560 0 0 0 0 0 60720 0 0
337 324005 0 81 0 7 0 60925 0 426475
338 325560 0 0 0 0 0 60720 0 0
339 325748 0 26 0 1 0 60896 0 60896
340 323385 0 125 0 13 0 59734 0 776542
341 315409 0 304 0 15 0 62969 0 944535
342 325560 0 0 0 0 0 60720 0 0
343 325560 0 0 0 0 0 60720 0 0
344 325560 0 0 0 0 0 60720 0 0
345 325560 0 0 0 0 0 60720 0 0
346 325560 0 0 0 0 0 60720 0 0
347 312275 0 119 0 6 0 59118 0 354708
348 325560 0 0 0 0 0 60720 0 0
349 325560 0 0 0 0 0 60720 0 0
350 325560 0 0 0 0 0 60720 0 0
351 320576 0 312 0 14 0 58598 0 820372
352 325246 0 60 0 10 0 61124 0 611240
353 332961 0 587 0 12 0 59595 0 715140
354 323010 0 135 0 2 0 62065 0 124130
355 325560 0 0 0 0 0 60720 0 0
356 325560 0 0 0 0 0 60720 0 0
357 345253 0 514 0 52 0 78780 0 4096560
358 325560 0 0 0 0 0 60720 0 0
359 325560 0 0 0 0 0 60720 0 0
360 325560 0 0 0 0 0 60720 0 0
361 325559 0 1 0 4 0 60722 0 242888
362 325560 0 0 0 0 0 60720 0 0
363 325560 0 0 0 0 0 60720 0 0
364 319634 1 58 58 3 3 61600 61600 184800
365 319951 0 180 0 11 0 59635 0 655985
366 325560 0 0 0 0 0 60720 0 0
367 325560 0 0 0 0 0 60720 0 0
368 325560 0 0 0 0 0 60720 0 0
369 325560 0 0 0 0 0 60720 0 0
370 325560 0 0 0 0 0 60720 0 0
371 325560 0 0 0 0 0 60720 0 0
372 318519 0 448 0 40 0 59781 0 2391240
373 343222 0 227 0 9 0 76644 0 689796
374 317234 0 174 0 1 0 64820 0 64820
375 325560 0 0 0 0 0 60720 0 0
376 325560 0 0 0 0 0 60720 0 0
377 314025 0 121 0 24 0 56178 0 1348272
378 320249 0 607 0 11 0 60436 0 664796
379 325560 0 0 0 0 0 60720 0 0
380 325560 0 0 0 0 0 60720 0 0
381 325560 0 0 0 0 0 60720 0 0
382 349365 0 530 0 60 0 73433 0 4405980
383 289197 0 571 0 80 0 41477 0 3318160
384 325560 0 0 0 0 0 60720 0 0
385 329245 0 78 0 16 0 62700 0 1003200
386 240869 0 2489 0 40 0 67804 0 2712160
387 327182 0 131 0 6 0 59661 0 357966
388 322876 0 923 0 8 0 58620 0 468960
389 323117 0 72 0 3 0 60398 0 181194
390 306351 0 572 0 16 0 58580 0 937280
391 335137 0 397 0 10 0 62710 0 627100
392 308271 0 450 0 8 0 59325 0 474600
393 301731 0 622 0 7 0 60950 0 426650
394 382409 0 694 0 8 0 68060 0 544480
395 279230 1 3425 3425 12 12 83620 83620 1003440
396 298731 0 562 0 13 0 58456 0 759928
397 243650 0 4917 0 42 0 52811 0 2218062
398 532682 1 1442 1442 118 118 121173 121173 14298414
399 319771 0 529 0 9 0 63870 0 574830
400 171493 1 2126 2126 138 138 21001 21001 2898138
401 347262 0 1061 0 5 0 70415 0 352075
402 343945 0 776 0 9 0 64230 0 578070
403 311874 0 611 0 8 0 59190 0 473520
404 302211 1 1526 1526 25 25 69351 69351 1733775
405 316708 0 592 0 7 0 64270 0 449890
406 333463 0 1182 0 13 0 70694 0 919022
407 344282 0 621 0 16 0 68005 0 1088080
408 319635 0 989 0 11 0 58930 0 648230
409 301186 0 438 0 11 0 58320 0 641520
410 300381 0 726 0 3 0 69980 0 209940
411 318765 0 1303 0 61 0 69863 0 4261643
412 286146 1 6341 6341 24 24 63255 63255 1518120
413 306844 1 1164 1164 17 17 57320 57320 974440
414 307705 1 3310 3310 33 33 75230 75230 2482590
415 312448 0 1366 0 7 0 79420 0 555940
416 299715 0 965 0 3 0 73490 0 220470
417 373399 0 3256 0 66 0 35250 0 2326500
418 299446 1 1135 1135 17 17 62285 62285 1058845
419 325586 0 1270 0 26 0 69206 0 1799356
420 291221 0 661 0 3 0 65920 0 197760
421 261173 0 1013 0 2 0 69770 0 139540
422 255027 0 2844 0 67 0 72683 0 4869761
423 -78375 1 11528 11528 70 70 -14545 -14545 -1018150
424 -58143 0 6526 0 26 0 55830 0 1451580
425 227033 0 2264 0 24 0 55174 0 1324176
426 235098 1 4461 4461 94 94 67038 67038 6301572
427 21267 0 3999 0 30 0 51252 0 1537560
428 238675 0 35624 0 223 0 157278 0 35072994
429 197687 0 9252 0 48 0 79510 0 3816480
430 418341 0 15236 0 90 0 77440 0 6969600
431 -297706 0 18073 0 180 0 27284 0 4911120
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) Group Costs GrCosts Trades GrTrades
1.638e+05 7.066e+04 -5.434e+00 3.857e+01 -5.819e+02 -7.036e+01
Dividends GrDiv TrDiv
2.704e+00 -1.643e+00 6.407e-03
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-743164 -18133 -2401 21979 3055410
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.638e+05 3.431e+04 4.774 2.50e-06 ***
Group 7.066e+04 5.587e+04 1.265 0.20671
Costs -5.434e+00 2.425e+00 -2.241 0.02553 *
GrCosts 3.857e+01 3.744e+00 10.302 < 2e-16 ***
Trades -5.819e+02 2.781e+02 -2.092 0.03703 *
GrTrades -7.036e+01 3.405e+02 -0.207 0.83636
Dividends 2.704e+00 4.765e-01 5.675 2.58e-08 ***
GrDiv -1.643e+00 7.088e-01 -2.318 0.02091 *
TrDiv 6.407e-03 2.115e-03 3.029 0.00261 **
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 208400 on 422 degrees of freedom
Multiple R-squared: 0.7871, Adjusted R-squared: 0.7831
F-statistic: 195 on 8 and 422 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 7.851251e-100 3.925626e-100
[2,] 1.0000000 2.191165e-113 1.095583e-113
[3,] 1.0000000 6.072675e-121 3.036338e-121
[4,] 1.0000000 2.671786e-125 1.335893e-125
[5,] 1.0000000 2.731368e-135 1.365684e-135
[6,] 1.0000000 1.148300e-137 5.741501e-138
[7,] 1.0000000 7.993861e-139 3.996931e-139
[8,] 1.0000000 7.659853e-146 3.829927e-146
[9,] 1.0000000 6.118071e-148 3.059036e-148
[10,] 1.0000000 8.791256e-154 4.395628e-154
[11,] 1.0000000 3.525189e-158 1.762594e-158
[12,] 1.0000000 4.274488e-159 2.137244e-159
[13,] 1.0000000 1.199647e-159 5.998233e-160
[14,] 1.0000000 6.364360e-159 3.182180e-159
[15,] 1.0000000 9.999843e-161 4.999921e-161
[16,] 1.0000000 4.209709e-164 2.104855e-164
[17,] 1.0000000 2.856631e-166 1.428315e-166
[18,] 1.0000000 5.112184e-179 2.556092e-179
[19,] 1.0000000 1.071187e-180 5.355933e-181
[20,] 1.0000000 2.749321e-184 1.374661e-184
[21,] 1.0000000 4.370172e-186 2.185086e-186
[22,] 1.0000000 2.256602e-192 1.128301e-192
[23,] 1.0000000 7.310422e-194 3.655211e-194
[24,] 1.0000000 7.134161e-193 3.567081e-193
[25,] 1.0000000 1.271116e-192 6.355578e-193
[26,] 1.0000000 7.610557e-200 3.805278e-200
[27,] 1.0000000 8.528374e-205 4.264187e-205
[28,] 1.0000000 4.029553e-204 2.014777e-204
[29,] 1.0000000 2.044370e-204 1.022185e-204
[30,] 1.0000000 2.069842e-206 1.034921e-206
[31,] 1.0000000 1.001833e-206 5.009163e-207
[32,] 1.0000000 7.168054e-207 3.584027e-207
[33,] 1.0000000 4.672260e-206 2.336130e-206
[34,] 1.0000000 1.134090e-205 5.670451e-206
[35,] 1.0000000 1.986613e-205 9.933063e-206
[36,] 1.0000000 2.011350e-204 1.005675e-204
[37,] 1.0000000 5.774158e-204 2.887079e-204
[38,] 1.0000000 2.333566e-203 1.166783e-203
[39,] 1.0000000 8.272180e-204 4.136090e-204
[40,] 1.0000000 7.636110e-205 3.818055e-205
[41,] 1.0000000 5.339899e-204 2.669949e-204
[42,] 1.0000000 2.070578e-203 1.035289e-203
[43,] 1.0000000 2.474863e-202 1.237432e-202
[44,] 1.0000000 2.854662e-207 1.427331e-207
[45,] 1.0000000 7.652691e-208 3.826346e-208
[46,] 1.0000000 2.359341e-212 1.179670e-212
[47,] 1.0000000 2.011848e-213 1.005924e-213
[48,] 1.0000000 1.527366e-212 7.636832e-213
[49,] 1.0000000 3.523830e-213 1.761915e-213
[50,] 1.0000000 1.206829e-214 6.034144e-215
[51,] 1.0000000 2.840292e-214 1.420146e-214
[52,] 1.0000000 1.118817e-213 5.594083e-214
[53,] 1.0000000 1.511618e-213 7.558088e-214
[54,] 1.0000000 3.172995e-213 1.586497e-213
[55,] 1.0000000 3.768960e-213 1.884480e-213
[56,] 1.0000000 3.646202e-214 1.823101e-214
[57,] 1.0000000 3.452764e-213 1.726382e-213
[58,] 1.0000000 1.532771e-212 7.663855e-213
[59,] 1.0000000 6.261774e-212 3.130887e-212
[60,] 1.0000000 1.366611e-211 6.833056e-212
[61,] 1.0000000 6.446728e-214 3.223364e-214
[62,] 1.0000000 1.523090e-213 7.615452e-214
[63,] 1.0000000 1.363732e-212 6.818661e-213
[64,] 1.0000000 1.055586e-211 5.277929e-212
[65,] 1.0000000 3.963795e-211 1.981897e-211
[66,] 1.0000000 1.178119e-210 5.890596e-211
[67,] 1.0000000 2.056105e-210 1.028052e-210
[68,] 1.0000000 2.275088e-209 1.137544e-209
[69,] 1.0000000 1.168665e-210 5.843323e-211
[70,] 1.0000000 8.805335e-210 4.402668e-210
[71,] 1.0000000 7.342751e-210 3.671376e-210
[72,] 1.0000000 3.768693e-210 1.884347e-210
[73,] 1.0000000 3.739645e-209 1.869822e-209
[74,] 1.0000000 3.589547e-208 1.794774e-208
[75,] 1.0000000 4.122509e-209 2.061254e-209
[76,] 1.0000000 3.914147e-209 1.957074e-209
[77,] 1.0000000 4.008049e-208 2.004024e-208
[78,] 1.0000000 1.199480e-207 5.997398e-208
[79,] 1.0000000 1.205720e-206 6.028602e-207
[80,] 1.0000000 1.288189e-205 6.440946e-206
[81,] 1.0000000 5.512733e-205 2.756367e-205
[82,] 1.0000000 7.561750e-206 3.780875e-206
[83,] 1.0000000 4.113166e-205 2.056583e-205
[84,] 1.0000000 5.673904e-206 2.836952e-206
[85,] 1.0000000 1.464725e-206 7.323624e-207
[86,] 1.0000000 6.771670e-206 3.385835e-206
[87,] 1.0000000 6.783791e-205 3.391896e-205
[88,] 1.0000000 5.970192e-204 2.985096e-204
[89,] 1.0000000 5.683514e-203 2.841757e-203
[90,] 1.0000000 5.157830e-202 2.578915e-202
[91,] 1.0000000 5.349010e-202 2.674505e-202
[92,] 1.0000000 2.878033e-201 1.439016e-201
[93,] 1.0000000 2.505846e-200 1.252923e-200
[94,] 1.0000000 1.052635e-199 5.263176e-200
[95,] 1.0000000 2.389309e-199 1.194655e-199
[96,] 1.0000000 2.079708e-200 1.039854e-200
[97,] 1.0000000 2.845376e-203 1.422688e-203
[98,] 1.0000000 2.353012e-202 1.176506e-202
[99,] 1.0000000 9.211686e-202 4.605843e-202
[100,] 1.0000000 4.993942e-202 2.496971e-202
[101,] 1.0000000 2.918974e-203 1.459487e-203
[102,] 1.0000000 6.757189e-205 3.378595e-205
[103,] 1.0000000 4.220176e-204 2.110088e-204
[104,] 1.0000000 1.434609e-203 7.173043e-204
[105,] 1.0000000 3.413438e-203 1.706719e-203
[106,] 1.0000000 3.329291e-202 1.664645e-202
[107,] 1.0000000 4.451413e-203 2.225707e-203
[108,] 1.0000000 8.239628e-204 4.119814e-204
[109,] 1.0000000 7.876436e-203 3.938218e-203
[110,] 1.0000000 4.385686e-204 2.192843e-204
[111,] 1.0000000 3.354392e-203 1.677196e-203
[112,] 1.0000000 1.351314e-202 6.756568e-203
[113,] 1.0000000 3.493354e-204 1.746677e-204
[114,] 1.0000000 3.420620e-203 1.710310e-203
[115,] 1.0000000 4.386036e-203 2.193018e-203
[116,] 1.0000000 4.941792e-202 2.470896e-202
[117,] 1.0000000 5.546052e-201 2.773026e-201
[118,] 1.0000000 6.158798e-200 3.079399e-200
[119,] 1.0000000 6.827413e-199 3.413707e-199
[120,] 1.0000000 7.536911e-198 3.768456e-198
[121,] 1.0000000 3.194940e-198 1.597470e-198
[122,] 1.0000000 3.521270e-197 1.760635e-197
[123,] 1.0000000 3.851471e-196 1.925736e-196
[124,] 1.0000000 4.200029e-195 2.100014e-195
[125,] 1.0000000 4.357628e-194 2.178814e-194
[126,] 1.0000000 4.711477e-193 2.355739e-193
[127,] 1.0000000 5.056614e-192 2.528307e-192
[128,] 1.0000000 5.398251e-191 2.699126e-191
[129,] 1.0000000 5.732416e-190 2.866208e-190
[130,] 1.0000000 6.013174e-189 3.006587e-189
[131,] 1.0000000 6.317827e-188 3.158913e-188
[132,] 1.0000000 5.250876e-187 2.625438e-187
[133,] 1.0000000 5.466402e-186 2.733201e-186
[134,] 1.0000000 4.527111e-185 2.263556e-185
[135,] 1.0000000 4.460852e-184 2.230426e-184
[136,] 1.0000000 4.571938e-183 2.285969e-183
[137,] 1.0000000 4.661299e-182 2.330649e-182
[138,] 1.0000000 4.727525e-181 2.363763e-181
[139,] 1.0000000 4.770023e-180 2.385012e-180
[140,] 1.0000000 4.760345e-179 2.380172e-179
[141,] 1.0000000 3.289288e-178 1.644644e-178
[142,] 1.0000000 3.271132e-177 1.635566e-177
[143,] 1.0000000 3.039135e-176 1.519567e-176
[144,] 1.0000000 2.992293e-175 1.496147e-175
[145,] 1.0000000 2.930747e-174 1.465373e-174
[146,] 1.0000000 2.855705e-173 1.427853e-173
[147,] 1.0000000 2.767403e-172 1.383701e-172
[148,] 1.0000000 2.232514e-171 1.116257e-171
[149,] 1.0000000 2.143563e-170 1.071782e-170
[150,] 1.0000000 2.047659e-169 1.023829e-169
[151,] 1.0000000 1.946088e-168 9.730439e-169
[152,] 1.0000000 1.840156e-167 9.200782e-168
[153,] 1.0000000 1.619920e-166 8.099599e-167
[154,] 1.0000000 1.515315e-165 7.576574e-166
[155,] 1.0000000 1.411462e-164 7.057310e-165
[156,] 1.0000000 1.308099e-163 6.540497e-164
[157,] 1.0000000 1.999802e-164 9.999011e-165
[158,] 1.0000000 1.864968e-163 9.324838e-164
[159,] 1.0000000 1.730530e-162 8.652649e-163
[160,] 1.0000000 1.597768e-161 7.988842e-162
[161,] 1.0000000 1.465500e-160 7.327501e-161
[162,] 1.0000000 1.339647e-159 6.698234e-160
[163,] 1.0000000 1.208928e-158 6.044638e-159
[164,] 1.0000000 1.094183e-157 5.470916e-158
[165,] 1.0000000 9.856067e-157 4.928033e-157
[166,] 1.0000000 8.832444e-156 4.416222e-156
[167,] 1.0000000 7.892415e-155 3.946208e-155
[168,] 1.0000000 7.002654e-154 3.501327e-154
[169,] 1.0000000 6.182615e-153 3.091308e-153
[170,] 1.0000000 1.699175e-152 8.495876e-153
[171,] 1.0000000 1.273459e-151 6.367296e-152
[172,] 1.0000000 8.111266e-151 4.055633e-151
[173,] 1.0000000 5.208978e-150 2.604489e-150
[174,] 1.0000000 4.530619e-149 2.265309e-149
[175,] 1.0000000 3.778162e-148 1.889081e-148
[176,] 1.0000000 2.040727e-147 1.020363e-147
[177,] 1.0000000 1.741004e-146 8.705019e-147
[178,] 1.0000000 1.228971e-145 6.144856e-146
[179,] 1.0000000 1.044676e-144 5.223381e-145
[180,] 1.0000000 8.320041e-144 4.160021e-144
[181,] 1.0000000 7.005587e-143 3.502793e-143
[182,] 1.0000000 5.752095e-142 2.876047e-142
[183,] 1.0000000 4.796909e-141 2.398455e-141
[184,] 1.0000000 3.980739e-140 1.990370e-140
[185,] 1.0000000 3.290844e-139 1.645422e-139
[186,] 1.0000000 2.390340e-138 1.195170e-138
[187,] 1.0000000 1.958266e-137 9.791331e-138
[188,] 1.0000000 1.596482e-136 7.982410e-137
[189,] 1.0000000 9.867813e-136 4.933906e-136
[190,] 1.0000000 7.973688e-135 3.986844e-135
[191,] 1.0000000 6.411728e-134 3.205864e-134
[192,] 1.0000000 5.130600e-133 2.565300e-133
[193,] 1.0000000 4.089697e-132 2.044848e-132
[194,] 1.0000000 2.919224e-131 1.459612e-131
[195,] 1.0000000 2.301929e-130 1.150965e-130
[196,] 1.0000000 1.796011e-129 8.980055e-130
[197,] 1.0000000 1.402682e-128 7.013411e-129
[198,] 1.0000000 2.338485e-128 1.169242e-128
[199,] 1.0000000 1.824807e-127 9.124033e-128
[200,] 1.0000000 1.415115e-126 7.075575e-127
[201,] 1.0000000 1.092154e-125 5.460771e-126
[202,] 1.0000000 8.122477e-125 4.061239e-125
[203,] 1.0000000 6.072738e-124 3.036369e-124
[204,] 1.0000000 4.567646e-123 2.283823e-123
[205,] 1.0000000 3.459292e-122 1.729646e-122
[206,] 1.0000000 2.607234e-121 1.303617e-121
[207,] 1.0000000 1.897200e-120 9.486001e-121
[208,] 1.0000000 1.416463e-119 7.082315e-120
[209,] 1.0000000 6.630758e-119 3.315379e-119
[210,] 1.0000000 4.210428e-118 2.105214e-118
[211,] 1.0000000 3.108270e-117 1.554135e-117
[212,] 1.0000000 2.283718e-116 1.141859e-116
[213,] 1.0000000 1.669771e-115 8.348853e-116
[214,] 1.0000000 1.214805e-114 6.074023e-115
[215,] 1.0000000 8.659823e-114 4.329911e-114
[216,] 1.0000000 6.241285e-113 3.120642e-113
[217,] 1.0000000 4.476565e-112 2.238283e-112
[218,] 1.0000000 3.070093e-111 1.535047e-111
[219,] 1.0000000 2.181458e-110 1.090729e-110
[220,] 1.0000000 1.527339e-109 7.636694e-110
[221,] 1.0000000 1.074895e-108 5.374477e-109
[222,] 1.0000000 7.528258e-108 3.764129e-108
[223,] 1.0000000 5.010585e-107 2.505293e-107
[224,] 1.0000000 3.476160e-106 1.738080e-106
[225,] 1.0000000 2.400047e-105 1.200023e-105
[226,] 1.0000000 1.648967e-104 8.244835e-105
[227,] 1.0000000 1.125397e-103 5.626986e-104
[228,] 1.0000000 7.657198e-103 3.828599e-103
[229,] 1.0000000 5.184612e-102 2.592306e-102
[230,] 1.0000000 3.493353e-101 1.746676e-101
[231,] 1.0000000 2.304431e-100 1.152215e-100
[232,] 1.0000000 1.537686e-99 7.688429e-100
[233,] 1.0000000 1.020333e-98 5.101666e-99
[234,] 1.0000000 4.464562e-98 2.232281e-98
[235,] 1.0000000 2.944635e-97 1.472318e-97
[236,] 1.0000000 1.932582e-96 9.662911e-97
[237,] 1.0000000 1.253793e-95 6.268963e-96
[238,] 1.0000000 8.147958e-95 4.073979e-95
[239,] 1.0000000 5.268838e-94 2.634419e-94
[240,] 1.0000000 3.390350e-93 1.695175e-93
[241,] 1.0000000 2.170645e-92 1.085322e-92
[242,] 1.0000000 1.346710e-91 6.733550e-92
[243,] 1.0000000 7.624759e-91 3.812380e-91
[244,] 1.0000000 2.223465e-90 1.111733e-90
[245,] 1.0000000 1.401791e-89 7.008956e-90
[246,] 1.0000000 7.100697e-89 3.550348e-89
[247,] 1.0000000 3.679878e-88 1.839939e-88
[248,] 1.0000000 2.293719e-87 1.146859e-87
[249,] 1.0000000 1.420145e-86 7.100723e-87
[250,] 1.0000000 8.762918e-86 4.381459e-86
[251,] 1.0000000 5.379695e-85 2.689847e-85
[252,] 1.0000000 3.275713e-84 1.637856e-84
[253,] 1.0000000 1.892728e-83 9.463641e-84
[254,] 1.0000000 1.144745e-82 5.723726e-83
[255,] 1.0000000 6.876357e-82 3.438178e-82
[256,] 1.0000000 4.116488e-81 2.058244e-81
[257,] 1.0000000 1.874660e-80 9.373301e-81
[258,] 1.0000000 1.112719e-79 5.563593e-80
[259,] 1.0000000 6.554998e-79 3.277499e-79
[260,] 1.0000000 3.850680e-78 1.925340e-78
[261,] 1.0000000 2.250268e-77 1.125134e-77
[262,] 1.0000000 1.308141e-76 6.540707e-77
[263,] 1.0000000 7.563932e-76 3.781966e-76
[264,] 1.0000000 4.344533e-75 2.172267e-75
[265,] 1.0000000 2.432268e-74 1.216134e-74
[266,] 1.0000000 1.351404e-73 6.757022e-74
[267,] 1.0000000 7.642411e-73 3.821205e-73
[268,] 1.0000000 3.979334e-72 1.989667e-72
[269,] 1.0000000 2.073574e-71 1.036787e-71
[270,] 1.0000000 1.157067e-70 5.785337e-71
[271,] 1.0000000 6.413274e-70 3.206637e-70
[272,] 1.0000000 3.424981e-69 1.712490e-69
[273,] 1.0000000 1.878607e-68 9.393035e-69
[274,] 1.0000000 9.962060e-68 4.981030e-68
[275,] 1.0000000 5.406997e-67 2.703499e-67
[276,] 1.0000000 2.922362e-66 1.461181e-66
[277,] 1.0000000 1.568347e-65 7.841735e-66
[278,] 1.0000000 8.383805e-65 4.191902e-65
[279,] 1.0000000 4.385742e-64 2.192871e-64
[280,] 1.0000000 2.318883e-63 1.159441e-63
[281,] 1.0000000 1.196050e-62 5.980249e-63
[282,] 1.0000000 6.092224e-62 3.046112e-62
[283,] 1.0000000 3.165327e-61 1.582663e-61
[284,] 1.0000000 1.637804e-60 8.189018e-61
[285,] 1.0000000 6.441852e-60 3.220926e-60
[286,] 1.0000000 3.298147e-59 1.649074e-59
[287,] 1.0000000 1.681981e-58 8.409904e-59
[288,] 1.0000000 8.513785e-58 4.256892e-58
[289,] 1.0000000 4.284748e-57 2.142374e-57
[290,] 1.0000000 2.147118e-56 1.073559e-56
[291,] 1.0000000 1.067749e-55 5.338743e-56
[292,] 1.0000000 5.281349e-55 2.640674e-55
[293,] 1.0000000 2.601522e-54 1.300761e-54
[294,] 1.0000000 1.036493e-53 5.182465e-54
[295,] 1.0000000 5.047364e-53 2.523682e-53
[296,] 1.0000000 2.447067e-52 1.223534e-52
[297,] 1.0000000 1.179226e-51 5.896132e-52
[298,] 1.0000000 5.062402e-51 2.531201e-51
[299,] 1.0000000 2.412903e-50 1.206452e-50
[300,] 1.0000000 1.143088e-49 5.715441e-50
[301,] 1.0000000 3.602734e-49 1.801367e-49
[302,] 1.0000000 1.692441e-48 8.462204e-49
[303,] 1.0000000 7.891565e-48 3.945783e-48
[304,] 1.0000000 3.653582e-47 1.826791e-47
[305,] 1.0000000 1.664452e-46 8.322259e-47
[306,] 1.0000000 7.589775e-46 3.794887e-46
[307,] 1.0000000 3.447869e-45 1.723934e-45
[308,] 1.0000000 1.557706e-44 7.788532e-45
[309,] 1.0000000 6.984510e-44 3.492255e-44
[310,] 1.0000000 3.111567e-43 1.555783e-43
[311,] 1.0000000 1.377199e-42 6.885994e-43
[312,] 1.0000000 6.055769e-42 3.027884e-42
[313,] 1.0000000 2.652684e-41 1.326342e-41
[314,] 1.0000000 1.151144e-40 5.755722e-41
[315,] 1.0000000 4.973864e-40 2.486932e-40
[316,] 1.0000000 2.135430e-39 1.067715e-39
[317,] 1.0000000 9.104121e-39 4.552060e-39
[318,] 1.0000000 3.845887e-38 1.922944e-38
[319,] 1.0000000 1.613490e-37 8.067452e-38
[320,] 1.0000000 6.722403e-37 3.361201e-37
[321,] 1.0000000 2.781298e-36 1.390649e-36
[322,] 1.0000000 1.142645e-35 5.713223e-36
[323,] 1.0000000 4.655590e-35 2.327795e-35
[324,] 1.0000000 1.886435e-34 9.432175e-35
[325,] 1.0000000 7.585461e-34 3.792731e-34
[326,] 1.0000000 3.033512e-33 1.516756e-33
[327,] 1.0000000 1.202116e-32 6.010580e-33
[328,] 1.0000000 4.727634e-32 2.363817e-32
[329,] 1.0000000 1.845013e-31 9.225067e-32
[330,] 1.0000000 7.120755e-31 3.560377e-31
[331,] 1.0000000 2.739196e-30 1.369598e-30
[332,] 1.0000000 1.045648e-29 5.228241e-30
[333,] 1.0000000 3.960789e-29 1.980395e-29
[334,] 1.0000000 1.488602e-28 7.443009e-29
[335,] 1.0000000 5.550629e-28 2.775315e-28
[336,] 1.0000000 2.058112e-27 1.029056e-27
[337,] 1.0000000 7.553423e-27 3.776711e-27
[338,] 1.0000000 2.749663e-26 1.374831e-26
[339,] 1.0000000 9.927456e-26 4.963728e-26
[340,] 1.0000000 3.531216e-25 1.765608e-25
[341,] 1.0000000 1.256723e-24 6.283614e-25
[342,] 1.0000000 4.127511e-24 2.063755e-24
[343,] 1.0000000 1.446441e-23 7.232207e-24
[344,] 1.0000000 5.008317e-23 2.504158e-23
[345,] 1.0000000 1.718906e-22 8.594531e-23
[346,] 1.0000000 5.164669e-22 2.582335e-22
[347,] 1.0000000 1.743900e-21 8.719502e-22
[348,] 1.0000000 5.834753e-21 2.917376e-21
[349,] 1.0000000 1.934160e-20 9.670800e-21
[350,] 1.0000000 6.363140e-20 3.181570e-20
[351,] 1.0000000 2.069828e-19 1.034914e-19
[352,] 1.0000000 6.668038e-19 3.334019e-19
[353,] 1.0000000 2.139356e-18 1.069678e-18
[354,] 1.0000000 6.782242e-18 3.391121e-18
[355,] 1.0000000 2.121598e-17 1.060799e-17
[356,] 1.0000000 6.569150e-17 3.284575e-17
[357,] 1.0000000 2.012988e-16 1.006494e-16
[358,] 1.0000000 6.103593e-16 3.051796e-16
[359,] 1.0000000 1.830904e-15 9.154522e-16
[360,] 1.0000000 5.432516e-15 2.716258e-15
[361,] 1.0000000 1.608125e-14 8.040626e-15
[362,] 1.0000000 4.629544e-14 2.314772e-14
[363,] 1.0000000 1.336757e-13 6.683787e-14
[364,] 1.0000000 3.797814e-13 1.898907e-13
[365,] 1.0000000 1.066256e-12 5.331282e-13
[366,] 1.0000000 2.972788e-12 1.486394e-12
[367,] 1.0000000 8.084023e-12 4.042012e-12
[368,] 1.0000000 2.187892e-11 1.093946e-11
[369,] 1.0000000 5.845773e-11 2.922886e-11
[370,] 1.0000000 1.541517e-10 7.707586e-11
[371,] 1.0000000 3.873563e-10 1.936781e-10
[372,] 1.0000000 9.260512e-10 4.630256e-10
[373,] 1.0000000 2.354000e-09 1.177000e-09
[374,] 1.0000000 5.951081e-09 2.975541e-09
[375,] 1.0000000 1.213660e-08 6.068299e-09
[376,] 1.0000000 2.948997e-08 1.474499e-08
[377,] 1.0000000 6.727029e-08 3.363515e-08
[378,] 0.9999999 1.603641e-07 8.018205e-08
[379,] 0.9999998 3.792423e-07 1.896212e-07
[380,] 0.9999996 8.630042e-07 4.315021e-07
[381,] 0.9999990 1.973840e-06 9.869201e-07
[382,] 0.9999978 4.469783e-06 2.234892e-06
[383,] 0.9999958 8.448100e-06 4.224050e-06
[384,] 0.9999914 1.722933e-05 8.614666e-06
[385,] 0.9999814 3.722305e-05 1.861153e-05
[386,] 0.9999656 6.872351e-05 3.436175e-05
[387,] 0.9999299 1.402851e-04 7.014253e-05
[388,] 0.9998565 2.869441e-04 1.434721e-04
[389,] 0.9997498 5.003118e-04 2.501559e-04
[390,] 0.9995169 9.661836e-04 4.830918e-04
[391,] 0.9991177 1.764696e-03 8.823479e-04
[392,] 0.9983704 3.259159e-03 1.629580e-03
[393,] 0.9969592 6.081585e-03 3.040793e-03
[394,] 0.9945000 1.099991e-02 5.499953e-03
[395,] 0.9902795 1.944102e-02 9.720508e-03
[396,] 0.9834526 3.309474e-02 1.654737e-02
[397,] 0.9740392 5.192155e-02 2.596078e-02
[398,] 0.9583871 8.322573e-02 4.161287e-02
[399,] 0.9330291 1.339418e-01 6.697091e-02
[400,] 0.8963089 2.073822e-01 1.036911e-01
[401,] 0.8463392 3.073216e-01 1.536608e-01
[402,] 0.7773996 4.452009e-01 2.226004e-01
[403,] 0.6896033 6.207933e-01 3.103967e-01
[404,] 0.5858268 8.283464e-01 4.141732e-01
[405,] 0.4683597 9.367194e-01 5.316403e-01
[406,] 0.7841718 4.316563e-01 2.158282e-01
[407,] 0.6533324 6.933352e-01 3.466676e-01
[408,] 0.4907134 9.814268e-01 5.092866e-01
> postscript(file="/var/www/html/freestat/rcomp/tmp/1frax1291330094.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> plot(x[,1], type='l', main='Actuals and Interpolation', ylab='value of Actuals and Interpolation (dots)', xlab='time or index')
> points(x[,1]-mysum$resid)
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/freestat/rcomp/tmp/2frax1291330094.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> plot(mysum$resid, type='b', pch=19, main='Residuals', ylab='value of Residuals', xlab='time or index')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/freestat/rcomp/tmp/3frax1291330094.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> hist(mysum$resid, main='Residual Histogram', xlab='values of Residuals')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/freestat/rcomp/tmp/4q0si1291330094.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> densityplot(~mysum$resid,col='black',main='Residual Density Plot', xlab='values of Residuals')
> dev.off()
null device
1
> postscript(file="/var/www/html/freestat/rcomp/tmp/5q0si1291330094.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
-335057.38899 3055410.06433 659410.54377 -405429.22905 -743164.38986
6 7 8 9 10
-89456.37662 957456.07911 358218.87563 -515485.62638 656837.66595
11 12 13 14 15
418629.65698 790336.11651 786396.36249 -496204.59220 -198517.34823
16 17 18 19 20
573684.82775 -252041.74304 -311194.45538 -428313.88915 -21431.49731
21 22 23 24 25
139295.64762 -335211.15220 114533.03501 -174367.24805 -66395.70200
26 27 28 29 30
301730.34172 -543727.48346 -174098.17192 593025.71736 341865.12044
31 32 33 34 35
252881.74560 -132561.70348 163074.51332 105181.03944 -175459.84008
36 37 38 39 40
26963.16687 306078.95930 161359.42712 138591.64288 -121232.08493
41 42 43 44 45
-315987.14980 21458.74625 104529.24113 -48445.33886 62799.44298
46 47 48 49 50
-127194.57210 118925.53513 -219278.61704 -14451.06244 -142981.86240
51 52 53 54 55
105768.72434 -34627.40188 -219322.05982 70651.42675 -281080.27993
56 57 58 59 60
124875.96679 232448.43092 85407.63172 -63984.83679 -220578.16036
61 62 63 64 65
-23272.87970 -78225.09783 -20749.17848 -216463.75382 68796.51174
66 67 68 69 70
23572.03863 -137137.24569 -268.85031 -84555.82611 -55112.87409
71 72 73 74 75
-89236.14587 153263.11954 -121884.06595 -46107.67637 -94572.99222
76 77 78 79 80
-73086.17103 -105602.61775 -85644.93657 37088.85778 21253.31059
81 82 83 84 85
29611.70914 -45462.15309 -134129.65691 36887.62176 -23125.20612
86 87 88 89 90
-116597.88333 70596.47129 -54748.11133 -88191.48526 -33468.48369
91 92 93 94 95
-20293.38577 -41947.40635 122154.90850 -228180.60334 119037.29664
96 97 98 99 100
-11477.33974 -1554.05235 49267.09848 -71063.71119 -26793.91038
101 102 103 104 105
-34155.89201 -116393.68974 -71182.11744 -63711.14507 39309.00263
106 107 108 109 110
-78901.81531 -233869.45225 107762.43360 -87569.45515 -104668.69729
111 112 113 114 115
82970.34063 -253759.38582 -175074.28537 -121209.46808 14252.74836
116 117 118 119 120
55220.38255 -26722.55853 65798.28707 -176082.65391 26570.83341
121 122 123 124 125
2477.23170 13431.45499 -36512.52694 80946.69568 -100397.40476
126 127 128 129 130
-201133.25705 -2401.22685 -987.01213 9809.10476 -2401.22685
131 132 133 134 135
-386.54690 85481.18719 -2401.22685 -3753.75972 -2401.22685
136 137 138 139 140
7599.65197 175.79713 -2401.22685 -2401.22685 -2401.22685
141 142 143 144 145
3951.31681 -2401.22685 17400.26177 -1675.86504 16223.04061
146 147 148 149 150
13349.51735 -2401.22685 -1714.53127 -2401.22685 -2564.61669
151 152 153 154 155
543.65038 22998.67080 -2401.22685 -14612.79492 -2327.64702
156 157 158 159 160
-2401.22685 -2401.22685 -2399.01202 15826.16411 -2401.22685
161 162 163 164 165
-2401.22685 -2401.22685 -2401.22685 11382.48822 -60.46908
166 167 168 169 170
-2401.22685 -2401.22685 110015.49676 -2401.22685 -2401.22685
171 172 173 174 175
-2401.22685 -4700.09393 -2401.22685 239.82022 -2401.22685
176 177 178 179 180
-1522.71355 -2401.22685 24853.40506 -2401.22685 -2401.22685
181 182 183 184 185
39696.20723 23643.86380 29378.05747 -12686.20265 29249.09785
186 187 188 189 190
11165.20288 33626.14463 -1350.71885 -32456.24329 -2401.22685
191 192 193 194 195
8980.09568 -2401.22685 3290.42860 -2401.22685 -2200.22282
196 197 198 199 200
26723.56888 -41404.00191 26723.56888 26723.56888 -13875.78642
201 202 203 204 205
26723.56888 26723.56888 26723.56888 25585.75864 -4271.33958
206 207 208 209 210
26784.69571 22287.18057 26723.56888 44234.67155 26723.56888
211 212 213 214 215
-2401.22685 -2401.22685 11362.80316 -14186.21899 -4637.33104
216 217 218 219 220
-2401.22685 -486.64678 6769.86053 -2401.22685 27804.36286
221 222 223 224 225
-27290.37324 -2401.22685 -1963.41995 -2401.22685 -1935.75325
226 227 228 229 230
4557.43857 -2401.22685 -2401.22685 -20062.40201 -2401.22685
231 232 233 234 235
-13551.83006 -2401.22685 -2401.22685 6727.90438 -2401.22685
236 237 238 239 240
-2333.83700 -2401.22685 339.50561 -2401.22685 -2401.22685
241 242 243 244 245
-2401.22685 3152.22517 -2401.22685 720.53012 -7170.17477
246 247 248 249 250
26723.56888 26723.56888 -8564.63380 26723.56888 26723.56888
251 252 253 254 255
3156.75624 26723.56888 -18278.38317 -17986.72143 -42817.86848
256 257 258 259 260
-2401.22685 26006.40867 32093.64563 26723.56888 -2401.22685
261 262 263 264 265
26723.56888 26723.56888 -9433.34743 10546.87598 26723.56888
266 267 268 269 270
-2401.22685 26723.56888 16171.00545 26723.56888 -14299.80626
271 272 273 274 275
26723.56888 26723.56888 26723.56888 -3108.46350 -2401.22685
276 277 278 279 280
8612.68765 8347.42108 -2401.22685 6172.76214 -20844.16653
281 282 283 284 285
26723.56888 -2401.22685 23272.10884 -2401.22685 -24024.75157
286 287 288 289 290
-2401.22685 26723.56888 -2486.74385 26723.56888 16783.25672
291 292 293 294 295
26723.56888 7464.49902 -6947.57562 -2401.22685 -1737.31668
296 297 298 299 300
51474.88847 -2401.22685 26723.56888 -2401.22685 -2401.22685
301 302 303 304 305
25184.63237 -9291.67440 -2401.22685 26723.56888 34355.02811
306 307 308 309 310
-666.00381 26723.56888 22464.83173 21670.74259 26723.56888
311 312 313 314 315
26723.56888 -26255.79903 26723.56888 29891.95209 -1414.38736
316 317 318 319 320
2902.88817 -643.29997 -9108.09383 -1238.86183 -2401.22685
321 322 323 324 325
-2401.22685 -2401.22685 -2401.22685 26723.56888 -2401.22685
326 327 328 329 330
-4812.37279 26723.56888 26723.56888 -2401.22685 -2401.22685
331 332 333 334 335
-2401.22685 -2401.22685 -2401.22685 -944.13962 -2022.08210
336 337 338 339 340
-2401.22685 -2729.53365 -2401.22685 -2356.13091 1358.77905
341 342 343 344 345
-14305.02958 -2401.22685 -2401.22685 -2401.22685 -2401.22685
346 347 348 349 350
-2401.22685 -9488.63760 -2401.22685 -2401.22685 -2401.22685
351 352 353 354 355
2939.06427 -1578.88762 13632.90069 -7486.21423 -2401.22685
356 357 358 359 360
-2401.22685 -24740.26924 -2401.22685 -2401.22685 -2401.22685
361 362 363 364 365
-1630.79340 -2401.22685 -2401.22685 18714.84729 -1899.99534
366 367 368 369 370
-2401.22685 -2401.22685 -2401.22685 -2401.22685 -2401.22685
371 372 373 374 375
-2401.22685 3486.88107 -25749.55270 -20702.21744 -2401.22685
376 377 378 379 380
-2401.22685 4330.91660 -1503.97547 -2401.22685 -2401.22685
381 382 383 384 385
-2401.22685 -3409.33390 41667.47187 -2401.22685 -763.76613
386 387 388 389 390
-86823.08487 3994.32689 7260.16022 -2997.41084 -9409.52872
391 392 393 394 395
5753.07660 -11857.94258 -22132.20075 39537.25471 -155995.44039
396 397 398 399 400
-17357.95088 -25973.91882 107289.31148 -12279.43491 -84213.96533
401 402 403 404 405
-496.61868 12242.61718 -7008.00459 -51150.14027 -16445.04507
406 407 408 409 410
-13369.73486 2334.62395 4136.64604 -15614.21169 -48274.99259
411 412 413 414 415
-18647.91267 -219572.24477 -22111.50353 -110588.42235 -58146.69165
416 417 418 419 420
-57201.32722 155507.03181 -34356.48016 -14820.25975 -46731.17993
421 422 423 424 425
-85486.26190 -82041.92819 -627186.38066 -331586.22377 -68145.58823
426 427 428 429 430
-197325.48460 -251752.72344 -251741.39555 -127326.89068 135682.64604
431
-363758.36870
> postscript(file="/var/www/html/freestat/rcomp/tmp/6q0si1291330094.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 -335057.38899 NA
1 3055410.06433 -335057.38899
2 659410.54377 3055410.06433
3 -405429.22905 659410.54377
4 -743164.38986 -405429.22905
5 -89456.37662 -743164.38986
6 957456.07911 -89456.37662
7 358218.87563 957456.07911
8 -515485.62638 358218.87563
9 656837.66595 -515485.62638
10 418629.65698 656837.66595
11 790336.11651 418629.65698
12 786396.36249 790336.11651
13 -496204.59220 786396.36249
14 -198517.34823 -496204.59220
15 573684.82775 -198517.34823
16 -252041.74304 573684.82775
17 -311194.45538 -252041.74304
18 -428313.88915 -311194.45538
19 -21431.49731 -428313.88915
20 139295.64762 -21431.49731
21 -335211.15220 139295.64762
22 114533.03501 -335211.15220
23 -174367.24805 114533.03501
24 -66395.70200 -174367.24805
25 301730.34172 -66395.70200
26 -543727.48346 301730.34172
27 -174098.17192 -543727.48346
28 593025.71736 -174098.17192
29 341865.12044 593025.71736
30 252881.74560 341865.12044
31 -132561.70348 252881.74560
32 163074.51332 -132561.70348
33 105181.03944 163074.51332
34 -175459.84008 105181.03944
35 26963.16687 -175459.84008
36 306078.95930 26963.16687
37 161359.42712 306078.95930
38 138591.64288 161359.42712
39 -121232.08493 138591.64288
40 -315987.14980 -121232.08493
41 21458.74625 -315987.14980
42 104529.24113 21458.74625
43 -48445.33886 104529.24113
44 62799.44298 -48445.33886
45 -127194.57210 62799.44298
46 118925.53513 -127194.57210
47 -219278.61704 118925.53513
48 -14451.06244 -219278.61704
49 -142981.86240 -14451.06244
50 105768.72434 -142981.86240
51 -34627.40188 105768.72434
52 -219322.05982 -34627.40188
53 70651.42675 -219322.05982
54 -281080.27993 70651.42675
55 124875.96679 -281080.27993
56 232448.43092 124875.96679
57 85407.63172 232448.43092
58 -63984.83679 85407.63172
59 -220578.16036 -63984.83679
60 -23272.87970 -220578.16036
61 -78225.09783 -23272.87970
62 -20749.17848 -78225.09783
63 -216463.75382 -20749.17848
64 68796.51174 -216463.75382
65 23572.03863 68796.51174
66 -137137.24569 23572.03863
67 -268.85031 -137137.24569
68 -84555.82611 -268.85031
69 -55112.87409 -84555.82611
70 -89236.14587 -55112.87409
71 153263.11954 -89236.14587
72 -121884.06595 153263.11954
73 -46107.67637 -121884.06595
74 -94572.99222 -46107.67637
75 -73086.17103 -94572.99222
76 -105602.61775 -73086.17103
77 -85644.93657 -105602.61775
78 37088.85778 -85644.93657
79 21253.31059 37088.85778
80 29611.70914 21253.31059
81 -45462.15309 29611.70914
82 -134129.65691 -45462.15309
83 36887.62176 -134129.65691
84 -23125.20612 36887.62176
85 -116597.88333 -23125.20612
86 70596.47129 -116597.88333
87 -54748.11133 70596.47129
88 -88191.48526 -54748.11133
89 -33468.48369 -88191.48526
90 -20293.38577 -33468.48369
91 -41947.40635 -20293.38577
92 122154.90850 -41947.40635
93 -228180.60334 122154.90850
94 119037.29664 -228180.60334
95 -11477.33974 119037.29664
96 -1554.05235 -11477.33974
97 49267.09848 -1554.05235
98 -71063.71119 49267.09848
99 -26793.91038 -71063.71119
100 -34155.89201 -26793.91038
101 -116393.68974 -34155.89201
102 -71182.11744 -116393.68974
103 -63711.14507 -71182.11744
104 39309.00263 -63711.14507
105 -78901.81531 39309.00263
106 -233869.45225 -78901.81531
107 107762.43360 -233869.45225
108 -87569.45515 107762.43360
109 -104668.69729 -87569.45515
110 82970.34063 -104668.69729
111 -253759.38582 82970.34063
112 -175074.28537 -253759.38582
113 -121209.46808 -175074.28537
114 14252.74836 -121209.46808
115 55220.38255 14252.74836
116 -26722.55853 55220.38255
117 65798.28707 -26722.55853
118 -176082.65391 65798.28707
119 26570.83341 -176082.65391
120 2477.23170 26570.83341
121 13431.45499 2477.23170
122 -36512.52694 13431.45499
123 80946.69568 -36512.52694
124 -100397.40476 80946.69568
125 -201133.25705 -100397.40476
126 -2401.22685 -201133.25705
127 -987.01213 -2401.22685
128 9809.10476 -987.01213
129 -2401.22685 9809.10476
130 -386.54690 -2401.22685
131 85481.18719 -386.54690
132 -2401.22685 85481.18719
133 -3753.75972 -2401.22685
134 -2401.22685 -3753.75972
135 7599.65197 -2401.22685
136 175.79713 7599.65197
137 -2401.22685 175.79713
138 -2401.22685 -2401.22685
139 -2401.22685 -2401.22685
140 3951.31681 -2401.22685
141 -2401.22685 3951.31681
142 17400.26177 -2401.22685
143 -1675.86504 17400.26177
144 16223.04061 -1675.86504
145 13349.51735 16223.04061
146 -2401.22685 13349.51735
147 -1714.53127 -2401.22685
148 -2401.22685 -1714.53127
149 -2564.61669 -2401.22685
150 543.65038 -2564.61669
151 22998.67080 543.65038
152 -2401.22685 22998.67080
153 -14612.79492 -2401.22685
154 -2327.64702 -14612.79492
155 -2401.22685 -2327.64702
156 -2401.22685 -2401.22685
157 -2399.01202 -2401.22685
158 15826.16411 -2399.01202
159 -2401.22685 15826.16411
160 -2401.22685 -2401.22685
161 -2401.22685 -2401.22685
162 -2401.22685 -2401.22685
163 11382.48822 -2401.22685
164 -60.46908 11382.48822
165 -2401.22685 -60.46908
166 -2401.22685 -2401.22685
167 110015.49676 -2401.22685
168 -2401.22685 110015.49676
169 -2401.22685 -2401.22685
170 -2401.22685 -2401.22685
171 -4700.09393 -2401.22685
172 -2401.22685 -4700.09393
173 239.82022 -2401.22685
174 -2401.22685 239.82022
175 -1522.71355 -2401.22685
176 -2401.22685 -1522.71355
177 24853.40506 -2401.22685
178 -2401.22685 24853.40506
179 -2401.22685 -2401.22685
180 39696.20723 -2401.22685
181 23643.86380 39696.20723
182 29378.05747 23643.86380
183 -12686.20265 29378.05747
184 29249.09785 -12686.20265
185 11165.20288 29249.09785
186 33626.14463 11165.20288
187 -1350.71885 33626.14463
188 -32456.24329 -1350.71885
189 -2401.22685 -32456.24329
190 8980.09568 -2401.22685
191 -2401.22685 8980.09568
192 3290.42860 -2401.22685
193 -2401.22685 3290.42860
194 -2200.22282 -2401.22685
195 26723.56888 -2200.22282
196 -41404.00191 26723.56888
197 26723.56888 -41404.00191
198 26723.56888 26723.56888
199 -13875.78642 26723.56888
200 26723.56888 -13875.78642
201 26723.56888 26723.56888
202 26723.56888 26723.56888
203 25585.75864 26723.56888
204 -4271.33958 25585.75864
205 26784.69571 -4271.33958
206 22287.18057 26784.69571
207 26723.56888 22287.18057
208 44234.67155 26723.56888
209 26723.56888 44234.67155
210 -2401.22685 26723.56888
211 -2401.22685 -2401.22685
212 11362.80316 -2401.22685
213 -14186.21899 11362.80316
214 -4637.33104 -14186.21899
215 -2401.22685 -4637.33104
216 -486.64678 -2401.22685
217 6769.86053 -486.64678
218 -2401.22685 6769.86053
219 27804.36286 -2401.22685
220 -27290.37324 27804.36286
221 -2401.22685 -27290.37324
222 -1963.41995 -2401.22685
223 -2401.22685 -1963.41995
224 -1935.75325 -2401.22685
225 4557.43857 -1935.75325
226 -2401.22685 4557.43857
227 -2401.22685 -2401.22685
228 -20062.40201 -2401.22685
229 -2401.22685 -20062.40201
230 -13551.83006 -2401.22685
231 -2401.22685 -13551.83006
232 -2401.22685 -2401.22685
233 6727.90438 -2401.22685
234 -2401.22685 6727.90438
235 -2333.83700 -2401.22685
236 -2401.22685 -2333.83700
237 339.50561 -2401.22685
238 -2401.22685 339.50561
239 -2401.22685 -2401.22685
240 -2401.22685 -2401.22685
241 3152.22517 -2401.22685
242 -2401.22685 3152.22517
243 720.53012 -2401.22685
244 -7170.17477 720.53012
245 26723.56888 -7170.17477
246 26723.56888 26723.56888
247 -8564.63380 26723.56888
248 26723.56888 -8564.63380
249 26723.56888 26723.56888
250 3156.75624 26723.56888
251 26723.56888 3156.75624
252 -18278.38317 26723.56888
253 -17986.72143 -18278.38317
254 -42817.86848 -17986.72143
255 -2401.22685 -42817.86848
256 26006.40867 -2401.22685
257 32093.64563 26006.40867
258 26723.56888 32093.64563
259 -2401.22685 26723.56888
260 26723.56888 -2401.22685
261 26723.56888 26723.56888
262 -9433.34743 26723.56888
263 10546.87598 -9433.34743
264 26723.56888 10546.87598
265 -2401.22685 26723.56888
266 26723.56888 -2401.22685
267 16171.00545 26723.56888
268 26723.56888 16171.00545
269 -14299.80626 26723.56888
270 26723.56888 -14299.80626
271 26723.56888 26723.56888
272 26723.56888 26723.56888
273 -3108.46350 26723.56888
274 -2401.22685 -3108.46350
275 8612.68765 -2401.22685
276 8347.42108 8612.68765
277 -2401.22685 8347.42108
278 6172.76214 -2401.22685
279 -20844.16653 6172.76214
280 26723.56888 -20844.16653
281 -2401.22685 26723.56888
282 23272.10884 -2401.22685
283 -2401.22685 23272.10884
284 -24024.75157 -2401.22685
285 -2401.22685 -24024.75157
286 26723.56888 -2401.22685
287 -2486.74385 26723.56888
288 26723.56888 -2486.74385
289 16783.25672 26723.56888
290 26723.56888 16783.25672
291 7464.49902 26723.56888
292 -6947.57562 7464.49902
293 -2401.22685 -6947.57562
294 -1737.31668 -2401.22685
295 51474.88847 -1737.31668
296 -2401.22685 51474.88847
297 26723.56888 -2401.22685
298 -2401.22685 26723.56888
299 -2401.22685 -2401.22685
300 25184.63237 -2401.22685
301 -9291.67440 25184.63237
302 -2401.22685 -9291.67440
303 26723.56888 -2401.22685
304 34355.02811 26723.56888
305 -666.00381 34355.02811
306 26723.56888 -666.00381
307 22464.83173 26723.56888
308 21670.74259 22464.83173
309 26723.56888 21670.74259
310 26723.56888 26723.56888
311 -26255.79903 26723.56888
312 26723.56888 -26255.79903
313 29891.95209 26723.56888
314 -1414.38736 29891.95209
315 2902.88817 -1414.38736
316 -643.29997 2902.88817
317 -9108.09383 -643.29997
318 -1238.86183 -9108.09383
319 -2401.22685 -1238.86183
320 -2401.22685 -2401.22685
321 -2401.22685 -2401.22685
322 -2401.22685 -2401.22685
323 26723.56888 -2401.22685
324 -2401.22685 26723.56888
325 -4812.37279 -2401.22685
326 26723.56888 -4812.37279
327 26723.56888 26723.56888
328 -2401.22685 26723.56888
329 -2401.22685 -2401.22685
330 -2401.22685 -2401.22685
331 -2401.22685 -2401.22685
332 -2401.22685 -2401.22685
333 -944.13962 -2401.22685
334 -2022.08210 -944.13962
335 -2401.22685 -2022.08210
336 -2729.53365 -2401.22685
337 -2401.22685 -2729.53365
338 -2356.13091 -2401.22685
339 1358.77905 -2356.13091
340 -14305.02958 1358.77905
341 -2401.22685 -14305.02958
342 -2401.22685 -2401.22685
343 -2401.22685 -2401.22685
344 -2401.22685 -2401.22685
345 -2401.22685 -2401.22685
346 -9488.63760 -2401.22685
347 -2401.22685 -9488.63760
348 -2401.22685 -2401.22685
349 -2401.22685 -2401.22685
350 2939.06427 -2401.22685
351 -1578.88762 2939.06427
352 13632.90069 -1578.88762
353 -7486.21423 13632.90069
354 -2401.22685 -7486.21423
355 -2401.22685 -2401.22685
356 -24740.26924 -2401.22685
357 -2401.22685 -24740.26924
358 -2401.22685 -2401.22685
359 -2401.22685 -2401.22685
360 -1630.79340 -2401.22685
361 -2401.22685 -1630.79340
362 -2401.22685 -2401.22685
363 18714.84729 -2401.22685
364 -1899.99534 18714.84729
365 -2401.22685 -1899.99534
366 -2401.22685 -2401.22685
367 -2401.22685 -2401.22685
368 -2401.22685 -2401.22685
369 -2401.22685 -2401.22685
370 -2401.22685 -2401.22685
371 3486.88107 -2401.22685
372 -25749.55270 3486.88107
373 -20702.21744 -25749.55270
374 -2401.22685 -20702.21744
375 -2401.22685 -2401.22685
376 4330.91660 -2401.22685
377 -1503.97547 4330.91660
378 -2401.22685 -1503.97547
379 -2401.22685 -2401.22685
380 -2401.22685 -2401.22685
381 -3409.33390 -2401.22685
382 41667.47187 -3409.33390
383 -2401.22685 41667.47187
384 -763.76613 -2401.22685
385 -86823.08487 -763.76613
386 3994.32689 -86823.08487
387 7260.16022 3994.32689
388 -2997.41084 7260.16022
389 -9409.52872 -2997.41084
390 5753.07660 -9409.52872
391 -11857.94258 5753.07660
392 -22132.20075 -11857.94258
393 39537.25471 -22132.20075
394 -155995.44039 39537.25471
395 -17357.95088 -155995.44039
396 -25973.91882 -17357.95088
397 107289.31148 -25973.91882
398 -12279.43491 107289.31148
399 -84213.96533 -12279.43491
400 -496.61868 -84213.96533
401 12242.61718 -496.61868
402 -7008.00459 12242.61718
403 -51150.14027 -7008.00459
404 -16445.04507 -51150.14027
405 -13369.73486 -16445.04507
406 2334.62395 -13369.73486
407 4136.64604 2334.62395
408 -15614.21169 4136.64604
409 -48274.99259 -15614.21169
410 -18647.91267 -48274.99259
411 -219572.24477 -18647.91267
412 -22111.50353 -219572.24477
413 -110588.42235 -22111.50353
414 -58146.69165 -110588.42235
415 -57201.32722 -58146.69165
416 155507.03181 -57201.32722
417 -34356.48016 155507.03181
418 -14820.25975 -34356.48016
419 -46731.17993 -14820.25975
420 -85486.26190 -46731.17993
421 -82041.92819 -85486.26190
422 -627186.38066 -82041.92819
423 -331586.22377 -627186.38066
424 -68145.58823 -331586.22377
425 -197325.48460 -68145.58823
426 -251752.72344 -197325.48460
427 -251741.39555 -251752.72344
428 -127326.89068 -251741.39555
429 135682.64604 -127326.89068
430 -363758.36870 135682.64604
431 NA -363758.36870
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] 3055410.06433 -335057.38899
[2,] 659410.54377 3055410.06433
[3,] -405429.22905 659410.54377
[4,] -743164.38986 -405429.22905
[5,] -89456.37662 -743164.38986
[6,] 957456.07911 -89456.37662
[7,] 358218.87563 957456.07911
[8,] -515485.62638 358218.87563
[9,] 656837.66595 -515485.62638
[10,] 418629.65698 656837.66595
[11,] 790336.11651 418629.65698
[12,] 786396.36249 790336.11651
[13,] -496204.59220 786396.36249
[14,] -198517.34823 -496204.59220
[15,] 573684.82775 -198517.34823
[16,] -252041.74304 573684.82775
[17,] -311194.45538 -252041.74304
[18,] -428313.88915 -311194.45538
[19,] -21431.49731 -428313.88915
[20,] 139295.64762 -21431.49731
[21,] -335211.15220 139295.64762
[22,] 114533.03501 -335211.15220
[23,] -174367.24805 114533.03501
[24,] -66395.70200 -174367.24805
[25,] 301730.34172 -66395.70200
[26,] -543727.48346 301730.34172
[27,] -174098.17192 -543727.48346
[28,] 593025.71736 -174098.17192
[29,] 341865.12044 593025.71736
[30,] 252881.74560 341865.12044
[31,] -132561.70348 252881.74560
[32,] 163074.51332 -132561.70348
[33,] 105181.03944 163074.51332
[34,] -175459.84008 105181.03944
[35,] 26963.16687 -175459.84008
[36,] 306078.95930 26963.16687
[37,] 161359.42712 306078.95930
[38,] 138591.64288 161359.42712
[39,] -121232.08493 138591.64288
[40,] -315987.14980 -121232.08493
[41,] 21458.74625 -315987.14980
[42,] 104529.24113 21458.74625
[43,] -48445.33886 104529.24113
[44,] 62799.44298 -48445.33886
[45,] -127194.57210 62799.44298
[46,] 118925.53513 -127194.57210
[47,] -219278.61704 118925.53513
[48,] -14451.06244 -219278.61704
[49,] -142981.86240 -14451.06244
[50,] 105768.72434 -142981.86240
[51,] -34627.40188 105768.72434
[52,] -219322.05982 -34627.40188
[53,] 70651.42675 -219322.05982
[54,] -281080.27993 70651.42675
[55,] 124875.96679 -281080.27993
[56,] 232448.43092 124875.96679
[57,] 85407.63172 232448.43092
[58,] -63984.83679 85407.63172
[59,] -220578.16036 -63984.83679
[60,] -23272.87970 -220578.16036
[61,] -78225.09783 -23272.87970
[62,] -20749.17848 -78225.09783
[63,] -216463.75382 -20749.17848
[64,] 68796.51174 -216463.75382
[65,] 23572.03863 68796.51174
[66,] -137137.24569 23572.03863
[67,] -268.85031 -137137.24569
[68,] -84555.82611 -268.85031
[69,] -55112.87409 -84555.82611
[70,] -89236.14587 -55112.87409
[71,] 153263.11954 -89236.14587
[72,] -121884.06595 153263.11954
[73,] -46107.67637 -121884.06595
[74,] -94572.99222 -46107.67637
[75,] -73086.17103 -94572.99222
[76,] -105602.61775 -73086.17103
[77,] -85644.93657 -105602.61775
[78,] 37088.85778 -85644.93657
[79,] 21253.31059 37088.85778
[80,] 29611.70914 21253.31059
[81,] -45462.15309 29611.70914
[82,] -134129.65691 -45462.15309
[83,] 36887.62176 -134129.65691
[84,] -23125.20612 36887.62176
[85,] -116597.88333 -23125.20612
[86,] 70596.47129 -116597.88333
[87,] -54748.11133 70596.47129
[88,] -88191.48526 -54748.11133
[89,] -33468.48369 -88191.48526
[90,] -20293.38577 -33468.48369
[91,] -41947.40635 -20293.38577
[92,] 122154.90850 -41947.40635
[93,] -228180.60334 122154.90850
[94,] 119037.29664 -228180.60334
[95,] -11477.33974 119037.29664
[96,] -1554.05235 -11477.33974
[97,] 49267.09848 -1554.05235
[98,] -71063.71119 49267.09848
[99,] -26793.91038 -71063.71119
[100,] -34155.89201 -26793.91038
[101,] -116393.68974 -34155.89201
[102,] -71182.11744 -116393.68974
[103,] -63711.14507 -71182.11744
[104,] 39309.00263 -63711.14507
[105,] -78901.81531 39309.00263
[106,] -233869.45225 -78901.81531
[107,] 107762.43360 -233869.45225
[108,] -87569.45515 107762.43360
[109,] -104668.69729 -87569.45515
[110,] 82970.34063 -104668.69729
[111,] -253759.38582 82970.34063
[112,] -175074.28537 -253759.38582
[113,] -121209.46808 -175074.28537
[114,] 14252.74836 -121209.46808
[115,] 55220.38255 14252.74836
[116,] -26722.55853 55220.38255
[117,] 65798.28707 -26722.55853
[118,] -176082.65391 65798.28707
[119,] 26570.83341 -176082.65391
[120,] 2477.23170 26570.83341
[121,] 13431.45499 2477.23170
[122,] -36512.52694 13431.45499
[123,] 80946.69568 -36512.52694
[124,] -100397.40476 80946.69568
[125,] -201133.25705 -100397.40476
[126,] -2401.22685 -201133.25705
[127,] -987.01213 -2401.22685
[128,] 9809.10476 -987.01213
[129,] -2401.22685 9809.10476
[130,] -386.54690 -2401.22685
[131,] 85481.18719 -386.54690
[132,] -2401.22685 85481.18719
[133,] -3753.75972 -2401.22685
[134,] -2401.22685 -3753.75972
[135,] 7599.65197 -2401.22685
[136,] 175.79713 7599.65197
[137,] -2401.22685 175.79713
[138,] -2401.22685 -2401.22685
[139,] -2401.22685 -2401.22685
[140,] 3951.31681 -2401.22685
[141,] -2401.22685 3951.31681
[142,] 17400.26177 -2401.22685
[143,] -1675.86504 17400.26177
[144,] 16223.04061 -1675.86504
[145,] 13349.51735 16223.04061
[146,] -2401.22685 13349.51735
[147,] -1714.53127 -2401.22685
[148,] -2401.22685 -1714.53127
[149,] -2564.61669 -2401.22685
[150,] 543.65038 -2564.61669
[151,] 22998.67080 543.65038
[152,] -2401.22685 22998.67080
[153,] -14612.79492 -2401.22685
[154,] -2327.64702 -14612.79492
[155,] -2401.22685 -2327.64702
[156,] -2401.22685 -2401.22685
[157,] -2399.01202 -2401.22685
[158,] 15826.16411 -2399.01202
[159,] -2401.22685 15826.16411
[160,] -2401.22685 -2401.22685
[161,] -2401.22685 -2401.22685
[162,] -2401.22685 -2401.22685
[163,] 11382.48822 -2401.22685
[164,] -60.46908 11382.48822
[165,] -2401.22685 -60.46908
[166,] -2401.22685 -2401.22685
[167,] 110015.49676 -2401.22685
[168,] -2401.22685 110015.49676
[169,] -2401.22685 -2401.22685
[170,] -2401.22685 -2401.22685
[171,] -4700.09393 -2401.22685
[172,] -2401.22685 -4700.09393
[173,] 239.82022 -2401.22685
[174,] -2401.22685 239.82022
[175,] -1522.71355 -2401.22685
[176,] -2401.22685 -1522.71355
[177,] 24853.40506 -2401.22685
[178,] -2401.22685 24853.40506
[179,] -2401.22685 -2401.22685
[180,] 39696.20723 -2401.22685
[181,] 23643.86380 39696.20723
[182,] 29378.05747 23643.86380
[183,] -12686.20265 29378.05747
[184,] 29249.09785 -12686.20265
[185,] 11165.20288 29249.09785
[186,] 33626.14463 11165.20288
[187,] -1350.71885 33626.14463
[188,] -32456.24329 -1350.71885
[189,] -2401.22685 -32456.24329
[190,] 8980.09568 -2401.22685
[191,] -2401.22685 8980.09568
[192,] 3290.42860 -2401.22685
[193,] -2401.22685 3290.42860
[194,] -2200.22282 -2401.22685
[195,] 26723.56888 -2200.22282
[196,] -41404.00191 26723.56888
[197,] 26723.56888 -41404.00191
[198,] 26723.56888 26723.56888
[199,] -13875.78642 26723.56888
[200,] 26723.56888 -13875.78642
[201,] 26723.56888 26723.56888
[202,] 26723.56888 26723.56888
[203,] 25585.75864 26723.56888
[204,] -4271.33958 25585.75864
[205,] 26784.69571 -4271.33958
[206,] 22287.18057 26784.69571
[207,] 26723.56888 22287.18057
[208,] 44234.67155 26723.56888
[209,] 26723.56888 44234.67155
[210,] -2401.22685 26723.56888
[211,] -2401.22685 -2401.22685
[212,] 11362.80316 -2401.22685
[213,] -14186.21899 11362.80316
[214,] -4637.33104 -14186.21899
[215,] -2401.22685 -4637.33104
[216,] -486.64678 -2401.22685
[217,] 6769.86053 -486.64678
[218,] -2401.22685 6769.86053
[219,] 27804.36286 -2401.22685
[220,] -27290.37324 27804.36286
[221,] -2401.22685 -27290.37324
[222,] -1963.41995 -2401.22685
[223,] -2401.22685 -1963.41995
[224,] -1935.75325 -2401.22685
[225,] 4557.43857 -1935.75325
[226,] -2401.22685 4557.43857
[227,] -2401.22685 -2401.22685
[228,] -20062.40201 -2401.22685
[229,] -2401.22685 -20062.40201
[230,] -13551.83006 -2401.22685
[231,] -2401.22685 -13551.83006
[232,] -2401.22685 -2401.22685
[233,] 6727.90438 -2401.22685
[234,] -2401.22685 6727.90438
[235,] -2333.83700 -2401.22685
[236,] -2401.22685 -2333.83700
[237,] 339.50561 -2401.22685
[238,] -2401.22685 339.50561
[239,] -2401.22685 -2401.22685
[240,] -2401.22685 -2401.22685
[241,] 3152.22517 -2401.22685
[242,] -2401.22685 3152.22517
[243,] 720.53012 -2401.22685
[244,] -7170.17477 720.53012
[245,] 26723.56888 -7170.17477
[246,] 26723.56888 26723.56888
[247,] -8564.63380 26723.56888
[248,] 26723.56888 -8564.63380
[249,] 26723.56888 26723.56888
[250,] 3156.75624 26723.56888
[251,] 26723.56888 3156.75624
[252,] -18278.38317 26723.56888
[253,] -17986.72143 -18278.38317
[254,] -42817.86848 -17986.72143
[255,] -2401.22685 -42817.86848
[256,] 26006.40867 -2401.22685
[257,] 32093.64563 26006.40867
[258,] 26723.56888 32093.64563
[259,] -2401.22685 26723.56888
[260,] 26723.56888 -2401.22685
[261,] 26723.56888 26723.56888
[262,] -9433.34743 26723.56888
[263,] 10546.87598 -9433.34743
[264,] 26723.56888 10546.87598
[265,] -2401.22685 26723.56888
[266,] 26723.56888 -2401.22685
[267,] 16171.00545 26723.56888
[268,] 26723.56888 16171.00545
[269,] -14299.80626 26723.56888
[270,] 26723.56888 -14299.80626
[271,] 26723.56888 26723.56888
[272,] 26723.56888 26723.56888
[273,] -3108.46350 26723.56888
[274,] -2401.22685 -3108.46350
[275,] 8612.68765 -2401.22685
[276,] 8347.42108 8612.68765
[277,] -2401.22685 8347.42108
[278,] 6172.76214 -2401.22685
[279,] -20844.16653 6172.76214
[280,] 26723.56888 -20844.16653
[281,] -2401.22685 26723.56888
[282,] 23272.10884 -2401.22685
[283,] -2401.22685 23272.10884
[284,] -24024.75157 -2401.22685
[285,] -2401.22685 -24024.75157
[286,] 26723.56888 -2401.22685
[287,] -2486.74385 26723.56888
[288,] 26723.56888 -2486.74385
[289,] 16783.25672 26723.56888
[290,] 26723.56888 16783.25672
[291,] 7464.49902 26723.56888
[292,] -6947.57562 7464.49902
[293,] -2401.22685 -6947.57562
[294,] -1737.31668 -2401.22685
[295,] 51474.88847 -1737.31668
[296,] -2401.22685 51474.88847
[297,] 26723.56888 -2401.22685
[298,] -2401.22685 26723.56888
[299,] -2401.22685 -2401.22685
[300,] 25184.63237 -2401.22685
[301,] -9291.67440 25184.63237
[302,] -2401.22685 -9291.67440
[303,] 26723.56888 -2401.22685
[304,] 34355.02811 26723.56888
[305,] -666.00381 34355.02811
[306,] 26723.56888 -666.00381
[307,] 22464.83173 26723.56888
[308,] 21670.74259 22464.83173
[309,] 26723.56888 21670.74259
[310,] 26723.56888 26723.56888
[311,] -26255.79903 26723.56888
[312,] 26723.56888 -26255.79903
[313,] 29891.95209 26723.56888
[314,] -1414.38736 29891.95209
[315,] 2902.88817 -1414.38736
[316,] -643.29997 2902.88817
[317,] -9108.09383 -643.29997
[318,] -1238.86183 -9108.09383
[319,] -2401.22685 -1238.86183
[320,] -2401.22685 -2401.22685
[321,] -2401.22685 -2401.22685
[322,] -2401.22685 -2401.22685
[323,] 26723.56888 -2401.22685
[324,] -2401.22685 26723.56888
[325,] -4812.37279 -2401.22685
[326,] 26723.56888 -4812.37279
[327,] 26723.56888 26723.56888
[328,] -2401.22685 26723.56888
[329,] -2401.22685 -2401.22685
[330,] -2401.22685 -2401.22685
[331,] -2401.22685 -2401.22685
[332,] -2401.22685 -2401.22685
[333,] -944.13962 -2401.22685
[334,] -2022.08210 -944.13962
[335,] -2401.22685 -2022.08210
[336,] -2729.53365 -2401.22685
[337,] -2401.22685 -2729.53365
[338,] -2356.13091 -2401.22685
[339,] 1358.77905 -2356.13091
[340,] -14305.02958 1358.77905
[341,] -2401.22685 -14305.02958
[342,] -2401.22685 -2401.22685
[343,] -2401.22685 -2401.22685
[344,] -2401.22685 -2401.22685
[345,] -2401.22685 -2401.22685
[346,] -9488.63760 -2401.22685
[347,] -2401.22685 -9488.63760
[348,] -2401.22685 -2401.22685
[349,] -2401.22685 -2401.22685
[350,] 2939.06427 -2401.22685
[351,] -1578.88762 2939.06427
[352,] 13632.90069 -1578.88762
[353,] -7486.21423 13632.90069
[354,] -2401.22685 -7486.21423
[355,] -2401.22685 -2401.22685
[356,] -24740.26924 -2401.22685
[357,] -2401.22685 -24740.26924
[358,] -2401.22685 -2401.22685
[359,] -2401.22685 -2401.22685
[360,] -1630.79340 -2401.22685
[361,] -2401.22685 -1630.79340
[362,] -2401.22685 -2401.22685
[363,] 18714.84729 -2401.22685
[364,] -1899.99534 18714.84729
[365,] -2401.22685 -1899.99534
[366,] -2401.22685 -2401.22685
[367,] -2401.22685 -2401.22685
[368,] -2401.22685 -2401.22685
[369,] -2401.22685 -2401.22685
[370,] -2401.22685 -2401.22685
[371,] 3486.88107 -2401.22685
[372,] -25749.55270 3486.88107
[373,] -20702.21744 -25749.55270
[374,] -2401.22685 -20702.21744
[375,] -2401.22685 -2401.22685
[376,] 4330.91660 -2401.22685
[377,] -1503.97547 4330.91660
[378,] -2401.22685 -1503.97547
[379,] -2401.22685 -2401.22685
[380,] -2401.22685 -2401.22685
[381,] -3409.33390 -2401.22685
[382,] 41667.47187 -3409.33390
[383,] -2401.22685 41667.47187
[384,] -763.76613 -2401.22685
[385,] -86823.08487 -763.76613
[386,] 3994.32689 -86823.08487
[387,] 7260.16022 3994.32689
[388,] -2997.41084 7260.16022
[389,] -9409.52872 -2997.41084
[390,] 5753.07660 -9409.52872
[391,] -11857.94258 5753.07660
[392,] -22132.20075 -11857.94258
[393,] 39537.25471 -22132.20075
[394,] -155995.44039 39537.25471
[395,] -17357.95088 -155995.44039
[396,] -25973.91882 -17357.95088
[397,] 107289.31148 -25973.91882
[398,] -12279.43491 107289.31148
[399,] -84213.96533 -12279.43491
[400,] -496.61868 -84213.96533
[401,] 12242.61718 -496.61868
[402,] -7008.00459 12242.61718
[403,] -51150.14027 -7008.00459
[404,] -16445.04507 -51150.14027
[405,] -13369.73486 -16445.04507
[406,] 2334.62395 -13369.73486
[407,] 4136.64604 2334.62395
[408,] -15614.21169 4136.64604
[409,] -48274.99259 -15614.21169
[410,] -18647.91267 -48274.99259
[411,] -219572.24477 -18647.91267
[412,] -22111.50353 -219572.24477
[413,] -110588.42235 -22111.50353
[414,] -58146.69165 -110588.42235
[415,] -57201.32722 -58146.69165
[416,] 155507.03181 -57201.32722
[417,] -34356.48016 155507.03181
[418,] -14820.25975 -34356.48016
[419,] -46731.17993 -14820.25975
[420,] -85486.26190 -46731.17993
[421,] -82041.92819 -85486.26190
[422,] -627186.38066 -82041.92819
[423,] -331586.22377 -627186.38066
[424,] -68145.58823 -331586.22377
[425,] -197325.48460 -68145.58823
[426,] -251752.72344 -197325.48460
[427,] -251741.39555 -251752.72344
[428,] -127326.89068 -251741.39555
[429,] 135682.64604 -127326.89068
[430,] -363758.36870 135682.64604
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 3055410.06433 -335057.38899
2 659410.54377 3055410.06433
3 -405429.22905 659410.54377
4 -743164.38986 -405429.22905
5 -89456.37662 -743164.38986
6 957456.07911 -89456.37662
7 358218.87563 957456.07911
8 -515485.62638 358218.87563
9 656837.66595 -515485.62638
10 418629.65698 656837.66595
11 790336.11651 418629.65698
12 786396.36249 790336.11651
13 -496204.59220 786396.36249
14 -198517.34823 -496204.59220
15 573684.82775 -198517.34823
16 -252041.74304 573684.82775
17 -311194.45538 -252041.74304
18 -428313.88915 -311194.45538
19 -21431.49731 -428313.88915
20 139295.64762 -21431.49731
21 -335211.15220 139295.64762
22 114533.03501 -335211.15220
23 -174367.24805 114533.03501
24 -66395.70200 -174367.24805
25 301730.34172 -66395.70200
26 -543727.48346 301730.34172
27 -174098.17192 -543727.48346
28 593025.71736 -174098.17192
29 341865.12044 593025.71736
30 252881.74560 341865.12044
31 -132561.70348 252881.74560
32 163074.51332 -132561.70348
33 105181.03944 163074.51332
34 -175459.84008 105181.03944
35 26963.16687 -175459.84008
36 306078.95930 26963.16687
37 161359.42712 306078.95930
38 138591.64288 161359.42712
39 -121232.08493 138591.64288
40 -315987.14980 -121232.08493
41 21458.74625 -315987.14980
42 104529.24113 21458.74625
43 -48445.33886 104529.24113
44 62799.44298 -48445.33886
45 -127194.57210 62799.44298
46 118925.53513 -127194.57210
47 -219278.61704 118925.53513
48 -14451.06244 -219278.61704
49 -142981.86240 -14451.06244
50 105768.72434 -142981.86240
51 -34627.40188 105768.72434
52 -219322.05982 -34627.40188
53 70651.42675 -219322.05982
54 -281080.27993 70651.42675
55 124875.96679 -281080.27993
56 232448.43092 124875.96679
57 85407.63172 232448.43092
58 -63984.83679 85407.63172
59 -220578.16036 -63984.83679
60 -23272.87970 -220578.16036
61 -78225.09783 -23272.87970
62 -20749.17848 -78225.09783
63 -216463.75382 -20749.17848
64 68796.51174 -216463.75382
65 23572.03863 68796.51174
66 -137137.24569 23572.03863
67 -268.85031 -137137.24569
68 -84555.82611 -268.85031
69 -55112.87409 -84555.82611
70 -89236.14587 -55112.87409
71 153263.11954 -89236.14587
72 -121884.06595 153263.11954
73 -46107.67637 -121884.06595
74 -94572.99222 -46107.67637
75 -73086.17103 -94572.99222
76 -105602.61775 -73086.17103
77 -85644.93657 -105602.61775
78 37088.85778 -85644.93657
79 21253.31059 37088.85778
80 29611.70914 21253.31059
81 -45462.15309 29611.70914
82 -134129.65691 -45462.15309
83 36887.62176 -134129.65691
84 -23125.20612 36887.62176
85 -116597.88333 -23125.20612
86 70596.47129 -116597.88333
87 -54748.11133 70596.47129
88 -88191.48526 -54748.11133
89 -33468.48369 -88191.48526
90 -20293.38577 -33468.48369
91 -41947.40635 -20293.38577
92 122154.90850 -41947.40635
93 -228180.60334 122154.90850
94 119037.29664 -228180.60334
95 -11477.33974 119037.29664
96 -1554.05235 -11477.33974
97 49267.09848 -1554.05235
98 -71063.71119 49267.09848
99 -26793.91038 -71063.71119
100 -34155.89201 -26793.91038
101 -116393.68974 -34155.89201
102 -71182.11744 -116393.68974
103 -63711.14507 -71182.11744
104 39309.00263 -63711.14507
105 -78901.81531 39309.00263
106 -233869.45225 -78901.81531
107 107762.43360 -233869.45225
108 -87569.45515 107762.43360
109 -104668.69729 -87569.45515
110 82970.34063 -104668.69729
111 -253759.38582 82970.34063
112 -175074.28537 -253759.38582
113 -121209.46808 -175074.28537
114 14252.74836 -121209.46808
115 55220.38255 14252.74836
116 -26722.55853 55220.38255
117 65798.28707 -26722.55853
118 -176082.65391 65798.28707
119 26570.83341 -176082.65391
120 2477.23170 26570.83341
121 13431.45499 2477.23170
122 -36512.52694 13431.45499
123 80946.69568 -36512.52694
124 -100397.40476 80946.69568
125 -201133.25705 -100397.40476
126 -2401.22685 -201133.25705
127 -987.01213 -2401.22685
128 9809.10476 -987.01213
129 -2401.22685 9809.10476
130 -386.54690 -2401.22685
131 85481.18719 -386.54690
132 -2401.22685 85481.18719
133 -3753.75972 -2401.22685
134 -2401.22685 -3753.75972
135 7599.65197 -2401.22685
136 175.79713 7599.65197
137 -2401.22685 175.79713
138 -2401.22685 -2401.22685
139 -2401.22685 -2401.22685
140 3951.31681 -2401.22685
141 -2401.22685 3951.31681
142 17400.26177 -2401.22685
143 -1675.86504 17400.26177
144 16223.04061 -1675.86504
145 13349.51735 16223.04061
146 -2401.22685 13349.51735
147 -1714.53127 -2401.22685
148 -2401.22685 -1714.53127
149 -2564.61669 -2401.22685
150 543.65038 -2564.61669
151 22998.67080 543.65038
152 -2401.22685 22998.67080
153 -14612.79492 -2401.22685
154 -2327.64702 -14612.79492
155 -2401.22685 -2327.64702
156 -2401.22685 -2401.22685
157 -2399.01202 -2401.22685
158 15826.16411 -2399.01202
159 -2401.22685 15826.16411
160 -2401.22685 -2401.22685
161 -2401.22685 -2401.22685
162 -2401.22685 -2401.22685
163 11382.48822 -2401.22685
164 -60.46908 11382.48822
165 -2401.22685 -60.46908
166 -2401.22685 -2401.22685
167 110015.49676 -2401.22685
168 -2401.22685 110015.49676
169 -2401.22685 -2401.22685
170 -2401.22685 -2401.22685
171 -4700.09393 -2401.22685
172 -2401.22685 -4700.09393
173 239.82022 -2401.22685
174 -2401.22685 239.82022
175 -1522.71355 -2401.22685
176 -2401.22685 -1522.71355
177 24853.40506 -2401.22685
178 -2401.22685 24853.40506
179 -2401.22685 -2401.22685
180 39696.20723 -2401.22685
181 23643.86380 39696.20723
182 29378.05747 23643.86380
183 -12686.20265 29378.05747
184 29249.09785 -12686.20265
185 11165.20288 29249.09785
186 33626.14463 11165.20288
187 -1350.71885 33626.14463
188 -32456.24329 -1350.71885
189 -2401.22685 -32456.24329
190 8980.09568 -2401.22685
191 -2401.22685 8980.09568
192 3290.42860 -2401.22685
193 -2401.22685 3290.42860
194 -2200.22282 -2401.22685
195 26723.56888 -2200.22282
196 -41404.00191 26723.56888
197 26723.56888 -41404.00191
198 26723.56888 26723.56888
199 -13875.78642 26723.56888
200 26723.56888 -13875.78642
201 26723.56888 26723.56888
202 26723.56888 26723.56888
203 25585.75864 26723.56888
204 -4271.33958 25585.75864
205 26784.69571 -4271.33958
206 22287.18057 26784.69571
207 26723.56888 22287.18057
208 44234.67155 26723.56888
209 26723.56888 44234.67155
210 -2401.22685 26723.56888
211 -2401.22685 -2401.22685
212 11362.80316 -2401.22685
213 -14186.21899 11362.80316
214 -4637.33104 -14186.21899
215 -2401.22685 -4637.33104
216 -486.64678 -2401.22685
217 6769.86053 -486.64678
218 -2401.22685 6769.86053
219 27804.36286 -2401.22685
220 -27290.37324 27804.36286
221 -2401.22685 -27290.37324
222 -1963.41995 -2401.22685
223 -2401.22685 -1963.41995
224 -1935.75325 -2401.22685
225 4557.43857 -1935.75325
226 -2401.22685 4557.43857
227 -2401.22685 -2401.22685
228 -20062.40201 -2401.22685
229 -2401.22685 -20062.40201
230 -13551.83006 -2401.22685
231 -2401.22685 -13551.83006
232 -2401.22685 -2401.22685
233 6727.90438 -2401.22685
234 -2401.22685 6727.90438
235 -2333.83700 -2401.22685
236 -2401.22685 -2333.83700
237 339.50561 -2401.22685
238 -2401.22685 339.50561
239 -2401.22685 -2401.22685
240 -2401.22685 -2401.22685
241 3152.22517 -2401.22685
242 -2401.22685 3152.22517
243 720.53012 -2401.22685
244 -7170.17477 720.53012
245 26723.56888 -7170.17477
246 26723.56888 26723.56888
247 -8564.63380 26723.56888
248 26723.56888 -8564.63380
249 26723.56888 26723.56888
250 3156.75624 26723.56888
251 26723.56888 3156.75624
252 -18278.38317 26723.56888
253 -17986.72143 -18278.38317
254 -42817.86848 -17986.72143
255 -2401.22685 -42817.86848
256 26006.40867 -2401.22685
257 32093.64563 26006.40867
258 26723.56888 32093.64563
259 -2401.22685 26723.56888
260 26723.56888 -2401.22685
261 26723.56888 26723.56888
262 -9433.34743 26723.56888
263 10546.87598 -9433.34743
264 26723.56888 10546.87598
265 -2401.22685 26723.56888
266 26723.56888 -2401.22685
267 16171.00545 26723.56888
268 26723.56888 16171.00545
269 -14299.80626 26723.56888
270 26723.56888 -14299.80626
271 26723.56888 26723.56888
272 26723.56888 26723.56888
273 -3108.46350 26723.56888
274 -2401.22685 -3108.46350
275 8612.68765 -2401.22685
276 8347.42108 8612.68765
277 -2401.22685 8347.42108
278 6172.76214 -2401.22685
279 -20844.16653 6172.76214
280 26723.56888 -20844.16653
281 -2401.22685 26723.56888
282 23272.10884 -2401.22685
283 -2401.22685 23272.10884
284 -24024.75157 -2401.22685
285 -2401.22685 -24024.75157
286 26723.56888 -2401.22685
287 -2486.74385 26723.56888
288 26723.56888 -2486.74385
289 16783.25672 26723.56888
290 26723.56888 16783.25672
291 7464.49902 26723.56888
292 -6947.57562 7464.49902
293 -2401.22685 -6947.57562
294 -1737.31668 -2401.22685
295 51474.88847 -1737.31668
296 -2401.22685 51474.88847
297 26723.56888 -2401.22685
298 -2401.22685 26723.56888
299 -2401.22685 -2401.22685
300 25184.63237 -2401.22685
301 -9291.67440 25184.63237
302 -2401.22685 -9291.67440
303 26723.56888 -2401.22685
304 34355.02811 26723.56888
305 -666.00381 34355.02811
306 26723.56888 -666.00381
307 22464.83173 26723.56888
308 21670.74259 22464.83173
309 26723.56888 21670.74259
310 26723.56888 26723.56888
311 -26255.79903 26723.56888
312 26723.56888 -26255.79903
313 29891.95209 26723.56888
314 -1414.38736 29891.95209
315 2902.88817 -1414.38736
316 -643.29997 2902.88817
317 -9108.09383 -643.29997
318 -1238.86183 -9108.09383
319 -2401.22685 -1238.86183
320 -2401.22685 -2401.22685
321 -2401.22685 -2401.22685
322 -2401.22685 -2401.22685
323 26723.56888 -2401.22685
324 -2401.22685 26723.56888
325 -4812.37279 -2401.22685
326 26723.56888 -4812.37279
327 26723.56888 26723.56888
328 -2401.22685 26723.56888
329 -2401.22685 -2401.22685
330 -2401.22685 -2401.22685
331 -2401.22685 -2401.22685
332 -2401.22685 -2401.22685
333 -944.13962 -2401.22685
334 -2022.08210 -944.13962
335 -2401.22685 -2022.08210
336 -2729.53365 -2401.22685
337 -2401.22685 -2729.53365
338 -2356.13091 -2401.22685
339 1358.77905 -2356.13091
340 -14305.02958 1358.77905
341 -2401.22685 -14305.02958
342 -2401.22685 -2401.22685
343 -2401.22685 -2401.22685
344 -2401.22685 -2401.22685
345 -2401.22685 -2401.22685
346 -9488.63760 -2401.22685
347 -2401.22685 -9488.63760
348 -2401.22685 -2401.22685
349 -2401.22685 -2401.22685
350 2939.06427 -2401.22685
351 -1578.88762 2939.06427
352 13632.90069 -1578.88762
353 -7486.21423 13632.90069
354 -2401.22685 -7486.21423
355 -2401.22685 -2401.22685
356 -24740.26924 -2401.22685
357 -2401.22685 -24740.26924
358 -2401.22685 -2401.22685
359 -2401.22685 -2401.22685
360 -1630.79340 -2401.22685
361 -2401.22685 -1630.79340
362 -2401.22685 -2401.22685
363 18714.84729 -2401.22685
364 -1899.99534 18714.84729
365 -2401.22685 -1899.99534
366 -2401.22685 -2401.22685
367 -2401.22685 -2401.22685
368 -2401.22685 -2401.22685
369 -2401.22685 -2401.22685
370 -2401.22685 -2401.22685
371 3486.88107 -2401.22685
372 -25749.55270 3486.88107
373 -20702.21744 -25749.55270
374 -2401.22685 -20702.21744
375 -2401.22685 -2401.22685
376 4330.91660 -2401.22685
377 -1503.97547 4330.91660
378 -2401.22685 -1503.97547
379 -2401.22685 -2401.22685
380 -2401.22685 -2401.22685
381 -3409.33390 -2401.22685
382 41667.47187 -3409.33390
383 -2401.22685 41667.47187
384 -763.76613 -2401.22685
385 -86823.08487 -763.76613
386 3994.32689 -86823.08487
387 7260.16022 3994.32689
388 -2997.41084 7260.16022
389 -9409.52872 -2997.41084
390 5753.07660 -9409.52872
391 -11857.94258 5753.07660
392 -22132.20075 -11857.94258
393 39537.25471 -22132.20075
394 -155995.44039 39537.25471
395 -17357.95088 -155995.44039
396 -25973.91882 -17357.95088
397 107289.31148 -25973.91882
398 -12279.43491 107289.31148
399 -84213.96533 -12279.43491
400 -496.61868 -84213.96533
401 12242.61718 -496.61868
402 -7008.00459 12242.61718
403 -51150.14027 -7008.00459
404 -16445.04507 -51150.14027
405 -13369.73486 -16445.04507
406 2334.62395 -13369.73486
407 4136.64604 2334.62395
408 -15614.21169 4136.64604
409 -48274.99259 -15614.21169
410 -18647.91267 -48274.99259
411 -219572.24477 -18647.91267
412 -22111.50353 -219572.24477
413 -110588.42235 -22111.50353
414 -58146.69165 -110588.42235
415 -57201.32722 -58146.69165
416 155507.03181 -57201.32722
417 -34356.48016 155507.03181
418 -14820.25975 -34356.48016
419 -46731.17993 -14820.25975
420 -85486.26190 -46731.17993
421 -82041.92819 -85486.26190
422 -627186.38066 -82041.92819
423 -331586.22377 -627186.38066
424 -68145.58823 -331586.22377
425 -197325.48460 -68145.58823
426 -251752.72344 -197325.48460
427 -251741.39555 -251752.72344
428 -127326.89068 -251741.39555
429 135682.64604 -127326.89068
430 -363758.36870 135682.64604
> plot(z,main=paste('Residual Lag plot, lowess, and regression line'), ylab='values of Residuals', xlab='lagged values of Residuals')
> lines(lowess(z))
> abline(lm(z))
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/freestat/rcomp/tmp/7i9rl1291330094.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> acf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Autocorrelation Function')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/freestat/rcomp/tmp/8tj861291330094.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> pacf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Partial Autocorrelation Function')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/freestat/rcomp/tmp/9tj861291330094.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> opar <- par(mfrow = c(2,2), oma = c(0, 0, 1.1, 0))
> plot(mylm, las = 1, sub='Residual Diagnostics')
> par(opar)
> dev.off()
null device
1
> if (n > n25) {
+ postscript(file="/var/www/html/freestat/rcomp/tmp/104a781291330094.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
+ plot(kp3:nmkm3,gqarr[,2], main='Goldfeld-Quandt test',ylab='2-sided p-value',xlab='breakpoint')
+ grid()
+ dev.off()
+ }
null device
1
>
> #Note: the /var/www/html/freestat/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/www/html/freestat/rcomp/createtable")
>
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Estimated Regression Equation', 1, TRUE)
> a<-table.row.end(a)
> myeq <- colnames(x)[1]
> myeq <- paste(myeq, '[t] = ', sep='')
> for (i in 1:k){
+ if (mysum$coefficients[i,1] > 0) myeq <- paste(myeq, '+', '')
+ myeq <- paste(myeq, mysum$coefficients[i,1], sep=' ')
+ if (rownames(mysum$coefficients)[i] != '(Intercept)') {
+ myeq <- paste(myeq, rownames(mysum$coefficients)[i], sep='')
+ if (rownames(mysum$coefficients)[i] != 't') myeq <- paste(myeq, '[t]', sep='')
+ }
+ }
> myeq <- paste(myeq, ' + e[t]')
> a<-table.row.start(a)
> a<-table.element(a, myeq)
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/www/html/freestat/rcomp/tmp/117a6e1291330094.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,hyperlink('http://www.xycoon.com/ols1.htm','Multiple Linear Regression - Ordinary Least Squares',''), 6, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Variable',header=TRUE)
> a<-table.element(a,'Parameter',header=TRUE)
> a<-table.element(a,'S.D.',header=TRUE)
> a<-table.element(a,'T-STAT
H0: parameter = 0',header=TRUE)
> a<-table.element(a,'2-tail p-value',header=TRUE)
> a<-table.element(a,'1-tail p-value',header=TRUE)
> a<-table.row.end(a)
> for (i in 1:k){
+ a<-table.row.start(a)
+ a<-table.element(a,rownames(mysum$coefficients)[i],header=TRUE)
+ a<-table.element(a,mysum$coefficients[i,1])
+ a<-table.element(a, round(mysum$coefficients[i,2],6))
+ a<-table.element(a, round(mysum$coefficients[i,3],4))
+ a<-table.element(a, round(mysum$coefficients[i,4],6))
+ a<-table.element(a, round(mysum$coefficients[i,4]/2,6))
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/www/html/freestat/rcomp/tmp/12st4k1291330094.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Regression Statistics', 2, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple R',1,TRUE)
> a<-table.element(a, sqrt(mysum$r.squared))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'R-squared',1,TRUE)
> a<-table.element(a, mysum$r.squared)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Adjusted R-squared',1,TRUE)
> a<-table.element(a, mysum$adj.r.squared)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (value)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[1])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (DF numerator)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[2])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (DF denominator)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[3])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'p-value',1,TRUE)
> a<-table.element(a, 1-pf(mysum$fstatistic[1],mysum$fstatistic[2],mysum$fstatistic[3]))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Residual Statistics', 2, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Residual Standard Deviation',1,TRUE)
> a<-table.element(a, mysum$sigma)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Sum Squared Residuals',1,TRUE)
> a<-table.element(a, sum(myerror*myerror))
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/www/html/freestat/rcomp/tmp/137l2b1291330094.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Actuals, Interpolation, and Residuals', 4, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Time or Index', 1, TRUE)
> a<-table.element(a, 'Actuals', 1, TRUE)
> a<-table.element(a, 'Interpolation
Forecast', 1, TRUE)
> a<-table.element(a, 'Residuals
Prediction Error', 1, TRUE)
> a<-table.row.end(a)
> for (i in 1:n) {
+ a<-table.row.start(a)
+ a<-table.element(a,i, 1, TRUE)
+ a<-table.element(a,x[i])
+ a<-table.element(a,x[i]-mysum$resid[i])
+ a<-table.element(a,mysum$resid[i])
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/www/html/freestat/rcomp/tmp/14s3jz1291330094.tab")
> if (n > n25) {
+ a<-table.start()
+ a<-table.row.start(a)
+ a<-table.element(a,'Goldfeld-Quandt test for Heteroskedasticity',4,TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'p-values',header=TRUE)
+ a<-table.element(a,'Alternative Hypothesis',3,header=TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'breakpoint index',header=TRUE)
+ a<-table.element(a,'greater',header=TRUE)
+ a<-table.element(a,'2-sided',header=TRUE)
+ a<-table.element(a,'less',header=TRUE)
+ a<-table.row.end(a)
+ for (mypoint in kp3:nmkm3) {
+ a<-table.row.start(a)
+ a<-table.element(a,mypoint,header=TRUE)
+ a<-table.element(a,gqarr[mypoint-kp3+1,1])
+ a<-table.element(a,gqarr[mypoint-kp3+1,2])
+ a<-table.element(a,gqarr[mypoint-kp3+1,3])
+ a<-table.row.end(a)
+ }
+ a<-table.end(a)
+ table.save(a,file="/var/www/html/freestat/rcomp/tmp/15v4h51291330094.tab")
+ a<-table.start()
+ a<-table.row.start(a)
+ a<-table.element(a,'Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity',4,TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'Description',header=TRUE)
+ a<-table.element(a,'# significant tests',header=TRUE)
+ a<-table.element(a,'% significant tests',header=TRUE)
+ a<-table.element(a,'OK/NOK',header=TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'1% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant1)
+ a<-table.element(a,numsignificant1/numgqtests)
+ if (numsignificant1/numgqtests < 0.01) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'5% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant5)
+ a<-table.element(a,numsignificant5/numgqtests)
+ if (numsignificant5/numgqtests < 0.05) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'10% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant10)
+ a<-table.element(a,numsignificant10/numgqtests)
+ if (numsignificant10/numgqtests < 0.1) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.end(a)
+ table.save(a,file="/var/www/html/freestat/rcomp/tmp/162nh81291330095.tab")
+ }
>
> try(system("convert tmp/1frax1291330094.ps tmp/1frax1291330094.png",intern=TRUE))
character(0)
> try(system("convert tmp/2frax1291330094.ps tmp/2frax1291330094.png",intern=TRUE))
character(0)
> try(system("convert tmp/3frax1291330094.ps tmp/3frax1291330094.png",intern=TRUE))
character(0)
> try(system("convert tmp/4q0si1291330094.ps tmp/4q0si1291330094.png",intern=TRUE))
character(0)
> try(system("convert tmp/5q0si1291330094.ps tmp/5q0si1291330094.png",intern=TRUE))
character(0)
> try(system("convert tmp/6q0si1291330094.ps tmp/6q0si1291330094.png",intern=TRUE))
character(0)
> try(system("convert tmp/7i9rl1291330094.ps tmp/7i9rl1291330094.png",intern=TRUE))
character(0)
> try(system("convert tmp/8tj861291330094.ps tmp/8tj861291330094.png",intern=TRUE))
character(0)
> try(system("convert tmp/9tj861291330094.ps tmp/9tj861291330094.png",intern=TRUE))
character(0)
> try(system("convert tmp/104a781291330094.ps tmp/104a781291330094.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
17.837 3.053 18.343