R version 2.8.0 (2008-10-20)
Copyright (C) 2008 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
Natural language support but running in an English locale
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> x <- array(list(1
+ ,162556
+ ,1081
+ ,807
+ ,213118
+ ,230380558
+ ,6282929
+ ,1
+ ,29790
+ ,309
+ ,444
+ ,81767
+ ,25266003
+ ,4324047
+ ,1
+ ,87550
+ ,458
+ ,412
+ ,153198
+ ,70164684
+ ,4108272
+ ,0
+ ,84738
+ ,588
+ ,428
+ ,-26007
+ ,-15292116
+ ,-1212617
+ ,1
+ ,54660
+ ,299
+ ,312
+ ,126942
+ ,37955658
+ ,1485329
+ ,1
+ ,42634
+ ,156
+ ,166
+ ,157214
+ ,24525384
+ ,1779876
+ ,0
+ ,40949
+ ,481
+ ,263
+ ,129352
+ ,62218312
+ ,1367203
+ ,1
+ ,42312
+ ,323
+ ,237
+ ,234817
+ ,75845891
+ ,2519076
+ ,1
+ ,37704
+ ,452
+ ,228
+ ,60448
+ ,27322496
+ ,912684
+ ,1
+ ,16275
+ ,109
+ ,129
+ ,47818
+ ,5212162
+ ,1443586
+ ,0
+ ,25830
+ ,115
+ ,104
+ ,245546
+ ,28237790
+ ,1220017
+ ,0
+ ,12679
+ ,110
+ ,122
+ ,48020
+ ,5282200
+ ,984885
+ ,1
+ ,18014
+ ,239
+ ,393
+ ,-1710
+ ,-408690
+ ,1457425
+ ,0
+ ,43556
+ ,247
+ ,190
+ ,32648
+ ,8064056
+ ,-572920
+ ,1
+ ,24524
+ ,497
+ ,275
+ ,95350
+ ,47388950
+ ,929144
+ ,0
+ ,6532
+ ,103
+ ,62
+ ,151352
+ ,15589256
+ ,1151176
+ ,0
+ ,7123
+ ,109
+ ,102
+ ,288170
+ ,31410530
+ ,790090
+ ,1
+ ,20813
+ ,502
+ ,255
+ ,114337
+ ,57397174
+ ,774497
+ ,1
+ ,37597
+ ,248
+ ,234
+ ,37884
+ ,9395232
+ ,990576
+ ,0
+ ,17821
+ ,373
+ ,277
+ ,122844
+ ,45820812
+ ,454195
+ ,1
+ ,12988
+ ,119
+ ,73
+ ,82340
+ ,9798460
+ ,876607
+ ,1
+ ,22330
+ ,84
+ ,67
+ ,79801
+ ,6703284
+ ,711969
+ ,0
+ ,13326
+ ,102
+ ,103
+ ,165548
+ ,16885896
+ ,702380
+ ,0
+ ,16189
+ ,295
+ ,290
+ ,116384
+ ,34333280
+ ,264449
+ ,0
+ ,7146
+ ,105
+ ,83
+ ,134028
+ ,14072940
+ ,450033
+ ,0
+ ,15824
+ ,64
+ ,56
+ ,63838
+ ,4085632
+ ,541063
+ ,1
+ ,26088
+ ,267
+ ,224
+ ,74996
+ ,20023932
+ ,588864
+ ,0
+ ,11326
+ ,129
+ ,64
+ ,31080
+ ,4009320
+ ,-37216
+ ,0
+ ,8568
+ ,37
+ ,34
+ ,32168
+ ,1190216
+ ,783310
+ ,0
+ ,14416
+ ,361
+ ,139
+ ,49857
+ ,17998377
+ ,467359
+ ,1
+ ,3369
+ ,28
+ ,26
+ ,87161
+ ,2440508
+ ,688779
+ ,1
+ ,11819
+ ,85
+ ,70
+ ,106113
+ ,9019605
+ ,608419
+ ,1
+ ,6620
+ ,44
+ ,40
+ ,80570
+ ,3545080
+ ,696348
+ ,1
+ ,4519
+ ,49
+ ,42
+ ,102129
+ ,5004321
+ ,597793
+ ,0
+ ,2220
+ ,22
+ ,12
+ ,301670
+ ,6636740
+ ,821730
+ ,0
+ ,18562
+ ,155
+ ,211
+ ,102313
+ ,15858515
+ ,377934
+ ,0
+ ,10327
+ ,91
+ ,74
+ ,88577
+ ,8060507
+ ,651939
+ ,1
+ ,5336
+ ,81
+ ,80
+ ,112477
+ ,9110637
+ ,697458
+ ,1
+ ,2365
+ ,79
+ ,83
+ ,191778
+ ,15150462
+ ,700368
+ ,0
+ ,4069
+ ,145
+ ,131
+ ,79804
+ ,11571580
+ ,225986
+ ,0
+ ,7710
+ ,816
+ ,187
+ ,128294
+ ,104687904
+ ,348695
+ ,0
+ ,13718
+ ,61
+ ,56
+ ,96448
+ ,5883328
+ ,373683
+ ,0
+ ,4525
+ ,226
+ ,89
+ ,93811
+ ,21201286
+ ,501709
+ ,0
+ ,6869
+ ,105
+ ,88
+ ,117520
+ ,12339600
+ ,413743
+ ,0
+ ,4628
+ ,62
+ ,39
+ ,69159
+ ,4287858
+ ,379825
+ ,1
+ ,3653
+ ,24
+ ,24
+ ,101792
+ ,2443008
+ ,336260
+ ,1
+ ,1265
+ ,26
+ ,19
+ ,210568
+ ,5474768
+ ,636765
+ ,1
+ ,7489
+ ,322
+ ,149
+ ,136996
+ ,44112712
+ ,481231
+ ,0
+ ,4901
+ ,84
+ ,58
+ ,121920
+ ,10241280
+ ,469107
+ ,0
+ ,2284
+ ,33
+ ,41
+ ,76403
+ ,2521299
+ ,211928
+ ,1
+ ,3160
+ ,108
+ ,90
+ ,108094
+ ,11674152
+ ,563925
+ ,1
+ ,4150
+ ,150
+ ,136
+ ,134759
+ ,20213850
+ ,511939
+ ,1
+ ,7285
+ ,115
+ ,97
+ ,188873
+ ,21720395
+ ,521016
+ ,1
+ ,1134
+ ,162
+ ,63
+ ,146216
+ ,23686992
+ ,543856
+ ,1
+ ,4658
+ ,158
+ ,114
+ ,156608
+ ,24744064
+ ,329304
+ ,0
+ ,2384
+ ,97
+ ,77
+ ,61348
+ ,5950756
+ ,423262
+ ,0
+ ,3748
+ ,9
+ ,6
+ ,50350
+ ,453150
+ ,509665
+ ,0
+ ,5371
+ ,66
+ ,47
+ ,87720
+ ,5789520
+ ,455881
+ ,0
+ ,1285
+ ,107
+ ,51
+ ,99489
+ ,10645323
+ ,367772
+ ,1
+ ,9327
+ ,101
+ ,85
+ ,87419
+ ,8829319
+ ,406339
+ ,1
+ ,5565
+ ,47
+ ,43
+ ,94355
+ ,4434685
+ ,493408
+ ,0
+ ,1528
+ ,38
+ ,32
+ ,60326
+ ,2292388
+ ,232942
+ ,1
+ ,3122
+ ,34
+ ,25
+ ,94670
+ ,3218780
+ ,416002
+ ,1
+ ,7317
+ ,84
+ ,74
+ ,82425
+ ,6923700
+ ,337430
+ ,0
+ ,2675
+ ,79
+ ,54
+ ,59017
+ ,4662343
+ ,361517
+ ,0
+ ,13253
+ ,947
+ ,251
+ ,90829
+ ,86015063
+ ,360962
+ ,0
+ ,880
+ ,74
+ ,15
+ ,80791
+ ,5978534
+ ,235561
+ ,1
+ ,2053
+ ,53
+ ,44
+ ,100423
+ ,5322419
+ ,408247
+ ,0
+ ,1424
+ ,94
+ ,73
+ ,131116
+ ,12324904
+ ,450296
+ ,1
+ ,4036
+ ,63
+ ,85
+ ,100269
+ ,6316947
+ ,418799
+ ,1
+ ,3045
+ ,58
+ ,49
+ ,27330
+ ,1585140
+ ,247405
+ ,0
+ ,5119
+ ,49
+ ,38
+ ,39039
+ ,1912911
+ ,378519
+ ,0
+ ,1431
+ ,34
+ ,35
+ ,106885
+ ,3634090
+ ,326638
+ ,0
+ ,554
+ ,11
+ ,9
+ ,79285
+ ,872135
+ ,328233
+ ,0
+ ,1975
+ ,35
+ ,34
+ ,118881
+ ,4160835
+ ,386225
+ ,1
+ ,1286
+ ,17
+ ,17
+ ,77623
+ ,1319591
+ ,283662
+ ,0
+ ,1012
+ ,47
+ ,29
+ ,114768
+ ,5394096
+ ,370225
+ ,0
+ ,810
+ ,43
+ ,11
+ ,74015
+ ,3182645
+ ,269236
+ ,0
+ ,1280
+ ,117
+ ,52
+ ,69465
+ ,8127405
+ ,365732
+ ,1
+ ,666
+ ,171
+ ,13
+ ,117869
+ ,20155599
+ ,420383
+ ,0
+ ,1380
+ ,26
+ ,29
+ ,60982
+ ,1585532
+ ,345811
+ ,1
+ ,4608
+ ,73
+ ,64
+ ,90131
+ ,6579563
+ ,431809
+ ,0
+ ,876
+ ,59
+ ,33
+ ,138971
+ ,8199289
+ ,418876
+ ,0
+ ,814
+ ,18
+ ,15
+ ,39625
+ ,713250
+ ,297476
+ ,0
+ ,514
+ ,15
+ ,15
+ ,102725
+ ,1540875
+ ,416776
+ ,1
+ ,5692
+ ,72
+ ,68
+ ,64239
+ ,4625208
+ ,357257
+ ,0
+ ,3642
+ ,86
+ ,100
+ ,90262
+ ,7762532
+ ,458343
+ ,0
+ ,540
+ ,14
+ ,13
+ ,103960
+ ,1455440
+ ,388386
+ ,0
+ ,2099
+ ,64
+ ,45
+ ,106611
+ ,6823104
+ ,358934
+ ,0
+ ,567
+ ,11
+ ,14
+ ,103345
+ ,1136795
+ ,407560
+ ,0
+ ,2001
+ ,52
+ ,36
+ ,95551
+ ,4968652
+ ,392558
+ ,1
+ ,2949
+ ,41
+ ,40
+ ,82903
+ ,3399023
+ ,373177
+ ,0
+ ,2253
+ ,99
+ ,68
+ ,63593
+ ,6295707
+ ,428370
+ ,1
+ ,6533
+ ,75
+ ,29
+ ,126910
+ ,9518250
+ ,369419
+ ,0
+ ,1889
+ ,45
+ ,43
+ ,37527
+ ,1688715
+ ,358649
+ ,1
+ ,3055
+ ,43
+ ,30
+ ,60247
+ ,2590621
+ ,376641
+ ,0
+ ,272
+ ,8
+ ,9
+ ,112995
+ ,903960
+ ,467427
+ ,1
+ ,1414
+ ,198
+ ,22
+ ,70184
+ ,13896432
+ ,364885
+ ,0
+ ,2564
+ ,22
+ ,19
+ ,130140
+ ,2863080
+ ,436230
+ ,1
+ ,1383
+ ,11
+ ,9
+ ,73221
+ ,805431
+ ,329118
+ ,1
+ ,1261
+ ,33
+ ,31
+ ,76114
+ ,2511762
+ ,317365
+ ,0
+ ,975
+ ,23
+ ,19
+ ,90534
+ ,2082282
+ ,286849
+ ,0
+ ,3366
+ ,80
+ ,55
+ ,108479
+ ,8678320
+ ,376685
+ ,0
+ ,576
+ ,18
+ ,8
+ ,113761
+ ,2047698
+ ,407198
+ ,0
+ ,1306
+ ,28
+ ,26
+ ,68696
+ ,1923488
+ ,377772
+ ,0
+ ,746
+ ,23
+ ,29
+ ,71561
+ ,1645903
+ ,271483
+ ,1
+ ,3192
+ ,60
+ ,48
+ ,59831
+ ,3589860
+ ,153661
+ ,1
+ ,2045
+ ,20
+ ,16
+ ,97890
+ ,1957800
+ ,513294
+ ,0
+ ,5477
+ ,59
+ ,45
+ ,101481
+ ,5987379
+ ,324881
+ ,1
+ ,1932
+ ,36
+ ,20
+ ,72954
+ ,2626344
+ ,264512
+ ,0
+ ,936
+ ,30
+ ,22
+ ,67939
+ ,2038170
+ ,420968
+ ,1
+ ,3437
+ ,47
+ ,33
+ ,48022
+ ,2257034
+ ,129302
+ ,0
+ ,5131
+ ,71
+ ,44
+ ,86111
+ ,6113881
+ ,191521
+ ,1
+ ,2397
+ ,14
+ ,13
+ ,74020
+ ,1036280
+ ,268673
+ ,1
+ ,1389
+ ,9
+ ,6
+ ,57530
+ ,517770
+ ,353179
+ ,0
+ ,1503
+ ,39
+ ,35
+ ,56364
+ ,2198196
+ ,354624
+ ,0
+ ,402
+ ,26
+ ,8
+ ,84990
+ ,2209740
+ ,363713
+ ,0
+ ,2239
+ ,21
+ ,17
+ ,88590
+ ,1860390
+ ,456657
+ ,1
+ ,2234
+ ,16
+ ,11
+ ,77200
+ ,1235200
+ ,211742
+ ,0
+ ,837
+ ,69
+ ,21
+ ,61262
+ ,4227078
+ ,338381
+ ,0
+ ,10579
+ ,92
+ ,92
+ ,110309
+ ,10148428
+ ,418530
+ ,0
+ ,875
+ ,14
+ ,12
+ ,67000
+ ,938000
+ ,351483
+ ,0
+ ,1395
+ ,103
+ ,108
+ ,93099
+ ,9589197
+ ,372928
+ ,1
+ ,1659
+ ,29
+ ,25
+ ,107577
+ ,3119733
+ ,485538
+ ,1
+ ,2647
+ ,37
+ ,17
+ ,62920
+ ,2328040
+ ,279268
+ ,1
+ ,3294
+ ,23
+ ,23
+ ,75832
+ ,1744136
+ ,219060
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,94
+ ,7
+ ,10
+ ,60793
+ ,425551
+ ,325314
+ ,0
+ ,422
+ ,28
+ ,23
+ ,57935
+ ,1622180
+ ,322046
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,34
+ ,8
+ ,7
+ ,60630
+ ,485040
+ ,325599
+ ,0
+ ,1558
+ ,63
+ ,25
+ ,55637
+ ,3505131
+ ,377028
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,43
+ ,3
+ ,20
+ ,60887
+ ,182661
+ ,323850
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,316
+ ,9
+ ,4
+ ,60505
+ ,544545
+ ,331514
+ ,0
+ ,115
+ ,13
+ ,10
+ ,60945
+ ,792285
+ ,325632
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,389
+ ,14
+ ,7
+ ,58990
+ ,825860
+ ,322265
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,1002
+ ,15
+ ,11
+ ,56750
+ ,851250
+ ,325906
+ ,0
+ ,36
+ ,3
+ ,4
+ ,60894
+ ,182682
+ ,325985
+ ,0
+ ,460
+ ,15
+ ,15
+ ,63346
+ ,950190
+ ,346145
+ ,0
+ ,309
+ ,11
+ ,9
+ ,56535
+ ,621885
+ ,325898
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,9
+ ,6
+ ,7
+ ,60835
+ ,365010
+ ,325356
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,14
+ ,1
+ ,0
+ ,61016
+ ,61016
+ ,325930
+ ,0
+ ,520
+ ,10
+ ,7
+ ,58650
+ ,586500
+ ,318020
+ ,0
+ ,1766
+ ,73
+ ,46
+ ,60438
+ ,4411974
+ ,326389
+ ,0
+ ,0
+ ,0
+ ,5
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,458
+ ,11
+ ,7
+ ,58625
+ ,644875
+ ,302925
+ ,0
+ ,20
+ ,3
+ ,2
+ ,60938
+ ,182814
+ ,325540
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,98
+ ,2
+ ,2
+ ,61490
+ ,122980
+ ,326736
+ ,0
+ ,405
+ ,7
+ ,5
+ ,60845
+ ,425915
+ ,340580
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,483
+ ,27
+ ,7
+ ,60830
+ ,1642410
+ ,331828
+ ,0
+ ,454
+ ,51
+ ,24
+ ,63261
+ ,3226311
+ ,323299
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,757
+ ,19
+ ,18
+ ,45689
+ ,868091
+ ,387722
+ ,0
+ ,0
+ ,0
+ ,4
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,36
+ ,4
+ ,3
+ ,61564
+ ,246256
+ ,324598
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,203
+ ,9
+ ,9
+ ,61938
+ ,557442
+ ,328726
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,90
+ ,8
+ ,6
+ ,60951
+ ,487608
+ ,325043
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,71
+ ,1
+ ,0
+ ,60745
+ ,60745
+ ,325806
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,972
+ ,34
+ ,19
+ ,71642
+ ,2435828
+ ,387732
+ ,1
+ ,531
+ ,10
+ ,11
+ ,71641
+ ,716410
+ ,349729
+ ,0
+ ,604
+ ,38
+ ,8
+ ,55792
+ ,2120096
+ ,332202
+ ,1
+ ,283
+ ,10
+ ,15
+ ,71873
+ ,718730
+ ,305442
+ ,1
+ ,23
+ ,5
+ ,5
+ ,62555
+ ,312775
+ ,329537
+ ,1
+ ,638
+ ,14
+ ,11
+ ,60370
+ ,845180
+ ,327055
+ ,1
+ ,699
+ ,16
+ ,10
+ ,64873
+ ,1037968
+ ,356245
+ ,0
+ ,149
+ ,5
+ ,6
+ ,62041
+ ,310205
+ ,328451
+ ,0
+ ,226
+ ,5
+ ,5
+ ,65745
+ ,328725
+ ,307062
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,275
+ ,4
+ ,7
+ ,59500
+ ,238000
+ ,331345
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,141
+ ,6
+ ,7
+ ,61630
+ ,369780
+ ,331824
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,28
+ ,2
+ ,3
+ ,60890
+ ,121780
+ ,325685
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,2566
+ ,80
+ ,74
+ ,113521
+ ,9081680
+ ,404480
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,472
+ ,20
+ ,19
+ ,80045
+ ,1600900
+ ,318314
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,203
+ ,27
+ ,12
+ ,50804
+ ,1371708
+ ,311807
+ ,1
+ ,496
+ ,17
+ ,12
+ ,87390
+ ,1485630
+ ,337724
+ ,1
+ ,10
+ ,2
+ ,5
+ ,61656
+ ,123312
+ ,326431
+ ,1
+ ,63
+ ,4
+ ,2
+ ,65688
+ ,262752
+ ,327556
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,1136
+ ,32
+ ,26
+ ,48522
+ ,1552704
+ ,356850
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,267
+ ,32
+ ,11
+ ,57640
+ ,1844480
+ ,322741
+ ,0
+ ,474
+ ,20
+ ,10
+ ,61977
+ ,1239540
+ ,310902
+ ,0
+ ,534
+ ,7
+ ,5
+ ,62620
+ ,438340
+ ,324295
+ ,0
+ ,0
+ ,0
+ ,2
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,15
+ ,8
+ ,6
+ ,60831
+ ,486648
+ ,326156
+ ,0
+ ,397
+ ,28
+ ,7
+ ,60646
+ ,1698088
+ ,326960
+ ,0
+ ,0
+ ,0
+ ,2
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,1061
+ ,20
+ ,22
+ ,56225
+ ,1124500
+ ,333411
+ ,0
+ ,288
+ ,4
+ ,3
+ ,60510
+ ,242040
+ ,297761
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,3
+ ,2
+ ,1
+ ,60698
+ ,121396
+ ,325536
+ ,0
+ ,0
+ ,0
+ ,6
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,20
+ ,2
+ ,1
+ ,60805
+ ,121610
+ ,325762
+ ,0
+ ,278
+ ,26
+ ,22
+ ,61404
+ ,1596504
+ ,327957
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,192
+ ,4
+ ,2
+ ,65276
+ ,261104
+ ,318521
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,317
+ ,9
+ ,7
+ ,63915
+ ,575235
+ ,319775
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,368
+ ,17
+ ,13
+ ,61686
+ ,1048662
+ ,332128
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,2
+ ,1
+ ,0
+ ,60743
+ ,60743
+ ,325486
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,53
+ ,6
+ ,6
+ ,60349
+ ,362094
+ ,325838
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,94
+ ,3
+ ,3
+ ,61360
+ ,184080
+ ,331767
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,24
+ ,8
+ ,7
+ ,59818
+ ,478544
+ ,324523
+ ,0
+ ,2332
+ ,4
+ ,2
+ ,72680
+ ,290720
+ ,339995
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,131
+ ,11
+ ,15
+ ,61808
+ ,679888
+ ,319582
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,206
+ ,9
+ ,9
+ ,53110
+ ,477990
+ ,307245
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,167
+ ,2
+ ,1
+ ,64245
+ ,128490
+ ,317967
+ ,0
+ ,622
+ ,73
+ ,38
+ ,73007
+ ,5329511
+ ,331488
+ ,0
+ ,885
+ ,85
+ ,49
+ ,82732
+ ,7032220
+ ,335452
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,365
+ ,8
+ ,6
+ ,54820
+ ,438560
+ ,334184
+ ,0
+ ,364
+ ,35
+ ,26
+ ,47705
+ ,1669675
+ ,313213
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,226
+ ,12
+ ,13
+ ,72835
+ ,874020
+ ,348678
+ ,0
+ ,307
+ ,15
+ ,10
+ ,58856
+ ,882840
+ ,328727
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,188
+ ,11
+ ,9
+ ,77655
+ ,854205
+ ,387978
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,138
+ ,6
+ ,26
+ ,69817
+ ,418902
+ ,336704
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,125
+ ,12
+ ,19
+ ,60798
+ ,729576
+ ,322076
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,282
+ ,30
+ ,12
+ ,62452
+ ,1873560
+ ,334272
+ ,0
+ ,335
+ ,33
+ ,23
+ ,64175
+ ,2117775
+ ,338197
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,813
+ ,82
+ ,24
+ ,67440
+ ,5530080
+ ,321024
+ ,0
+ ,176
+ ,28
+ ,8
+ ,68136
+ ,1907808
+ ,322145
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,249
+ ,72
+ ,26
+ ,56726
+ ,4084272
+ ,323351
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,333
+ ,13
+ ,9
+ ,70811
+ ,920543
+ ,327748
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,30
+ ,4
+ ,3
+ ,62045
+ ,248180
+ ,328157
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,249
+ ,62
+ ,13
+ ,54323
+ ,3368026
+ ,311594
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,165
+ ,24
+ ,12
+ ,62841
+ ,1508184
+ ,335962
+ ,0
+ ,453
+ ,21
+ ,19
+ ,81125
+ ,1703625
+ ,372426
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,53
+ ,14
+ ,10
+ ,59506
+ ,833084
+ ,319844
+ ,1
+ ,382
+ ,21
+ ,9
+ ,59365
+ ,1246665
+ ,355822
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,9
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,30
+ ,4
+ ,4
+ ,60798
+ ,243192
+ ,324047
+ ,0
+ ,290
+ ,2
+ ,1
+ ,58790
+ ,117580
+ ,311464
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,366
+ ,53
+ ,14
+ ,61808
+ ,3275824
+ ,353417
+ ,0
+ ,2
+ ,9
+ ,12
+ ,60735
+ ,546615
+ ,325590
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,209
+ ,13
+ ,19
+ ,64016
+ ,832208
+ ,328576
+ ,0
+ ,384
+ ,22
+ ,17
+ ,54683
+ ,1203026
+ ,326126
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,365
+ ,83
+ ,32
+ ,87192
+ ,7236936
+ ,369376
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,49
+ ,8
+ ,14
+ ,64107
+ ,512856
+ ,332013
+ ,0
+ ,3
+ ,4
+ ,8
+ ,60761
+ ,243044
+ ,325871
+ ,0
+ ,133
+ ,14
+ ,4
+ ,65990
+ ,923860
+ ,342165
+ ,0
+ ,32
+ ,1
+ ,0
+ ,59988
+ ,59988
+ ,324967
+ ,0
+ ,368
+ ,17
+ ,20
+ ,61167
+ ,1039839
+ ,314832
+ ,0
+ ,1
+ ,6
+ ,5
+ ,60719
+ ,364314
+ ,325557
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,22
+ ,2
+ ,1
+ ,60722
+ ,121444
+ ,322649
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,96
+ ,5
+ ,4
+ ,60379
+ ,301895
+ ,324598
+ ,0
+ ,1
+ ,2
+ ,1
+ ,60727
+ ,121454
+ ,325567
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,81
+ ,7
+ ,4
+ ,60925
+ ,426475
+ ,324005
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,26
+ ,1
+ ,1
+ ,60896
+ ,60896
+ ,325748
+ ,0
+ ,125
+ ,13
+ ,10
+ ,59734
+ ,776542
+ ,323385
+ ,0
+ ,304
+ ,15
+ ,12
+ ,62969
+ ,944535
+ ,315409
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,119
+ ,6
+ ,3
+ ,59118
+ ,354708
+ ,312275
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,312
+ ,14
+ ,3
+ ,58598
+ ,820372
+ ,320576
+ ,0
+ ,60
+ ,10
+ ,7
+ ,61124
+ ,611240
+ ,325246
+ ,0
+ ,587
+ ,12
+ ,10
+ ,59595
+ ,715140
+ ,332961
+ ,0
+ ,135
+ ,2
+ ,1
+ ,62065
+ ,124130
+ ,323010
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,514
+ ,52
+ ,15
+ ,78780
+ ,4096560
+ ,345253
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,1
+ ,4
+ ,4
+ ,60722
+ ,242888
+ ,325559
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,58
+ ,3
+ ,2
+ ,61600
+ ,184800
+ ,319634
+ ,0
+ ,180
+ ,11
+ ,9
+ ,59635
+ ,655985
+ ,319951
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,448
+ ,40
+ ,7
+ ,59781
+ ,2391240
+ ,318519
+ ,0
+ ,227
+ ,9
+ ,7
+ ,76644
+ ,689796
+ ,343222
+ ,0
+ ,174
+ ,1
+ ,3
+ ,64820
+ ,64820
+ ,317234
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,121
+ ,24
+ ,11
+ ,56178
+ ,1348272
+ ,314025
+ ,0
+ ,607
+ ,11
+ ,7
+ ,60436
+ ,664796
+ ,320249
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,530
+ ,60
+ ,18
+ ,73433
+ ,4405980
+ ,349365
+ ,0
+ ,571
+ ,80
+ ,14
+ ,41477
+ ,3318160
+ ,289197
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,78
+ ,16
+ ,12
+ ,62700
+ ,1003200
+ ,329245
+ ,0
+ ,2489
+ ,40
+ ,29
+ ,67804
+ ,2712160
+ ,240869
+ ,0
+ ,131
+ ,6
+ ,3
+ ,59661
+ ,357966
+ ,327182
+ ,0
+ ,923
+ ,8
+ ,6
+ ,58620
+ ,468960
+ ,322876
+ ,0
+ ,72
+ ,3
+ ,3
+ ,60398
+ ,181194
+ ,323117
+ ,0
+ ,572
+ ,16
+ ,8
+ ,58580
+ ,937280
+ ,306351
+ ,0
+ ,397
+ ,10
+ ,10
+ ,62710
+ ,627100
+ ,335137
+ ,0
+ ,450
+ ,8
+ ,6
+ ,59325
+ ,474600
+ ,308271
+ ,0
+ ,622
+ ,7
+ ,8
+ ,60950
+ ,426650
+ ,301731
+ ,0
+ ,694
+ ,8
+ ,6
+ ,68060
+ ,544480
+ ,382409
+ ,1
+ ,3425
+ ,12
+ ,9
+ ,83620
+ ,1003440
+ ,279230
+ ,0
+ ,562
+ ,13
+ ,8
+ ,58456
+ ,759928
+ ,298731
+ ,0
+ ,4917
+ ,42
+ ,26
+ ,52811
+ ,2218062
+ ,243650
+ ,1
+ ,1442
+ ,118
+ ,239
+ ,121173
+ ,14298414
+ ,532682
+ ,0
+ ,529
+ ,9
+ ,7
+ ,63870
+ ,574830
+ ,319771
+ ,1
+ ,2126
+ ,138
+ ,41
+ ,21001
+ ,2898138
+ ,171493
+ ,0
+ ,1061
+ ,5
+ ,3
+ ,70415
+ ,352075
+ ,347262
+ ,0
+ ,776
+ ,9
+ ,8
+ ,64230
+ ,578070
+ ,343945
+ ,0
+ ,611
+ ,8
+ ,6
+ ,59190
+ ,473520
+ ,311874
+ ,1
+ ,1526
+ ,25
+ ,21
+ ,69351
+ ,1733775
+ ,302211
+ ,0
+ ,592
+ ,7
+ ,7
+ ,64270
+ ,449890
+ ,316708
+ ,0
+ ,1182
+ ,13
+ ,11
+ ,70694
+ ,919022
+ ,333463
+ ,0
+ ,621
+ ,16
+ ,11
+ ,68005
+ ,1088080
+ ,344282
+ ,0
+ ,989
+ ,11
+ ,12
+ ,58930
+ ,648230
+ ,319635
+ ,0
+ ,438
+ ,11
+ ,9
+ ,58320
+ ,641520
+ ,301186
+ ,0
+ ,726
+ ,3
+ ,3
+ ,69980
+ ,209940
+ ,300381
+ ,0
+ ,1303
+ ,61
+ ,57
+ ,69863
+ ,4261643
+ ,318765
+ ,1
+ ,6341
+ ,24
+ ,16
+ ,63255
+ ,1518120
+ ,286146
+ ,1
+ ,1164
+ ,17
+ ,15
+ ,57320
+ ,974440
+ ,306844
+ ,1
+ ,3310
+ ,33
+ ,32
+ ,75230
+ ,2482590
+ ,307705
+ ,0
+ ,1366
+ ,7
+ ,5
+ ,79420
+ ,555940
+ ,312448
+ ,0
+ ,965
+ ,3
+ ,2
+ ,73490
+ ,220470
+ ,299715
+ ,0
+ ,3256
+ ,66
+ ,23
+ ,35250
+ ,2326500
+ ,373399
+ ,1
+ ,1135
+ ,17
+ ,20
+ ,62285
+ ,1058845
+ ,299446
+ ,0
+ ,1270
+ ,26
+ ,24
+ ,69206
+ ,1799356
+ ,325586
+ ,0
+ ,661
+ ,3
+ ,1
+ ,65920
+ ,197760
+ ,291221
+ ,0
+ ,1013
+ ,2
+ ,1
+ ,69770
+ ,139540
+ ,261173
+ ,0
+ ,2844
+ ,67
+ ,74
+ ,72683
+ ,4869761
+ ,255027
+ ,1
+ ,11528
+ ,70
+ ,68
+ ,-14545
+ ,-1018150
+ ,-78375
+ ,0
+ ,6526
+ ,26
+ ,20
+ ,55830
+ ,1451580
+ ,-58143
+ ,0
+ ,2264
+ ,24
+ ,20
+ ,55174
+ ,1324176
+ ,227033
+ ,1
+ ,4461
+ ,94
+ ,79
+ ,67038
+ ,6301572
+ ,235098
+ ,0
+ ,3999
+ ,30
+ ,21
+ ,51252
+ ,1537560
+ ,21267
+ ,0
+ ,35624
+ ,223
+ ,244
+ ,157278
+ ,35072994
+ ,238675
+ ,0
+ ,9252
+ ,48
+ ,32
+ ,79510
+ ,3816480
+ ,197687
+ ,0
+ ,15236
+ ,90
+ ,86
+ ,77440
+ ,6969600
+ ,418341
+ ,0
+ ,18073
+ ,180
+ ,69
+ ,27284
+ ,4911120
+ ,-297706)
+ ,dim=c(7
+ ,431)
+ ,dimnames=list(c('Group'
+ ,'Costs'
+ ,'Trades'
+ ,'Orders'
+ ,'Dividends'
+ ,'TrDiv'
+ ,'Wealth')
+ ,1:431))
> y <- array(NA,dim=c(7,431),dimnames=list(c('Group','Costs','Trades','Orders','Dividends','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 = '7'
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Prof. Dr. P. Wessa
> #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #Technical description: Write here your technical program description (don't use hard returns!)
> library(lattice)
> library(lmtest)
Loading required package: zoo
Attaching package: 'zoo'
The following object(s) are masked from package:base :
as.Date.numeric
> n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test
> par1 <- as.numeric(par1)
> x <- t(y)
> k <- length(x[1,])
> n <- length(x[,1])
> x1 <- cbind(x[,par1], x[,1:k!=par1])
> mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1])
> colnames(x1) <- mycolnames #colnames(x)[par1]
> x <- x1
> if (par3 == 'First Differences'){
+ x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep='')))
+ for (i in 1:n-1) {
+ for (j in 1:k) {
+ x2[i,j] <- x[i+1,j] - x[i,j]
+ }
+ }
+ x <- x2
+ }
> if (par2 == 'Include Monthly Dummies'){
+ x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep ='')))
+ for (i in 1:11){
+ x2[seq(i,n,12),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> if (par2 == 'Include Quarterly Dummies'){
+ x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep ='')))
+ for (i in 1:3){
+ x2[seq(i,n,4),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> k <- length(x[1,])
> if (par3 == 'Linear Trend'){
+ x <- cbind(x, c(1:n))
+ colnames(x)[k+1] <- 't'
+ }
> x
Wealth Group Costs Trades Orders Dividends TrDiv
1 6282929 1 162556 1081 807 213118 230380558
2 4324047 1 29790 309 444 81767 25266003
3 4108272 1 87550 458 412 153198 70164684
4 -1212617 0 84738 588 428 -26007 -15292116
5 1485329 1 54660 299 312 126942 37955658
6 1779876 1 42634 156 166 157214 24525384
7 1367203 0 40949 481 263 129352 62218312
8 2519076 1 42312 323 237 234817 75845891
9 912684 1 37704 452 228 60448 27322496
10 1443586 1 16275 109 129 47818 5212162
11 1220017 0 25830 115 104 245546 28237790
12 984885 0 12679 110 122 48020 5282200
13 1457425 1 18014 239 393 -1710 -408690
14 -572920 0 43556 247 190 32648 8064056
15 929144 1 24524 497 275 95350 47388950
16 1151176 0 6532 103 62 151352 15589256
17 790090 0 7123 109 102 288170 31410530
18 774497 1 20813 502 255 114337 57397174
19 990576 1 37597 248 234 37884 9395232
20 454195 0 17821 373 277 122844 45820812
21 876607 1 12988 119 73 82340 9798460
22 711969 1 22330 84 67 79801 6703284
23 702380 0 13326 102 103 165548 16885896
24 264449 0 16189 295 290 116384 34333280
25 450033 0 7146 105 83 134028 14072940
26 541063 0 15824 64 56 63838 4085632
27 588864 1 26088 267 224 74996 20023932
28 -37216 0 11326 129 64 31080 4009320
29 783310 0 8568 37 34 32168 1190216
30 467359 0 14416 361 139 49857 17998377
31 688779 1 3369 28 26 87161 2440508
32 608419 1 11819 85 70 106113 9019605
33 696348 1 6620 44 40 80570 3545080
34 597793 1 4519 49 42 102129 5004321
35 821730 0 2220 22 12 301670 6636740
36 377934 0 18562 155 211 102313 15858515
37 651939 0 10327 91 74 88577 8060507
38 697458 1 5336 81 80 112477 9110637
39 700368 1 2365 79 83 191778 15150462
40 225986 0 4069 145 131 79804 11571580
41 348695 0 7710 816 187 128294 104687904
42 373683 0 13718 61 56 96448 5883328
43 501709 0 4525 226 89 93811 21201286
44 413743 0 6869 105 88 117520 12339600
45 379825 0 4628 62 39 69159 4287858
46 336260 1 3653 24 24 101792 2443008
47 636765 1 1265 26 19 210568 5474768
48 481231 1 7489 322 149 136996 44112712
49 469107 0 4901 84 58 121920 10241280
50 211928 0 2284 33 41 76403 2521299
51 563925 1 3160 108 90 108094 11674152
52 511939 1 4150 150 136 134759 20213850
53 521016 1 7285 115 97 188873 21720395
54 543856 1 1134 162 63 146216 23686992
55 329304 1 4658 158 114 156608 24744064
56 423262 0 2384 97 77 61348 5950756
57 509665 0 3748 9 6 50350 453150
58 455881 0 5371 66 47 87720 5789520
59 367772 0 1285 107 51 99489 10645323
60 406339 1 9327 101 85 87419 8829319
61 493408 1 5565 47 43 94355 4434685
62 232942 0 1528 38 32 60326 2292388
63 416002 1 3122 34 25 94670 3218780
64 337430 1 7317 84 74 82425 6923700
65 361517 0 2675 79 54 59017 4662343
66 360962 0 13253 947 251 90829 86015063
67 235561 0 880 74 15 80791 5978534
68 408247 1 2053 53 44 100423 5322419
69 450296 0 1424 94 73 131116 12324904
70 418799 1 4036 63 85 100269 6316947
71 247405 1 3045 58 49 27330 1585140
72 378519 0 5119 49 38 39039 1912911
73 326638 0 1431 34 35 106885 3634090
74 328233 0 554 11 9 79285 872135
75 386225 0 1975 35 34 118881 4160835
76 283662 1 1286 17 17 77623 1319591
77 370225 0 1012 47 29 114768 5394096
78 269236 0 810 43 11 74015 3182645
79 365732 0 1280 117 52 69465 8127405
80 420383 1 666 171 13 117869 20155599
81 345811 0 1380 26 29 60982 1585532
82 431809 1 4608 73 64 90131 6579563
83 418876 0 876 59 33 138971 8199289
84 297476 0 814 18 15 39625 713250
85 416776 0 514 15 15 102725 1540875
86 357257 1 5692 72 68 64239 4625208
87 458343 0 3642 86 100 90262 7762532
88 388386 0 540 14 13 103960 1455440
89 358934 0 2099 64 45 106611 6823104
90 407560 0 567 11 14 103345 1136795
91 392558 0 2001 52 36 95551 4968652
92 373177 1 2949 41 40 82903 3399023
93 428370 0 2253 99 68 63593 6295707
94 369419 1 6533 75 29 126910 9518250
95 358649 0 1889 45 43 37527 1688715
96 376641 1 3055 43 30 60247 2590621
97 467427 0 272 8 9 112995 903960
98 364885 1 1414 198 22 70184 13896432
99 436230 0 2564 22 19 130140 2863080
100 329118 1 1383 11 9 73221 805431
101 317365 1 1261 33 31 76114 2511762
102 286849 0 975 23 19 90534 2082282
103 376685 0 3366 80 55 108479 8678320
104 407198 0 576 18 8 113761 2047698
105 377772 0 1306 28 26 68696 1923488
106 271483 0 746 23 29 71561 1645903
107 153661 1 3192 60 48 59831 3589860
108 513294 1 2045 20 16 97890 1957800
109 324881 0 5477 59 45 101481 5987379
110 264512 1 1932 36 20 72954 2626344
111 420968 0 936 30 22 67939 2038170
112 129302 1 3437 47 33 48022 2257034
113 191521 0 5131 71 44 86111 6113881
114 268673 1 2397 14 13 74020 1036280
115 353179 1 1389 9 6 57530 517770
116 354624 0 1503 39 35 56364 2198196
117 363713 0 402 26 8 84990 2209740
118 456657 0 2239 21 17 88590 1860390
119 211742 1 2234 16 11 77200 1235200
120 338381 0 837 69 21 61262 4227078
121 418530 0 10579 92 92 110309 10148428
122 351483 0 875 14 12 67000 938000
123 372928 0 1395 103 108 93099 9589197
124 485538 1 1659 29 25 107577 3119733
125 279268 1 2647 37 17 62920 2328040
126 219060 1 3294 23 23 75832 1744136
127 325560 0 0 0 0 60720 0
128 325314 0 94 7 10 60793 425551
129 322046 0 422 28 23 57935 1622180
130 325560 0 0 0 0 60720 0
131 325599 0 34 8 7 60630 485040
132 377028 0 1558 63 25 55637 3505131
133 325560 0 0 0 0 60720 0
134 323850 0 43 3 20 60887 182661
135 325560 0 0 0 0 60720 0
136 331514 0 316 9 4 60505 544545
137 325632 0 115 13 10 60945 792285
138 325560 0 0 0 0 60720 0
139 325560 0 0 0 0 60720 0
140 325560 0 0 0 0 60720 0
141 322265 0 389 14 7 58990 825860
142 325560 0 0 0 0 60720 0
143 325906 0 1002 15 11 56750 851250
144 325985 0 36 3 4 60894 182682
145 346145 0 460 15 15 63346 950190
146 325898 0 309 11 9 56535 621885
147 325560 0 0 0 0 60720 0
148 325356 0 9 6 7 60835 365010
149 325560 0 0 0 0 60720 0
150 325930 0 14 1 0 61016 61016
151 318020 0 520 10 7 58650 586500
152 326389 0 1766 73 46 60438 4411974
153 325560 0 0 0 5 60720 0
154 302925 0 458 11 7 58625 644875
155 325540 0 20 3 2 60938 182814
156 325560 0 0 0 0 60720 0
157 325560 0 0 0 0 60720 0
158 326736 0 98 2 2 61490 122980
159 340580 0 405 7 5 60845 425915
160 325560 0 0 0 0 60720 0
161 325560 0 0 0 0 60720 0
162 325560 0 0 0 0 60720 0
163 325560 0 0 0 0 60720 0
164 331828 0 483 27 7 60830 1642410
165 323299 0 454 51 24 63261 3226311
166 325560 0 0 0 0 60720 0
167 325560 0 0 0 0 60720 0
168 387722 0 757 19 18 45689 868091
169 325560 0 0 0 4 60720 0
170 325560 0 0 0 0 60720 0
171 325560 0 0 0 0 60720 0
172 324598 0 36 4 3 61564 246256
173 325560 0 0 0 0 60720 0
174 328726 0 203 9 9 61938 557442
175 325560 0 0 0 0 60720 0
176 325043 0 90 8 6 60951 487608
177 325560 0 0 0 0 60720 0
178 325806 1 71 1 0 60745 60745
179 325560 0 0 0 0 60720 0
180 325560 0 0 0 0 60720 0
181 387732 0 972 34 19 71642 2435828
182 349729 1 531 10 11 71641 716410
183 332202 0 604 38 8 55792 2120096
184 305442 1 283 10 15 71873 718730
185 329537 1 23 5 5 62555 312775
186 327055 1 638 14 11 60370 845180
187 356245 1 699 16 10 64873 1037968
188 328451 0 149 5 6 62041 310205
189 307062 0 226 5 5 65745 328725
190 325560 0 0 0 0 60720 0
191 331345 0 275 4 7 59500 238000
192 325560 0 0 0 0 60720 0
193 331824 0 141 6 7 61630 369780
194 325560 0 0 0 0 60720 0
195 325685 0 28 2 3 60890 121780
196 325560 1 0 0 0 60720 0
197 404480 1 2566 80 74 113521 9081680
198 325560 1 0 0 0 60720 0
199 325560 1 0 0 0 60720 0
200 318314 1 472 20 19 80045 1600900
201 325560 1 0 0 0 60720 0
202 325560 1 0 0 0 60720 0
203 325560 1 0 0 0 60720 0
204 311807 1 203 27 12 50804 1371708
205 337724 1 496 17 12 87390 1485630
206 326431 1 10 2 5 61656 123312
207 327556 1 63 4 2 65688 262752
208 325560 1 0 0 0 60720 0
209 356850 1 1136 32 26 48522 1552704
210 325560 1 0 0 0 60720 0
211 325560 0 0 0 0 60720 0
212 325560 0 0 0 0 60720 0
213 322741 0 267 32 11 57640 1844480
214 310902 0 474 20 10 61977 1239540
215 324295 0 534 7 5 62620 438340
216 325560 0 0 0 2 60720 0
217 326156 0 15 8 6 60831 486648
218 326960 0 397 28 7 60646 1698088
219 325560 0 0 0 2 60720 0
220 333411 0 1061 20 22 56225 1124500
221 297761 0 288 4 3 60510 242040
222 325560 0 0 0 0 60720 0
223 325536 0 3 2 1 60698 121396
224 325560 0 0 0 6 60720 0
225 325762 0 20 2 1 60805 121610
226 327957 0 278 26 22 61404 1596504
227 325560 0 0 0 0 60720 0
228 325560 0 0 0 0 60720 0
229 318521 0 192 4 2 65276 261104
230 325560 0 0 0 0 60720 0
231 319775 0 317 9 7 63915 575235
232 325560 0 0 0 0 60720 0
233 325560 0 0 0 0 60720 0
234 332128 0 368 17 13 61686 1048662
235 325560 0 0 0 0 60720 0
236 325486 0 2 1 0 60743 60743
237 325560 0 0 0 0 60720 0
238 325838 0 53 6 6 60349 362094
239 325560 0 0 0 0 60720 0
240 325560 0 0 0 0 60720 0
241 325560 0 0 0 0 60720 0
242 331767 0 94 3 3 61360 184080
243 325560 0 0 0 0 60720 0
244 324523 0 24 8 7 59818 478544
245 339995 0 2332 4 2 72680 290720
246 325560 1 0 0 0 60720 0
247 325560 1 0 0 0 60720 0
248 319582 0 131 11 15 61808 679888
249 325560 1 0 0 0 60720 0
250 325560 1 0 0 0 60720 0
251 307245 0 206 9 9 53110 477990
252 325560 1 0 0 0 60720 0
253 317967 0 167 2 1 64245 128490
254 331488 0 622 73 38 73007 5329511
255 335452 0 885 85 49 82732 7032220
256 325560 0 0 0 0 60720 0
257 334184 0 365 8 6 54820 438560
258 313213 0 364 35 26 47705 1669675
259 325560 1 0 0 0 60720 0
260 325560 0 0 0 0 60720 0
261 325560 1 0 0 0 60720 0
262 325560 1 0 0 0 60720 0
263 348678 0 226 12 13 72835 874020
264 328727 0 307 15 10 58856 882840
265 325560 1 0 0 0 60720 0
266 325560 0 0 0 0 60720 0
267 325560 1 0 0 0 60720 0
268 387978 0 188 11 9 77655 854205
269 325560 1 0 0 0 60720 0
270 336704 0 138 6 26 69817 418902
271 325560 1 0 0 0 60720 0
272 325560 1 0 0 0 60720 0
273 325560 1 0 0 0 60720 0
274 322076 0 125 12 19 60798 729576
275 325560 0 0 0 0 60720 0
276 334272 0 282 30 12 62452 1873560
277 338197 0 335 33 23 64175 2117775
278 325560 0 0 0 0 60720 0
279 321024 1 813 82 24 67440 5530080
280 322145 0 176 28 8 68136 1907808
281 325560 1 0 0 0 60720 0
282 325560 0 0 0 0 60720 0
283 323351 0 249 72 26 56726 4084272
284 325560 0 0 0 0 60720 0
285 327748 0 333 13 9 70811 920543
286 325560 0 0 0 0 60720 0
287 325560 1 0 0 0 60720 0
288 328157 0 30 4 3 62045 248180
289 325560 1 0 0 0 60720 0
290 311594 0 249 62 13 54323 3368026
291 325560 1 0 0 0 60720 0
292 335962 0 165 24 12 62841 1508184
293 372426 0 453 21 19 81125 1703625
294 325560 0 0 0 0 60720 0
295 319844 0 53 14 10 59506 833084
296 355822 1 382 21 9 59365 1246665
297 325560 0 0 0 0 60720 0
298 325560 1 0 0 0 60720 0
299 325560 0 0 0 0 60720 0
300 325560 0 0 0 9 60720 0
301 324047 1 30 4 4 60798 243192
302 311464 0 290 2 1 58790 117580
303 325560 0 0 0 0 60720 0
304 325560 1 0 0 0 60720 0
305 353417 0 366 53 14 61808 3275824
306 325590 0 2 9 12 60735 546615
307 325560 1 0 0 0 60720 0
308 328576 1 209 13 19 64016 832208
309 326126 0 384 22 17 54683 1203026
310 325560 1 0 0 0 60720 0
311 325560 1 0 0 0 60720 0
312 369376 0 365 83 32 87192 7236936
313 325560 1 0 0 0 60720 0
314 332013 1 49 8 14 64107 512856
315 325871 0 3 4 8 60761 243044
316 342165 0 133 14 4 65990 923860
317 324967 0 32 1 0 59988 59988
318 314832 0 368 17 20 61167 1039839
319 325557 0 1 6 5 60719 364314
320 325560 0 0 0 0 60720 0
321 325560 0 0 0 0 60720 0
322 325560 0 0 0 0 60720 0
323 325560 0 0 0 0 60720 0
324 325560 1 0 0 0 60720 0
325 325560 0 0 0 0 60720 0
326 322649 0 22 2 1 60722 121444
327 325560 1 0 0 0 60720 0
328 325560 1 0 0 0 60720 0
329 325560 0 0 0 0 60720 0
330 325560 0 0 0 0 60720 0
331 325560 0 0 0 0 60720 0
332 325560 0 0 0 0 60720 0
333 325560 0 0 0 0 60720 0
334 324598 0 96 5 4 60379 301895
335 325567 0 1 2 1 60727 121454
336 325560 0 0 0 0 60720 0
337 324005 0 81 7 4 60925 426475
338 325560 0 0 0 0 60720 0
339 325748 0 26 1 1 60896 60896
340 323385 0 125 13 10 59734 776542
341 315409 0 304 15 12 62969 944535
342 325560 0 0 0 0 60720 0
343 325560 0 0 0 0 60720 0
344 325560 0 0 0 0 60720 0
345 325560 0 0 0 0 60720 0
346 325560 0 0 0 0 60720 0
347 312275 0 119 6 3 59118 354708
348 325560 0 0 0 0 60720 0
349 325560 0 0 0 0 60720 0
350 325560 0 0 0 0 60720 0
351 320576 0 312 14 3 58598 820372
352 325246 0 60 10 7 61124 611240
353 332961 0 587 12 10 59595 715140
354 323010 0 135 2 1 62065 124130
355 325560 0 0 0 0 60720 0
356 325560 0 0 0 0 60720 0
357 345253 0 514 52 15 78780 4096560
358 325560 0 0 0 0 60720 0
359 325560 0 0 0 0 60720 0
360 325560 0 0 0 0 60720 0
361 325559 0 1 4 4 60722 242888
362 325560 0 0 0 0 60720 0
363 325560 0 0 0 0 60720 0
364 319634 1 58 3 2 61600 184800
365 319951 0 180 11 9 59635 655985
366 325560 0 0 0 0 60720 0
367 325560 0 0 0 0 60720 0
368 325560 0 0 0 0 60720 0
369 325560 0 0 0 0 60720 0
370 325560 0 0 0 0 60720 0
371 325560 0 0 0 0 60720 0
372 318519 0 448 40 7 59781 2391240
373 343222 0 227 9 7 76644 689796
374 317234 0 174 1 3 64820 64820
375 325560 0 0 0 0 60720 0
376 325560 0 0 0 0 60720 0
377 314025 0 121 24 11 56178 1348272
378 320249 0 607 11 7 60436 664796
379 325560 0 0 0 0 60720 0
380 325560 0 0 0 0 60720 0
381 325560 0 0 0 0 60720 0
382 349365 0 530 60 18 73433 4405980
383 289197 0 571 80 14 41477 3318160
384 325560 0 0 0 0 60720 0
385 329245 0 78 16 12 62700 1003200
386 240869 0 2489 40 29 67804 2712160
387 327182 0 131 6 3 59661 357966
388 322876 0 923 8 6 58620 468960
389 323117 0 72 3 3 60398 181194
390 306351 0 572 16 8 58580 937280
391 335137 0 397 10 10 62710 627100
392 308271 0 450 8 6 59325 474600
393 301731 0 622 7 8 60950 426650
394 382409 0 694 8 6 68060 544480
395 279230 1 3425 12 9 83620 1003440
396 298731 0 562 13 8 58456 759928
397 243650 0 4917 42 26 52811 2218062
398 532682 1 1442 118 239 121173 14298414
399 319771 0 529 9 7 63870 574830
400 171493 1 2126 138 41 21001 2898138
401 347262 0 1061 5 3 70415 352075
402 343945 0 776 9 8 64230 578070
403 311874 0 611 8 6 59190 473520
404 302211 1 1526 25 21 69351 1733775
405 316708 0 592 7 7 64270 449890
406 333463 0 1182 13 11 70694 919022
407 344282 0 621 16 11 68005 1088080
408 319635 0 989 11 12 58930 648230
409 301186 0 438 11 9 58320 641520
410 300381 0 726 3 3 69980 209940
411 318765 0 1303 61 57 69863 4261643
412 286146 1 6341 24 16 63255 1518120
413 306844 1 1164 17 15 57320 974440
414 307705 1 3310 33 32 75230 2482590
415 312448 0 1366 7 5 79420 555940
416 299715 0 965 3 2 73490 220470
417 373399 0 3256 66 23 35250 2326500
418 299446 1 1135 17 20 62285 1058845
419 325586 0 1270 26 24 69206 1799356
420 291221 0 661 3 1 65920 197760
421 261173 0 1013 2 1 69770 139540
422 255027 0 2844 67 74 72683 4869761
423 -78375 1 11528 70 68 -14545 -1018150
424 -58143 0 6526 26 20 55830 1451580
425 227033 0 2264 24 20 55174 1324176
426 235098 1 4461 94 79 67038 6301572
427 21267 0 3999 30 21 51252 1537560
428 238675 0 35624 223 244 157278 35072994
429 197687 0 9252 48 32 79510 3816480
430 418341 0 15236 90 86 77440 6969600
431 -297706 0 18073 180 69 27284 4911120
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) Group Costs Trades Orders Dividends
3.354e+05 5.511e+04 5.927e+00 -4.315e+03 4.170e+03 -5.567e-01
TrDiv
2.875e-02
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-1284038 -31164 18420 23944 2557838
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 3.354e+05 3.252e+04 10.314 < 2e-16 ***
Group 5.511e+04 2.528e+04 2.180 0.02979 *
Costs 5.927e+00 2.091e+00 2.834 0.00482 **
Trades -4.315e+03 2.752e+02 -15.680 < 2e-16 ***
Orders 4.170e+03 4.287e+02 9.726 < 2e-16 ***
Dividends -5.567e-01 4.374e-01 -1.273 0.20380
TrDiv 2.875e-02 1.767e-03 16.273 < 2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 228400 on 424 degrees of freedom
Multiple R-squared: 0.7431, Adjusted R-squared: 0.7395
F-statistic: 204.4 on 6 and 424 DF, p-value: < 2.2e-16
> if (n > n25) {
+ kp3 <- k + 3
+ nmkm3 <- n - k - 3
+ gqarr <- array(NA, dim=c(nmkm3-kp3+1,3))
+ numgqtests <- 0
+ numsignificant1 <- 0
+ numsignificant5 <- 0
+ numsignificant10 <- 0
+ for (mypoint in kp3:nmkm3) {
+ j <- 0
+ numgqtests <- numgqtests + 1
+ for (myalt in c('greater', 'two.sided', 'less')) {
+ j <- j + 1
+ gqarr[mypoint-kp3+1,j] <- gqtest(mylm, point=mypoint, alternative=myalt)$p.value
+ }
+ if (gqarr[mypoint-kp3+1,2] < 0.01) numsignificant1 <- numsignificant1 + 1
+ if (gqarr[mypoint-kp3+1,2] < 0.05) numsignificant5 <- numsignificant5 + 1
+ if (gqarr[mypoint-kp3+1,2] < 0.10) numsignificant10 <- numsignificant10 + 1
+ }
+ gqarr
+ }
[,1] [,2] [,3]
[1,] 1.0000000 6.002798e-36 3.001399e-36
[2,] 1.0000000 2.824924e-37 1.412462e-37
[3,] 1.0000000 1.595201e-38 7.976005e-39
[4,] 1.0000000 1.467970e-62 7.339852e-63
[5,] 1.0000000 2.653330e-76 1.326665e-76
[6,] 1.0000000 1.195983e-80 5.979917e-81
[7,] 1.0000000 2.246340e-88 1.123170e-88
[8,] 1.0000000 5.632362e-94 2.816181e-94
[9,] 1.0000000 3.075679e-97 1.537839e-97
[10,] 1.0000000 5.127542e-104 2.563771e-104
[11,] 1.0000000 6.942451e-108 3.471226e-108
[12,] 1.0000000 3.627250e-113 1.813625e-113
[13,] 1.0000000 8.517863e-117 4.258931e-117
[14,] 1.0000000 5.687746e-117 2.843873e-117
[15,] 1.0000000 1.112129e-125 5.560644e-126
[16,] 1.0000000 9.704956e-125 4.852478e-125
[17,] 1.0000000 4.084167e-127 2.042084e-127
[18,] 1.0000000 2.740490e-133 1.370245e-133
[19,] 1.0000000 1.323807e-134 6.619034e-135
[20,] 1.0000000 9.832949e-149 4.916474e-149
[21,] 1.0000000 2.362293e-159 1.181146e-159
[22,] 1.0000000 7.670821e-164 3.835410e-164
[23,] 1.0000000 7.720804e-168 3.860402e-168
[24,] 1.0000000 3.601832e-175 1.800916e-175
[25,] 1.0000000 1.510464e-177 7.552322e-178
[26,] 1.0000000 8.266993e-178 4.133496e-178
[27,] 1.0000000 1.566351e-180 7.831755e-181
[28,] 1.0000000 2.428257e-188 1.214128e-188
[29,] 1.0000000 3.309918e-194 1.654959e-194
[30,] 1.0000000 1.556186e-195 7.780932e-196
[31,] 1.0000000 1.292405e-196 6.462024e-197
[32,] 1.0000000 7.452755e-199 3.726378e-199
[33,] 1.0000000 2.325587e-199 1.162793e-199
[34,] 1.0000000 2.289329e-200 1.144665e-200
[35,] 1.0000000 2.268201e-199 1.134101e-199
[36,] 1.0000000 4.373332e-199 2.186666e-199
[37,] 1.0000000 3.439152e-199 1.719576e-199
[38,] 1.0000000 1.426529e-198 7.132645e-199
[39,] 1.0000000 9.397522e-200 4.698761e-200
[40,] 1.0000000 4.364562e-199 2.182281e-199
[41,] 1.0000000 1.730090e-199 8.650450e-200
[42,] 1.0000000 9.464975e-201 4.732487e-201
[43,] 1.0000000 6.352046e-202 3.176023e-202
[44,] 1.0000000 7.366366e-203 3.683183e-203
[45,] 1.0000000 6.049162e-202 3.024581e-202
[46,] 1.0000000 2.334303e-208 1.167151e-208
[47,] 1.0000000 2.294666e-208 1.147333e-208
[48,] 1.0000000 7.090789e-214 3.545395e-214
[49,] 1.0000000 6.656088e-215 3.328044e-215
[50,] 1.0000000 4.429154e-214 2.214577e-214
[51,] 1.0000000 2.545440e-215 1.272720e-215
[52,] 1.0000000 1.111902e-217 5.559511e-218
[53,] 1.0000000 2.703935e-217 1.351968e-217
[54,] 1.0000000 5.974204e-217 2.987102e-217
[55,] 1.0000000 7.994260e-217 3.997130e-217
[56,] 1.0000000 2.594442e-216 1.297221e-216
[57,] 1.0000000 1.454873e-219 7.274365e-220
[58,] 1.0000000 6.123754e-220 3.061877e-220
[59,] 1.0000000 4.614312e-219 2.307156e-219
[60,] 1.0000000 2.513912e-218 1.256956e-218
[61,] 1.0000000 5.599035e-218 2.799518e-218
[62,] 1.0000000 3.015624e-217 1.507812e-217
[63,] 1.0000000 1.260718e-219 6.303592e-220
[64,] 1.0000000 3.505790e-219 1.752895e-219
[65,] 1.0000000 3.397266e-218 1.698633e-218
[66,] 1.0000000 3.137009e-217 1.568505e-217
[67,] 1.0000000 1.194170e-216 5.970851e-217
[68,] 1.0000000 4.991704e-216 2.495852e-216
[69,] 1.0000000 1.469281e-215 7.346403e-216
[70,] 1.0000000 9.691553e-215 4.845776e-215
[71,] 1.0000000 7.731077e-215 3.865538e-215
[72,] 1.0000000 6.471011e-214 3.235505e-214
[73,] 1.0000000 4.484858e-214 2.242429e-214
[74,] 1.0000000 4.768882e-214 2.384441e-214
[75,] 1.0000000 4.877578e-213 2.438789e-213
[76,] 1.0000000 3.971417e-212 1.985708e-212
[77,] 1.0000000 8.543522e-213 4.271761e-213
[78,] 1.0000000 9.944656e-213 4.972328e-213
[79,] 1.0000000 9.597108e-212 4.798554e-212
[80,] 1.0000000 3.929530e-211 1.964765e-211
[81,] 1.0000000 3.507644e-210 1.753822e-210
[82,] 1.0000000 3.620712e-209 1.810356e-209
[83,] 1.0000000 1.668600e-208 8.343002e-209
[84,] 1.0000000 1.980987e-209 9.904936e-210
[85,] 1.0000000 1.479132e-208 7.395658e-209
[86,] 1.0000000 2.510024e-209 1.255012e-209
[87,] 1.0000000 7.771629e-210 3.885815e-210
[88,] 1.0000000 2.226946e-209 1.113473e-209
[89,] 1.0000000 7.982705e-209 3.991352e-209
[90,] 1.0000000 6.597699e-208 3.298850e-208
[91,] 1.0000000 6.318079e-207 3.159040e-207
[92,] 1.0000000 4.705927e-206 2.352964e-206
[93,] 1.0000000 6.327862e-206 3.163931e-206
[94,] 1.0000000 4.403444e-205 2.201722e-205
[95,] 1.0000000 3.590693e-204 1.795346e-204
[96,] 1.0000000 1.417922e-203 7.089608e-204
[97,] 1.0000000 3.653709e-203 1.826854e-203
[98,] 1.0000000 2.092618e-204 1.046309e-204
[99,] 1.0000000 6.477906e-207 3.238953e-207
[100,] 1.0000000 6.875801e-206 3.437901e-206
[101,] 1.0000000 2.510838e-205 1.255419e-205
[102,] 1.0000000 1.055984e-205 5.279918e-206
[103,] 1.0000000 7.895318e-207 3.947659e-207
[104,] 1.0000000 4.953772e-208 2.476886e-208
[105,] 1.0000000 2.710585e-207 1.355292e-207
[106,] 1.0000000 8.919880e-207 4.459940e-207
[107,] 1.0000000 2.636587e-206 1.318293e-206
[108,] 1.0000000 2.469933e-205 1.234967e-205
[109,] 1.0000000 9.358925e-207 4.679463e-207
[110,] 1.0000000 1.396874e-207 6.984371e-208
[111,] 1.0000000 1.016909e-206 5.084546e-207
[112,] 1.0000000 2.280916e-208 1.140458e-208
[113,] 1.0000000 1.660128e-207 8.300642e-208
[114,] 1.0000000 5.227016e-207 2.613508e-207
[115,] 1.0000000 2.302634e-208 1.151317e-208
[116,] 1.0000000 2.451444e-207 1.225722e-207
[117,] 1.0000000 2.666786e-207 1.333393e-207
[118,] 1.0000000 2.910706e-206 1.455353e-206
[119,] 1.0000000 3.280003e-205 1.640001e-205
[120,] 1.0000000 3.693380e-204 1.846690e-204
[121,] 1.0000000 3.990237e-203 1.995119e-203
[122,] 1.0000000 4.360552e-202 2.180276e-202
[123,] 1.0000000 1.544943e-202 7.724715e-203
[124,] 1.0000000 1.677675e-201 8.388377e-202
[125,] 1.0000000 1.868701e-200 9.343507e-201
[126,] 1.0000000 2.010425e-199 1.005212e-199
[127,] 1.0000000 2.030963e-198 1.015481e-198
[128,] 1.0000000 2.200412e-197 1.100206e-197
[129,] 1.0000000 2.341951e-196 1.170975e-196
[130,] 1.0000000 2.483175e-195 1.241587e-195
[131,] 1.0000000 2.622741e-194 1.311370e-194
[132,] 1.0000000 2.741870e-193 1.370935e-193
[133,] 1.0000000 2.873026e-192 1.436513e-192
[134,] 1.0000000 2.521973e-191 1.260986e-191
[135,] 1.0000000 2.647161e-190 1.323581e-190
[136,] 1.0000000 2.216655e-189 1.108328e-189
[137,] 1.0000000 2.246036e-188 1.123018e-188
[138,] 1.0000000 2.306315e-187 1.153157e-187
[139,] 1.0000000 2.386137e-186 1.193068e-186
[140,] 1.0000000 2.428429e-185 1.214214e-185
[141,] 1.0000000 2.454408e-184 1.227204e-184
[142,] 1.0000000 2.498241e-183 1.249120e-183
[143,] 1.0000000 1.846172e-182 9.230859e-183
[144,] 1.0000000 1.874281e-181 9.371405e-182
[145,] 1.0000000 1.766246e-180 8.831229e-181
[146,] 1.0000000 1.752562e-179 8.762812e-180
[147,] 1.0000000 1.729442e-178 8.647208e-179
[148,] 1.0000000 1.699032e-177 8.495161e-178
[149,] 1.0000000 1.664753e-176 8.323765e-177
[150,] 1.0000000 1.350890e-175 6.754449e-176
[151,] 1.0000000 1.310423e-174 6.552114e-175
[152,] 1.0000000 1.265427e-173 6.327134e-174
[153,] 1.0000000 1.216441e-172 6.082205e-173
[154,] 1.0000000 1.164041e-171 5.820206e-172
[155,] 1.0000000 9.803631e-171 4.901816e-171
[156,] 1.0000000 9.153033e-170 4.576516e-170
[157,] 1.0000000 8.648284e-169 4.324142e-169
[158,] 1.0000000 8.133850e-168 4.066925e-168
[159,] 1.0000000 2.607154e-168 1.303577e-168
[160,] 1.0000000 2.485401e-167 1.242701e-167
[161,] 1.0000000 2.337313e-166 1.168657e-166
[162,] 1.0000000 2.187902e-165 1.093951e-165
[163,] 1.0000000 2.041177e-164 1.020588e-164
[164,] 1.0000000 1.892998e-163 9.464989e-164
[165,] 1.0000000 1.746163e-162 8.730814e-163
[166,] 1.0000000 1.604275e-161 8.021376e-162
[167,] 1.0000000 1.472332e-160 7.361659e-161
[168,] 1.0000000 1.340108e-159 6.700541e-160
[169,] 1.0000000 1.184561e-158 5.922804e-159
[170,] 1.0000000 1.067843e-157 5.339216e-158
[171,] 1.0000000 9.581788e-157 4.790894e-157
[172,] 1.0000000 2.085174e-156 1.042587e-156
[173,] 1.0000000 1.594143e-155 7.970716e-156
[174,] 1.0000000 9.994779e-155 4.997389e-155
[175,] 1.0000000 5.513397e-154 2.756699e-154
[176,] 1.0000000 4.752891e-153 2.376446e-153
[177,] 1.0000000 3.883291e-152 1.941645e-152
[178,] 1.0000000 2.119812e-151 1.059906e-151
[179,] 1.0000000 1.853598e-150 9.267988e-151
[180,] 1.0000000 1.407746e-149 7.038730e-150
[181,] 1.0000000 1.220464e-148 6.102320e-149
[182,] 1.0000000 1.011774e-147 5.058870e-148
[183,] 1.0000000 8.692031e-147 4.346015e-147
[184,] 1.0000000 7.321939e-146 3.660970e-146
[185,] 1.0000000 6.232571e-145 3.116286e-145
[186,] 1.0000000 5.299201e-144 2.649601e-144
[187,] 1.0000000 4.385964e-143 2.192982e-143
[188,] 1.0000000 2.157530e-142 1.078765e-142
[189,] 1.0000000 1.779330e-141 8.896650e-142
[190,] 1.0000000 1.462099e-140 7.310494e-141
[191,] 1.0000000 8.171960e-140 4.085980e-140
[192,] 1.0000000 6.674793e-139 3.337396e-139
[193,] 1.0000000 5.430631e-138 2.715316e-138
[194,] 1.0000000 4.400840e-137 2.200420e-137
[195,] 1.0000000 3.558390e-136 1.779195e-136
[196,] 1.0000000 2.442963e-135 1.221482e-135
[197,] 1.0000000 1.929561e-134 9.647805e-135
[198,] 1.0000000 1.527932e-133 7.639658e-134
[199,] 1.0000000 1.213008e-132 6.065040e-133
[200,] 1.0000000 1.783082e-132 8.915412e-133
[201,] 1.0000000 1.418570e-131 7.092849e-132
[202,] 1.0000000 1.128402e-130 5.642008e-131
[203,] 1.0000000 8.935136e-130 4.467568e-130
[204,] 1.0000000 6.712429e-129 3.356214e-129
[205,] 1.0000000 5.185937e-128 2.592968e-128
[206,] 1.0000000 4.013507e-127 2.006754e-127
[207,] 1.0000000 3.132851e-126 1.566426e-126
[208,] 1.0000000 2.430512e-125 1.215256e-125
[209,] 1.0000000 1.759180e-124 8.795901e-125
[210,] 1.0000000 1.354653e-123 6.773263e-124
[211,] 1.0000000 7.173323e-123 3.586661e-123
[212,] 1.0000000 4.729795e-122 2.364897e-122
[213,] 1.0000000 3.588150e-121 1.794075e-121
[214,] 1.0000000 2.708732e-120 1.354366e-120
[215,] 1.0000000 2.049289e-119 1.024644e-119
[216,] 1.0000000 1.532110e-118 7.660552e-119
[217,] 1.0000000 1.135773e-117 5.678867e-118
[218,] 1.0000000 8.417690e-117 4.208845e-117
[219,] 1.0000000 6.209841e-116 3.104920e-116
[220,] 1.0000000 4.453989e-115 2.226995e-115
[221,] 1.0000000 3.255831e-114 1.627916e-114
[222,] 1.0000000 2.365443e-113 1.182721e-113
[223,] 1.0000000 1.713044e-112 8.565219e-113
[224,] 1.0000000 1.234800e-111 6.174000e-112
[225,] 1.0000000 8.511066e-111 4.255533e-111
[226,] 1.0000000 6.078732e-110 3.039366e-110
[227,] 1.0000000 4.317244e-109 2.158622e-109
[228,] 1.0000000 3.054662e-108 1.527331e-108
[229,] 1.0000000 2.155014e-107 1.077507e-107
[230,] 1.0000000 1.510410e-106 7.552051e-107
[231,] 1.0000000 1.053633e-105 5.268163e-106
[232,] 1.0000000 7.315254e-105 3.657627e-105
[233,] 1.0000000 4.967653e-104 2.483826e-104
[234,] 1.0000000 3.416655e-103 1.708327e-103
[235,] 1.0000000 2.345060e-102 1.172530e-102
[236,] 1.0000000 9.635195e-102 4.817598e-102
[237,] 1.0000000 6.502954e-101 3.251477e-101
[238,] 1.0000000 4.369610e-100 2.184805e-100
[239,] 1.0000000 2.928495e-99 1.464248e-99
[240,] 1.0000000 1.949721e-98 9.748605e-99
[241,] 1.0000000 1.292247e-97 6.461236e-98
[242,] 1.0000000 8.585482e-97 4.292741e-97
[243,] 1.0000000 5.637233e-96 2.818617e-96
[244,] 1.0000000 3.649198e-95 1.824599e-95
[245,] 1.0000000 2.210509e-94 1.105254e-94
[246,] 1.0000000 8.656227e-94 4.328113e-94
[247,] 1.0000000 5.626883e-93 2.813441e-93
[248,] 1.0000000 3.095274e-92 1.547637e-92
[249,] 1.0000000 1.809838e-91 9.049189e-92
[250,] 1.0000000 1.155237e-90 5.776184e-91
[251,] 1.0000000 7.373300e-90 3.686650e-90
[252,] 1.0000000 4.661453e-89 2.330727e-89
[253,] 1.0000000 2.933174e-88 1.466587e-88
[254,] 1.0000000 1.813141e-87 9.065706e-88
[255,] 1.0000000 1.091623e-86 5.458116e-87
[256,] 1.0000000 6.771084e-86 3.385542e-86
[257,] 1.0000000 4.195986e-85 2.097993e-85
[258,] 1.0000000 2.577269e-84 1.288634e-84
[259,] 1.0000000 1.004948e-83 5.024741e-84
[260,] 1.0000000 6.129898e-83 3.064949e-83
[261,] 1.0000000 3.681522e-82 1.840761e-82
[262,] 1.0000000 2.223619e-81 1.111810e-81
[263,] 1.0000000 1.336425e-80 6.682127e-81
[264,] 1.0000000 7.992157e-80 3.996078e-80
[265,] 1.0000000 4.770203e-79 2.385101e-79
[266,] 1.0000000 2.832675e-78 1.416338e-78
[267,] 1.0000000 1.606535e-77 8.032676e-78
[268,] 1.0000000 9.158524e-77 4.579262e-77
[269,] 1.0000000 5.357910e-76 2.678955e-76
[270,] 1.0000000 2.908977e-75 1.454488e-75
[271,] 1.0000000 1.599946e-74 7.999729e-75
[272,] 1.0000000 9.184860e-74 4.592430e-74
[273,] 1.0000000 5.267360e-73 2.633680e-73
[274,] 1.0000000 2.809018e-72 1.404509e-72
[275,] 1.0000000 1.595049e-71 7.975247e-72
[276,] 1.0000000 8.959682e-71 4.479841e-71
[277,] 1.0000000 5.034532e-70 2.517266e-70
[278,] 1.0000000 2.801642e-69 1.400821e-69
[279,] 1.0000000 1.555724e-68 7.778620e-69
[280,] 1.0000000 8.565057e-68 4.282529e-68
[281,] 1.0000000 4.406702e-67 2.203351e-67
[282,] 1.0000000 2.401065e-66 1.200533e-66
[283,] 1.0000000 1.271454e-65 6.357270e-66
[284,] 1.0000000 6.229709e-65 3.114854e-65
[285,] 1.0000000 3.358146e-64 1.679073e-64
[286,] 1.0000000 1.805286e-63 9.026428e-64
[287,] 1.0000000 7.432713e-63 3.716357e-63
[288,] 1.0000000 3.949071e-62 1.974536e-62
[289,] 1.0000000 2.081487e-61 1.040744e-61
[290,] 1.0000000 1.093651e-60 5.468254e-61
[291,] 1.0000000 5.724377e-60 2.862188e-60
[292,] 1.0000000 2.953599e-59 1.476799e-59
[293,] 1.0000000 1.525815e-58 7.629074e-59
[294,] 1.0000000 7.837902e-58 3.918951e-58
[295,] 1.0000000 3.993526e-57 1.996763e-57
[296,] 1.0000000 1.537549e-56 7.687743e-57
[297,] 1.0000000 7.801712e-56 3.900856e-56
[298,] 1.0000000 3.916165e-55 1.958082e-55
[299,] 1.0000000 1.915961e-54 9.579807e-55
[300,] 1.0000000 8.885438e-54 4.442719e-54
[301,] 1.0000000 4.385613e-53 2.192806e-53
[302,] 1.0000000 2.151387e-52 1.075694e-52
[303,] 1.0000000 8.878233e-52 4.439116e-52
[304,] 1.0000000 4.309386e-51 2.154693e-51
[305,] 1.0000000 2.058336e-50 1.029168e-50
[306,] 1.0000000 9.918849e-50 4.959425e-50
[307,] 1.0000000 4.570749e-49 2.285375e-49
[308,] 1.0000000 2.169033e-48 1.084516e-48
[309,] 1.0000000 1.022335e-47 5.111677e-48
[310,] 1.0000000 4.804484e-47 2.402242e-47
[311,] 1.0000000 2.240942e-46 1.120471e-46
[312,] 1.0000000 1.038664e-45 5.193319e-46
[313,] 1.0000000 4.783678e-45 2.391839e-45
[314,] 1.0000000 2.189137e-44 1.094569e-44
[315,] 1.0000000 9.927348e-44 4.963674e-44
[316,] 1.0000000 4.485069e-43 2.242534e-43
[317,] 1.0000000 2.016889e-42 1.008444e-42
[318,] 1.0000000 8.965039e-42 4.482520e-42
[319,] 1.0000000 3.955446e-41 1.977723e-41
[320,] 1.0000000 1.740893e-40 8.704467e-41
[321,] 1.0000000 7.610872e-40 3.805436e-40
[322,] 1.0000000 3.304914e-39 1.652457e-39
[323,] 1.0000000 1.425367e-38 7.126836e-39
[324,] 1.0000000 6.105360e-38 3.052680e-38
[325,] 1.0000000 2.595886e-37 1.297943e-37
[326,] 1.0000000 1.096565e-36 5.482826e-37
[327,] 1.0000000 4.599661e-36 2.299831e-36
[328,] 1.0000000 1.917371e-35 9.586853e-36
[329,] 1.0000000 7.928914e-35 3.964457e-35
[330,] 1.0000000 3.254633e-34 1.627316e-34
[331,] 1.0000000 1.326541e-33 6.632707e-34
[332,] 1.0000000 5.362006e-33 2.681003e-33
[333,] 1.0000000 2.153727e-32 1.076864e-32
[334,] 1.0000000 8.586370e-32 4.293185e-32
[335,] 1.0000000 3.397456e-31 1.698728e-31
[336,] 1.0000000 1.334113e-30 6.670565e-31
[337,] 1.0000000 5.198691e-30 2.599346e-30
[338,] 1.0000000 2.009276e-29 1.004638e-29
[339,] 1.0000000 7.709729e-29 3.854865e-29
[340,] 1.0000000 2.934955e-28 1.467478e-28
[341,] 1.0000000 1.108387e-27 5.541934e-28
[342,] 1.0000000 4.123429e-27 2.061714e-27
[343,] 1.0000000 1.533604e-26 7.668019e-27
[344,] 1.0000000 5.318953e-26 2.659476e-26
[345,] 1.0000000 1.948710e-25 9.743552e-26
[346,] 1.0000000 7.066181e-25 3.533090e-25
[347,] 1.0000000 2.540468e-24 1.270234e-24
[348,] 1.0000000 8.778946e-24 4.389473e-24
[349,] 1.0000000 3.104047e-23 1.552024e-23
[350,] 1.0000000 1.087860e-22 5.439302e-23
[351,] 1.0000000 3.778587e-22 1.889294e-22
[352,] 1.0000000 1.301628e-21 6.508142e-22
[353,] 1.0000000 4.439317e-21 2.219658e-21
[354,] 1.0000000 1.500047e-20 7.500233e-21
[355,] 1.0000000 4.963784e-20 2.481892e-20
[356,] 1.0000000 1.649745e-19 8.248726e-20
[357,] 1.0000000 5.419878e-19 2.709939e-19
[358,] 1.0000000 1.763185e-18 8.815926e-19
[359,] 1.0000000 5.679108e-18 2.839554e-18
[360,] 1.0000000 1.810806e-17 9.054031e-18
[361,] 1.0000000 5.714874e-17 2.857437e-17
[362,] 1.0000000 1.784904e-16 8.924522e-17
[363,] 1.0000000 5.465338e-16 2.732669e-16
[364,] 1.0000000 1.675998e-15 8.379989e-16
[365,] 1.0000000 5.084768e-15 2.542384e-15
[366,] 1.0000000 1.521808e-14 7.609042e-15
[367,] 1.0000000 4.503430e-14 2.251715e-14
[368,] 1.0000000 1.320324e-13 6.601618e-14
[369,] 1.0000000 3.789784e-13 1.894892e-13
[370,] 1.0000000 1.082890e-12 5.414452e-13
[371,] 1.0000000 3.056873e-12 1.528436e-12
[372,] 1.0000000 8.522921e-12 4.261461e-12
[373,] 1.0000000 2.332347e-11 1.166173e-11
[374,] 1.0000000 5.600820e-11 2.800410e-11
[375,] 1.0000000 1.507013e-10 7.535066e-11
[376,] 1.0000000 3.996184e-10 1.998092e-10
[377,] 1.0000000 9.117798e-10 4.558899e-10
[378,] 1.0000000 2.348823e-09 1.174412e-09
[379,] 1.0000000 5.804373e-09 2.902187e-09
[380,] 1.0000000 1.466598e-08 7.332992e-09
[381,] 1.0000000 3.665646e-08 1.832823e-08
[382,] 1.0000000 8.777207e-08 4.388603e-08
[383,] 0.9999999 2.130467e-07 1.065234e-07
[384,] 0.9999997 5.097429e-07 2.548714e-07
[385,] 0.9999995 9.650960e-07 4.825480e-07
[386,] 0.9999990 2.098654e-06 1.049327e-06
[387,] 0.9999976 4.823996e-06 2.411998e-06
[388,] 0.9999949 1.017128e-05 5.085642e-06
[389,] 0.9999932 1.369560e-05 6.847801e-06
[390,] 0.9999849 3.021501e-05 1.510750e-05
[391,] 0.9999750 5.006669e-05 2.503334e-05
[392,] 0.9999484 1.032830e-04 5.164149e-05
[393,] 0.9998980 2.040522e-04 1.020261e-04
[394,] 0.9997933 4.133326e-04 2.066663e-04
[395,] 0.9995787 8.426896e-04 4.213448e-04
[396,] 0.9991633 1.673403e-03 8.367016e-04
[397,] 0.9983938 3.212325e-03 1.606163e-03
[398,] 0.9970859 5.828218e-03 2.914109e-03
[399,] 0.9948777 1.024464e-02 5.122321e-03
[400,] 0.9909115 1.817702e-02 9.088511e-03
[401,] 0.9837571 3.248589e-02 1.624295e-02
[402,] 0.9719053 5.618932e-02 2.809466e-02
[403,] 0.9572517 8.549669e-02 4.274834e-02
[404,] 0.9325706 1.348587e-01 6.742936e-02
[405,] 0.8936672 2.126657e-01 1.063328e-01
[406,] 0.8372058 3.255884e-01 1.627942e-01
[407,] 0.7608735 4.782530e-01 2.391265e-01
[408,] 0.9735282 5.294361e-02 2.647180e-02
[409,] 0.9434052 1.131895e-01 5.659476e-02
[410,] 0.9051731 1.896538e-01 9.482689e-02
[411,] 0.8616411 2.767177e-01 1.383589e-01
[412,] 0.7500628 4.998745e-01 2.499372e-01
> postscript(file="/var/www/html/freestat/rcomp/tmp/198301291358215.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/2kh231291358215.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/3kh231291358215.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/4kh231291358215.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/5kh231291358215.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
-277397.7292 2557837.5403 1524849.5589 -872756.8463 -260723.9355
6 7 8 9 10
499912.6898 50803.3342 233096.1568 546213.5869 765744.5259
11 12 13 14 15
118815.3978 415065.5112 363390.6666 -1106733.5215 81428.2065
16 17 18 19 20
598941.0692 -285286.3107 -223449.8448 222456.0823 -781630.3797
21 22 23 24 25
382247.6700 123835.2468 -94773.1470 -1025742.4420 -150822.8397
26 27 28 29 30
72550.9570 -272316.0688 -248027.7430 398664.0991 534699.5439
31 32 33 34 35
269029.6805 22432.4533 232558.4609 129733.1617 495152.7680
36 37 38 39 40
-677577.9129 156925.7774 91863.3885 -38270.7135 -342467.4789
41 42 43 44 45
-230150.3517 -128825.6111 186076.1321 -165671.9290 37071.3337
46 47 48 49 50
-86017.8086 231502.1696 -377790.7035 -1374.9665 -195568.0272
51 52 53 54 55
-30126.9217 -329273.4078 -340355.4879 -16799.7902 -506760.5780
56 57 58 59 60
34203.7748 180846.5192 59780.1536 23037.7797 -163328.4059
61 62 63 64 65
18400.9273 -113341.5416 9576.2456 -195793.0162 24724.8221
66 67 68 69 70
563636.9711 24732.4057 -46375.2413 -73770.5888 -204071.2810
71 72 73 74 75
-145607.0407 32454.6517 -61495.9718 18524.9487 -5115.8936
76 77 78 79 80
-106754.6017 19465.0654 18368.4564 115684.9807 195574.6939
81 82 83 84 85
-18171.0815 -76936.2201 36834.4012 -26100.2435 93365.5534
86 87 88 89 90
-137127.4478 -117531.4463 71990.4549 -37266.1893 82710.6706
91 92 93 94 95
29855.4842 -76301.2867 77576.4644 -60184.6443 -774.7735
96 97 98 99 100
-12507.9392 164297.9901 368028.3966 91437.6969 -42069.6757
101 102 103 104 105
-97365.9034 -43810.1211 -51991.9092 117123.1199 29944.0083
106 107 108 109 110
-97531.6477 -266971.0776 128424.8455 -91738.3138 -100438.5848
111 112 113 114 115
96922.0217 -254576.2834 -179297.6440 -118452.4114 -14628.1591
116 117 118 119 120
796.5017 88509.4743 123514.9171 -161398.0342 120703.3278
121 122 123 124 125
-196652.6095 31573.9512 -200582.8122 76243.4774 -70105.2464
126 127 128 129 130
-195592.4510 23944.4897 -549.6192 -5360.9040 23944.4897
131 132 133 134 135
15113.5220 130137.6496 23944.4897 -53630.8243 23944.4897
136 137 138 139 140
34400.4677 15071.0163 23944.4897 23944.4897 23944.4897
141 142 143 144 145
24850.1705 23944.4897 10516.6532 15264.9363 18116.0680
146 147 148 149 150
12172.3994 23944.4897 9954.9738 23944.4897 26956.4369
151 152 153 154 155
9263.7127 10444.8738 3095.7455 -2841.6848 23274.9586
156 157 158 159 160
23944.4897 23944.4897 21721.8609 33740.4414 23944.4897
161 162 163 164 165
23944.4897 23944.4897 23944.4897 67491.0308 47608.3635
166 167 168 169 170
23944.4897 23944.4897 55212.5717 7265.4943 23944.4897
171 172 173 174 175
23944.4897 20907.2606 23944.4897 11860.3540 23944.4897
176 177 178 179 180
18500.2461 23944.4897 -28757.0335 23944.4897 23944.4897
181 182 183 184 185
83867.6908 -27383.2329 93898.3554 -86816.9886 -34571.1742
186 187 188 189 190
-43410.3807 -4819.4832 14322.6612 -1823.4228 23944.4897
191 192 193 194 195
8647.1684 23944.4897 15946.0910 23944.4897 16616.4565
196 197 198 199 200
-31164.0924 -162580.2171 -31164.0924 -31164.0924 -69414.3865
201 202 203 204 205
-31164.0924 -31164.0924 -31164.0924 -24625.9356 -26498.6333
206 207 208 209 210
-45596.5627 -25412.0680 -31164.0924 -28390.8870 -31164.0924
211 212 213 214 215
23944.4897 23944.4897 56991.7930 16129.6807 17321.8019
216 217 218 219 220
15604.9920 20018.5381 65743.9018 15604.9920 -14771.6417
221 222 223 224 225
-7888.7766 23944.4897 24859.2736 -1074.0034 25037.9354
226 227 228 229 230
-386.2961 23944.4897 23944.4897 19715.0099 23944.4897
231 232 233 234 235
11162.2202 23944.4897 23944.4897 17857.3750 23944.4897
236 237 238 239 240
26439.4266 23944.4897 14159.2402 23944.4897 23944.4897
241 242 243 244 245
23944.4897 25092.1689 23944.4897 13831.5253 31776.3908
246 247 248 249 250
-31164.0924 -31164.0924 -16839.4265 -31164.0924 -31164.0924
251 252 253 254 255
-12268.5111 -31164.0924 18088.9535 36297.4716 1066.2229
256 257 258 259 260
23944.4897 24008.5743 -3218.5157 -31164.0924 23944.4897
261 262 263 264 265
-31164.0924 -31164.0924 24904.4705 21890.4878 -31164.0924
266 267 268 269 270
23944.4897 -31164.0924 80047.1915 -31164.0924 -55236.0503
271 272 273 274 275
-31164.0924 -31164.0924 -31164.0924 -28665.2822 23944.4897
276 277 278 279 280
57477.8478 22102.3950 23944.4897 57933.9111 56208.6106
281 282 283 284 285
-31164.0924 23944.4897 102833.7758 23944.4897 21869.4117
286 287 288 289 290
23944.4897 -31164.0924 24714.2721 -31164.0924 121394.6428
291 292 293 294 295
-31164.0924 44696.2917 41880.9308 23944.4897 11990.8023
296 297 298 299 300
13311.6374 23944.4897 -31164.0924 23944.4897 -13583.2499
301 302 303 304 305
-39224.8440 8133.8845 23944.4897 -31164.0924 126341.8884
306 307 308 309 310
-2952.0544 -31164.0924 -74616.6713 8317.1496 -31164.0924
311 312 313 314 315
-31164.0924 96925.5331 -31164.0924 -61722.3414 1172.4179
316 317 318 319 320
59855.7592 25344.0301 -28662.1095 18498.3189 23944.4897
321 322 323 324 325
23944.4897 23944.4897 23944.4897 -31164.0924 23944.4897
326 327 328 329 330
21871.6492 -31164.0924 -31164.0924 23944.4897 23944.4897
331 332 333 334 335
23944.4897 23944.4897 23944.4897 18436.9685 24916.6033
336 337 338 339 340
23944.4897 23283.8362 23944.4897 22470.2156 12543.2403
341 342 343 344 345
766.5812 23944.4897 23944.4897 23944.4897 23944.4897
346 347 348 349 350
23944.4897 12241.4031 23944.4897 23944.4897 23944.4897
351 352 353 354 355
40236.0807 19881.8725 16754.7417 22233.3625 23944.4897
356 357 358 359 360
23944.4897 94666.5242 23944.4897 23944.4897 23944.4897
361 362 363 364 365
17534.0406 23944.4897 23944.4897 -37653.4017 7735.2309
366 367 368 369 370
23944.4897 23944.4897 23944.4897 23944.4897 23944.4897
371 372 373 374 375
23944.4897 88363.2415 38934.9047 6811.2441 23944.4897
376 377 378 379 380
23944.4897 28078.4082 14034.6561 23944.4897 23944.4897
381 382 383 384 385
23944.4897 108817.1919 164863.0184 23944.4897 18419.9700
386 387 388 389 390
-97878.2616 27285.8956 10634.9196 16119.9896 8879.0068
391 392 393 394 395
15693.2771 -936.5007 -18866.5525 74608.9461 -99649.4056
396 397 398 399 400
-6594.9425 -82488.7376 -697508.8098 9888.3789 121170.3199
401 402 403 404 405
43695.7911 28536.0182 1668.2186 -88304.8133 1638.0231
406 407 408 409 410
14191.8476 34921.3543 -10054.1141 -12874.9790 -5984.2293
411 412 413 414 415
-82507.1098 -113565.0487 -75888.7341 -122992.5639 6515.1212
416 417 418 419 420
-2245.6856 160267.5784 -103626.4598 -18465.1137 -8329.4936
421 422 423 424 425
-40960.8889 -216291.7420 -497573.8832 -414112.1634 -109008.0932
426 427 428 429 430
-249580.5640 -311658.1463 -1284038.3256 -184370.2558 -134950.8637
431
-377355.7461
> postscript(file="/var/www/html/freestat/rcomp/tmp/6icqu1291358215.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 -277397.7292 NA
1 2557837.5403 -277397.7292
2 1524849.5589 2557837.5403
3 -872756.8463 1524849.5589
4 -260723.9355 -872756.8463
5 499912.6898 -260723.9355
6 50803.3342 499912.6898
7 233096.1568 50803.3342
8 546213.5869 233096.1568
9 765744.5259 546213.5869
10 118815.3978 765744.5259
11 415065.5112 118815.3978
12 363390.6666 415065.5112
13 -1106733.5215 363390.6666
14 81428.2065 -1106733.5215
15 598941.0692 81428.2065
16 -285286.3107 598941.0692
17 -223449.8448 -285286.3107
18 222456.0823 -223449.8448
19 -781630.3797 222456.0823
20 382247.6700 -781630.3797
21 123835.2468 382247.6700
22 -94773.1470 123835.2468
23 -1025742.4420 -94773.1470
24 -150822.8397 -1025742.4420
25 72550.9570 -150822.8397
26 -272316.0688 72550.9570
27 -248027.7430 -272316.0688
28 398664.0991 -248027.7430
29 534699.5439 398664.0991
30 269029.6805 534699.5439
31 22432.4533 269029.6805
32 232558.4609 22432.4533
33 129733.1617 232558.4609
34 495152.7680 129733.1617
35 -677577.9129 495152.7680
36 156925.7774 -677577.9129
37 91863.3885 156925.7774
38 -38270.7135 91863.3885
39 -342467.4789 -38270.7135
40 -230150.3517 -342467.4789
41 -128825.6111 -230150.3517
42 186076.1321 -128825.6111
43 -165671.9290 186076.1321
44 37071.3337 -165671.9290
45 -86017.8086 37071.3337
46 231502.1696 -86017.8086
47 -377790.7035 231502.1696
48 -1374.9665 -377790.7035
49 -195568.0272 -1374.9665
50 -30126.9217 -195568.0272
51 -329273.4078 -30126.9217
52 -340355.4879 -329273.4078
53 -16799.7902 -340355.4879
54 -506760.5780 -16799.7902
55 34203.7748 -506760.5780
56 180846.5192 34203.7748
57 59780.1536 180846.5192
58 23037.7797 59780.1536
59 -163328.4059 23037.7797
60 18400.9273 -163328.4059
61 -113341.5416 18400.9273
62 9576.2456 -113341.5416
63 -195793.0162 9576.2456
64 24724.8221 -195793.0162
65 563636.9711 24724.8221
66 24732.4057 563636.9711
67 -46375.2413 24732.4057
68 -73770.5888 -46375.2413
69 -204071.2810 -73770.5888
70 -145607.0407 -204071.2810
71 32454.6517 -145607.0407
72 -61495.9718 32454.6517
73 18524.9487 -61495.9718
74 -5115.8936 18524.9487
75 -106754.6017 -5115.8936
76 19465.0654 -106754.6017
77 18368.4564 19465.0654
78 115684.9807 18368.4564
79 195574.6939 115684.9807
80 -18171.0815 195574.6939
81 -76936.2201 -18171.0815
82 36834.4012 -76936.2201
83 -26100.2435 36834.4012
84 93365.5534 -26100.2435
85 -137127.4478 93365.5534
86 -117531.4463 -137127.4478
87 71990.4549 -117531.4463
88 -37266.1893 71990.4549
89 82710.6706 -37266.1893
90 29855.4842 82710.6706
91 -76301.2867 29855.4842
92 77576.4644 -76301.2867
93 -60184.6443 77576.4644
94 -774.7735 -60184.6443
95 -12507.9392 -774.7735
96 164297.9901 -12507.9392
97 368028.3966 164297.9901
98 91437.6969 368028.3966
99 -42069.6757 91437.6969
100 -97365.9034 -42069.6757
101 -43810.1211 -97365.9034
102 -51991.9092 -43810.1211
103 117123.1199 -51991.9092
104 29944.0083 117123.1199
105 -97531.6477 29944.0083
106 -266971.0776 -97531.6477
107 128424.8455 -266971.0776
108 -91738.3138 128424.8455
109 -100438.5848 -91738.3138
110 96922.0217 -100438.5848
111 -254576.2834 96922.0217
112 -179297.6440 -254576.2834
113 -118452.4114 -179297.6440
114 -14628.1591 -118452.4114
115 796.5017 -14628.1591
116 88509.4743 796.5017
117 123514.9171 88509.4743
118 -161398.0342 123514.9171
119 120703.3278 -161398.0342
120 -196652.6095 120703.3278
121 31573.9512 -196652.6095
122 -200582.8122 31573.9512
123 76243.4774 -200582.8122
124 -70105.2464 76243.4774
125 -195592.4510 -70105.2464
126 23944.4897 -195592.4510
127 -549.6192 23944.4897
128 -5360.9040 -549.6192
129 23944.4897 -5360.9040
130 15113.5220 23944.4897
131 130137.6496 15113.5220
132 23944.4897 130137.6496
133 -53630.8243 23944.4897
134 23944.4897 -53630.8243
135 34400.4677 23944.4897
136 15071.0163 34400.4677
137 23944.4897 15071.0163
138 23944.4897 23944.4897
139 23944.4897 23944.4897
140 24850.1705 23944.4897
141 23944.4897 24850.1705
142 10516.6532 23944.4897
143 15264.9363 10516.6532
144 18116.0680 15264.9363
145 12172.3994 18116.0680
146 23944.4897 12172.3994
147 9954.9738 23944.4897
148 23944.4897 9954.9738
149 26956.4369 23944.4897
150 9263.7127 26956.4369
151 10444.8738 9263.7127
152 3095.7455 10444.8738
153 -2841.6848 3095.7455
154 23274.9586 -2841.6848
155 23944.4897 23274.9586
156 23944.4897 23944.4897
157 21721.8609 23944.4897
158 33740.4414 21721.8609
159 23944.4897 33740.4414
160 23944.4897 23944.4897
161 23944.4897 23944.4897
162 23944.4897 23944.4897
163 67491.0308 23944.4897
164 47608.3635 67491.0308
165 23944.4897 47608.3635
166 23944.4897 23944.4897
167 55212.5717 23944.4897
168 7265.4943 55212.5717
169 23944.4897 7265.4943
170 23944.4897 23944.4897
171 20907.2606 23944.4897
172 23944.4897 20907.2606
173 11860.3540 23944.4897
174 23944.4897 11860.3540
175 18500.2461 23944.4897
176 23944.4897 18500.2461
177 -28757.0335 23944.4897
178 23944.4897 -28757.0335
179 23944.4897 23944.4897
180 83867.6908 23944.4897
181 -27383.2329 83867.6908
182 93898.3554 -27383.2329
183 -86816.9886 93898.3554
184 -34571.1742 -86816.9886
185 -43410.3807 -34571.1742
186 -4819.4832 -43410.3807
187 14322.6612 -4819.4832
188 -1823.4228 14322.6612
189 23944.4897 -1823.4228
190 8647.1684 23944.4897
191 23944.4897 8647.1684
192 15946.0910 23944.4897
193 23944.4897 15946.0910
194 16616.4565 23944.4897
195 -31164.0924 16616.4565
196 -162580.2171 -31164.0924
197 -31164.0924 -162580.2171
198 -31164.0924 -31164.0924
199 -69414.3865 -31164.0924
200 -31164.0924 -69414.3865
201 -31164.0924 -31164.0924
202 -31164.0924 -31164.0924
203 -24625.9356 -31164.0924
204 -26498.6333 -24625.9356
205 -45596.5627 -26498.6333
206 -25412.0680 -45596.5627
207 -31164.0924 -25412.0680
208 -28390.8870 -31164.0924
209 -31164.0924 -28390.8870
210 23944.4897 -31164.0924
211 23944.4897 23944.4897
212 56991.7930 23944.4897
213 16129.6807 56991.7930
214 17321.8019 16129.6807
215 15604.9920 17321.8019
216 20018.5381 15604.9920
217 65743.9018 20018.5381
218 15604.9920 65743.9018
219 -14771.6417 15604.9920
220 -7888.7766 -14771.6417
221 23944.4897 -7888.7766
222 24859.2736 23944.4897
223 -1074.0034 24859.2736
224 25037.9354 -1074.0034
225 -386.2961 25037.9354
226 23944.4897 -386.2961
227 23944.4897 23944.4897
228 19715.0099 23944.4897
229 23944.4897 19715.0099
230 11162.2202 23944.4897
231 23944.4897 11162.2202
232 23944.4897 23944.4897
233 17857.3750 23944.4897
234 23944.4897 17857.3750
235 26439.4266 23944.4897
236 23944.4897 26439.4266
237 14159.2402 23944.4897
238 23944.4897 14159.2402
239 23944.4897 23944.4897
240 23944.4897 23944.4897
241 25092.1689 23944.4897
242 23944.4897 25092.1689
243 13831.5253 23944.4897
244 31776.3908 13831.5253
245 -31164.0924 31776.3908
246 -31164.0924 -31164.0924
247 -16839.4265 -31164.0924
248 -31164.0924 -16839.4265
249 -31164.0924 -31164.0924
250 -12268.5111 -31164.0924
251 -31164.0924 -12268.5111
252 18088.9535 -31164.0924
253 36297.4716 18088.9535
254 1066.2229 36297.4716
255 23944.4897 1066.2229
256 24008.5743 23944.4897
257 -3218.5157 24008.5743
258 -31164.0924 -3218.5157
259 23944.4897 -31164.0924
260 -31164.0924 23944.4897
261 -31164.0924 -31164.0924
262 24904.4705 -31164.0924
263 21890.4878 24904.4705
264 -31164.0924 21890.4878
265 23944.4897 -31164.0924
266 -31164.0924 23944.4897
267 80047.1915 -31164.0924
268 -31164.0924 80047.1915
269 -55236.0503 -31164.0924
270 -31164.0924 -55236.0503
271 -31164.0924 -31164.0924
272 -31164.0924 -31164.0924
273 -28665.2822 -31164.0924
274 23944.4897 -28665.2822
275 57477.8478 23944.4897
276 22102.3950 57477.8478
277 23944.4897 22102.3950
278 57933.9111 23944.4897
279 56208.6106 57933.9111
280 -31164.0924 56208.6106
281 23944.4897 -31164.0924
282 102833.7758 23944.4897
283 23944.4897 102833.7758
284 21869.4117 23944.4897
285 23944.4897 21869.4117
286 -31164.0924 23944.4897
287 24714.2721 -31164.0924
288 -31164.0924 24714.2721
289 121394.6428 -31164.0924
290 -31164.0924 121394.6428
291 44696.2917 -31164.0924
292 41880.9308 44696.2917
293 23944.4897 41880.9308
294 11990.8023 23944.4897
295 13311.6374 11990.8023
296 23944.4897 13311.6374
297 -31164.0924 23944.4897
298 23944.4897 -31164.0924
299 -13583.2499 23944.4897
300 -39224.8440 -13583.2499
301 8133.8845 -39224.8440
302 23944.4897 8133.8845
303 -31164.0924 23944.4897
304 126341.8884 -31164.0924
305 -2952.0544 126341.8884
306 -31164.0924 -2952.0544
307 -74616.6713 -31164.0924
308 8317.1496 -74616.6713
309 -31164.0924 8317.1496
310 -31164.0924 -31164.0924
311 96925.5331 -31164.0924
312 -31164.0924 96925.5331
313 -61722.3414 -31164.0924
314 1172.4179 -61722.3414
315 59855.7592 1172.4179
316 25344.0301 59855.7592
317 -28662.1095 25344.0301
318 18498.3189 -28662.1095
319 23944.4897 18498.3189
320 23944.4897 23944.4897
321 23944.4897 23944.4897
322 23944.4897 23944.4897
323 -31164.0924 23944.4897
324 23944.4897 -31164.0924
325 21871.6492 23944.4897
326 -31164.0924 21871.6492
327 -31164.0924 -31164.0924
328 23944.4897 -31164.0924
329 23944.4897 23944.4897
330 23944.4897 23944.4897
331 23944.4897 23944.4897
332 23944.4897 23944.4897
333 18436.9685 23944.4897
334 24916.6033 18436.9685
335 23944.4897 24916.6033
336 23283.8362 23944.4897
337 23944.4897 23283.8362
338 22470.2156 23944.4897
339 12543.2403 22470.2156
340 766.5812 12543.2403
341 23944.4897 766.5812
342 23944.4897 23944.4897
343 23944.4897 23944.4897
344 23944.4897 23944.4897
345 23944.4897 23944.4897
346 12241.4031 23944.4897
347 23944.4897 12241.4031
348 23944.4897 23944.4897
349 23944.4897 23944.4897
350 40236.0807 23944.4897
351 19881.8725 40236.0807
352 16754.7417 19881.8725
353 22233.3625 16754.7417
354 23944.4897 22233.3625
355 23944.4897 23944.4897
356 94666.5242 23944.4897
357 23944.4897 94666.5242
358 23944.4897 23944.4897
359 23944.4897 23944.4897
360 17534.0406 23944.4897
361 23944.4897 17534.0406
362 23944.4897 23944.4897
363 -37653.4017 23944.4897
364 7735.2309 -37653.4017
365 23944.4897 7735.2309
366 23944.4897 23944.4897
367 23944.4897 23944.4897
368 23944.4897 23944.4897
369 23944.4897 23944.4897
370 23944.4897 23944.4897
371 88363.2415 23944.4897
372 38934.9047 88363.2415
373 6811.2441 38934.9047
374 23944.4897 6811.2441
375 23944.4897 23944.4897
376 28078.4082 23944.4897
377 14034.6561 28078.4082
378 23944.4897 14034.6561
379 23944.4897 23944.4897
380 23944.4897 23944.4897
381 108817.1919 23944.4897
382 164863.0184 108817.1919
383 23944.4897 164863.0184
384 18419.9700 23944.4897
385 -97878.2616 18419.9700
386 27285.8956 -97878.2616
387 10634.9196 27285.8956
388 16119.9896 10634.9196
389 8879.0068 16119.9896
390 15693.2771 8879.0068
391 -936.5007 15693.2771
392 -18866.5525 -936.5007
393 74608.9461 -18866.5525
394 -99649.4056 74608.9461
395 -6594.9425 -99649.4056
396 -82488.7376 -6594.9425
397 -697508.8098 -82488.7376
398 9888.3789 -697508.8098
399 121170.3199 9888.3789
400 43695.7911 121170.3199
401 28536.0182 43695.7911
402 1668.2186 28536.0182
403 -88304.8133 1668.2186
404 1638.0231 -88304.8133
405 14191.8476 1638.0231
406 34921.3543 14191.8476
407 -10054.1141 34921.3543
408 -12874.9790 -10054.1141
409 -5984.2293 -12874.9790
410 -82507.1098 -5984.2293
411 -113565.0487 -82507.1098
412 -75888.7341 -113565.0487
413 -122992.5639 -75888.7341
414 6515.1212 -122992.5639
415 -2245.6856 6515.1212
416 160267.5784 -2245.6856
417 -103626.4598 160267.5784
418 -18465.1137 -103626.4598
419 -8329.4936 -18465.1137
420 -40960.8889 -8329.4936
421 -216291.7420 -40960.8889
422 -497573.8832 -216291.7420
423 -414112.1634 -497573.8832
424 -109008.0932 -414112.1634
425 -249580.5640 -109008.0932
426 -311658.1463 -249580.5640
427 -1284038.3256 -311658.1463
428 -184370.2558 -1284038.3256
429 -134950.8637 -184370.2558
430 -377355.7461 -134950.8637
431 NA -377355.7461
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] 2557837.5403 -277397.7292
[2,] 1524849.5589 2557837.5403
[3,] -872756.8463 1524849.5589
[4,] -260723.9355 -872756.8463
[5,] 499912.6898 -260723.9355
[6,] 50803.3342 499912.6898
[7,] 233096.1568 50803.3342
[8,] 546213.5869 233096.1568
[9,] 765744.5259 546213.5869
[10,] 118815.3978 765744.5259
[11,] 415065.5112 118815.3978
[12,] 363390.6666 415065.5112
[13,] -1106733.5215 363390.6666
[14,] 81428.2065 -1106733.5215
[15,] 598941.0692 81428.2065
[16,] -285286.3107 598941.0692
[17,] -223449.8448 -285286.3107
[18,] 222456.0823 -223449.8448
[19,] -781630.3797 222456.0823
[20,] 382247.6700 -781630.3797
[21,] 123835.2468 382247.6700
[22,] -94773.1470 123835.2468
[23,] -1025742.4420 -94773.1470
[24,] -150822.8397 -1025742.4420
[25,] 72550.9570 -150822.8397
[26,] -272316.0688 72550.9570
[27,] -248027.7430 -272316.0688
[28,] 398664.0991 -248027.7430
[29,] 534699.5439 398664.0991
[30,] 269029.6805 534699.5439
[31,] 22432.4533 269029.6805
[32,] 232558.4609 22432.4533
[33,] 129733.1617 232558.4609
[34,] 495152.7680 129733.1617
[35,] -677577.9129 495152.7680
[36,] 156925.7774 -677577.9129
[37,] 91863.3885 156925.7774
[38,] -38270.7135 91863.3885
[39,] -342467.4789 -38270.7135
[40,] -230150.3517 -342467.4789
[41,] -128825.6111 -230150.3517
[42,] 186076.1321 -128825.6111
[43,] -165671.9290 186076.1321
[44,] 37071.3337 -165671.9290
[45,] -86017.8086 37071.3337
[46,] 231502.1696 -86017.8086
[47,] -377790.7035 231502.1696
[48,] -1374.9665 -377790.7035
[49,] -195568.0272 -1374.9665
[50,] -30126.9217 -195568.0272
[51,] -329273.4078 -30126.9217
[52,] -340355.4879 -329273.4078
[53,] -16799.7902 -340355.4879
[54,] -506760.5780 -16799.7902
[55,] 34203.7748 -506760.5780
[56,] 180846.5192 34203.7748
[57,] 59780.1536 180846.5192
[58,] 23037.7797 59780.1536
[59,] -163328.4059 23037.7797
[60,] 18400.9273 -163328.4059
[61,] -113341.5416 18400.9273
[62,] 9576.2456 -113341.5416
[63,] -195793.0162 9576.2456
[64,] 24724.8221 -195793.0162
[65,] 563636.9711 24724.8221
[66,] 24732.4057 563636.9711
[67,] -46375.2413 24732.4057
[68,] -73770.5888 -46375.2413
[69,] -204071.2810 -73770.5888
[70,] -145607.0407 -204071.2810
[71,] 32454.6517 -145607.0407
[72,] -61495.9718 32454.6517
[73,] 18524.9487 -61495.9718
[74,] -5115.8936 18524.9487
[75,] -106754.6017 -5115.8936
[76,] 19465.0654 -106754.6017
[77,] 18368.4564 19465.0654
[78,] 115684.9807 18368.4564
[79,] 195574.6939 115684.9807
[80,] -18171.0815 195574.6939
[81,] -76936.2201 -18171.0815
[82,] 36834.4012 -76936.2201
[83,] -26100.2435 36834.4012
[84,] 93365.5534 -26100.2435
[85,] -137127.4478 93365.5534
[86,] -117531.4463 -137127.4478
[87,] 71990.4549 -117531.4463
[88,] -37266.1893 71990.4549
[89,] 82710.6706 -37266.1893
[90,] 29855.4842 82710.6706
[91,] -76301.2867 29855.4842
[92,] 77576.4644 -76301.2867
[93,] -60184.6443 77576.4644
[94,] -774.7735 -60184.6443
[95,] -12507.9392 -774.7735
[96,] 164297.9901 -12507.9392
[97,] 368028.3966 164297.9901
[98,] 91437.6969 368028.3966
[99,] -42069.6757 91437.6969
[100,] -97365.9034 -42069.6757
[101,] -43810.1211 -97365.9034
[102,] -51991.9092 -43810.1211
[103,] 117123.1199 -51991.9092
[104,] 29944.0083 117123.1199
[105,] -97531.6477 29944.0083
[106,] -266971.0776 -97531.6477
[107,] 128424.8455 -266971.0776
[108,] -91738.3138 128424.8455
[109,] -100438.5848 -91738.3138
[110,] 96922.0217 -100438.5848
[111,] -254576.2834 96922.0217
[112,] -179297.6440 -254576.2834
[113,] -118452.4114 -179297.6440
[114,] -14628.1591 -118452.4114
[115,] 796.5017 -14628.1591
[116,] 88509.4743 796.5017
[117,] 123514.9171 88509.4743
[118,] -161398.0342 123514.9171
[119,] 120703.3278 -161398.0342
[120,] -196652.6095 120703.3278
[121,] 31573.9512 -196652.6095
[122,] -200582.8122 31573.9512
[123,] 76243.4774 -200582.8122
[124,] -70105.2464 76243.4774
[125,] -195592.4510 -70105.2464
[126,] 23944.4897 -195592.4510
[127,] -549.6192 23944.4897
[128,] -5360.9040 -549.6192
[129,] 23944.4897 -5360.9040
[130,] 15113.5220 23944.4897
[131,] 130137.6496 15113.5220
[132,] 23944.4897 130137.6496
[133,] -53630.8243 23944.4897
[134,] 23944.4897 -53630.8243
[135,] 34400.4677 23944.4897
[136,] 15071.0163 34400.4677
[137,] 23944.4897 15071.0163
[138,] 23944.4897 23944.4897
[139,] 23944.4897 23944.4897
[140,] 24850.1705 23944.4897
[141,] 23944.4897 24850.1705
[142,] 10516.6532 23944.4897
[143,] 15264.9363 10516.6532
[144,] 18116.0680 15264.9363
[145,] 12172.3994 18116.0680
[146,] 23944.4897 12172.3994
[147,] 9954.9738 23944.4897
[148,] 23944.4897 9954.9738
[149,] 26956.4369 23944.4897
[150,] 9263.7127 26956.4369
[151,] 10444.8738 9263.7127
[152,] 3095.7455 10444.8738
[153,] -2841.6848 3095.7455
[154,] 23274.9586 -2841.6848
[155,] 23944.4897 23274.9586
[156,] 23944.4897 23944.4897
[157,] 21721.8609 23944.4897
[158,] 33740.4414 21721.8609
[159,] 23944.4897 33740.4414
[160,] 23944.4897 23944.4897
[161,] 23944.4897 23944.4897
[162,] 23944.4897 23944.4897
[163,] 67491.0308 23944.4897
[164,] 47608.3635 67491.0308
[165,] 23944.4897 47608.3635
[166,] 23944.4897 23944.4897
[167,] 55212.5717 23944.4897
[168,] 7265.4943 55212.5717
[169,] 23944.4897 7265.4943
[170,] 23944.4897 23944.4897
[171,] 20907.2606 23944.4897
[172,] 23944.4897 20907.2606
[173,] 11860.3540 23944.4897
[174,] 23944.4897 11860.3540
[175,] 18500.2461 23944.4897
[176,] 23944.4897 18500.2461
[177,] -28757.0335 23944.4897
[178,] 23944.4897 -28757.0335
[179,] 23944.4897 23944.4897
[180,] 83867.6908 23944.4897
[181,] -27383.2329 83867.6908
[182,] 93898.3554 -27383.2329
[183,] -86816.9886 93898.3554
[184,] -34571.1742 -86816.9886
[185,] -43410.3807 -34571.1742
[186,] -4819.4832 -43410.3807
[187,] 14322.6612 -4819.4832
[188,] -1823.4228 14322.6612
[189,] 23944.4897 -1823.4228
[190,] 8647.1684 23944.4897
[191,] 23944.4897 8647.1684
[192,] 15946.0910 23944.4897
[193,] 23944.4897 15946.0910
[194,] 16616.4565 23944.4897
[195,] -31164.0924 16616.4565
[196,] -162580.2171 -31164.0924
[197,] -31164.0924 -162580.2171
[198,] -31164.0924 -31164.0924
[199,] -69414.3865 -31164.0924
[200,] -31164.0924 -69414.3865
[201,] -31164.0924 -31164.0924
[202,] -31164.0924 -31164.0924
[203,] -24625.9356 -31164.0924
[204,] -26498.6333 -24625.9356
[205,] -45596.5627 -26498.6333
[206,] -25412.0680 -45596.5627
[207,] -31164.0924 -25412.0680
[208,] -28390.8870 -31164.0924
[209,] -31164.0924 -28390.8870
[210,] 23944.4897 -31164.0924
[211,] 23944.4897 23944.4897
[212,] 56991.7930 23944.4897
[213,] 16129.6807 56991.7930
[214,] 17321.8019 16129.6807
[215,] 15604.9920 17321.8019
[216,] 20018.5381 15604.9920
[217,] 65743.9018 20018.5381
[218,] 15604.9920 65743.9018
[219,] -14771.6417 15604.9920
[220,] -7888.7766 -14771.6417
[221,] 23944.4897 -7888.7766
[222,] 24859.2736 23944.4897
[223,] -1074.0034 24859.2736
[224,] 25037.9354 -1074.0034
[225,] -386.2961 25037.9354
[226,] 23944.4897 -386.2961
[227,] 23944.4897 23944.4897
[228,] 19715.0099 23944.4897
[229,] 23944.4897 19715.0099
[230,] 11162.2202 23944.4897
[231,] 23944.4897 11162.2202
[232,] 23944.4897 23944.4897
[233,] 17857.3750 23944.4897
[234,] 23944.4897 17857.3750
[235,] 26439.4266 23944.4897
[236,] 23944.4897 26439.4266
[237,] 14159.2402 23944.4897
[238,] 23944.4897 14159.2402
[239,] 23944.4897 23944.4897
[240,] 23944.4897 23944.4897
[241,] 25092.1689 23944.4897
[242,] 23944.4897 25092.1689
[243,] 13831.5253 23944.4897
[244,] 31776.3908 13831.5253
[245,] -31164.0924 31776.3908
[246,] -31164.0924 -31164.0924
[247,] -16839.4265 -31164.0924
[248,] -31164.0924 -16839.4265
[249,] -31164.0924 -31164.0924
[250,] -12268.5111 -31164.0924
[251,] -31164.0924 -12268.5111
[252,] 18088.9535 -31164.0924
[253,] 36297.4716 18088.9535
[254,] 1066.2229 36297.4716
[255,] 23944.4897 1066.2229
[256,] 24008.5743 23944.4897
[257,] -3218.5157 24008.5743
[258,] -31164.0924 -3218.5157
[259,] 23944.4897 -31164.0924
[260,] -31164.0924 23944.4897
[261,] -31164.0924 -31164.0924
[262,] 24904.4705 -31164.0924
[263,] 21890.4878 24904.4705
[264,] -31164.0924 21890.4878
[265,] 23944.4897 -31164.0924
[266,] -31164.0924 23944.4897
[267,] 80047.1915 -31164.0924
[268,] -31164.0924 80047.1915
[269,] -55236.0503 -31164.0924
[270,] -31164.0924 -55236.0503
[271,] -31164.0924 -31164.0924
[272,] -31164.0924 -31164.0924
[273,] -28665.2822 -31164.0924
[274,] 23944.4897 -28665.2822
[275,] 57477.8478 23944.4897
[276,] 22102.3950 57477.8478
[277,] 23944.4897 22102.3950
[278,] 57933.9111 23944.4897
[279,] 56208.6106 57933.9111
[280,] -31164.0924 56208.6106
[281,] 23944.4897 -31164.0924
[282,] 102833.7758 23944.4897
[283,] 23944.4897 102833.7758
[284,] 21869.4117 23944.4897
[285,] 23944.4897 21869.4117
[286,] -31164.0924 23944.4897
[287,] 24714.2721 -31164.0924
[288,] -31164.0924 24714.2721
[289,] 121394.6428 -31164.0924
[290,] -31164.0924 121394.6428
[291,] 44696.2917 -31164.0924
[292,] 41880.9308 44696.2917
[293,] 23944.4897 41880.9308
[294,] 11990.8023 23944.4897
[295,] 13311.6374 11990.8023
[296,] 23944.4897 13311.6374
[297,] -31164.0924 23944.4897
[298,] 23944.4897 -31164.0924
[299,] -13583.2499 23944.4897
[300,] -39224.8440 -13583.2499
[301,] 8133.8845 -39224.8440
[302,] 23944.4897 8133.8845
[303,] -31164.0924 23944.4897
[304,] 126341.8884 -31164.0924
[305,] -2952.0544 126341.8884
[306,] -31164.0924 -2952.0544
[307,] -74616.6713 -31164.0924
[308,] 8317.1496 -74616.6713
[309,] -31164.0924 8317.1496
[310,] -31164.0924 -31164.0924
[311,] 96925.5331 -31164.0924
[312,] -31164.0924 96925.5331
[313,] -61722.3414 -31164.0924
[314,] 1172.4179 -61722.3414
[315,] 59855.7592 1172.4179
[316,] 25344.0301 59855.7592
[317,] -28662.1095 25344.0301
[318,] 18498.3189 -28662.1095
[319,] 23944.4897 18498.3189
[320,] 23944.4897 23944.4897
[321,] 23944.4897 23944.4897
[322,] 23944.4897 23944.4897
[323,] -31164.0924 23944.4897
[324,] 23944.4897 -31164.0924
[325,] 21871.6492 23944.4897
[326,] -31164.0924 21871.6492
[327,] -31164.0924 -31164.0924
[328,] 23944.4897 -31164.0924
[329,] 23944.4897 23944.4897
[330,] 23944.4897 23944.4897
[331,] 23944.4897 23944.4897
[332,] 23944.4897 23944.4897
[333,] 18436.9685 23944.4897
[334,] 24916.6033 18436.9685
[335,] 23944.4897 24916.6033
[336,] 23283.8362 23944.4897
[337,] 23944.4897 23283.8362
[338,] 22470.2156 23944.4897
[339,] 12543.2403 22470.2156
[340,] 766.5812 12543.2403
[341,] 23944.4897 766.5812
[342,] 23944.4897 23944.4897
[343,] 23944.4897 23944.4897
[344,] 23944.4897 23944.4897
[345,] 23944.4897 23944.4897
[346,] 12241.4031 23944.4897
[347,] 23944.4897 12241.4031
[348,] 23944.4897 23944.4897
[349,] 23944.4897 23944.4897
[350,] 40236.0807 23944.4897
[351,] 19881.8725 40236.0807
[352,] 16754.7417 19881.8725
[353,] 22233.3625 16754.7417
[354,] 23944.4897 22233.3625
[355,] 23944.4897 23944.4897
[356,] 94666.5242 23944.4897
[357,] 23944.4897 94666.5242
[358,] 23944.4897 23944.4897
[359,] 23944.4897 23944.4897
[360,] 17534.0406 23944.4897
[361,] 23944.4897 17534.0406
[362,] 23944.4897 23944.4897
[363,] -37653.4017 23944.4897
[364,] 7735.2309 -37653.4017
[365,] 23944.4897 7735.2309
[366,] 23944.4897 23944.4897
[367,] 23944.4897 23944.4897
[368,] 23944.4897 23944.4897
[369,] 23944.4897 23944.4897
[370,] 23944.4897 23944.4897
[371,] 88363.2415 23944.4897
[372,] 38934.9047 88363.2415
[373,] 6811.2441 38934.9047
[374,] 23944.4897 6811.2441
[375,] 23944.4897 23944.4897
[376,] 28078.4082 23944.4897
[377,] 14034.6561 28078.4082
[378,] 23944.4897 14034.6561
[379,] 23944.4897 23944.4897
[380,] 23944.4897 23944.4897
[381,] 108817.1919 23944.4897
[382,] 164863.0184 108817.1919
[383,] 23944.4897 164863.0184
[384,] 18419.9700 23944.4897
[385,] -97878.2616 18419.9700
[386,] 27285.8956 -97878.2616
[387,] 10634.9196 27285.8956
[388,] 16119.9896 10634.9196
[389,] 8879.0068 16119.9896
[390,] 15693.2771 8879.0068
[391,] -936.5007 15693.2771
[392,] -18866.5525 -936.5007
[393,] 74608.9461 -18866.5525
[394,] -99649.4056 74608.9461
[395,] -6594.9425 -99649.4056
[396,] -82488.7376 -6594.9425
[397,] -697508.8098 -82488.7376
[398,] 9888.3789 -697508.8098
[399,] 121170.3199 9888.3789
[400,] 43695.7911 121170.3199
[401,] 28536.0182 43695.7911
[402,] 1668.2186 28536.0182
[403,] -88304.8133 1668.2186
[404,] 1638.0231 -88304.8133
[405,] 14191.8476 1638.0231
[406,] 34921.3543 14191.8476
[407,] -10054.1141 34921.3543
[408,] -12874.9790 -10054.1141
[409,] -5984.2293 -12874.9790
[410,] -82507.1098 -5984.2293
[411,] -113565.0487 -82507.1098
[412,] -75888.7341 -113565.0487
[413,] -122992.5639 -75888.7341
[414,] 6515.1212 -122992.5639
[415,] -2245.6856 6515.1212
[416,] 160267.5784 -2245.6856
[417,] -103626.4598 160267.5784
[418,] -18465.1137 -103626.4598
[419,] -8329.4936 -18465.1137
[420,] -40960.8889 -8329.4936
[421,] -216291.7420 -40960.8889
[422,] -497573.8832 -216291.7420
[423,] -414112.1634 -497573.8832
[424,] -109008.0932 -414112.1634
[425,] -249580.5640 -109008.0932
[426,] -311658.1463 -249580.5640
[427,] -1284038.3256 -311658.1463
[428,] -184370.2558 -1284038.3256
[429,] -134950.8637 -184370.2558
[430,] -377355.7461 -134950.8637
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 2557837.5403 -277397.7292
2 1524849.5589 2557837.5403
3 -872756.8463 1524849.5589
4 -260723.9355 -872756.8463
5 499912.6898 -260723.9355
6 50803.3342 499912.6898
7 233096.1568 50803.3342
8 546213.5869 233096.1568
9 765744.5259 546213.5869
10 118815.3978 765744.5259
11 415065.5112 118815.3978
12 363390.6666 415065.5112
13 -1106733.5215 363390.6666
14 81428.2065 -1106733.5215
15 598941.0692 81428.2065
16 -285286.3107 598941.0692
17 -223449.8448 -285286.3107
18 222456.0823 -223449.8448
19 -781630.3797 222456.0823
20 382247.6700 -781630.3797
21 123835.2468 382247.6700
22 -94773.1470 123835.2468
23 -1025742.4420 -94773.1470
24 -150822.8397 -1025742.4420
25 72550.9570 -150822.8397
26 -272316.0688 72550.9570
27 -248027.7430 -272316.0688
28 398664.0991 -248027.7430
29 534699.5439 398664.0991
30 269029.6805 534699.5439
31 22432.4533 269029.6805
32 232558.4609 22432.4533
33 129733.1617 232558.4609
34 495152.7680 129733.1617
35 -677577.9129 495152.7680
36 156925.7774 -677577.9129
37 91863.3885 156925.7774
38 -38270.7135 91863.3885
39 -342467.4789 -38270.7135
40 -230150.3517 -342467.4789
41 -128825.6111 -230150.3517
42 186076.1321 -128825.6111
43 -165671.9290 186076.1321
44 37071.3337 -165671.9290
45 -86017.8086 37071.3337
46 231502.1696 -86017.8086
47 -377790.7035 231502.1696
48 -1374.9665 -377790.7035
49 -195568.0272 -1374.9665
50 -30126.9217 -195568.0272
51 -329273.4078 -30126.9217
52 -340355.4879 -329273.4078
53 -16799.7902 -340355.4879
54 -506760.5780 -16799.7902
55 34203.7748 -506760.5780
56 180846.5192 34203.7748
57 59780.1536 180846.5192
58 23037.7797 59780.1536
59 -163328.4059 23037.7797
60 18400.9273 -163328.4059
61 -113341.5416 18400.9273
62 9576.2456 -113341.5416
63 -195793.0162 9576.2456
64 24724.8221 -195793.0162
65 563636.9711 24724.8221
66 24732.4057 563636.9711
67 -46375.2413 24732.4057
68 -73770.5888 -46375.2413
69 -204071.2810 -73770.5888
70 -145607.0407 -204071.2810
71 32454.6517 -145607.0407
72 -61495.9718 32454.6517
73 18524.9487 -61495.9718
74 -5115.8936 18524.9487
75 -106754.6017 -5115.8936
76 19465.0654 -106754.6017
77 18368.4564 19465.0654
78 115684.9807 18368.4564
79 195574.6939 115684.9807
80 -18171.0815 195574.6939
81 -76936.2201 -18171.0815
82 36834.4012 -76936.2201
83 -26100.2435 36834.4012
84 93365.5534 -26100.2435
85 -137127.4478 93365.5534
86 -117531.4463 -137127.4478
87 71990.4549 -117531.4463
88 -37266.1893 71990.4549
89 82710.6706 -37266.1893
90 29855.4842 82710.6706
91 -76301.2867 29855.4842
92 77576.4644 -76301.2867
93 -60184.6443 77576.4644
94 -774.7735 -60184.6443
95 -12507.9392 -774.7735
96 164297.9901 -12507.9392
97 368028.3966 164297.9901
98 91437.6969 368028.3966
99 -42069.6757 91437.6969
100 -97365.9034 -42069.6757
101 -43810.1211 -97365.9034
102 -51991.9092 -43810.1211
103 117123.1199 -51991.9092
104 29944.0083 117123.1199
105 -97531.6477 29944.0083
106 -266971.0776 -97531.6477
107 128424.8455 -266971.0776
108 -91738.3138 128424.8455
109 -100438.5848 -91738.3138
110 96922.0217 -100438.5848
111 -254576.2834 96922.0217
112 -179297.6440 -254576.2834
113 -118452.4114 -179297.6440
114 -14628.1591 -118452.4114
115 796.5017 -14628.1591
116 88509.4743 796.5017
117 123514.9171 88509.4743
118 -161398.0342 123514.9171
119 120703.3278 -161398.0342
120 -196652.6095 120703.3278
121 31573.9512 -196652.6095
122 -200582.8122 31573.9512
123 76243.4774 -200582.8122
124 -70105.2464 76243.4774
125 -195592.4510 -70105.2464
126 23944.4897 -195592.4510
127 -549.6192 23944.4897
128 -5360.9040 -549.6192
129 23944.4897 -5360.9040
130 15113.5220 23944.4897
131 130137.6496 15113.5220
132 23944.4897 130137.6496
133 -53630.8243 23944.4897
134 23944.4897 -53630.8243
135 34400.4677 23944.4897
136 15071.0163 34400.4677
137 23944.4897 15071.0163
138 23944.4897 23944.4897
139 23944.4897 23944.4897
140 24850.1705 23944.4897
141 23944.4897 24850.1705
142 10516.6532 23944.4897
143 15264.9363 10516.6532
144 18116.0680 15264.9363
145 12172.3994 18116.0680
146 23944.4897 12172.3994
147 9954.9738 23944.4897
148 23944.4897 9954.9738
149 26956.4369 23944.4897
150 9263.7127 26956.4369
151 10444.8738 9263.7127
152 3095.7455 10444.8738
153 -2841.6848 3095.7455
154 23274.9586 -2841.6848
155 23944.4897 23274.9586
156 23944.4897 23944.4897
157 21721.8609 23944.4897
158 33740.4414 21721.8609
159 23944.4897 33740.4414
160 23944.4897 23944.4897
161 23944.4897 23944.4897
162 23944.4897 23944.4897
163 67491.0308 23944.4897
164 47608.3635 67491.0308
165 23944.4897 47608.3635
166 23944.4897 23944.4897
167 55212.5717 23944.4897
168 7265.4943 55212.5717
169 23944.4897 7265.4943
170 23944.4897 23944.4897
171 20907.2606 23944.4897
172 23944.4897 20907.2606
173 11860.3540 23944.4897
174 23944.4897 11860.3540
175 18500.2461 23944.4897
176 23944.4897 18500.2461
177 -28757.0335 23944.4897
178 23944.4897 -28757.0335
179 23944.4897 23944.4897
180 83867.6908 23944.4897
181 -27383.2329 83867.6908
182 93898.3554 -27383.2329
183 -86816.9886 93898.3554
184 -34571.1742 -86816.9886
185 -43410.3807 -34571.1742
186 -4819.4832 -43410.3807
187 14322.6612 -4819.4832
188 -1823.4228 14322.6612
189 23944.4897 -1823.4228
190 8647.1684 23944.4897
191 23944.4897 8647.1684
192 15946.0910 23944.4897
193 23944.4897 15946.0910
194 16616.4565 23944.4897
195 -31164.0924 16616.4565
196 -162580.2171 -31164.0924
197 -31164.0924 -162580.2171
198 -31164.0924 -31164.0924
199 -69414.3865 -31164.0924
200 -31164.0924 -69414.3865
201 -31164.0924 -31164.0924
202 -31164.0924 -31164.0924
203 -24625.9356 -31164.0924
204 -26498.6333 -24625.9356
205 -45596.5627 -26498.6333
206 -25412.0680 -45596.5627
207 -31164.0924 -25412.0680
208 -28390.8870 -31164.0924
209 -31164.0924 -28390.8870
210 23944.4897 -31164.0924
211 23944.4897 23944.4897
212 56991.7930 23944.4897
213 16129.6807 56991.7930
214 17321.8019 16129.6807
215 15604.9920 17321.8019
216 20018.5381 15604.9920
217 65743.9018 20018.5381
218 15604.9920 65743.9018
219 -14771.6417 15604.9920
220 -7888.7766 -14771.6417
221 23944.4897 -7888.7766
222 24859.2736 23944.4897
223 -1074.0034 24859.2736
224 25037.9354 -1074.0034
225 -386.2961 25037.9354
226 23944.4897 -386.2961
227 23944.4897 23944.4897
228 19715.0099 23944.4897
229 23944.4897 19715.0099
230 11162.2202 23944.4897
231 23944.4897 11162.2202
232 23944.4897 23944.4897
233 17857.3750 23944.4897
234 23944.4897 17857.3750
235 26439.4266 23944.4897
236 23944.4897 26439.4266
237 14159.2402 23944.4897
238 23944.4897 14159.2402
239 23944.4897 23944.4897
240 23944.4897 23944.4897
241 25092.1689 23944.4897
242 23944.4897 25092.1689
243 13831.5253 23944.4897
244 31776.3908 13831.5253
245 -31164.0924 31776.3908
246 -31164.0924 -31164.0924
247 -16839.4265 -31164.0924
248 -31164.0924 -16839.4265
249 -31164.0924 -31164.0924
250 -12268.5111 -31164.0924
251 -31164.0924 -12268.5111
252 18088.9535 -31164.0924
253 36297.4716 18088.9535
254 1066.2229 36297.4716
255 23944.4897 1066.2229
256 24008.5743 23944.4897
257 -3218.5157 24008.5743
258 -31164.0924 -3218.5157
259 23944.4897 -31164.0924
260 -31164.0924 23944.4897
261 -31164.0924 -31164.0924
262 24904.4705 -31164.0924
263 21890.4878 24904.4705
264 -31164.0924 21890.4878
265 23944.4897 -31164.0924
266 -31164.0924 23944.4897
267 80047.1915 -31164.0924
268 -31164.0924 80047.1915
269 -55236.0503 -31164.0924
270 -31164.0924 -55236.0503
271 -31164.0924 -31164.0924
272 -31164.0924 -31164.0924
273 -28665.2822 -31164.0924
274 23944.4897 -28665.2822
275 57477.8478 23944.4897
276 22102.3950 57477.8478
277 23944.4897 22102.3950
278 57933.9111 23944.4897
279 56208.6106 57933.9111
280 -31164.0924 56208.6106
281 23944.4897 -31164.0924
282 102833.7758 23944.4897
283 23944.4897 102833.7758
284 21869.4117 23944.4897
285 23944.4897 21869.4117
286 -31164.0924 23944.4897
287 24714.2721 -31164.0924
288 -31164.0924 24714.2721
289 121394.6428 -31164.0924
290 -31164.0924 121394.6428
291 44696.2917 -31164.0924
292 41880.9308 44696.2917
293 23944.4897 41880.9308
294 11990.8023 23944.4897
295 13311.6374 11990.8023
296 23944.4897 13311.6374
297 -31164.0924 23944.4897
298 23944.4897 -31164.0924
299 -13583.2499 23944.4897
300 -39224.8440 -13583.2499
301 8133.8845 -39224.8440
302 23944.4897 8133.8845
303 -31164.0924 23944.4897
304 126341.8884 -31164.0924
305 -2952.0544 126341.8884
306 -31164.0924 -2952.0544
307 -74616.6713 -31164.0924
308 8317.1496 -74616.6713
309 -31164.0924 8317.1496
310 -31164.0924 -31164.0924
311 96925.5331 -31164.0924
312 -31164.0924 96925.5331
313 -61722.3414 -31164.0924
314 1172.4179 -61722.3414
315 59855.7592 1172.4179
316 25344.0301 59855.7592
317 -28662.1095 25344.0301
318 18498.3189 -28662.1095
319 23944.4897 18498.3189
320 23944.4897 23944.4897
321 23944.4897 23944.4897
322 23944.4897 23944.4897
323 -31164.0924 23944.4897
324 23944.4897 -31164.0924
325 21871.6492 23944.4897
326 -31164.0924 21871.6492
327 -31164.0924 -31164.0924
328 23944.4897 -31164.0924
329 23944.4897 23944.4897
330 23944.4897 23944.4897
331 23944.4897 23944.4897
332 23944.4897 23944.4897
333 18436.9685 23944.4897
334 24916.6033 18436.9685
335 23944.4897 24916.6033
336 23283.8362 23944.4897
337 23944.4897 23283.8362
338 22470.2156 23944.4897
339 12543.2403 22470.2156
340 766.5812 12543.2403
341 23944.4897 766.5812
342 23944.4897 23944.4897
343 23944.4897 23944.4897
344 23944.4897 23944.4897
345 23944.4897 23944.4897
346 12241.4031 23944.4897
347 23944.4897 12241.4031
348 23944.4897 23944.4897
349 23944.4897 23944.4897
350 40236.0807 23944.4897
351 19881.8725 40236.0807
352 16754.7417 19881.8725
353 22233.3625 16754.7417
354 23944.4897 22233.3625
355 23944.4897 23944.4897
356 94666.5242 23944.4897
357 23944.4897 94666.5242
358 23944.4897 23944.4897
359 23944.4897 23944.4897
360 17534.0406 23944.4897
361 23944.4897 17534.0406
362 23944.4897 23944.4897
363 -37653.4017 23944.4897
364 7735.2309 -37653.4017
365 23944.4897 7735.2309
366 23944.4897 23944.4897
367 23944.4897 23944.4897
368 23944.4897 23944.4897
369 23944.4897 23944.4897
370 23944.4897 23944.4897
371 88363.2415 23944.4897
372 38934.9047 88363.2415
373 6811.2441 38934.9047
374 23944.4897 6811.2441
375 23944.4897 23944.4897
376 28078.4082 23944.4897
377 14034.6561 28078.4082
378 23944.4897 14034.6561
379 23944.4897 23944.4897
380 23944.4897 23944.4897
381 108817.1919 23944.4897
382 164863.0184 108817.1919
383 23944.4897 164863.0184
384 18419.9700 23944.4897
385 -97878.2616 18419.9700
386 27285.8956 -97878.2616
387 10634.9196 27285.8956
388 16119.9896 10634.9196
389 8879.0068 16119.9896
390 15693.2771 8879.0068
391 -936.5007 15693.2771
392 -18866.5525 -936.5007
393 74608.9461 -18866.5525
394 -99649.4056 74608.9461
395 -6594.9425 -99649.4056
396 -82488.7376 -6594.9425
397 -697508.8098 -82488.7376
398 9888.3789 -697508.8098
399 121170.3199 9888.3789
400 43695.7911 121170.3199
401 28536.0182 43695.7911
402 1668.2186 28536.0182
403 -88304.8133 1668.2186
404 1638.0231 -88304.8133
405 14191.8476 1638.0231
406 34921.3543 14191.8476
407 -10054.1141 34921.3543
408 -12874.9790 -10054.1141
409 -5984.2293 -12874.9790
410 -82507.1098 -5984.2293
411 -113565.0487 -82507.1098
412 -75888.7341 -113565.0487
413 -122992.5639 -75888.7341
414 6515.1212 -122992.5639
415 -2245.6856 6515.1212
416 160267.5784 -2245.6856
417 -103626.4598 160267.5784
418 -18465.1137 -103626.4598
419 -8329.4936 -18465.1137
420 -40960.8889 -8329.4936
421 -216291.7420 -40960.8889
422 -497573.8832 -216291.7420
423 -414112.1634 -497573.8832
424 -109008.0932 -414112.1634
425 -249580.5640 -109008.0932
426 -311658.1463 -249580.5640
427 -1284038.3256 -311658.1463
428 -184370.2558 -1284038.3256
429 -134950.8637 -184370.2558
430 -377355.7461 -134950.8637
> 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/7ni0r1291358215.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/8ni0r1291358215.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/9ni0r1291358215.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/10gric1291358215.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/111ayh1291358215.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/125sf51291358215.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/13jkue1291358215.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/14m2t21291358215.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/1583981291358215.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/16blqw1291358215.tab")
+ }
>
> try(system("convert tmp/198301291358215.ps tmp/198301291358215.png",intern=TRUE))
character(0)
> try(system("convert tmp/2kh231291358215.ps tmp/2kh231291358215.png",intern=TRUE))
character(0)
> try(system("convert tmp/3kh231291358215.ps tmp/3kh231291358215.png",intern=TRUE))
character(0)
> try(system("convert tmp/4kh231291358215.ps tmp/4kh231291358215.png",intern=TRUE))
character(0)
> try(system("convert tmp/5kh231291358215.ps tmp/5kh231291358215.png",intern=TRUE))
character(0)
> try(system("convert tmp/6icqu1291358215.ps tmp/6icqu1291358215.png",intern=TRUE))
character(0)
> try(system("convert tmp/7ni0r1291358215.ps tmp/7ni0r1291358215.png",intern=TRUE))
character(0)
> try(system("convert tmp/8ni0r1291358215.ps tmp/8ni0r1291358215.png",intern=TRUE))
character(0)
> try(system("convert tmp/9ni0r1291358215.ps tmp/9ni0r1291358215.png",intern=TRUE))
character(0)
> try(system("convert tmp/10gric1291358215.ps tmp/10gric1291358215.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
15.740 3.059 16.373