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(162556
+ ,1081
+ ,213118
+ ,6282154
+ ,29790
+ ,309
+ ,81767
+ ,4321023
+ ,87550
+ ,458
+ ,153198
+ ,4111912
+ ,84738
+ ,588
+ ,-26007
+ ,223193
+ ,54660
+ ,302
+ ,126942
+ ,1491348
+ ,42634
+ ,156
+ ,157214
+ ,1629616
+ ,40949
+ ,481
+ ,129352
+ ,1398893
+ ,45187
+ ,353
+ ,234817
+ ,1926517
+ ,37704
+ ,452
+ ,60448
+ ,983660
+ ,16275
+ ,109
+ ,47818
+ ,1443586
+ ,25830
+ ,115
+ ,245546
+ ,1073089
+ ,12679
+ ,110
+ ,48020
+ ,984885
+ ,18014
+ ,239
+ ,-1710
+ ,1405225
+ ,43556
+ ,247
+ ,32648
+ ,227132
+ ,24811
+ ,505
+ ,95350
+ ,929118
+ ,6575
+ ,159
+ ,151352
+ ,1071292
+ ,7123
+ ,109
+ ,288170
+ ,638830
+ ,21950
+ ,519
+ ,114337
+ ,856956
+ ,37597
+ ,248
+ ,37884
+ ,992426
+ ,17821
+ ,373
+ ,122844
+ ,444477
+ ,12988
+ ,119
+ ,82340
+ ,857217
+ ,22330
+ ,84
+ ,79801
+ ,711969
+ ,13326
+ ,102
+ ,165548
+ ,702380
+ ,16189
+ ,295
+ ,116384
+ ,358589
+ ,7146
+ ,105
+ ,134028
+ ,297978
+ ,15824
+ ,64
+ ,63838
+ ,585715
+ ,27664
+ ,282
+ ,74996
+ ,657954
+ ,11920
+ ,182
+ ,31080
+ ,209458
+ ,8568
+ ,37
+ ,32168
+ ,786690
+ ,14416
+ ,361
+ ,49857
+ ,439798
+ ,3369
+ ,28
+ ,87161
+ ,688779
+ ,11819
+ ,85
+ ,106113
+ ,574339
+ ,6984
+ ,45
+ ,80570
+ ,741409
+ ,4519
+ ,49
+ ,102129
+ ,597793
+ ,2220
+ ,22
+ ,301670
+ ,644190
+ ,18562
+ ,155
+ ,102313
+ ,377934
+ ,10327
+ ,91
+ ,88577
+ ,640273
+ ,5336
+ ,81
+ ,112477
+ ,697458
+ ,2365
+ ,79
+ ,191778
+ ,550608
+ ,4069
+ ,145
+ ,79804
+ ,207393
+ ,8636
+ ,855
+ ,128294
+ ,301607
+ ,13718
+ ,61
+ ,96448
+ ,345783
+ ,4525
+ ,226
+ ,93811
+ ,501749
+ ,6869
+ ,105
+ ,117520
+ ,379983
+ ,4628
+ ,62
+ ,69159
+ ,387475
+ ,3689
+ ,25
+ ,101792
+ ,377305
+ ,4891
+ ,217
+ ,210568
+ ,370837
+ ,7489
+ ,322
+ ,136996
+ ,430866
+ ,4901
+ ,84
+ ,121920
+ ,469107
+ ,2284
+ ,33
+ ,76403
+ ,194493
+ ,3160
+ ,108
+ ,108094
+ ,530670
+ ,4150
+ ,150
+ ,134759
+ ,518365
+ ,7285
+ ,115
+ ,188873
+ ,491303
+ ,1134
+ ,162
+ ,146216
+ ,527021
+ ,4658
+ ,158
+ ,156608
+ ,233773
+ ,2384
+ ,97
+ ,61348
+ ,405972
+ ,3748
+ ,9
+ ,50350
+ ,652925
+ ,5371
+ ,66
+ ,87720
+ ,446211
+ ,1285
+ ,107
+ ,99489
+ ,341340
+ ,9327
+ ,101
+ ,87419
+ ,387699
+ ,5565
+ ,47
+ ,94355
+ ,493408
+ ,1528
+ ,38
+ ,60326
+ ,146494
+ ,3122
+ ,34
+ ,94670
+ ,414462
+ ,7561
+ ,87
+ ,82425
+ ,364304
+ ,2675
+ ,79
+ ,59017
+ ,355178
+ ,13253
+ ,947
+ ,90829
+ ,357760
+ ,880
+ ,74
+ ,80791
+ ,261216
+ ,2053
+ ,53
+ ,100423
+ ,397144
+ ,1424
+ ,94
+ ,131116
+ ,374943
+ ,4036
+ ,63
+ ,100269
+ ,424898
+ ,3045
+ ,58
+ ,27330
+ ,202055
+ ,5119
+ ,49
+ ,39039
+ ,378525
+ ,1431
+ ,34
+ ,106885
+ ,310768
+ ,554
+ ,11
+ ,79285
+ ,325738
+ ,1975
+ ,35
+ ,118881
+ ,394510
+ ,1765
+ ,20
+ ,77623
+ ,247060
+ ,1012
+ ,47
+ ,114768
+ ,368078
+ ,810
+ ,43
+ ,74015
+ ,236761
+ ,1280
+ ,117
+ ,69465
+ ,312378
+ ,666
+ ,171
+ ,117869
+ ,339836
+ ,1380
+ ,26
+ ,60982
+ ,347385
+ ,4677
+ ,75
+ ,90131
+ ,426280
+ ,876
+ ,59
+ ,138971
+ ,352850
+ ,814
+ ,18
+ ,39625
+ ,301881
+ ,514
+ ,15
+ ,102725
+ ,377516
+ ,5692
+ ,72
+ ,64239
+ ,357312
+ ,3642
+ ,86
+ ,90262
+ ,458343
+ ,540
+ ,14
+ ,103960
+ ,354228
+ ,2099
+ ,64
+ ,106611
+ ,308636
+ ,567
+ ,11
+ ,103345
+ ,386212
+ ,2001
+ ,52
+ ,95551
+ ,393343
+ ,2949
+ ,41
+ ,82903
+ ,378509
+ ,2253
+ ,99
+ ,63593
+ ,452469
+ ,6533
+ ,75
+ ,126910
+ ,364839
+ ,1889
+ ,45
+ ,37527
+ ,358649
+ ,3055
+ ,43
+ ,60247
+ ,376641
+ ,272
+ ,8
+ ,112995
+ ,429112
+ ,1414
+ ,198
+ ,70184
+ ,330546
+ ,2564
+ ,22
+ ,130140
+ ,403560
+ ,1383
+ ,11
+ ,73221
+ ,317892
+ ,1261
+ ,33
+ ,76114
+ ,307528
+ ,975
+ ,23
+ ,90534
+ ,235133
+ ,3366
+ ,80
+ ,108479
+ ,299243
+ ,576
+ ,18
+ ,113761
+ ,314073
+ ,1686
+ ,40
+ ,68696
+ ,368186
+ ,746
+ ,23
+ ,71561
+ ,269661
+ ,3192
+ ,60
+ ,59831
+ ,125390
+ ,2045
+ ,20
+ ,97890
+ ,510834
+ ,5702
+ ,61
+ ,101481
+ ,321896
+ ,1932
+ ,36
+ ,72954
+ ,249898
+ ,936
+ ,30
+ ,67939
+ ,408881
+ ,3437
+ ,47
+ ,48022
+ ,158492
+ ,5131
+ ,71
+ ,86111
+ ,292154
+ ,2397
+ ,14
+ ,74020
+ ,289513
+ ,1389
+ ,9
+ ,57530
+ ,378049
+ ,1503
+ ,39
+ ,56364
+ ,343466
+ ,402
+ ,26
+ ,84990
+ ,332743
+ ,2239
+ ,21
+ ,88590
+ ,442882
+ ,2234
+ ,16
+ ,77200
+ ,214215
+ ,837
+ ,69
+ ,61262
+ ,315688
+ ,10579
+ ,92
+ ,110309
+ ,375195
+ ,875
+ ,14
+ ,67000
+ ,334280
+ ,1585
+ ,107
+ ,93099
+ ,355864
+ ,1659
+ ,29
+ ,107577
+ ,480382
+ ,2647
+ ,37
+ ,62920
+ ,353058
+ ,3294
+ ,23
+ ,75832
+ ,217193
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,94
+ ,7
+ ,60793
+ ,314533
+ ,422
+ ,28
+ ,57935
+ ,318056
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,34
+ ,8
+ ,60630
+ ,314353
+ ,1558
+ ,63
+ ,55637
+ ,369448
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,43
+ ,3
+ ,60887
+ ,312846
+ ,645
+ ,5
+ ,60720
+ ,312075
+ ,316
+ ,9
+ ,60505
+ ,315009
+ ,115
+ ,13
+ ,60945
+ ,318903
+ ,5
+ ,2
+ ,60720
+ ,314887
+ ,897
+ ,5
+ ,60720
+ ,314913
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,389
+ ,14
+ ,58990
+ ,325506
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1002
+ ,15
+ ,56750
+ ,298568
+ ,36
+ ,3
+ ,60894
+ ,315834
+ ,460
+ ,15
+ ,63346
+ ,329784
+ ,309
+ ,11
+ ,56535
+ ,312878
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,9
+ ,6
+ ,60835
+ ,314987
+ ,271
+ ,2
+ ,60720
+ ,325249
+ ,14
+ ,1
+ ,61016
+ ,315877
+ ,520
+ ,10
+ ,58650
+ ,291650
+ ,1766
+ ,73
+ ,60438
+ ,305959
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,458
+ ,11
+ ,58625
+ ,297765
+ ,20
+ ,3
+ ,60938
+ ,315245
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,98
+ ,2
+ ,61490
+ ,315236
+ ,405
+ ,7
+ ,60845
+ ,336425
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,483
+ ,27
+ ,60830
+ ,306268
+ ,454
+ ,51
+ ,63261
+ ,302187
+ ,47
+ ,3
+ ,60720
+ ,314882
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,757
+ ,19
+ ,45689
+ ,382712
+ ,4655
+ ,393
+ ,60720
+ ,341570
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,36
+ ,4
+ ,61564
+ ,312412
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,203
+ ,9
+ ,61938
+ ,309596
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,126
+ ,10
+ ,60951
+ ,315547
+ ,400
+ ,152
+ ,60720
+ ,313267
+ ,71
+ ,1
+ ,60745
+ ,316176
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,972
+ ,34
+ ,71642
+ ,359335
+ ,531
+ ,10
+ ,71641
+ ,330068
+ ,2461
+ ,57
+ ,55792
+ ,314289
+ ,378
+ ,52
+ ,71873
+ ,297413
+ ,23
+ ,5
+ ,62555
+ ,314806
+ ,638
+ ,14
+ ,60370
+ ,333210
+ ,2300
+ ,29
+ ,64873
+ ,352108
+ ,149
+ ,5
+ ,62041
+ ,313332
+ ,226
+ ,5
+ ,65745
+ ,291787
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,275
+ ,4
+ ,59500
+ ,318745
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,141
+ ,6
+ ,61630
+ ,315366
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,28
+ ,2
+ ,60890
+ ,315688
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,4980
+ ,91
+ ,113521
+ ,409642
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,472
+ ,20
+ ,80045
+ ,269587
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,203
+ ,27
+ ,50804
+ ,300962
+ ,496
+ ,17
+ ,87390
+ ,325479
+ ,10
+ ,2
+ ,61656
+ ,316155
+ ,63
+ ,4
+ ,65688
+ ,318574
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1136
+ ,32
+ ,48522
+ ,343613
+ ,265
+ ,31
+ ,60720
+ ,306948
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,267
+ ,32
+ ,57640
+ ,330059
+ ,474
+ ,20
+ ,61977
+ ,288985
+ ,534
+ ,7
+ ,62620
+ ,304485
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,15
+ ,8
+ ,60831
+ ,315688
+ ,397
+ ,28
+ ,60646
+ ,317736
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1866
+ ,29
+ ,56225
+ ,322331
+ ,288
+ ,4
+ ,60510
+ ,296656
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,3
+ ,2
+ ,60698
+ ,315354
+ ,468
+ ,21
+ ,60720
+ ,312161
+ ,20
+ ,2
+ ,60805
+ ,315576
+ ,278
+ ,26
+ ,61404
+ ,314922
+ ,61
+ ,14
+ ,60720
+ ,314551
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,192
+ ,4
+ ,65276
+ ,312339
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,317
+ ,9
+ ,63915
+ ,298700
+ ,738
+ ,10
+ ,60720
+ ,321376
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,368
+ ,17
+ ,61686
+ ,303230
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,2
+ ,1
+ ,60743
+ ,315487
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,53
+ ,6
+ ,60349
+ ,315793
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,94
+ ,3
+ ,61360
+ ,312887
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,24
+ ,8
+ ,59818
+ ,315637
+ ,2332
+ ,4
+ ,72680
+ ,324385
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,131
+ ,11
+ ,61808
+ ,308989
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,206
+ ,9
+ ,53110
+ ,296702
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,167
+ ,2
+ ,64245
+ ,307322
+ ,622
+ ,73
+ ,73007
+ ,304376
+ ,2328
+ ,94
+ ,82732
+ ,253588
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,365
+ ,8
+ ,54820
+ ,309560
+ ,364
+ ,35
+ ,47705
+ ,298466
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,226
+ ,12
+ ,72835
+ ,343929
+ ,307
+ ,15
+ ,58856
+ ,331955
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,188
+ ,11
+ ,77655
+ ,381180
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,138
+ ,6
+ ,69817
+ ,331420
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,125
+ ,12
+ ,60798
+ ,310201
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,282
+ ,30
+ ,62452
+ ,320016
+ ,335
+ ,33
+ ,64175
+ ,320398
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1324
+ ,117
+ ,67440
+ ,291841
+ ,176
+ ,28
+ ,68136
+ ,310670
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,249
+ ,72
+ ,56726
+ ,313491
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,333
+ ,13
+ ,70811
+ ,331323
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,601
+ ,6
+ ,60720
+ ,319210
+ ,30
+ ,4
+ ,62045
+ ,318098
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,249
+ ,62
+ ,54323
+ ,292754
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,165
+ ,24
+ ,62841
+ ,325176
+ ,453
+ ,21
+ ,81125
+ ,365959
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,53
+ ,14
+ ,59506
+ ,302409
+ ,382
+ ,21
+ ,59365
+ ,340968
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,30
+ ,4
+ ,60798
+ ,313164
+ ,290
+ ,2
+ ,58790
+ ,301164
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,366
+ ,53
+ ,61808
+ ,344425
+ ,2
+ ,9
+ ,60735
+ ,315394
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,209
+ ,13
+ ,64016
+ ,316647
+ ,384
+ ,22
+ ,54683
+ ,309836
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,365
+ ,83
+ ,87192
+ ,346611
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,49
+ ,8
+ ,64107
+ ,322031
+ ,3
+ ,4
+ ,60761
+ ,315656
+ ,133
+ ,14
+ ,65990
+ ,339445
+ ,32
+ ,1
+ ,59988
+ ,314964
+ ,368
+ ,17
+ ,61167
+ ,297141
+ ,1
+ ,6
+ ,60719
+ ,315372
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,22
+ ,2
+ ,60722
+ ,312502
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,96
+ ,5
+ ,60379
+ ,313729
+ ,1
+ ,2
+ ,60727
+ ,315388
+ ,314
+ ,5
+ ,60720
+ ,315371
+ ,844
+ ,78
+ ,60925
+ ,296139
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,26
+ ,1
+ ,60896
+ ,313880
+ ,125
+ ,13
+ ,59734
+ ,317698
+ ,304
+ ,15
+ ,62969
+ ,295580
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,621
+ ,48
+ ,60720
+ ,308256
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,119
+ ,6
+ ,59118
+ ,303677
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1595
+ ,17
+ ,60720
+ ,319369
+ ,312
+ ,14
+ ,58598
+ ,318690
+ ,60
+ ,10
+ ,61124
+ ,314049
+ ,587
+ ,12
+ ,59595
+ ,325699
+ ,135
+ ,2
+ ,62065
+ ,314210
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,514
+ ,52
+ ,78780
+ ,322378
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,4
+ ,60722
+ ,315398
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1763
+ ,24
+ ,61600
+ ,308336
+ ,180
+ ,11
+ ,59635
+ ,316386
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,218
+ ,21
+ ,60720
+ ,315553
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,448
+ ,40
+ ,59781
+ ,323361
+ ,227
+ ,9
+ ,76644
+ ,336639
+ ,174
+ ,1
+ ,64820
+ ,307424
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,121
+ ,24
+ ,56178
+ ,295370
+ ,607
+ ,11
+ ,60436
+ ,322340
+ ,2212
+ ,14
+ ,60720
+ ,319864
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,530
+ ,60
+ ,73433
+ ,317291
+ ,571
+ ,80
+ ,41477
+ ,280398
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,78
+ ,16
+ ,62700
+ ,317330
+ ,2489
+ ,40
+ ,67804
+ ,238125
+ ,131
+ ,6
+ ,59661
+ ,327071
+ ,923
+ ,8
+ ,58620
+ ,309038
+ ,72
+ ,3
+ ,60398
+ ,314210
+ ,572
+ ,16
+ ,58580
+ ,307930
+ ,397
+ ,10
+ ,62710
+ ,322327
+ ,450
+ ,8
+ ,59325
+ ,292136
+ ,622
+ ,7
+ ,60950
+ ,263276
+ ,694
+ ,8
+ ,68060
+ ,367655
+ ,3425
+ ,12
+ ,83620
+ ,283910
+ ,562
+ ,13
+ ,58456
+ ,283587
+ ,4917
+ ,42
+ ,52811
+ ,243650
+ ,1442
+ ,118
+ ,121173
+ ,438493
+ ,529
+ ,9
+ ,63870
+ ,296261
+ ,2126
+ ,138
+ ,21001
+ ,230621
+ ,1061
+ ,5
+ ,70415
+ ,304252
+ ,776
+ ,9
+ ,64230
+ ,333505
+ ,611
+ ,8
+ ,59190
+ ,296919
+ ,1526
+ ,25
+ ,69351
+ ,278990
+ ,592
+ ,7
+ ,64270
+ ,276898
+ ,1182
+ ,13
+ ,70694
+ ,327007
+ ,621
+ ,16
+ ,68005
+ ,317046
+ ,989
+ ,11
+ ,58930
+ ,304555
+ ,438
+ ,11
+ ,58320
+ ,298096
+ ,726
+ ,3
+ ,69980
+ ,231861
+ ,1303
+ ,61
+ ,69863
+ ,309422
+ ,7419
+ ,29
+ ,63255
+ ,286963
+ ,1164
+ ,17
+ ,57320
+ ,269753
+ ,3310
+ ,33
+ ,75230
+ ,448243
+ ,1920
+ ,15
+ ,79420
+ ,165404
+ ,965
+ ,3
+ ,73490
+ ,204325
+ ,3256
+ ,66
+ ,35250
+ ,407159
+ ,1135
+ ,17
+ ,62285
+ ,290476
+ ,1270
+ ,26
+ ,69206
+ ,275311
+ ,661
+ ,3
+ ,65920
+ ,246541
+ ,1013
+ ,2
+ ,69770
+ ,253468
+ ,2844
+ ,67
+ ,72683
+ ,240897
+ ,11528
+ ,70
+ ,-14545
+ ,-83265
+ ,6526
+ ,26
+ ,55830
+ ,-42143
+ ,2264
+ ,24
+ ,55174
+ ,272713
+ ,5109
+ ,97
+ ,67038
+ ,215362
+ ,3999
+ ,30
+ ,51252
+ ,42754
+ ,35624
+ ,223
+ ,157278
+ ,306275
+ ,9252
+ ,48
+ ,79510
+ ,253537
+ ,15236
+ ,90
+ ,77440
+ ,372631
+ ,18073
+ ,180
+ ,27284
+ ,-7170)
+ ,dim=c(4
+ ,431)
+ ,dimnames=list(c('Costs'
+ ,'Trades'
+ ,'Dividends'
+ ,'Wealth
')
+ ,1:431))
> y <- array(NA,dim=c(4,431),dimnames=list(c('Costs','Trades','Dividends','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 = '4'
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Prof. Dr. P. Wessa
> #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #Technical description: Write here your technical program description (don't use hard returns!)
> library(lattice)
> library(lmtest)
Loading required package: zoo
Attaching package: 'zoo'
The following object(s) are masked from package:base :
as.Date.numeric
> n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test
> par1 <- as.numeric(par1)
> x <- t(y)
> k <- length(x[1,])
> n <- length(x[,1])
> x1 <- cbind(x[,par1], x[,1:k!=par1])
> mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1])
> colnames(x1) <- mycolnames #colnames(x)[par1]
> x <- x1
> if (par3 == 'First Differences'){
+ x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep='')))
+ for (i in 1:n-1) {
+ for (j in 1:k) {
+ x2[i,j] <- x[i+1,j] - x[i,j]
+ }
+ }
+ x <- x2
+ }
> if (par2 == 'Include Monthly Dummies'){
+ x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep ='')))
+ for (i in 1:11){
+ x2[seq(i,n,12),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> if (par2 == 'Include Quarterly Dummies'){
+ x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep ='')))
+ for (i in 1:3){
+ x2[seq(i,n,4),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> k <- length(x[1,])
> if (par3 == 'Linear Trend'){
+ x <- cbind(x, c(1:n))
+ colnames(x)[k+1] <- 't'
+ }
> x
Wealth\r Costs Trades Dividends
1 6282154 162556 1081 213118
2 4321023 29790 309 81767
3 4111912 87550 458 153198
4 223193 84738 588 -26007
5 1491348 54660 302 126942
6 1629616 42634 156 157214
7 1398893 40949 481 129352
8 1926517 45187 353 234817
9 983660 37704 452 60448
10 1443586 16275 109 47818
11 1073089 25830 115 245546
12 984885 12679 110 48020
13 1405225 18014 239 -1710
14 227132 43556 247 32648
15 929118 24811 505 95350
16 1071292 6575 159 151352
17 638830 7123 109 288170
18 856956 21950 519 114337
19 992426 37597 248 37884
20 444477 17821 373 122844
21 857217 12988 119 82340
22 711969 22330 84 79801
23 702380 13326 102 165548
24 358589 16189 295 116384
25 297978 7146 105 134028
26 585715 15824 64 63838
27 657954 27664 282 74996
28 209458 11920 182 31080
29 786690 8568 37 32168
30 439798 14416 361 49857
31 688779 3369 28 87161
32 574339 11819 85 106113
33 741409 6984 45 80570
34 597793 4519 49 102129
35 644190 2220 22 301670
36 377934 18562 155 102313
37 640273 10327 91 88577
38 697458 5336 81 112477
39 550608 2365 79 191778
40 207393 4069 145 79804
41 301607 8636 855 128294
42 345783 13718 61 96448
43 501749 4525 226 93811
44 379983 6869 105 117520
45 387475 4628 62 69159
46 377305 3689 25 101792
47 370837 4891 217 210568
48 430866 7489 322 136996
49 469107 4901 84 121920
50 194493 2284 33 76403
51 530670 3160 108 108094
52 518365 4150 150 134759
53 491303 7285 115 188873
54 527021 1134 162 146216
55 233773 4658 158 156608
56 405972 2384 97 61348
57 652925 3748 9 50350
58 446211 5371 66 87720
59 341340 1285 107 99489
60 387699 9327 101 87419
61 493408 5565 47 94355
62 146494 1528 38 60326
63 414462 3122 34 94670
64 364304 7561 87 82425
65 355178 2675 79 59017
66 357760 13253 947 90829
67 261216 880 74 80791
68 397144 2053 53 100423
69 374943 1424 94 131116
70 424898 4036 63 100269
71 202055 3045 58 27330
72 378525 5119 49 39039
73 310768 1431 34 106885
74 325738 554 11 79285
75 394510 1975 35 118881
76 247060 1765 20 77623
77 368078 1012 47 114768
78 236761 810 43 74015
79 312378 1280 117 69465
80 339836 666 171 117869
81 347385 1380 26 60982
82 426280 4677 75 90131
83 352850 876 59 138971
84 301881 814 18 39625
85 377516 514 15 102725
86 357312 5692 72 64239
87 458343 3642 86 90262
88 354228 540 14 103960
89 308636 2099 64 106611
90 386212 567 11 103345
91 393343 2001 52 95551
92 378509 2949 41 82903
93 452469 2253 99 63593
94 364839 6533 75 126910
95 358649 1889 45 37527
96 376641 3055 43 60247
97 429112 272 8 112995
98 330546 1414 198 70184
99 403560 2564 22 130140
100 317892 1383 11 73221
101 307528 1261 33 76114
102 235133 975 23 90534
103 299243 3366 80 108479
104 314073 576 18 113761
105 368186 1686 40 68696
106 269661 746 23 71561
107 125390 3192 60 59831
108 510834 2045 20 97890
109 321896 5702 61 101481
110 249898 1932 36 72954
111 408881 936 30 67939
112 158492 3437 47 48022
113 292154 5131 71 86111
114 289513 2397 14 74020
115 378049 1389 9 57530
116 343466 1503 39 56364
117 332743 402 26 84990
118 442882 2239 21 88590
119 214215 2234 16 77200
120 315688 837 69 61262
121 375195 10579 92 110309
122 334280 875 14 67000
123 355864 1585 107 93099
124 480382 1659 29 107577
125 353058 2647 37 62920
126 217193 3294 23 75832
127 315380 0 0 60720
128 314533 94 7 60793
129 318056 422 28 57935
130 315380 0 0 60720
131 314353 34 8 60630
132 369448 1558 63 55637
133 315380 0 0 60720
134 312846 43 3 60887
135 312075 645 5 60720
136 315009 316 9 60505
137 318903 115 13 60945
138 314887 5 2 60720
139 314913 897 5 60720
140 315380 0 0 60720
141 325506 389 14 58990
142 315380 0 0 60720
143 298568 1002 15 56750
144 315834 36 3 60894
145 329784 460 15 63346
146 312878 309 11 56535
147 315380 0 0 60720
148 314987 9 6 60835
149 325249 271 2 60720
150 315877 14 1 61016
151 291650 520 10 58650
152 305959 1766 73 60438
153 315380 0 0 60720
154 297765 458 11 58625
155 315245 20 3 60938
156 315380 0 0 60720
157 315380 0 0 60720
158 315236 98 2 61490
159 336425 405 7 60845
160 315380 0 0 60720
161 315380 0 0 60720
162 315380 0 0 60720
163 315380 0 0 60720
164 306268 483 27 60830
165 302187 454 51 63261
166 314882 47 3 60720
167 315380 0 0 60720
168 382712 757 19 45689
169 341570 4655 393 60720
170 315380 0 0 60720
171 315380 0 0 60720
172 312412 36 4 61564
173 315380 0 0 60720
174 309596 203 9 61938
175 315380 0 0 60720
176 315547 126 10 60951
177 313267 400 152 60720
178 316176 71 1 60745
179 315380 0 0 60720
180 315380 0 0 60720
181 359335 972 34 71642
182 330068 531 10 71641
183 314289 2461 57 55792
184 297413 378 52 71873
185 314806 23 5 62555
186 333210 638 14 60370
187 352108 2300 29 64873
188 313332 149 5 62041
189 291787 226 5 65745
190 315380 0 0 60720
191 318745 275 4 59500
192 315380 0 0 60720
193 315366 141 6 61630
194 315380 0 0 60720
195 315688 28 2 60890
196 315380 0 0 60720
197 409642 4980 91 113521
198 315380 0 0 60720
199 315380 0 0 60720
200 269587 472 20 80045
201 315380 0 0 60720
202 315380 0 0 60720
203 315380 0 0 60720
204 300962 203 27 50804
205 325479 496 17 87390
206 316155 10 2 61656
207 318574 63 4 65688
208 315380 0 0 60720
209 343613 1136 32 48522
210 306948 265 31 60720
211 315380 0 0 60720
212 315380 0 0 60720
213 330059 267 32 57640
214 288985 474 20 61977
215 304485 534 7 62620
216 315380 0 0 60720
217 315688 15 8 60831
218 317736 397 28 60646
219 315380 0 0 60720
220 322331 1866 29 56225
221 296656 288 4 60510
222 315380 0 0 60720
223 315354 3 2 60698
224 312161 468 21 60720
225 315576 20 2 60805
226 314922 278 26 61404
227 314551 61 14 60720
228 315380 0 0 60720
229 312339 192 4 65276
230 315380 0 0 60720
231 298700 317 9 63915
232 321376 738 10 60720
233 315380 0 0 60720
234 303230 368 17 61686
235 315380 0 0 60720
236 315487 2 1 60743
237 315380 0 0 60720
238 315793 53 6 60349
239 315380 0 0 60720
240 315380 0 0 60720
241 315380 0 0 60720
242 312887 94 3 61360
243 315380 0 0 60720
244 315637 24 8 59818
245 324385 2332 4 72680
246 315380 0 0 60720
247 315380 0 0 60720
248 308989 131 11 61808
249 315380 0 0 60720
250 315380 0 0 60720
251 296702 206 9 53110
252 315380 0 0 60720
253 307322 167 2 64245
254 304376 622 73 73007
255 253588 2328 94 82732
256 315380 0 0 60720
257 309560 365 8 54820
258 298466 364 35 47705
259 315380 0 0 60720
260 315380 0 0 60720
261 315380 0 0 60720
262 315380 0 0 60720
263 343929 226 12 72835
264 331955 307 15 58856
265 315380 0 0 60720
266 315380 0 0 60720
267 315380 0 0 60720
268 381180 188 11 77655
269 315380 0 0 60720
270 331420 138 6 69817
271 315380 0 0 60720
272 315380 0 0 60720
273 315380 0 0 60720
274 310201 125 12 60798
275 315380 0 0 60720
276 320016 282 30 62452
277 320398 335 33 64175
278 315380 0 0 60720
279 291841 1324 117 67440
280 310670 176 28 68136
281 315380 0 0 60720
282 315380 0 0 60720
283 313491 249 72 56726
284 315380 0 0 60720
285 331323 333 13 70811
286 315380 0 0 60720
287 319210 601 6 60720
288 318098 30 4 62045
289 315380 0 0 60720
290 292754 249 62 54323
291 315380 0 0 60720
292 325176 165 24 62841
293 365959 453 21 81125
294 315380 0 0 60720
295 302409 53 14 59506
296 340968 382 21 59365
297 315380 0 0 60720
298 315380 0 0 60720
299 315380 0 0 60720
300 315380 0 0 60720
301 313164 30 4 60798
302 301164 290 2 58790
303 315380 0 0 60720
304 315380 0 0 60720
305 344425 366 53 61808
306 315394 2 9 60735
307 315380 0 0 60720
308 316647 209 13 64016
309 309836 384 22 54683
310 315380 0 0 60720
311 315380 0 0 60720
312 346611 365 83 87192
313 315380 0 0 60720
314 322031 49 8 64107
315 315656 3 4 60761
316 339445 133 14 65990
317 314964 32 1 59988
318 297141 368 17 61167
319 315372 1 6 60719
320 315380 0 0 60720
321 315380 0 0 60720
322 315380 0 0 60720
323 315380 0 0 60720
324 315380 0 0 60720
325 315380 0 0 60720
326 312502 22 2 60722
327 315380 0 0 60720
328 315380 0 0 60720
329 315380 0 0 60720
330 315380 0 0 60720
331 315380 0 0 60720
332 315380 0 0 60720
333 315380 0 0 60720
334 313729 96 5 60379
335 315388 1 2 60727
336 315371 314 5 60720
337 296139 844 78 60925
338 315380 0 0 60720
339 313880 26 1 60896
340 317698 125 13 59734
341 295580 304 15 62969
342 315380 0 0 60720
343 315380 0 0 60720
344 315380 0 0 60720
345 308256 621 48 60720
346 315380 0 0 60720
347 303677 119 6 59118
348 315380 0 0 60720
349 315380 0 0 60720
350 319369 1595 17 60720
351 318690 312 14 58598
352 314049 60 10 61124
353 325699 587 12 59595
354 314210 135 2 62065
355 315380 0 0 60720
356 315380 0 0 60720
357 322378 514 52 78780
358 315380 0 0 60720
359 315380 0 0 60720
360 315380 0 0 60720
361 315398 1 4 60722
362 315380 0 0 60720
363 315380 0 0 60720
364 308336 1763 24 61600
365 316386 180 11 59635
366 315380 0 0 60720
367 315380 0 0 60720
368 315380 0 0 60720
369 315380 0 0 60720
370 315553 218 21 60720
371 315380 0 0 60720
372 323361 448 40 59781
373 336639 227 9 76644
374 307424 174 1 64820
375 315380 0 0 60720
376 315380 0 0 60720
377 295370 121 24 56178
378 322340 607 11 60436
379 319864 2212 14 60720
380 315380 0 0 60720
381 315380 0 0 60720
382 317291 530 60 73433
383 280398 571 80 41477
384 315380 0 0 60720
385 317330 78 16 62700
386 238125 2489 40 67804
387 327071 131 6 59661
388 309038 923 8 58620
389 314210 72 3 60398
390 307930 572 16 58580
391 322327 397 10 62710
392 292136 450 8 59325
393 263276 622 7 60950
394 367655 694 8 68060
395 283910 3425 12 83620
396 283587 562 13 58456
397 243650 4917 42 52811
398 438493 1442 118 121173
399 296261 529 9 63870
400 230621 2126 138 21001
401 304252 1061 5 70415
402 333505 776 9 64230
403 296919 611 8 59190
404 278990 1526 25 69351
405 276898 592 7 64270
406 327007 1182 13 70694
407 317046 621 16 68005
408 304555 989 11 58930
409 298096 438 11 58320
410 231861 726 3 69980
411 309422 1303 61 69863
412 286963 7419 29 63255
413 269753 1164 17 57320
414 448243 3310 33 75230
415 165404 1920 15 79420
416 204325 965 3 73490
417 407159 3256 66 35250
418 290476 1135 17 62285
419 275311 1270 26 69206
420 246541 661 3 65920
421 253468 1013 2 69770
422 240897 2844 67 72683
423 -83265 11528 70 -14545
424 -42143 6526 26 55830
425 272713 2264 24 55174
426 215362 5109 97 67038
427 42754 3999 30 51252
428 306275 35624 223 157278
429 253537 9252 48 79510
430 372631 15236 90 77440
431 -7170 18073 180 27284
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) Costs Trades Dividends
159337.317 30.122 -404.577 2.000
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-2198710 -4394 34580 34580 3225801
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.593e+05 3.072e+04 5.186 3.32e-07 ***
Costs 3.012e+01 1.576e+00 19.114 < 2e-16 ***
Trades -4.046e+02 1.669e+02 -2.424 0.0158 *
Dividends 2.000e+00 4.103e-01 4.876 1.53e-06 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 254200 on 427 degrees of freedom
Multiple R-squared: 0.6537, Adjusted R-squared: 0.6512
F-statistic: 268.7 on 3 and 427 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 9.875609e-93 4.937805e-93
[2,] 1.0000000 6.214917e-110 3.107459e-110
[3,] 1.0000000 4.427687e-110 2.213844e-110
[4,] 1.0000000 2.773469e-123 1.386734e-123
[5,] 1.0000000 6.380867e-130 3.190434e-130
[6,] 1.0000000 5.015475e-134 2.507738e-134
[7,] 1.0000000 1.772949e-153 8.864743e-154
[8,] 1.0000000 1.107371e-160 5.536855e-161
[9,] 1.0000000 9.758343e-165 4.879171e-165
[10,] 1.0000000 1.633175e-172 8.165877e-173
[11,] 1.0000000 3.984119e-175 1.992060e-175
[12,] 1.0000000 6.034183e-178 3.017092e-178
[13,] 1.0000000 3.954370e-183 1.977185e-183
[14,] 1.0000000 2.586062e-183 1.293031e-183
[15,] 1.0000000 7.113834e-189 3.556917e-189
[16,] 1.0000000 1.139728e-190 5.698639e-191
[17,] 1.0000000 2.389425e-191 1.194713e-191
[18,] 1.0000000 2.283247e-191 1.141623e-191
[19,] 1.0000000 2.697203e-191 1.348602e-191
[20,] 1.0000000 5.183394e-192 2.591697e-192
[21,] 1.0000000 1.198780e-194 5.993899e-195
[22,] 1.0000000 9.930032e-194 4.965016e-194
[23,] 1.0000000 1.866795e-204 9.333976e-205
[24,] 1.0000000 1.217792e-204 6.088962e-205
[25,] 1.0000000 4.154427e-209 2.077214e-209
[26,] 1.0000000 3.191013e-210 1.595506e-210
[27,] 1.0000000 2.590078e-218 1.295039e-218
[28,] 1.0000000 2.928903e-220 1.464452e-220
[29,] 1.0000000 6.016902e-220 3.008451e-220
[30,] 1.0000000 7.060014e-220 3.530007e-220
[31,] 1.0000000 4.325746e-225 2.162873e-225
[32,] 1.0000000 5.370087e-231 2.685043e-231
[33,] 1.0000000 5.293847e-230 2.646923e-230
[34,] 1.0000000 7.734743e-230 3.867372e-230
[35,] 1.0000000 1.359907e-229 6.799535e-230
[36,] 1.0000000 4.296520e-229 2.148260e-229
[37,] 1.0000000 1.600355e-229 8.001777e-230
[38,] 1.0000000 1.782341e-228 8.911703e-229
[39,] 1.0000000 7.587601e-228 3.793800e-228
[40,] 1.0000000 9.591632e-227 4.795816e-227
[41,] 1.0000000 5.476649e-228 2.738324e-228
[42,] 1.0000000 6.657804e-227 3.328902e-227
[43,] 1.0000000 3.104703e-226 1.552351e-226
[44,] 1.0000000 2.190141e-226 1.095071e-226
[45,] 1.0000000 4.576708e-227 2.288354e-227
[46,] 1.0000000 9.139549e-227 4.569775e-227
[47,] 1.0000000 6.981042e-226 3.490521e-226
[48,] 1.0000000 2.494295e-225 1.247148e-225
[49,] 1.0000000 6.357183e-228 3.178591e-228
[50,] 1.0000000 1.013824e-227 5.069121e-228
[51,] 1.0000000 1.108036e-239 5.540181e-240
[52,] 1.0000000 6.622115e-240 3.311057e-240
[53,] 1.0000000 6.961421e-239 3.480710e-239
[54,] 1.0000000 1.452938e-238 7.264689e-239
[55,] 1.0000000 5.350273e-240 2.675137e-240
[56,] 1.0000000 2.628926e-241 1.314463e-241
[57,] 1.0000000 1.262599e-240 6.312994e-241
[58,] 1.0000000 5.672154e-240 2.836077e-240
[59,] 1.0000000 3.047876e-239 1.523938e-239
[60,] 1.0000000 3.783268e-238 1.891634e-238
[61,] 1.0000000 1.046125e-237 5.230625e-238
[62,] 1.0000000 1.109933e-236 5.549664e-237
[63,] 1.0000000 7.987114e-236 3.993557e-236
[64,] 1.0000000 2.963754e-235 1.481877e-235
[65,] 1.0000000 2.974734e-234 1.487367e-234
[66,] 1.0000000 3.048998e-235 1.524499e-235
[67,] 1.0000000 1.135725e-234 5.678627e-235
[68,] 1.0000000 1.387209e-233 6.936044e-234
[69,] 1.0000000 1.802432e-232 9.012162e-233
[70,] 1.0000000 4.573414e-232 2.286707e-232
[71,] 1.0000000 4.359277e-231 2.179639e-231
[72,] 1.0000000 5.873097e-231 2.936548e-231
[73,] 1.0000000 6.755034e-230 3.377517e-230
[74,] 1.0000000 1.169430e-229 5.847150e-230
[75,] 1.0000000 8.288844e-229 4.144422e-229
[76,] 1.0000000 1.217363e-228 6.086817e-229
[77,] 1.0000000 8.132851e-229 4.066425e-229
[78,] 1.0000000 6.296147e-228 3.148074e-228
[79,] 1.0000000 7.385872e-227 3.692936e-227
[80,] 1.0000000 1.583384e-226 7.916921e-227
[81,] 1.0000000 6.919061e-227 3.459531e-227
[82,] 1.0000000 5.805052e-226 2.902526e-226
[83,] 1.0000000 1.156444e-225 5.782221e-226
[84,] 1.0000000 1.378359e-224 6.891795e-225
[85,] 1.0000000 1.430831e-223 7.154154e-224
[86,] 1.0000000 1.003700e-222 5.018501e-223
[87,] 1.0000000 2.856121e-224 1.428060e-224
[88,] 1.0000000 2.159437e-223 1.079718e-223
[89,] 1.0000000 7.728519e-224 3.864259e-224
[90,] 1.0000000 8.297363e-224 4.148681e-224
[91,] 1.0000000 9.551591e-223 4.775795e-223
[92,] 1.0000000 1.062431e-221 5.312154e-222
[93,] 1.0000000 9.105336e-221 4.552668e-221
[94,] 1.0000000 1.071484e-219 5.357421e-220
[95,] 1.0000000 1.042746e-218 5.213732e-219
[96,] 1.0000000 7.487027e-220 3.743514e-220
[97,] 1.0000000 5.823251e-220 2.911625e-220
[98,] 1.0000000 9.018533e-221 4.509266e-221
[99,] 1.0000000 4.792609e-220 2.396305e-220
[100,] 1.0000000 1.773933e-219 8.869667e-220
[101,] 1.0000000 9.041595e-222 4.520797e-222
[102,] 1.0000000 5.532620e-223 2.766310e-223
[103,] 1.0000000 3.713151e-222 1.856576e-222
[104,] 1.0000000 7.065588e-222 3.532794e-222
[105,] 1.0000000 6.374601e-222 3.187300e-222
[106,] 1.0000000 3.840291e-222 1.920146e-222
[107,] 1.0000000 2.653352e-221 1.326676e-221
[108,] 1.0000000 2.305012e-220 1.152506e-220
[109,] 1.0000000 1.950256e-220 9.751279e-221
[110,] 1.0000000 8.449564e-220 4.224782e-220
[111,] 1.0000000 6.552957e-219 3.276479e-219
[112,] 1.0000000 6.635318e-219 3.317659e-219
[113,] 1.0000000 6.463719e-220 3.231859e-220
[114,] 1.0000000 7.504225e-219 3.752112e-219
[115,] 1.0000000 3.615875e-218 1.807938e-218
[116,] 1.0000000 3.915776e-217 1.957888e-217
[117,] 1.0000000 3.478489e-216 1.739245e-216
[118,] 1.0000000 6.400805e-216 3.200402e-216
[119,] 1.0000000 2.201908e-215 1.100954e-215
[120,] 1.0000000 6.604227e-216 3.302114e-216
[121,] 1.0000000 7.549642e-215 3.774821e-215
[122,] 1.0000000 8.636219e-214 4.318109e-214
[123,] 1.0000000 9.041269e-213 4.520635e-213
[124,] 1.0000000 1.022146e-211 5.110732e-212
[125,] 1.0000000 1.154369e-210 5.771843e-211
[126,] 1.0000000 1.257251e-210 6.286256e-211
[127,] 1.0000000 1.416033e-209 7.080166e-210
[128,] 1.0000000 1.604869e-208 8.024343e-209
[129,] 1.0000000 1.808183e-207 9.040913e-208
[130,] 1.0000000 2.003049e-206 1.001525e-206
[131,] 1.0000000 2.180228e-205 1.090114e-205
[132,] 1.0000000 2.409455e-204 1.204728e-204
[133,] 1.0000000 2.575562e-203 1.287781e-203
[134,] 1.0000000 2.815797e-202 1.407898e-202
[135,] 1.0000000 2.599239e-201 1.299620e-201
[136,] 1.0000000 2.818858e-200 1.409429e-200
[137,] 1.0000000 3.105819e-199 1.552909e-199
[138,] 1.0000000 3.336564e-198 1.668282e-198
[139,] 1.0000000 3.284344e-197 1.642172e-197
[140,] 1.0000000 3.274545e-196 1.637273e-196
[141,] 1.0000000 3.474687e-195 1.737344e-195
[142,] 1.0000000 3.680834e-194 1.840417e-194
[143,] 1.0000000 3.505095e-193 1.752548e-193
[144,] 1.0000000 3.674480e-192 1.837240e-192
[145,] 1.0000000 3.829009e-191 1.914504e-191
[146,] 1.0000000 4.041399e-190 2.020699e-190
[147,] 1.0000000 4.175814e-189 2.087907e-189
[148,] 1.0000000 4.380274e-188 2.190137e-188
[149,] 1.0000000 4.497351e-187 2.248676e-187
[150,] 1.0000000 4.585045e-186 2.292523e-186
[151,] 1.0000000 4.653728e-185 2.326864e-185
[152,] 1.0000000 4.735252e-184 2.367626e-184
[153,] 1.0000000 3.569242e-183 1.784621e-183
[154,] 1.0000000 3.577261e-182 1.788630e-182
[155,] 1.0000000 3.569299e-181 1.784650e-181
[156,] 1.0000000 3.545443e-180 1.772722e-180
[157,] 1.0000000 3.505994e-179 1.752997e-179
[158,] 1.0000000 3.513381e-178 1.756691e-178
[159,] 1.0000000 3.277433e-177 1.638717e-177
[160,] 1.0000000 3.203972e-176 1.601986e-176
[161,] 1.0000000 3.112640e-175 1.556320e-175
[162,] 1.0000000 3.607545e-176 1.803772e-176
[163,] 1.0000000 1.341746e-175 6.708729e-176
[164,] 1.0000000 1.316268e-174 6.581342e-175
[165,] 1.0000000 1.285459e-173 6.427295e-174
[166,] 1.0000000 1.264426e-172 6.322130e-173
[167,] 1.0000000 1.223736e-171 6.118678e-172
[168,] 1.0000000 1.196037e-170 5.980186e-171
[169,] 1.0000000 1.147125e-169 5.735624e-170
[170,] 1.0000000 1.092220e-168 5.461100e-169
[171,] 1.0000000 9.880837e-168 4.940418e-168
[172,] 1.0000000 9.309706e-167 4.654853e-167
[173,] 1.0000000 8.774024e-166 4.387012e-166
[174,] 1.0000000 8.231790e-165 4.115895e-165
[175,] 1.0000000 5.620470e-164 2.810235e-164
[176,] 1.0000000 5.322883e-163 2.661441e-163
[177,] 1.0000000 3.053001e-162 1.526500e-162
[178,] 1.0000000 1.923482e-161 9.617409e-162
[179,] 1.0000000 1.791673e-160 8.958367e-161
[180,] 1.0000000 1.240046e-159 6.200232e-160
[181,] 1.0000000 4.710561e-159 2.355280e-159
[182,] 1.0000000 4.353831e-158 2.176916e-158
[183,] 1.0000000 3.141849e-157 1.570925e-157
[184,] 1.0000000 2.844780e-156 1.422390e-156
[185,] 1.0000000 2.422479e-155 1.211240e-155
[186,] 1.0000000 2.173839e-154 1.086920e-154
[187,] 1.0000000 1.951216e-153 9.756079e-154
[188,] 1.0000000 1.735142e-152 8.675712e-153
[189,] 1.0000000 1.535140e-151 7.675702e-152
[190,] 1.0000000 1.352829e-150 6.764143e-151
[191,] 1.0000000 1.085638e-149 5.428188e-150
[192,] 1.0000000 9.484256e-149 4.742128e-149
[193,] 1.0000000 8.248111e-148 4.124056e-148
[194,] 1.0000000 8.572229e-148 4.286114e-148
[195,] 1.0000000 7.463794e-147 3.731897e-147
[196,] 1.0000000 6.469330e-146 3.234665e-146
[197,] 1.0000000 5.582023e-145 2.791012e-145
[198,] 1.0000000 4.336191e-144 2.168095e-144
[199,] 1.0000000 1.913698e-143 9.568490e-144
[200,] 1.0000000 1.643158e-142 8.215791e-143
[201,] 1.0000000 1.410649e-141 7.053244e-142
[202,] 1.0000000 1.194740e-140 5.973702e-141
[203,] 1.0000000 1.708399e-140 8.541995e-141
[204,] 1.0000000 1.461851e-139 7.309254e-140
[205,] 1.0000000 1.230659e-138 6.153293e-139
[206,] 1.0000000 1.031322e-137 5.156612e-138
[207,] 1.0000000 6.314833e-137 3.157417e-137
[208,] 1.0000000 4.820284e-136 2.410142e-136
[209,] 1.0000000 4.036397e-135 2.018198e-135
[210,] 1.0000000 3.328507e-134 1.664253e-134
[211,] 1.0000000 2.725239e-133 1.362620e-133
[212,] 1.0000000 2.139781e-132 1.069890e-132
[213,] 1.0000000 1.740816e-131 8.704079e-132
[214,] 1.0000000 8.436053e-131 4.218026e-131
[215,] 1.0000000 6.768944e-130 3.384472e-130
[216,] 1.0000000 5.444065e-129 2.722032e-129
[217,] 1.0000000 4.355243e-128 2.177622e-128
[218,] 1.0000000 3.455383e-127 1.727691e-127
[219,] 1.0000000 2.738613e-126 1.369306e-126
[220,] 1.0000000 2.152252e-125 1.076126e-125
[221,] 1.0000000 1.688557e-124 8.442784e-125
[222,] 1.0000000 1.321247e-123 6.606237e-124
[223,] 1.0000000 1.032621e-122 5.163105e-123
[224,] 1.0000000 8.007525e-122 4.003763e-122
[225,] 1.0000000 5.887986e-121 2.943993e-121
[226,] 1.0000000 4.140144e-120 2.070072e-120
[227,] 1.0000000 3.168814e-119 1.584407e-119
[228,] 1.0000000 2.440999e-118 1.220499e-118
[229,] 1.0000000 1.851262e-117 9.256310e-118
[230,] 1.0000000 1.396814e-116 6.984068e-117
[231,] 1.0000000 1.049519e-115 5.247597e-116
[232,] 1.0000000 7.777962e-115 3.888981e-115
[233,] 1.0000000 5.789815e-114 2.894907e-114
[234,] 1.0000000 4.289704e-113 2.144852e-113
[235,] 1.0000000 3.163373e-112 1.581687e-112
[236,] 1.0000000 2.345837e-111 1.172918e-111
[237,] 1.0000000 1.713720e-110 8.568598e-111
[238,] 1.0000000 1.227663e-109 6.138316e-110
[239,] 1.0000000 8.887127e-109 4.443563e-109
[240,] 1.0000000 6.400576e-108 3.200288e-108
[241,] 1.0000000 4.587926e-107 2.293963e-107
[242,] 1.0000000 3.321162e-106 1.660581e-106
[243,] 1.0000000 2.358307e-105 1.179153e-105
[244,] 1.0000000 1.666627e-104 8.333135e-105
[245,] 1.0000000 1.155720e-103 5.778598e-104
[246,] 1.0000000 8.089828e-103 4.044914e-103
[247,] 1.0000000 5.631279e-102 2.815640e-102
[248,] 1.0000000 3.210088e-101 1.605044e-101
[249,] 1.0000000 2.432782e-101 1.216391e-101
[250,] 1.0000000 1.692684e-100 8.463420e-101
[251,] 1.0000000 1.066985e-99 5.334926e-100
[252,] 1.0000000 6.157806e-99 3.078903e-99
[253,] 1.0000000 4.227291e-98 2.113646e-98
[254,] 1.0000000 2.887882e-97 1.443941e-97
[255,] 1.0000000 1.963235e-96 9.816177e-97
[256,] 1.0000000 1.328114e-95 6.640571e-96
[257,] 1.0000000 9.029732e-95 4.514866e-95
[258,] 1.0000000 4.800266e-94 2.400133e-94
[259,] 1.0000000 3.202844e-93 1.601422e-93
[260,] 1.0000000 2.126437e-92 1.063219e-92
[261,] 1.0000000 1.404784e-91 7.023921e-92
[262,] 1.0000000 7.637153e-91 3.818576e-91
[263,] 1.0000000 4.998906e-90 2.499453e-90
[264,] 1.0000000 3.317701e-89 1.658851e-89
[265,] 1.0000000 2.150197e-88 1.075099e-88
[266,] 1.0000000 1.386521e-87 6.932604e-88
[267,] 1.0000000 8.895584e-87 4.447792e-87
[268,] 1.0000000 5.769418e-86 2.884709e-86
[269,] 1.0000000 3.664388e-85 1.832194e-85
[270,] 1.0000000 2.309817e-84 1.154908e-84
[271,] 1.0000000 1.469383e-83 7.346913e-84
[272,] 1.0000000 9.193088e-83 4.596544e-83
[273,] 1.0000000 4.780492e-82 2.390246e-82
[274,] 1.0000000 2.846048e-81 1.423024e-81
[275,] 1.0000000 1.755835e-80 8.779173e-81
[276,] 1.0000000 1.077606e-79 5.388029e-80
[277,] 1.0000000 6.430182e-79 3.215091e-79
[278,] 1.0000000 3.907637e-78 1.953819e-78
[279,] 1.0000000 2.413419e-77 1.206709e-77
[280,] 1.0000000 1.451536e-76 7.257680e-77
[281,] 1.0000000 8.238338e-76 4.119169e-76
[282,] 1.0000000 4.922932e-75 2.461466e-75
[283,] 1.0000000 2.914186e-74 1.457093e-74
[284,] 1.0000000 1.749147e-73 8.745736e-74
[285,] 1.0000000 1.024604e-72 5.123021e-73
[286,] 1.0000000 5.881046e-72 2.940523e-72
[287,] 1.0000000 3.444738e-71 1.722369e-71
[288,] 1.0000000 1.985892e-70 9.929459e-71
[289,] 1.0000000 1.164554e-69 5.822770e-70
[290,] 1.0000000 4.858391e-69 2.429195e-69
[291,] 1.0000000 2.760285e-68 1.380143e-68
[292,] 1.0000000 1.559553e-67 7.797767e-68
[293,] 1.0000000 8.762331e-67 4.381165e-67
[294,] 1.0000000 4.895535e-66 2.447767e-66
[295,] 1.0000000 2.747026e-65 1.373513e-65
[296,] 1.0000000 1.549794e-64 7.748968e-65
[297,] 1.0000000 8.513971e-64 4.256986e-64
[298,] 1.0000000 4.650528e-63 2.325264e-63
[299,] 1.0000000 2.050656e-62 1.025328e-62
[300,] 1.0000000 1.113027e-61 5.565137e-62
[301,] 1.0000000 5.987148e-61 2.993574e-61
[302,] 1.0000000 3.272556e-60 1.636278e-60
[303,] 1.0000000 1.605672e-59 8.028362e-60
[304,] 1.0000000 8.491598e-59 4.245799e-59
[305,] 1.0000000 4.464165e-58 2.232082e-58
[306,] 1.0000000 1.905455e-57 9.527277e-58
[307,] 1.0000000 9.915633e-57 4.957817e-57
[308,] 1.0000000 5.195221e-56 2.597611e-56
[309,] 1.0000000 2.675413e-55 1.337706e-55
[310,] 1.0000000 1.298822e-54 6.494109e-55
[311,] 1.0000000 6.543833e-54 3.271916e-54
[312,] 1.0000000 3.391570e-53 1.695785e-53
[313,] 1.0000000 1.708983e-52 8.544914e-53
[314,] 1.0000000 8.527669e-52 4.263835e-52
[315,] 1.0000000 4.228446e-51 2.114223e-51
[316,] 1.0000000 2.083387e-50 1.041693e-50
[317,] 1.0000000 1.019948e-49 5.099738e-50
[318,] 1.0000000 4.961176e-49 2.480588e-49
[319,] 1.0000000 2.397560e-48 1.198780e-48
[320,] 1.0000000 1.165094e-47 5.825471e-48
[321,] 1.0000000 5.557371e-47 2.778685e-47
[322,] 1.0000000 2.633245e-46 1.316622e-46
[323,] 1.0000000 1.239377e-45 6.196887e-46
[324,] 1.0000000 5.794059e-45 2.897030e-45
[325,] 1.0000000 2.690326e-44 1.345163e-44
[326,] 1.0000000 1.240637e-43 6.203183e-44
[327,] 1.0000000 5.681673e-43 2.840836e-43
[328,] 1.0000000 2.592458e-42 1.296229e-42
[329,] 1.0000000 1.172813e-41 5.864067e-42
[330,] 1.0000000 5.207948e-41 2.603974e-41
[331,] 1.0000000 2.372217e-40 1.186109e-40
[332,] 1.0000000 1.049442e-39 5.247209e-40
[333,] 1.0000000 4.648204e-39 2.324102e-39
[334,] 1.0000000 1.991062e-38 9.955308e-39
[335,] 1.0000000 8.885298e-38 4.442649e-38
[336,] 1.0000000 3.822618e-37 1.911309e-37
[337,] 1.0000000 1.632186e-36 8.160930e-37
[338,] 1.0000000 6.916143e-36 3.458071e-36
[339,] 1.0000000 3.017103e-35 1.508552e-35
[340,] 1.0000000 1.259771e-34 6.298855e-35
[341,] 1.0000000 5.363560e-34 2.681780e-34
[342,] 1.0000000 2.205094e-33 1.102547e-33
[343,] 1.0000000 8.993246e-33 4.496623e-33
[344,] 1.0000000 3.282821e-32 1.641410e-32
[345,] 1.0000000 1.252869e-31 6.264343e-32
[346,] 1.0000000 5.081274e-31 2.540637e-31
[347,] 1.0000000 1.803898e-30 9.019490e-31
[348,] 1.0000000 7.196742e-30 3.598371e-30
[349,] 1.0000000 2.796494e-29 1.398247e-29
[350,] 1.0000000 1.077198e-28 5.385992e-29
[351,] 1.0000000 3.915609e-28 1.957804e-28
[352,] 1.0000000 1.486898e-27 7.434490e-28
[353,] 1.0000000 5.595423e-27 2.797712e-27
[354,] 1.0000000 2.086393e-26 1.043196e-26
[355,] 1.0000000 7.752073e-26 3.876036e-26
[356,] 1.0000000 2.837219e-25 1.418610e-25
[357,] 1.0000000 1.028468e-24 5.142340e-25
[358,] 1.0000000 3.672638e-24 1.836319e-24
[359,] 1.0000000 1.291124e-23 6.455618e-24
[360,] 1.0000000 4.548003e-23 2.274001e-23
[361,] 1.0000000 1.585700e-22 7.928500e-23
[362,] 1.0000000 5.471218e-22 2.735609e-22
[363,] 1.0000000 1.867755e-21 9.338776e-22
[364,] 1.0000000 6.464886e-21 3.232443e-21
[365,] 1.0000000 2.161510e-20 1.080755e-20
[366,] 1.0000000 7.127188e-20 3.563594e-20
[367,] 1.0000000 2.517037e-19 1.258518e-19
[368,] 1.0000000 8.712157e-19 4.356079e-19
[369,] 1.0000000 2.797258e-18 1.398629e-18
[370,] 1.0000000 8.873494e-18 4.436747e-18
[371,] 1.0000000 2.951247e-17 1.475623e-17
[372,] 1.0000000 8.645670e-17 4.322835e-17
[373,] 1.0000000 2.275641e-16 1.137821e-16
[374,] 1.0000000 6.877425e-16 3.438713e-16
[375,] 1.0000000 2.050173e-15 1.025086e-15
[376,] 1.0000000 6.521844e-15 3.260922e-15
[377,] 1.0000000 1.983022e-14 9.915108e-15
[378,] 1.0000000 5.728861e-14 2.864430e-14
[379,] 1.0000000 1.707799e-13 8.538995e-14
[380,] 1.0000000 4.468576e-13 2.234288e-13
[381,] 1.0000000 1.156656e-12 5.783280e-13
[382,] 1.0000000 3.095127e-12 1.547563e-12
[383,] 1.0000000 8.453329e-12 4.226665e-12
[384,] 1.0000000 2.298339e-11 1.149169e-11
[385,] 1.0000000 6.040574e-11 3.020287e-11
[386,] 1.0000000 1.696765e-10 8.483827e-11
[387,] 1.0000000 4.916376e-10 2.458188e-10
[388,] 1.0000000 9.655498e-10 4.827749e-10
[389,] 1.0000000 2.666488e-09 1.333244e-09
[390,] 1.0000000 7.261224e-09 3.630612e-09
[391,] 1.0000000 1.882463e-08 9.412317e-09
[392,] 1.0000000 4.058986e-08 2.029493e-08
[393,] 0.9999999 1.067754e-07 5.338771e-08
[394,] 0.9999999 2.798909e-07 1.399455e-07
[395,] 0.9999996 7.160232e-07 3.580116e-07
[396,] 0.9999992 1.529035e-06 7.645177e-07
[397,] 0.9999982 3.539000e-06 1.769500e-06
[398,] 0.9999956 8.791658e-06 4.395829e-06
[399,] 0.9999894 2.129955e-05 1.064977e-05
[400,] 0.9999764 4.714256e-05 2.357128e-05
[401,] 0.9999475 1.050786e-04 5.253932e-05
[402,] 0.9998952 2.095866e-04 1.047933e-04
[403,] 0.9997876 4.248413e-04 2.124206e-04
[404,] 0.9995456 9.087123e-04 4.543562e-04
[405,] 0.9990200 1.960069e-03 9.800345e-04
[406,] 0.9984653 3.069301e-03 1.534651e-03
[407,] 0.9969555 6.088941e-03 3.044471e-03
[408,] 0.9973612 5.277579e-03 2.638790e-03
[409,] 0.9958550 8.289994e-03 4.144997e-03
[410,] 0.9920121 1.597584e-02 7.987918e-03
[411,] 0.9976614 4.677224e-03 2.338612e-03
[412,] 0.9951718 9.656323e-03 4.828162e-03
[413,] 0.9889162 2.216767e-02 1.108384e-02
[414,] 0.9750855 4.982890e-02 2.491445e-02
[415,] 0.9470880 1.058241e-01 5.291205e-02
[416,] 0.8920542 2.158916e-01 1.079458e-01
[417,] 0.8048247 3.903505e-01 1.951753e-01
[418,] 0.8344600 3.310800e-01 1.655400e-01
> postscript(file="/var/www/html/freestat/rcomp/tmp/1ec8y1291211846.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/263811291211846.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/363811291211846.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/463811291211846.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/5hvpl1291211846.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
1237333.7716 3225800.6189 1194235.0281 -2198709.5189 -446207.8858
6 7 8 9 10
-65315.8652 -58061.4111 79150.5670 -249447.6804 742457.9366
11 12 13 14 15
-308956.7059 392076.2392 803384.0750 -1209577.7725 35999.2506
16 17 18 19 20
475470.1181 -267414.2259 17699.4160 -274855.9885 -346491.1382
21 22 23 24 25
190089.0283 -245640.1815 -148253.8127 -401854.6120 -302236.3074
26 27 28 29 30
-152079.9239 -370608.0313 -297472.3304 319888.6245 -107458.5634
31 32 33 34 35
264934.4893 -118886.7400 228735.6187 117862.7404 -176569.1359
36 37 38 39 40
-482482.7536 29495.2712 185164.6824 -31633.2784 -175484.6937
41 42 43 44 45
-28586.2323 -395020.7371 109887.2197 -178865.3985 -24526.8599
46 47 48 49 50
-86659.7032 -269247.8351 -97824.0452 -47758.9998 -173126.2485
51 52 53 54 55
103613.4092 25139.9739 -218762.7914 106580.6851 -315223.4826
56 57 58 59 60
91349.0492 283612.8987 -23682.0837 -12429.1879 -186594.5055
61 62 63 64 65
-3288.1749 -164170.1642 -14535.6869 -152468.2407 29170.0161
66 67 68 69 70
658.3290 -56301.8915 -3474.4077 -51538.3527 -31098.5552
71 72 73 74 75
-80208.5263 6724.8692 -91727.8629 -4435.9069 -47964.1058
76 77 78 79 80
-112625.8705 -32306.1501 -77635.7304 22864.3015 -6161.5256
81 82 83 84 85
35011.7284 -23890.0975 -86997.5856 46042.0802 3276.6093
86 87 88 89 90
-72852.1034 43537.4967 -23669.5973 -101295.9632 7517.6045
91 92 93 94 95
3632.1632 -18907.1836 138110.3578 -214809.1834 85549.2598
96 97 98 99 100
22161.4187 38786.3001 68328.4234 -84437.6455 -25122.8216
101 102 103 104 105
-28698.3104 -125369.5171 -146117.0270 -82896.3085 36828.6130
106 107 108 109 110
-45990.5534 -225506.3353 102172.4638 -187517.4506 -99005.2886
111 112 113 114 115
97583.6413 -181421.3633 -165268.3173 -84430.1190 65432.0980
116 117 118 119 120
41884.9294 1804.1853 47384.8091 -160370.2516 36507.6963
121 122 123 124 125
-286240.9429 20225.1998 5840.5691 67611.1743 3093.7535
126 127 128 129 130
-183753.0522 34580.2165 33587.7584 41443.9094 34580.2165
131 132 133 134 135
35945.7172 77374.3357 34580.2165 31630.6387 13869.3958
136 137 138 139 140
28761.9295 39448.6018 34745.7604 9116.6456 34580.2165
141 142 143 144 145
42113.4674 34580.2165 1596.0722 34815.4903 35943.7690
146 147 148 149 150
35592.4067 34580.2165 36113.5379 37095.3019 34467.9757
151 152 153 154 155
3373.2999 2061.9484 34580.2165 11810.4517 34620.4264
156 157 158 159 160
34580.2165 34580.2165 30753.1273 46007.7895 34580.2165
161 162 163 164 165
34580.2165 34580.2165 34580.2165 21622.8177 23262.3055
166 167 168 169 170
33880.2124 34580.2165 116864.3690 79550.9637 34580.2165
171 172 173 174 175
34580.2165 30457.8195 34580.2165 23886.1881 34580.2165
176 177 178 179 180
34535.5262 81914.1147 33592.1205 34580.2165 34580.2165
181 182 183 184 185
41164.1864 15473.1509 -7722.3720 3954.9700 31665.6161
186 187 188 189 190
39556.5727 5452.7572 27424.4312 -3849.3352 34580.2165
191 192 193 194 195
33720.4193 34580.2165 30926.1365 34580.2165 34513.8910
196 197 198 199 200
34580.2165 -89970.4929 34580.2165 34580.2165 -55995.9889
201 202 203 204 205
34580.2165 34580.2165 34580.2165 44806.6942 -16733.3662
206 207 208 209 210
33990.8040 27556.9989 34580.2165 65941.5748 30707.7678
211 212 213 214 215
34580.2165 34580.2165 60324.2404 -515.5477 6631.3916
216 217 218 219 220
34580.2165 37450.9613 36453.9569 34580.2165 6047.9156
221 222 223 224 225
9219.4592 34580.2165 35317.0126 25760.2266 34812.8986
226 227 228 229 230
34899.0433 37577.8514 34580.2165 18260.4102 34580.2165
231 232 233 234 235
5601.5458 22391.9327 34580.2165 16290.7635 34580.2165
236 237 238 239 240
34985.5410 34580.2165 36566.3486 34580.2165 34580.2165
241 242 243 244 245
34580.2165 29189.2405 34580.2165 39155.2379 -48965.4619
246 247 248 249 250
34580.2165 34580.2165 26517.1760 34580.2165 34580.2165
251 252 253 254 255
28561.0884 34580.2165 15249.6882 9795.8940 -103337.7607
256 257 258 259 260
34580.2165 32804.4768 46896.8111 34580.2165 34580.2165
261 262 263 264 265
34580.2165 34580.2165 36942.0807 51705.1001 34580.2165
266 267 268 269 270
34580.2165 34580.2165 65291.3572 34580.2165 30693.4734
271 272 273 274 275
34580.2165 34580.2165 34580.2165 30334.8589 34580.2165
276 277 278 279 280
39394.4755 35947.1018 34580.2165 5052.6817 21062.1530
281 282 283 284 285
34580.2165 34580.2165 62309.8561 34580.2165 25566.3500
286 287 288 289 290
34580.2165 22734.3419 35362.3736 34580.2165 42333.9750
291 292 293 294 295
34580.2165 44873.1463 39192.5070 34580.2165 28105.2768
296 297 298 299 300
59868.2221 34580.2165 34580.2165 34580.2165 34580.2165
301 302 303 304 305
32922.8350 16298.6976 34580.2165 34580.2165 71866.7354
306 307 308 309 310
38145.1602 34580.2165 28217.9952 38446.2874 34580.2165
311 312 313 314 315
34580.2165 35442.7760 34580.2165 36216.6003 36302.1434
316 317 318 319 320
49761.1157 35069.1596 11239.9555 36971.5570 34580.2165
321 322 323 324 325
34580.2165 34580.2165 34580.2165 34580.2165 34580.2165
326 327 328 329 330
31844.6853 34580.2165 34580.2165 34580.2165 34580.2165
331 332 333 334 335
34580.2165 34580.2165 34580.2165 32742.5135 35353.2460
336 337 338 339 340
27135.7859 21063.1600 34580.2165 32349.5558 40364.8296
341 342 343 344 345
7192.9443 34580.2165 34580.2165 34580.2165 28170.1365
346 347 348 349 350
34580.2165 24924.7505 34580.2165 34580.2165 -2597.6032
351 352 353 354 355
38401.0084 34679.5158 34322.9287 27462.3996 34580.2165
356 357 358 359 360
34580.2165 11006.8190 34580.2165 34580.2165 34580.2165
361 362 363 364 365
36182.4019 34580.2165 34580.2165 -17619.3897 36785.0008
366 367 368 369 370
34580.2165 34580.2165 34580.2165 34580.2165 36682.7327
371 372 373 374 375
34580.2165 47127.9780 20788.8182 13586.0443 34580.2165
376 377 378 379 380
34580.2165 29720.9807 28274.6000 -21901.6234 34580.2165
381 382 383 384 385
34580.2165 19370.4612 53257.8226 34580.2165 36693.1984
386 387 388 389 390
-115636.0426 46871.0853 7872.9811 33099.2810 20654.4427
391 392 393 394 395
29633.8067 3808.4378 -33887.7287 54504.4319 -141011.2663
396 397 398 399 400
-4353.0222 -152446.6167 41068.9784 -3133.3067 21066.1182
401 402 403 404 405
-25877.9536 25950.4200 4011.8418 -54926.7608 -26003.2963
406 407 408 409 410
-4089.2056 9440.9762 1995.5438 13354.0051 -88117.0685
411 412 413 414 415
-4236.9658 -210650.2878 -32429.7526 52064.9884 -204568.3343
416 417 418 419 420
-129873.5310 105933.4123 -20765.0510 -50199.8918 -63357.6347
421 422 423 424 425
-75139.5888 -122393.5869 -532433.2423 -499218.3037 -55479.1464
426 427 428 429 430
-192725.5730 -327426.9457 -1150522.8360 -324119.0097 -364142.2036
431
-692656.8969
> postscript(file="/var/www/html/freestat/rcomp/tmp/6hvpl1291211846.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 1237333.7716 NA
1 3225800.6189 1237333.7716
2 1194235.0281 3225800.6189
3 -2198709.5189 1194235.0281
4 -446207.8858 -2198709.5189
5 -65315.8652 -446207.8858
6 -58061.4111 -65315.8652
7 79150.5670 -58061.4111
8 -249447.6804 79150.5670
9 742457.9366 -249447.6804
10 -308956.7059 742457.9366
11 392076.2392 -308956.7059
12 803384.0750 392076.2392
13 -1209577.7725 803384.0750
14 35999.2506 -1209577.7725
15 475470.1181 35999.2506
16 -267414.2259 475470.1181
17 17699.4160 -267414.2259
18 -274855.9885 17699.4160
19 -346491.1382 -274855.9885
20 190089.0283 -346491.1382
21 -245640.1815 190089.0283
22 -148253.8127 -245640.1815
23 -401854.6120 -148253.8127
24 -302236.3074 -401854.6120
25 -152079.9239 -302236.3074
26 -370608.0313 -152079.9239
27 -297472.3304 -370608.0313
28 319888.6245 -297472.3304
29 -107458.5634 319888.6245
30 264934.4893 -107458.5634
31 -118886.7400 264934.4893
32 228735.6187 -118886.7400
33 117862.7404 228735.6187
34 -176569.1359 117862.7404
35 -482482.7536 -176569.1359
36 29495.2712 -482482.7536
37 185164.6824 29495.2712
38 -31633.2784 185164.6824
39 -175484.6937 -31633.2784
40 -28586.2323 -175484.6937
41 -395020.7371 -28586.2323
42 109887.2197 -395020.7371
43 -178865.3985 109887.2197
44 -24526.8599 -178865.3985
45 -86659.7032 -24526.8599
46 -269247.8351 -86659.7032
47 -97824.0452 -269247.8351
48 -47758.9998 -97824.0452
49 -173126.2485 -47758.9998
50 103613.4092 -173126.2485
51 25139.9739 103613.4092
52 -218762.7914 25139.9739
53 106580.6851 -218762.7914
54 -315223.4826 106580.6851
55 91349.0492 -315223.4826
56 283612.8987 91349.0492
57 -23682.0837 283612.8987
58 -12429.1879 -23682.0837
59 -186594.5055 -12429.1879
60 -3288.1749 -186594.5055
61 -164170.1642 -3288.1749
62 -14535.6869 -164170.1642
63 -152468.2407 -14535.6869
64 29170.0161 -152468.2407
65 658.3290 29170.0161
66 -56301.8915 658.3290
67 -3474.4077 -56301.8915
68 -51538.3527 -3474.4077
69 -31098.5552 -51538.3527
70 -80208.5263 -31098.5552
71 6724.8692 -80208.5263
72 -91727.8629 6724.8692
73 -4435.9069 -91727.8629
74 -47964.1058 -4435.9069
75 -112625.8705 -47964.1058
76 -32306.1501 -112625.8705
77 -77635.7304 -32306.1501
78 22864.3015 -77635.7304
79 -6161.5256 22864.3015
80 35011.7284 -6161.5256
81 -23890.0975 35011.7284
82 -86997.5856 -23890.0975
83 46042.0802 -86997.5856
84 3276.6093 46042.0802
85 -72852.1034 3276.6093
86 43537.4967 -72852.1034
87 -23669.5973 43537.4967
88 -101295.9632 -23669.5973
89 7517.6045 -101295.9632
90 3632.1632 7517.6045
91 -18907.1836 3632.1632
92 138110.3578 -18907.1836
93 -214809.1834 138110.3578
94 85549.2598 -214809.1834
95 22161.4187 85549.2598
96 38786.3001 22161.4187
97 68328.4234 38786.3001
98 -84437.6455 68328.4234
99 -25122.8216 -84437.6455
100 -28698.3104 -25122.8216
101 -125369.5171 -28698.3104
102 -146117.0270 -125369.5171
103 -82896.3085 -146117.0270
104 36828.6130 -82896.3085
105 -45990.5534 36828.6130
106 -225506.3353 -45990.5534
107 102172.4638 -225506.3353
108 -187517.4506 102172.4638
109 -99005.2886 -187517.4506
110 97583.6413 -99005.2886
111 -181421.3633 97583.6413
112 -165268.3173 -181421.3633
113 -84430.1190 -165268.3173
114 65432.0980 -84430.1190
115 41884.9294 65432.0980
116 1804.1853 41884.9294
117 47384.8091 1804.1853
118 -160370.2516 47384.8091
119 36507.6963 -160370.2516
120 -286240.9429 36507.6963
121 20225.1998 -286240.9429
122 5840.5691 20225.1998
123 67611.1743 5840.5691
124 3093.7535 67611.1743
125 -183753.0522 3093.7535
126 34580.2165 -183753.0522
127 33587.7584 34580.2165
128 41443.9094 33587.7584
129 34580.2165 41443.9094
130 35945.7172 34580.2165
131 77374.3357 35945.7172
132 34580.2165 77374.3357
133 31630.6387 34580.2165
134 13869.3958 31630.6387
135 28761.9295 13869.3958
136 39448.6018 28761.9295
137 34745.7604 39448.6018
138 9116.6456 34745.7604
139 34580.2165 9116.6456
140 42113.4674 34580.2165
141 34580.2165 42113.4674
142 1596.0722 34580.2165
143 34815.4903 1596.0722
144 35943.7690 34815.4903
145 35592.4067 35943.7690
146 34580.2165 35592.4067
147 36113.5379 34580.2165
148 37095.3019 36113.5379
149 34467.9757 37095.3019
150 3373.2999 34467.9757
151 2061.9484 3373.2999
152 34580.2165 2061.9484
153 11810.4517 34580.2165
154 34620.4264 11810.4517
155 34580.2165 34620.4264
156 34580.2165 34580.2165
157 30753.1273 34580.2165
158 46007.7895 30753.1273
159 34580.2165 46007.7895
160 34580.2165 34580.2165
161 34580.2165 34580.2165
162 34580.2165 34580.2165
163 21622.8177 34580.2165
164 23262.3055 21622.8177
165 33880.2124 23262.3055
166 34580.2165 33880.2124
167 116864.3690 34580.2165
168 79550.9637 116864.3690
169 34580.2165 79550.9637
170 34580.2165 34580.2165
171 30457.8195 34580.2165
172 34580.2165 30457.8195
173 23886.1881 34580.2165
174 34580.2165 23886.1881
175 34535.5262 34580.2165
176 81914.1147 34535.5262
177 33592.1205 81914.1147
178 34580.2165 33592.1205
179 34580.2165 34580.2165
180 41164.1864 34580.2165
181 15473.1509 41164.1864
182 -7722.3720 15473.1509
183 3954.9700 -7722.3720
184 31665.6161 3954.9700
185 39556.5727 31665.6161
186 5452.7572 39556.5727
187 27424.4312 5452.7572
188 -3849.3352 27424.4312
189 34580.2165 -3849.3352
190 33720.4193 34580.2165
191 34580.2165 33720.4193
192 30926.1365 34580.2165
193 34580.2165 30926.1365
194 34513.8910 34580.2165
195 34580.2165 34513.8910
196 -89970.4929 34580.2165
197 34580.2165 -89970.4929
198 34580.2165 34580.2165
199 -55995.9889 34580.2165
200 34580.2165 -55995.9889
201 34580.2165 34580.2165
202 34580.2165 34580.2165
203 44806.6942 34580.2165
204 -16733.3662 44806.6942
205 33990.8040 -16733.3662
206 27556.9989 33990.8040
207 34580.2165 27556.9989
208 65941.5748 34580.2165
209 30707.7678 65941.5748
210 34580.2165 30707.7678
211 34580.2165 34580.2165
212 60324.2404 34580.2165
213 -515.5477 60324.2404
214 6631.3916 -515.5477
215 34580.2165 6631.3916
216 37450.9613 34580.2165
217 36453.9569 37450.9613
218 34580.2165 36453.9569
219 6047.9156 34580.2165
220 9219.4592 6047.9156
221 34580.2165 9219.4592
222 35317.0126 34580.2165
223 25760.2266 35317.0126
224 34812.8986 25760.2266
225 34899.0433 34812.8986
226 37577.8514 34899.0433
227 34580.2165 37577.8514
228 18260.4102 34580.2165
229 34580.2165 18260.4102
230 5601.5458 34580.2165
231 22391.9327 5601.5458
232 34580.2165 22391.9327
233 16290.7635 34580.2165
234 34580.2165 16290.7635
235 34985.5410 34580.2165
236 34580.2165 34985.5410
237 36566.3486 34580.2165
238 34580.2165 36566.3486
239 34580.2165 34580.2165
240 34580.2165 34580.2165
241 29189.2405 34580.2165
242 34580.2165 29189.2405
243 39155.2379 34580.2165
244 -48965.4619 39155.2379
245 34580.2165 -48965.4619
246 34580.2165 34580.2165
247 26517.1760 34580.2165
248 34580.2165 26517.1760
249 34580.2165 34580.2165
250 28561.0884 34580.2165
251 34580.2165 28561.0884
252 15249.6882 34580.2165
253 9795.8940 15249.6882
254 -103337.7607 9795.8940
255 34580.2165 -103337.7607
256 32804.4768 34580.2165
257 46896.8111 32804.4768
258 34580.2165 46896.8111
259 34580.2165 34580.2165
260 34580.2165 34580.2165
261 34580.2165 34580.2165
262 36942.0807 34580.2165
263 51705.1001 36942.0807
264 34580.2165 51705.1001
265 34580.2165 34580.2165
266 34580.2165 34580.2165
267 65291.3572 34580.2165
268 34580.2165 65291.3572
269 30693.4734 34580.2165
270 34580.2165 30693.4734
271 34580.2165 34580.2165
272 34580.2165 34580.2165
273 30334.8589 34580.2165
274 34580.2165 30334.8589
275 39394.4755 34580.2165
276 35947.1018 39394.4755
277 34580.2165 35947.1018
278 5052.6817 34580.2165
279 21062.1530 5052.6817
280 34580.2165 21062.1530
281 34580.2165 34580.2165
282 62309.8561 34580.2165
283 34580.2165 62309.8561
284 25566.3500 34580.2165
285 34580.2165 25566.3500
286 22734.3419 34580.2165
287 35362.3736 22734.3419
288 34580.2165 35362.3736
289 42333.9750 34580.2165
290 34580.2165 42333.9750
291 44873.1463 34580.2165
292 39192.5070 44873.1463
293 34580.2165 39192.5070
294 28105.2768 34580.2165
295 59868.2221 28105.2768
296 34580.2165 59868.2221
297 34580.2165 34580.2165
298 34580.2165 34580.2165
299 34580.2165 34580.2165
300 32922.8350 34580.2165
301 16298.6976 32922.8350
302 34580.2165 16298.6976
303 34580.2165 34580.2165
304 71866.7354 34580.2165
305 38145.1602 71866.7354
306 34580.2165 38145.1602
307 28217.9952 34580.2165
308 38446.2874 28217.9952
309 34580.2165 38446.2874
310 34580.2165 34580.2165
311 35442.7760 34580.2165
312 34580.2165 35442.7760
313 36216.6003 34580.2165
314 36302.1434 36216.6003
315 49761.1157 36302.1434
316 35069.1596 49761.1157
317 11239.9555 35069.1596
318 36971.5570 11239.9555
319 34580.2165 36971.5570
320 34580.2165 34580.2165
321 34580.2165 34580.2165
322 34580.2165 34580.2165
323 34580.2165 34580.2165
324 34580.2165 34580.2165
325 31844.6853 34580.2165
326 34580.2165 31844.6853
327 34580.2165 34580.2165
328 34580.2165 34580.2165
329 34580.2165 34580.2165
330 34580.2165 34580.2165
331 34580.2165 34580.2165
332 34580.2165 34580.2165
333 32742.5135 34580.2165
334 35353.2460 32742.5135
335 27135.7859 35353.2460
336 21063.1600 27135.7859
337 34580.2165 21063.1600
338 32349.5558 34580.2165
339 40364.8296 32349.5558
340 7192.9443 40364.8296
341 34580.2165 7192.9443
342 34580.2165 34580.2165
343 34580.2165 34580.2165
344 28170.1365 34580.2165
345 34580.2165 28170.1365
346 24924.7505 34580.2165
347 34580.2165 24924.7505
348 34580.2165 34580.2165
349 -2597.6032 34580.2165
350 38401.0084 -2597.6032
351 34679.5158 38401.0084
352 34322.9287 34679.5158
353 27462.3996 34322.9287
354 34580.2165 27462.3996
355 34580.2165 34580.2165
356 11006.8190 34580.2165
357 34580.2165 11006.8190
358 34580.2165 34580.2165
359 34580.2165 34580.2165
360 36182.4019 34580.2165
361 34580.2165 36182.4019
362 34580.2165 34580.2165
363 -17619.3897 34580.2165
364 36785.0008 -17619.3897
365 34580.2165 36785.0008
366 34580.2165 34580.2165
367 34580.2165 34580.2165
368 34580.2165 34580.2165
369 36682.7327 34580.2165
370 34580.2165 36682.7327
371 47127.9780 34580.2165
372 20788.8182 47127.9780
373 13586.0443 20788.8182
374 34580.2165 13586.0443
375 34580.2165 34580.2165
376 29720.9807 34580.2165
377 28274.6000 29720.9807
378 -21901.6234 28274.6000
379 34580.2165 -21901.6234
380 34580.2165 34580.2165
381 19370.4612 34580.2165
382 53257.8226 19370.4612
383 34580.2165 53257.8226
384 36693.1984 34580.2165
385 -115636.0426 36693.1984
386 46871.0853 -115636.0426
387 7872.9811 46871.0853
388 33099.2810 7872.9811
389 20654.4427 33099.2810
390 29633.8067 20654.4427
391 3808.4378 29633.8067
392 -33887.7287 3808.4378
393 54504.4319 -33887.7287
394 -141011.2663 54504.4319
395 -4353.0222 -141011.2663
396 -152446.6167 -4353.0222
397 41068.9784 -152446.6167
398 -3133.3067 41068.9784
399 21066.1182 -3133.3067
400 -25877.9536 21066.1182
401 25950.4200 -25877.9536
402 4011.8418 25950.4200
403 -54926.7608 4011.8418
404 -26003.2963 -54926.7608
405 -4089.2056 -26003.2963
406 9440.9762 -4089.2056
407 1995.5438 9440.9762
408 13354.0051 1995.5438
409 -88117.0685 13354.0051
410 -4236.9658 -88117.0685
411 -210650.2878 -4236.9658
412 -32429.7526 -210650.2878
413 52064.9884 -32429.7526
414 -204568.3343 52064.9884
415 -129873.5310 -204568.3343
416 105933.4123 -129873.5310
417 -20765.0510 105933.4123
418 -50199.8918 -20765.0510
419 -63357.6347 -50199.8918
420 -75139.5888 -63357.6347
421 -122393.5869 -75139.5888
422 -532433.2423 -122393.5869
423 -499218.3037 -532433.2423
424 -55479.1464 -499218.3037
425 -192725.5730 -55479.1464
426 -327426.9457 -192725.5730
427 -1150522.8360 -327426.9457
428 -324119.0097 -1150522.8360
429 -364142.2036 -324119.0097
430 -692656.8969 -364142.2036
431 NA -692656.8969
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] 3225800.6189 1237333.7716
[2,] 1194235.0281 3225800.6189
[3,] -2198709.5189 1194235.0281
[4,] -446207.8858 -2198709.5189
[5,] -65315.8652 -446207.8858
[6,] -58061.4111 -65315.8652
[7,] 79150.5670 -58061.4111
[8,] -249447.6804 79150.5670
[9,] 742457.9366 -249447.6804
[10,] -308956.7059 742457.9366
[11,] 392076.2392 -308956.7059
[12,] 803384.0750 392076.2392
[13,] -1209577.7725 803384.0750
[14,] 35999.2506 -1209577.7725
[15,] 475470.1181 35999.2506
[16,] -267414.2259 475470.1181
[17,] 17699.4160 -267414.2259
[18,] -274855.9885 17699.4160
[19,] -346491.1382 -274855.9885
[20,] 190089.0283 -346491.1382
[21,] -245640.1815 190089.0283
[22,] -148253.8127 -245640.1815
[23,] -401854.6120 -148253.8127
[24,] -302236.3074 -401854.6120
[25,] -152079.9239 -302236.3074
[26,] -370608.0313 -152079.9239
[27,] -297472.3304 -370608.0313
[28,] 319888.6245 -297472.3304
[29,] -107458.5634 319888.6245
[30,] 264934.4893 -107458.5634
[31,] -118886.7400 264934.4893
[32,] 228735.6187 -118886.7400
[33,] 117862.7404 228735.6187
[34,] -176569.1359 117862.7404
[35,] -482482.7536 -176569.1359
[36,] 29495.2712 -482482.7536
[37,] 185164.6824 29495.2712
[38,] -31633.2784 185164.6824
[39,] -175484.6937 -31633.2784
[40,] -28586.2323 -175484.6937
[41,] -395020.7371 -28586.2323
[42,] 109887.2197 -395020.7371
[43,] -178865.3985 109887.2197
[44,] -24526.8599 -178865.3985
[45,] -86659.7032 -24526.8599
[46,] -269247.8351 -86659.7032
[47,] -97824.0452 -269247.8351
[48,] -47758.9998 -97824.0452
[49,] -173126.2485 -47758.9998
[50,] 103613.4092 -173126.2485
[51,] 25139.9739 103613.4092
[52,] -218762.7914 25139.9739
[53,] 106580.6851 -218762.7914
[54,] -315223.4826 106580.6851
[55,] 91349.0492 -315223.4826
[56,] 283612.8987 91349.0492
[57,] -23682.0837 283612.8987
[58,] -12429.1879 -23682.0837
[59,] -186594.5055 -12429.1879
[60,] -3288.1749 -186594.5055
[61,] -164170.1642 -3288.1749
[62,] -14535.6869 -164170.1642
[63,] -152468.2407 -14535.6869
[64,] 29170.0161 -152468.2407
[65,] 658.3290 29170.0161
[66,] -56301.8915 658.3290
[67,] -3474.4077 -56301.8915
[68,] -51538.3527 -3474.4077
[69,] -31098.5552 -51538.3527
[70,] -80208.5263 -31098.5552
[71,] 6724.8692 -80208.5263
[72,] -91727.8629 6724.8692
[73,] -4435.9069 -91727.8629
[74,] -47964.1058 -4435.9069
[75,] -112625.8705 -47964.1058
[76,] -32306.1501 -112625.8705
[77,] -77635.7304 -32306.1501
[78,] 22864.3015 -77635.7304
[79,] -6161.5256 22864.3015
[80,] 35011.7284 -6161.5256
[81,] -23890.0975 35011.7284
[82,] -86997.5856 -23890.0975
[83,] 46042.0802 -86997.5856
[84,] 3276.6093 46042.0802
[85,] -72852.1034 3276.6093
[86,] 43537.4967 -72852.1034
[87,] -23669.5973 43537.4967
[88,] -101295.9632 -23669.5973
[89,] 7517.6045 -101295.9632
[90,] 3632.1632 7517.6045
[91,] -18907.1836 3632.1632
[92,] 138110.3578 -18907.1836
[93,] -214809.1834 138110.3578
[94,] 85549.2598 -214809.1834
[95,] 22161.4187 85549.2598
[96,] 38786.3001 22161.4187
[97,] 68328.4234 38786.3001
[98,] -84437.6455 68328.4234
[99,] -25122.8216 -84437.6455
[100,] -28698.3104 -25122.8216
[101,] -125369.5171 -28698.3104
[102,] -146117.0270 -125369.5171
[103,] -82896.3085 -146117.0270
[104,] 36828.6130 -82896.3085
[105,] -45990.5534 36828.6130
[106,] -225506.3353 -45990.5534
[107,] 102172.4638 -225506.3353
[108,] -187517.4506 102172.4638
[109,] -99005.2886 -187517.4506
[110,] 97583.6413 -99005.2886
[111,] -181421.3633 97583.6413
[112,] -165268.3173 -181421.3633
[113,] -84430.1190 -165268.3173
[114,] 65432.0980 -84430.1190
[115,] 41884.9294 65432.0980
[116,] 1804.1853 41884.9294
[117,] 47384.8091 1804.1853
[118,] -160370.2516 47384.8091
[119,] 36507.6963 -160370.2516
[120,] -286240.9429 36507.6963
[121,] 20225.1998 -286240.9429
[122,] 5840.5691 20225.1998
[123,] 67611.1743 5840.5691
[124,] 3093.7535 67611.1743
[125,] -183753.0522 3093.7535
[126,] 34580.2165 -183753.0522
[127,] 33587.7584 34580.2165
[128,] 41443.9094 33587.7584
[129,] 34580.2165 41443.9094
[130,] 35945.7172 34580.2165
[131,] 77374.3357 35945.7172
[132,] 34580.2165 77374.3357
[133,] 31630.6387 34580.2165
[134,] 13869.3958 31630.6387
[135,] 28761.9295 13869.3958
[136,] 39448.6018 28761.9295
[137,] 34745.7604 39448.6018
[138,] 9116.6456 34745.7604
[139,] 34580.2165 9116.6456
[140,] 42113.4674 34580.2165
[141,] 34580.2165 42113.4674
[142,] 1596.0722 34580.2165
[143,] 34815.4903 1596.0722
[144,] 35943.7690 34815.4903
[145,] 35592.4067 35943.7690
[146,] 34580.2165 35592.4067
[147,] 36113.5379 34580.2165
[148,] 37095.3019 36113.5379
[149,] 34467.9757 37095.3019
[150,] 3373.2999 34467.9757
[151,] 2061.9484 3373.2999
[152,] 34580.2165 2061.9484
[153,] 11810.4517 34580.2165
[154,] 34620.4264 11810.4517
[155,] 34580.2165 34620.4264
[156,] 34580.2165 34580.2165
[157,] 30753.1273 34580.2165
[158,] 46007.7895 30753.1273
[159,] 34580.2165 46007.7895
[160,] 34580.2165 34580.2165
[161,] 34580.2165 34580.2165
[162,] 34580.2165 34580.2165
[163,] 21622.8177 34580.2165
[164,] 23262.3055 21622.8177
[165,] 33880.2124 23262.3055
[166,] 34580.2165 33880.2124
[167,] 116864.3690 34580.2165
[168,] 79550.9637 116864.3690
[169,] 34580.2165 79550.9637
[170,] 34580.2165 34580.2165
[171,] 30457.8195 34580.2165
[172,] 34580.2165 30457.8195
[173,] 23886.1881 34580.2165
[174,] 34580.2165 23886.1881
[175,] 34535.5262 34580.2165
[176,] 81914.1147 34535.5262
[177,] 33592.1205 81914.1147
[178,] 34580.2165 33592.1205
[179,] 34580.2165 34580.2165
[180,] 41164.1864 34580.2165
[181,] 15473.1509 41164.1864
[182,] -7722.3720 15473.1509
[183,] 3954.9700 -7722.3720
[184,] 31665.6161 3954.9700
[185,] 39556.5727 31665.6161
[186,] 5452.7572 39556.5727
[187,] 27424.4312 5452.7572
[188,] -3849.3352 27424.4312
[189,] 34580.2165 -3849.3352
[190,] 33720.4193 34580.2165
[191,] 34580.2165 33720.4193
[192,] 30926.1365 34580.2165
[193,] 34580.2165 30926.1365
[194,] 34513.8910 34580.2165
[195,] 34580.2165 34513.8910
[196,] -89970.4929 34580.2165
[197,] 34580.2165 -89970.4929
[198,] 34580.2165 34580.2165
[199,] -55995.9889 34580.2165
[200,] 34580.2165 -55995.9889
[201,] 34580.2165 34580.2165
[202,] 34580.2165 34580.2165
[203,] 44806.6942 34580.2165
[204,] -16733.3662 44806.6942
[205,] 33990.8040 -16733.3662
[206,] 27556.9989 33990.8040
[207,] 34580.2165 27556.9989
[208,] 65941.5748 34580.2165
[209,] 30707.7678 65941.5748
[210,] 34580.2165 30707.7678
[211,] 34580.2165 34580.2165
[212,] 60324.2404 34580.2165
[213,] -515.5477 60324.2404
[214,] 6631.3916 -515.5477
[215,] 34580.2165 6631.3916
[216,] 37450.9613 34580.2165
[217,] 36453.9569 37450.9613
[218,] 34580.2165 36453.9569
[219,] 6047.9156 34580.2165
[220,] 9219.4592 6047.9156
[221,] 34580.2165 9219.4592
[222,] 35317.0126 34580.2165
[223,] 25760.2266 35317.0126
[224,] 34812.8986 25760.2266
[225,] 34899.0433 34812.8986
[226,] 37577.8514 34899.0433
[227,] 34580.2165 37577.8514
[228,] 18260.4102 34580.2165
[229,] 34580.2165 18260.4102
[230,] 5601.5458 34580.2165
[231,] 22391.9327 5601.5458
[232,] 34580.2165 22391.9327
[233,] 16290.7635 34580.2165
[234,] 34580.2165 16290.7635
[235,] 34985.5410 34580.2165
[236,] 34580.2165 34985.5410
[237,] 36566.3486 34580.2165
[238,] 34580.2165 36566.3486
[239,] 34580.2165 34580.2165
[240,] 34580.2165 34580.2165
[241,] 29189.2405 34580.2165
[242,] 34580.2165 29189.2405
[243,] 39155.2379 34580.2165
[244,] -48965.4619 39155.2379
[245,] 34580.2165 -48965.4619
[246,] 34580.2165 34580.2165
[247,] 26517.1760 34580.2165
[248,] 34580.2165 26517.1760
[249,] 34580.2165 34580.2165
[250,] 28561.0884 34580.2165
[251,] 34580.2165 28561.0884
[252,] 15249.6882 34580.2165
[253,] 9795.8940 15249.6882
[254,] -103337.7607 9795.8940
[255,] 34580.2165 -103337.7607
[256,] 32804.4768 34580.2165
[257,] 46896.8111 32804.4768
[258,] 34580.2165 46896.8111
[259,] 34580.2165 34580.2165
[260,] 34580.2165 34580.2165
[261,] 34580.2165 34580.2165
[262,] 36942.0807 34580.2165
[263,] 51705.1001 36942.0807
[264,] 34580.2165 51705.1001
[265,] 34580.2165 34580.2165
[266,] 34580.2165 34580.2165
[267,] 65291.3572 34580.2165
[268,] 34580.2165 65291.3572
[269,] 30693.4734 34580.2165
[270,] 34580.2165 30693.4734
[271,] 34580.2165 34580.2165
[272,] 34580.2165 34580.2165
[273,] 30334.8589 34580.2165
[274,] 34580.2165 30334.8589
[275,] 39394.4755 34580.2165
[276,] 35947.1018 39394.4755
[277,] 34580.2165 35947.1018
[278,] 5052.6817 34580.2165
[279,] 21062.1530 5052.6817
[280,] 34580.2165 21062.1530
[281,] 34580.2165 34580.2165
[282,] 62309.8561 34580.2165
[283,] 34580.2165 62309.8561
[284,] 25566.3500 34580.2165
[285,] 34580.2165 25566.3500
[286,] 22734.3419 34580.2165
[287,] 35362.3736 22734.3419
[288,] 34580.2165 35362.3736
[289,] 42333.9750 34580.2165
[290,] 34580.2165 42333.9750
[291,] 44873.1463 34580.2165
[292,] 39192.5070 44873.1463
[293,] 34580.2165 39192.5070
[294,] 28105.2768 34580.2165
[295,] 59868.2221 28105.2768
[296,] 34580.2165 59868.2221
[297,] 34580.2165 34580.2165
[298,] 34580.2165 34580.2165
[299,] 34580.2165 34580.2165
[300,] 32922.8350 34580.2165
[301,] 16298.6976 32922.8350
[302,] 34580.2165 16298.6976
[303,] 34580.2165 34580.2165
[304,] 71866.7354 34580.2165
[305,] 38145.1602 71866.7354
[306,] 34580.2165 38145.1602
[307,] 28217.9952 34580.2165
[308,] 38446.2874 28217.9952
[309,] 34580.2165 38446.2874
[310,] 34580.2165 34580.2165
[311,] 35442.7760 34580.2165
[312,] 34580.2165 35442.7760
[313,] 36216.6003 34580.2165
[314,] 36302.1434 36216.6003
[315,] 49761.1157 36302.1434
[316,] 35069.1596 49761.1157
[317,] 11239.9555 35069.1596
[318,] 36971.5570 11239.9555
[319,] 34580.2165 36971.5570
[320,] 34580.2165 34580.2165
[321,] 34580.2165 34580.2165
[322,] 34580.2165 34580.2165
[323,] 34580.2165 34580.2165
[324,] 34580.2165 34580.2165
[325,] 31844.6853 34580.2165
[326,] 34580.2165 31844.6853
[327,] 34580.2165 34580.2165
[328,] 34580.2165 34580.2165
[329,] 34580.2165 34580.2165
[330,] 34580.2165 34580.2165
[331,] 34580.2165 34580.2165
[332,] 34580.2165 34580.2165
[333,] 32742.5135 34580.2165
[334,] 35353.2460 32742.5135
[335,] 27135.7859 35353.2460
[336,] 21063.1600 27135.7859
[337,] 34580.2165 21063.1600
[338,] 32349.5558 34580.2165
[339,] 40364.8296 32349.5558
[340,] 7192.9443 40364.8296
[341,] 34580.2165 7192.9443
[342,] 34580.2165 34580.2165
[343,] 34580.2165 34580.2165
[344,] 28170.1365 34580.2165
[345,] 34580.2165 28170.1365
[346,] 24924.7505 34580.2165
[347,] 34580.2165 24924.7505
[348,] 34580.2165 34580.2165
[349,] -2597.6032 34580.2165
[350,] 38401.0084 -2597.6032
[351,] 34679.5158 38401.0084
[352,] 34322.9287 34679.5158
[353,] 27462.3996 34322.9287
[354,] 34580.2165 27462.3996
[355,] 34580.2165 34580.2165
[356,] 11006.8190 34580.2165
[357,] 34580.2165 11006.8190
[358,] 34580.2165 34580.2165
[359,] 34580.2165 34580.2165
[360,] 36182.4019 34580.2165
[361,] 34580.2165 36182.4019
[362,] 34580.2165 34580.2165
[363,] -17619.3897 34580.2165
[364,] 36785.0008 -17619.3897
[365,] 34580.2165 36785.0008
[366,] 34580.2165 34580.2165
[367,] 34580.2165 34580.2165
[368,] 34580.2165 34580.2165
[369,] 36682.7327 34580.2165
[370,] 34580.2165 36682.7327
[371,] 47127.9780 34580.2165
[372,] 20788.8182 47127.9780
[373,] 13586.0443 20788.8182
[374,] 34580.2165 13586.0443
[375,] 34580.2165 34580.2165
[376,] 29720.9807 34580.2165
[377,] 28274.6000 29720.9807
[378,] -21901.6234 28274.6000
[379,] 34580.2165 -21901.6234
[380,] 34580.2165 34580.2165
[381,] 19370.4612 34580.2165
[382,] 53257.8226 19370.4612
[383,] 34580.2165 53257.8226
[384,] 36693.1984 34580.2165
[385,] -115636.0426 36693.1984
[386,] 46871.0853 -115636.0426
[387,] 7872.9811 46871.0853
[388,] 33099.2810 7872.9811
[389,] 20654.4427 33099.2810
[390,] 29633.8067 20654.4427
[391,] 3808.4378 29633.8067
[392,] -33887.7287 3808.4378
[393,] 54504.4319 -33887.7287
[394,] -141011.2663 54504.4319
[395,] -4353.0222 -141011.2663
[396,] -152446.6167 -4353.0222
[397,] 41068.9784 -152446.6167
[398,] -3133.3067 41068.9784
[399,] 21066.1182 -3133.3067
[400,] -25877.9536 21066.1182
[401,] 25950.4200 -25877.9536
[402,] 4011.8418 25950.4200
[403,] -54926.7608 4011.8418
[404,] -26003.2963 -54926.7608
[405,] -4089.2056 -26003.2963
[406,] 9440.9762 -4089.2056
[407,] 1995.5438 9440.9762
[408,] 13354.0051 1995.5438
[409,] -88117.0685 13354.0051
[410,] -4236.9658 -88117.0685
[411,] -210650.2878 -4236.9658
[412,] -32429.7526 -210650.2878
[413,] 52064.9884 -32429.7526
[414,] -204568.3343 52064.9884
[415,] -129873.5310 -204568.3343
[416,] 105933.4123 -129873.5310
[417,] -20765.0510 105933.4123
[418,] -50199.8918 -20765.0510
[419,] -63357.6347 -50199.8918
[420,] -75139.5888 -63357.6347
[421,] -122393.5869 -75139.5888
[422,] -532433.2423 -122393.5869
[423,] -499218.3037 -532433.2423
[424,] -55479.1464 -499218.3037
[425,] -192725.5730 -55479.1464
[426,] -327426.9457 -192725.5730
[427,] -1150522.8360 -327426.9457
[428,] -324119.0097 -1150522.8360
[429,] -364142.2036 -324119.0097
[430,] -692656.8969 -364142.2036
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 3225800.6189 1237333.7716
2 1194235.0281 3225800.6189
3 -2198709.5189 1194235.0281
4 -446207.8858 -2198709.5189
5 -65315.8652 -446207.8858
6 -58061.4111 -65315.8652
7 79150.5670 -58061.4111
8 -249447.6804 79150.5670
9 742457.9366 -249447.6804
10 -308956.7059 742457.9366
11 392076.2392 -308956.7059
12 803384.0750 392076.2392
13 -1209577.7725 803384.0750
14 35999.2506 -1209577.7725
15 475470.1181 35999.2506
16 -267414.2259 475470.1181
17 17699.4160 -267414.2259
18 -274855.9885 17699.4160
19 -346491.1382 -274855.9885
20 190089.0283 -346491.1382
21 -245640.1815 190089.0283
22 -148253.8127 -245640.1815
23 -401854.6120 -148253.8127
24 -302236.3074 -401854.6120
25 -152079.9239 -302236.3074
26 -370608.0313 -152079.9239
27 -297472.3304 -370608.0313
28 319888.6245 -297472.3304
29 -107458.5634 319888.6245
30 264934.4893 -107458.5634
31 -118886.7400 264934.4893
32 228735.6187 -118886.7400
33 117862.7404 228735.6187
34 -176569.1359 117862.7404
35 -482482.7536 -176569.1359
36 29495.2712 -482482.7536
37 185164.6824 29495.2712
38 -31633.2784 185164.6824
39 -175484.6937 -31633.2784
40 -28586.2323 -175484.6937
41 -395020.7371 -28586.2323
42 109887.2197 -395020.7371
43 -178865.3985 109887.2197
44 -24526.8599 -178865.3985
45 -86659.7032 -24526.8599
46 -269247.8351 -86659.7032
47 -97824.0452 -269247.8351
48 -47758.9998 -97824.0452
49 -173126.2485 -47758.9998
50 103613.4092 -173126.2485
51 25139.9739 103613.4092
52 -218762.7914 25139.9739
53 106580.6851 -218762.7914
54 -315223.4826 106580.6851
55 91349.0492 -315223.4826
56 283612.8987 91349.0492
57 -23682.0837 283612.8987
58 -12429.1879 -23682.0837
59 -186594.5055 -12429.1879
60 -3288.1749 -186594.5055
61 -164170.1642 -3288.1749
62 -14535.6869 -164170.1642
63 -152468.2407 -14535.6869
64 29170.0161 -152468.2407
65 658.3290 29170.0161
66 -56301.8915 658.3290
67 -3474.4077 -56301.8915
68 -51538.3527 -3474.4077
69 -31098.5552 -51538.3527
70 -80208.5263 -31098.5552
71 6724.8692 -80208.5263
72 -91727.8629 6724.8692
73 -4435.9069 -91727.8629
74 -47964.1058 -4435.9069
75 -112625.8705 -47964.1058
76 -32306.1501 -112625.8705
77 -77635.7304 -32306.1501
78 22864.3015 -77635.7304
79 -6161.5256 22864.3015
80 35011.7284 -6161.5256
81 -23890.0975 35011.7284
82 -86997.5856 -23890.0975
83 46042.0802 -86997.5856
84 3276.6093 46042.0802
85 -72852.1034 3276.6093
86 43537.4967 -72852.1034
87 -23669.5973 43537.4967
88 -101295.9632 -23669.5973
89 7517.6045 -101295.9632
90 3632.1632 7517.6045
91 -18907.1836 3632.1632
92 138110.3578 -18907.1836
93 -214809.1834 138110.3578
94 85549.2598 -214809.1834
95 22161.4187 85549.2598
96 38786.3001 22161.4187
97 68328.4234 38786.3001
98 -84437.6455 68328.4234
99 -25122.8216 -84437.6455
100 -28698.3104 -25122.8216
101 -125369.5171 -28698.3104
102 -146117.0270 -125369.5171
103 -82896.3085 -146117.0270
104 36828.6130 -82896.3085
105 -45990.5534 36828.6130
106 -225506.3353 -45990.5534
107 102172.4638 -225506.3353
108 -187517.4506 102172.4638
109 -99005.2886 -187517.4506
110 97583.6413 -99005.2886
111 -181421.3633 97583.6413
112 -165268.3173 -181421.3633
113 -84430.1190 -165268.3173
114 65432.0980 -84430.1190
115 41884.9294 65432.0980
116 1804.1853 41884.9294
117 47384.8091 1804.1853
118 -160370.2516 47384.8091
119 36507.6963 -160370.2516
120 -286240.9429 36507.6963
121 20225.1998 -286240.9429
122 5840.5691 20225.1998
123 67611.1743 5840.5691
124 3093.7535 67611.1743
125 -183753.0522 3093.7535
126 34580.2165 -183753.0522
127 33587.7584 34580.2165
128 41443.9094 33587.7584
129 34580.2165 41443.9094
130 35945.7172 34580.2165
131 77374.3357 35945.7172
132 34580.2165 77374.3357
133 31630.6387 34580.2165
134 13869.3958 31630.6387
135 28761.9295 13869.3958
136 39448.6018 28761.9295
137 34745.7604 39448.6018
138 9116.6456 34745.7604
139 34580.2165 9116.6456
140 42113.4674 34580.2165
141 34580.2165 42113.4674
142 1596.0722 34580.2165
143 34815.4903 1596.0722
144 35943.7690 34815.4903
145 35592.4067 35943.7690
146 34580.2165 35592.4067
147 36113.5379 34580.2165
148 37095.3019 36113.5379
149 34467.9757 37095.3019
150 3373.2999 34467.9757
151 2061.9484 3373.2999
152 34580.2165 2061.9484
153 11810.4517 34580.2165
154 34620.4264 11810.4517
155 34580.2165 34620.4264
156 34580.2165 34580.2165
157 30753.1273 34580.2165
158 46007.7895 30753.1273
159 34580.2165 46007.7895
160 34580.2165 34580.2165
161 34580.2165 34580.2165
162 34580.2165 34580.2165
163 21622.8177 34580.2165
164 23262.3055 21622.8177
165 33880.2124 23262.3055
166 34580.2165 33880.2124
167 116864.3690 34580.2165
168 79550.9637 116864.3690
169 34580.2165 79550.9637
170 34580.2165 34580.2165
171 30457.8195 34580.2165
172 34580.2165 30457.8195
173 23886.1881 34580.2165
174 34580.2165 23886.1881
175 34535.5262 34580.2165
176 81914.1147 34535.5262
177 33592.1205 81914.1147
178 34580.2165 33592.1205
179 34580.2165 34580.2165
180 41164.1864 34580.2165
181 15473.1509 41164.1864
182 -7722.3720 15473.1509
183 3954.9700 -7722.3720
184 31665.6161 3954.9700
185 39556.5727 31665.6161
186 5452.7572 39556.5727
187 27424.4312 5452.7572
188 -3849.3352 27424.4312
189 34580.2165 -3849.3352
190 33720.4193 34580.2165
191 34580.2165 33720.4193
192 30926.1365 34580.2165
193 34580.2165 30926.1365
194 34513.8910 34580.2165
195 34580.2165 34513.8910
196 -89970.4929 34580.2165
197 34580.2165 -89970.4929
198 34580.2165 34580.2165
199 -55995.9889 34580.2165
200 34580.2165 -55995.9889
201 34580.2165 34580.2165
202 34580.2165 34580.2165
203 44806.6942 34580.2165
204 -16733.3662 44806.6942
205 33990.8040 -16733.3662
206 27556.9989 33990.8040
207 34580.2165 27556.9989
208 65941.5748 34580.2165
209 30707.7678 65941.5748
210 34580.2165 30707.7678
211 34580.2165 34580.2165
212 60324.2404 34580.2165
213 -515.5477 60324.2404
214 6631.3916 -515.5477
215 34580.2165 6631.3916
216 37450.9613 34580.2165
217 36453.9569 37450.9613
218 34580.2165 36453.9569
219 6047.9156 34580.2165
220 9219.4592 6047.9156
221 34580.2165 9219.4592
222 35317.0126 34580.2165
223 25760.2266 35317.0126
224 34812.8986 25760.2266
225 34899.0433 34812.8986
226 37577.8514 34899.0433
227 34580.2165 37577.8514
228 18260.4102 34580.2165
229 34580.2165 18260.4102
230 5601.5458 34580.2165
231 22391.9327 5601.5458
232 34580.2165 22391.9327
233 16290.7635 34580.2165
234 34580.2165 16290.7635
235 34985.5410 34580.2165
236 34580.2165 34985.5410
237 36566.3486 34580.2165
238 34580.2165 36566.3486
239 34580.2165 34580.2165
240 34580.2165 34580.2165
241 29189.2405 34580.2165
242 34580.2165 29189.2405
243 39155.2379 34580.2165
244 -48965.4619 39155.2379
245 34580.2165 -48965.4619
246 34580.2165 34580.2165
247 26517.1760 34580.2165
248 34580.2165 26517.1760
249 34580.2165 34580.2165
250 28561.0884 34580.2165
251 34580.2165 28561.0884
252 15249.6882 34580.2165
253 9795.8940 15249.6882
254 -103337.7607 9795.8940
255 34580.2165 -103337.7607
256 32804.4768 34580.2165
257 46896.8111 32804.4768
258 34580.2165 46896.8111
259 34580.2165 34580.2165
260 34580.2165 34580.2165
261 34580.2165 34580.2165
262 36942.0807 34580.2165
263 51705.1001 36942.0807
264 34580.2165 51705.1001
265 34580.2165 34580.2165
266 34580.2165 34580.2165
267 65291.3572 34580.2165
268 34580.2165 65291.3572
269 30693.4734 34580.2165
270 34580.2165 30693.4734
271 34580.2165 34580.2165
272 34580.2165 34580.2165
273 30334.8589 34580.2165
274 34580.2165 30334.8589
275 39394.4755 34580.2165
276 35947.1018 39394.4755
277 34580.2165 35947.1018
278 5052.6817 34580.2165
279 21062.1530 5052.6817
280 34580.2165 21062.1530
281 34580.2165 34580.2165
282 62309.8561 34580.2165
283 34580.2165 62309.8561
284 25566.3500 34580.2165
285 34580.2165 25566.3500
286 22734.3419 34580.2165
287 35362.3736 22734.3419
288 34580.2165 35362.3736
289 42333.9750 34580.2165
290 34580.2165 42333.9750
291 44873.1463 34580.2165
292 39192.5070 44873.1463
293 34580.2165 39192.5070
294 28105.2768 34580.2165
295 59868.2221 28105.2768
296 34580.2165 59868.2221
297 34580.2165 34580.2165
298 34580.2165 34580.2165
299 34580.2165 34580.2165
300 32922.8350 34580.2165
301 16298.6976 32922.8350
302 34580.2165 16298.6976
303 34580.2165 34580.2165
304 71866.7354 34580.2165
305 38145.1602 71866.7354
306 34580.2165 38145.1602
307 28217.9952 34580.2165
308 38446.2874 28217.9952
309 34580.2165 38446.2874
310 34580.2165 34580.2165
311 35442.7760 34580.2165
312 34580.2165 35442.7760
313 36216.6003 34580.2165
314 36302.1434 36216.6003
315 49761.1157 36302.1434
316 35069.1596 49761.1157
317 11239.9555 35069.1596
318 36971.5570 11239.9555
319 34580.2165 36971.5570
320 34580.2165 34580.2165
321 34580.2165 34580.2165
322 34580.2165 34580.2165
323 34580.2165 34580.2165
324 34580.2165 34580.2165
325 31844.6853 34580.2165
326 34580.2165 31844.6853
327 34580.2165 34580.2165
328 34580.2165 34580.2165
329 34580.2165 34580.2165
330 34580.2165 34580.2165
331 34580.2165 34580.2165
332 34580.2165 34580.2165
333 32742.5135 34580.2165
334 35353.2460 32742.5135
335 27135.7859 35353.2460
336 21063.1600 27135.7859
337 34580.2165 21063.1600
338 32349.5558 34580.2165
339 40364.8296 32349.5558
340 7192.9443 40364.8296
341 34580.2165 7192.9443
342 34580.2165 34580.2165
343 34580.2165 34580.2165
344 28170.1365 34580.2165
345 34580.2165 28170.1365
346 24924.7505 34580.2165
347 34580.2165 24924.7505
348 34580.2165 34580.2165
349 -2597.6032 34580.2165
350 38401.0084 -2597.6032
351 34679.5158 38401.0084
352 34322.9287 34679.5158
353 27462.3996 34322.9287
354 34580.2165 27462.3996
355 34580.2165 34580.2165
356 11006.8190 34580.2165
357 34580.2165 11006.8190
358 34580.2165 34580.2165
359 34580.2165 34580.2165
360 36182.4019 34580.2165
361 34580.2165 36182.4019
362 34580.2165 34580.2165
363 -17619.3897 34580.2165
364 36785.0008 -17619.3897
365 34580.2165 36785.0008
366 34580.2165 34580.2165
367 34580.2165 34580.2165
368 34580.2165 34580.2165
369 36682.7327 34580.2165
370 34580.2165 36682.7327
371 47127.9780 34580.2165
372 20788.8182 47127.9780
373 13586.0443 20788.8182
374 34580.2165 13586.0443
375 34580.2165 34580.2165
376 29720.9807 34580.2165
377 28274.6000 29720.9807
378 -21901.6234 28274.6000
379 34580.2165 -21901.6234
380 34580.2165 34580.2165
381 19370.4612 34580.2165
382 53257.8226 19370.4612
383 34580.2165 53257.8226
384 36693.1984 34580.2165
385 -115636.0426 36693.1984
386 46871.0853 -115636.0426
387 7872.9811 46871.0853
388 33099.2810 7872.9811
389 20654.4427 33099.2810
390 29633.8067 20654.4427
391 3808.4378 29633.8067
392 -33887.7287 3808.4378
393 54504.4319 -33887.7287
394 -141011.2663 54504.4319
395 -4353.0222 -141011.2663
396 -152446.6167 -4353.0222
397 41068.9784 -152446.6167
398 -3133.3067 41068.9784
399 21066.1182 -3133.3067
400 -25877.9536 21066.1182
401 25950.4200 -25877.9536
402 4011.8418 25950.4200
403 -54926.7608 4011.8418
404 -26003.2963 -54926.7608
405 -4089.2056 -26003.2963
406 9440.9762 -4089.2056
407 1995.5438 9440.9762
408 13354.0051 1995.5438
409 -88117.0685 13354.0051
410 -4236.9658 -88117.0685
411 -210650.2878 -4236.9658
412 -32429.7526 -210650.2878
413 52064.9884 -32429.7526
414 -204568.3343 52064.9884
415 -129873.5310 -204568.3343
416 105933.4123 -129873.5310
417 -20765.0510 105933.4123
418 -50199.8918 -20765.0510
419 -63357.6347 -50199.8918
420 -75139.5888 -63357.6347
421 -122393.5869 -75139.5888
422 -532433.2423 -122393.5869
423 -499218.3037 -532433.2423
424 -55479.1464 -499218.3037
425 -192725.5730 -55479.1464
426 -327426.9457 -192725.5730
427 -1150522.8360 -327426.9457
428 -324119.0097 -1150522.8360
429 -364142.2036 -324119.0097
430 -692656.8969 -364142.2036
> 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/7s46o1291211846.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/8s46o1291211846.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/9kv5r1291211846.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/10kv5r1291211846.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/116vmx1291211846.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/12rwkl1291211846.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/13yxhf1291211846.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/1496z01291211846.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/15cpfo1291211846.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/168hve1291211846.tab")
+ }
>
> try(system("convert tmp/1ec8y1291211846.ps tmp/1ec8y1291211846.png",intern=TRUE))
character(0)
> try(system("convert tmp/263811291211846.ps tmp/263811291211846.png",intern=TRUE))
character(0)
> try(system("convert tmp/363811291211846.ps tmp/363811291211846.png",intern=TRUE))
character(0)
> try(system("convert tmp/463811291211846.ps tmp/463811291211846.png",intern=TRUE))
character(0)
> try(system("convert tmp/5hvpl1291211846.ps tmp/5hvpl1291211846.png",intern=TRUE))
character(0)
> try(system("convert tmp/6hvpl1291211846.ps tmp/6hvpl1291211846.png",intern=TRUE))
character(0)
> try(system("convert tmp/7s46o1291211846.ps tmp/7s46o1291211846.png",intern=TRUE))
character(0)
> try(system("convert tmp/8s46o1291211846.ps tmp/8s46o1291211846.png",intern=TRUE))
character(0)
> try(system("convert tmp/9kv5r1291211846.ps tmp/9kv5r1291211846.png",intern=TRUE))
character(0)
> try(system("convert tmp/10kv5r1291211846.ps tmp/10kv5r1291211846.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
13.352 3.110 13.872