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
+ ,807
+ ,213118
+ ,6282154
+ ,29790
+ ,444
+ ,81767
+ ,4321023
+ ,87550
+ ,412
+ ,153198
+ ,4111912
+ ,84738
+ ,428
+ ,-26007
+ ,223193
+ ,54660
+ ,315
+ ,126942
+ ,1491348
+ ,42634
+ ,168
+ ,157214
+ ,1629616
+ ,40949
+ ,263
+ ,129352
+ ,1398893
+ ,45187
+ ,267
+ ,234817
+ ,1926517
+ ,37704
+ ,228
+ ,60448
+ ,983660
+ ,16275
+ ,129
+ ,47818
+ ,1443586
+ ,25830
+ ,104
+ ,245546
+ ,1073089
+ ,12679
+ ,122
+ ,48020
+ ,984885
+ ,18014
+ ,393
+ ,-1710
+ ,1405225
+ ,43556
+ ,190
+ ,32648
+ ,227132
+ ,24811
+ ,280
+ ,95350
+ ,929118
+ ,6575
+ ,63
+ ,151352
+ ,1071292
+ ,7123
+ ,102
+ ,288170
+ ,638830
+ ,21950
+ ,265
+ ,114337
+ ,856956
+ ,37597
+ ,234
+ ,37884
+ ,992426
+ ,17821
+ ,277
+ ,122844
+ ,444477
+ ,12988
+ ,73
+ ,82340
+ ,857217
+ ,22330
+ ,67
+ ,79801
+ ,711969
+ ,13326
+ ,103
+ ,165548
+ ,702380
+ ,16189
+ ,290
+ ,116384
+ ,358589
+ ,7146
+ ,83
+ ,134028
+ ,297978
+ ,15824
+ ,56
+ ,63838
+ ,585715
+ ,27664
+ ,236
+ ,74996
+ ,657954
+ ,11920
+ ,73
+ ,31080
+ ,209458
+ ,8568
+ ,34
+ ,32168
+ ,786690
+ ,14416
+ ,139
+ ,49857
+ ,439798
+ ,3369
+ ,26
+ ,87161
+ ,688779
+ ,11819
+ ,70
+ ,106113
+ ,574339
+ ,6984
+ ,40
+ ,80570
+ ,741409
+ ,4519
+ ,42
+ ,102129
+ ,597793
+ ,2220
+ ,12
+ ,301670
+ ,644190
+ ,18562
+ ,211
+ ,102313
+ ,377934
+ ,10327
+ ,74
+ ,88577
+ ,640273
+ ,5336
+ ,80
+ ,112477
+ ,697458
+ ,2365
+ ,83
+ ,191778
+ ,550608
+ ,4069
+ ,131
+ ,79804
+ ,207393
+ ,8636
+ ,203
+ ,128294
+ ,301607
+ ,13718
+ ,56
+ ,96448
+ ,345783
+ ,4525
+ ,89
+ ,93811
+ ,501749
+ ,6869
+ ,88
+ ,117520
+ ,379983
+ ,4628
+ ,39
+ ,69159
+ ,387475
+ ,3689
+ ,25
+ ,101792
+ ,377305
+ ,4891
+ ,49
+ ,210568
+ ,370837
+ ,7489
+ ,149
+ ,136996
+ ,430866
+ ,4901
+ ,58
+ ,121920
+ ,469107
+ ,2284
+ ,41
+ ,76403
+ ,194493
+ ,3160
+ ,90
+ ,108094
+ ,530670
+ ,4150
+ ,136
+ ,134759
+ ,518365
+ ,7285
+ ,97
+ ,188873
+ ,491303
+ ,1134
+ ,63
+ ,146216
+ ,527021
+ ,4658
+ ,114
+ ,156608
+ ,233773
+ ,2384
+ ,77
+ ,61348
+ ,405972
+ ,3748
+ ,6
+ ,50350
+ ,652925
+ ,5371
+ ,47
+ ,87720
+ ,446211
+ ,1285
+ ,51
+ ,99489
+ ,341340
+ ,9327
+ ,85
+ ,87419
+ ,387699
+ ,5565
+ ,43
+ ,94355
+ ,493408
+ ,1528
+ ,32
+ ,60326
+ ,146494
+ ,3122
+ ,25
+ ,94670
+ ,414462
+ ,7561
+ ,77
+ ,82425
+ ,364304
+ ,2675
+ ,54
+ ,59017
+ ,355178
+ ,13253
+ ,251
+ ,90829
+ ,357760
+ ,880
+ ,15
+ ,80791
+ ,261216
+ ,2053
+ ,44
+ ,100423
+ ,397144
+ ,1424
+ ,73
+ ,131116
+ ,374943
+ ,4036
+ ,85
+ ,100269
+ ,424898
+ ,3045
+ ,49
+ ,27330
+ ,202055
+ ,5119
+ ,38
+ ,39039
+ ,378525
+ ,1431
+ ,35
+ ,106885
+ ,310768
+ ,554
+ ,9
+ ,79285
+ ,325738
+ ,1975
+ ,34
+ ,118881
+ ,394510
+ ,1765
+ ,20
+ ,77623
+ ,247060
+ ,1012
+ ,29
+ ,114768
+ ,368078
+ ,810
+ ,11
+ ,74015
+ ,236761
+ ,1280
+ ,52
+ ,69465
+ ,312378
+ ,666
+ ,13
+ ,117869
+ ,339836
+ ,1380
+ ,29
+ ,60982
+ ,347385
+ ,4677
+ ,66
+ ,90131
+ ,426280
+ ,876
+ ,33
+ ,138971
+ ,352850
+ ,814
+ ,15
+ ,39625
+ ,301881
+ ,514
+ ,15
+ ,102725
+ ,377516
+ ,5692
+ ,68
+ ,64239
+ ,357312
+ ,3642
+ ,100
+ ,90262
+ ,458343
+ ,540
+ ,13
+ ,103960
+ ,354228
+ ,2099
+ ,45
+ ,106611
+ ,308636
+ ,567
+ ,14
+ ,103345
+ ,386212
+ ,2001
+ ,36
+ ,95551
+ ,393343
+ ,2949
+ ,40
+ ,82903
+ ,378509
+ ,2253
+ ,68
+ ,63593
+ ,452469
+ ,6533
+ ,29
+ ,126910
+ ,364839
+ ,1889
+ ,43
+ ,37527
+ ,358649
+ ,3055
+ ,30
+ ,60247
+ ,376641
+ ,272
+ ,9
+ ,112995
+ ,429112
+ ,1414
+ ,22
+ ,70184
+ ,330546
+ ,2564
+ ,19
+ ,130140
+ ,403560
+ ,1383
+ ,9
+ ,73221
+ ,317892
+ ,1261
+ ,31
+ ,76114
+ ,307528
+ ,975
+ ,19
+ ,90534
+ ,235133
+ ,3366
+ ,55
+ ,108479
+ ,299243
+ ,576
+ ,8
+ ,113761
+ ,314073
+ ,1686
+ ,28
+ ,68696
+ ,368186
+ ,746
+ ,29
+ ,71561
+ ,269661
+ ,3192
+ ,48
+ ,59831
+ ,125390
+ ,2045
+ ,16
+ ,97890
+ ,510834
+ ,5702
+ ,47
+ ,101481
+ ,321896
+ ,1932
+ ,20
+ ,72954
+ ,249898
+ ,936
+ ,22
+ ,67939
+ ,408881
+ ,3437
+ ,33
+ ,48022
+ ,158492
+ ,5131
+ ,44
+ ,86111
+ ,292154
+ ,2397
+ ,13
+ ,74020
+ ,289513
+ ,1389
+ ,6
+ ,57530
+ ,378049
+ ,1503
+ ,35
+ ,56364
+ ,343466
+ ,402
+ ,8
+ ,84990
+ ,332743
+ ,2239
+ ,17
+ ,88590
+ ,442882
+ ,2234
+ ,11
+ ,77200
+ ,214215
+ ,837
+ ,21
+ ,61262
+ ,315688
+ ,10579
+ ,92
+ ,110309
+ ,375195
+ ,875
+ ,12
+ ,67000
+ ,334280
+ ,1585
+ ,112
+ ,93099
+ ,355864
+ ,1659
+ ,25
+ ,107577
+ ,480382
+ ,2647
+ ,17
+ ,62920
+ ,353058
+ ,3294
+ ,23
+ ,75832
+ ,217193
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,94
+ ,10
+ ,60793
+ ,314533
+ ,422
+ ,23
+ ,57935
+ ,318056
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,34
+ ,7
+ ,60630
+ ,314353
+ ,1558
+ ,25
+ ,55637
+ ,369448
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,43
+ ,20
+ ,60887
+ ,312846
+ ,645
+ ,4
+ ,60720
+ ,312075
+ ,316
+ ,4
+ ,60505
+ ,315009
+ ,115
+ ,10
+ ,60945
+ ,318903
+ ,5
+ ,1
+ ,60720
+ ,314887
+ ,897
+ ,4
+ ,60720
+ ,314913
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,389
+ ,8
+ ,58990
+ ,325506
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1002
+ ,11
+ ,56750
+ ,298568
+ ,36
+ ,4
+ ,60894
+ ,315834
+ ,460
+ ,15
+ ,63346
+ ,329784
+ ,309
+ ,9
+ ,56535
+ ,312878
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,9
+ ,7
+ ,60835
+ ,314987
+ ,271
+ ,2
+ ,60720
+ ,325249
+ ,14
+ ,0
+ ,61016
+ ,315877
+ ,520
+ ,7
+ ,58650
+ ,291650
+ ,1766
+ ,46
+ ,60438
+ ,305959
+ ,0
+ ,5
+ ,60720
+ ,315380
+ ,458
+ ,7
+ ,58625
+ ,297765
+ ,20
+ ,2
+ ,60938
+ ,315245
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,98
+ ,2
+ ,61490
+ ,315236
+ ,405
+ ,5
+ ,60845
+ ,336425
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,483
+ ,7
+ ,60830
+ ,306268
+ ,454
+ ,24
+ ,63261
+ ,302187
+ ,47
+ ,1
+ ,60720
+ ,314882
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,757
+ ,18
+ ,45689
+ ,382712
+ ,4655
+ ,55
+ ,60720
+ ,341570
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,36
+ ,3
+ ,61564
+ ,312412
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,203
+ ,9
+ ,61938
+ ,309596
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,126
+ ,8
+ ,60951
+ ,315547
+ ,400
+ ,113
+ ,60720
+ ,313267
+ ,71
+ ,0
+ ,60745
+ ,316176
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,972
+ ,19
+ ,71642
+ ,359335
+ ,531
+ ,11
+ ,71641
+ ,330068
+ ,2461
+ ,25
+ ,55792
+ ,314289
+ ,378
+ ,16
+ ,71873
+ ,297413
+ ,23
+ ,5
+ ,62555
+ ,314806
+ ,638
+ ,11
+ ,60370
+ ,333210
+ ,2300
+ ,23
+ ,64873
+ ,352108
+ ,149
+ ,6
+ ,62041
+ ,313332
+ ,226
+ ,5
+ ,65745
+ ,291787
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,275
+ ,7
+ ,59500
+ ,318745
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,141
+ ,7
+ ,61630
+ ,315366
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,28
+ ,3
+ ,60890
+ ,315688
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,4980
+ ,89
+ ,113521
+ ,409642
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,472
+ ,19
+ ,80045
+ ,269587
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,203
+ ,12
+ ,50804
+ ,300962
+ ,496
+ ,12
+ ,87390
+ ,325479
+ ,10
+ ,5
+ ,61656
+ ,316155
+ ,63
+ ,2
+ ,65688
+ ,318574
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1136
+ ,26
+ ,48522
+ ,343613
+ ,265
+ ,3
+ ,60720
+ ,306948
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,267
+ ,11
+ ,57640
+ ,330059
+ ,474
+ ,10
+ ,61977
+ ,288985
+ ,534
+ ,5
+ ,62620
+ ,304485
+ ,0
+ ,2
+ ,60720
+ ,315380
+ ,15
+ ,6
+ ,60831
+ ,315688
+ ,397
+ ,7
+ ,60646
+ ,317736
+ ,0
+ ,2
+ ,60720
+ ,315380
+ ,1866
+ ,28
+ ,56225
+ ,322331
+ ,288
+ ,3
+ ,60510
+ ,296656
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,3
+ ,1
+ ,60698
+ ,315354
+ ,468
+ ,20
+ ,60720
+ ,312161
+ ,20
+ ,1
+ ,60805
+ ,315576
+ ,278
+ ,22
+ ,61404
+ ,314922
+ ,61
+ ,9
+ ,60720
+ ,314551
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,192
+ ,2
+ ,65276
+ ,312339
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,317
+ ,7
+ ,63915
+ ,298700
+ ,738
+ ,9
+ ,60720
+ ,321376
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,368
+ ,13
+ ,61686
+ ,303230
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,2
+ ,0
+ ,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
+ ,7
+ ,59818
+ ,315637
+ ,2332
+ ,2
+ ,72680
+ ,324385
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,131
+ ,15
+ ,61808
+ ,308989
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,206
+ ,9
+ ,53110
+ ,296702
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,167
+ ,1
+ ,64245
+ ,307322
+ ,622
+ ,38
+ ,73007
+ ,304376
+ ,2328
+ ,57
+ ,82732
+ ,253588
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,365
+ ,7
+ ,54820
+ ,309560
+ ,364
+ ,26
+ ,47705
+ ,298466
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,226
+ ,13
+ ,72835
+ ,343929
+ ,307
+ ,10
+ ,58856
+ ,331955
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,188
+ ,9
+ ,77655
+ ,381180
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,138
+ ,26
+ ,69817
+ ,331420
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,125
+ ,19
+ ,60798
+ ,310201
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,282
+ ,12
+ ,62452
+ ,320016
+ ,335
+ ,23
+ ,64175
+ ,320398
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1324
+ ,29
+ ,67440
+ ,291841
+ ,176
+ ,8
+ ,68136
+ ,310670
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,249
+ ,26
+ ,56726
+ ,313491
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,333
+ ,9
+ ,70811
+ ,331323
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,601
+ ,5
+ ,60720
+ ,319210
+ ,30
+ ,3
+ ,62045
+ ,318098
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,249
+ ,13
+ ,54323
+ ,292754
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,165
+ ,12
+ ,62841
+ ,325176
+ ,453
+ ,19
+ ,81125
+ ,365959
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,53
+ ,10
+ ,59506
+ ,302409
+ ,382
+ ,9
+ ,59365
+ ,340968
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,9
+ ,60720
+ ,315380
+ ,30
+ ,4
+ ,60798
+ ,313164
+ ,290
+ ,1
+ ,58790
+ ,301164
+ ,0
+ ,1
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,366
+ ,14
+ ,61808
+ ,344425
+ ,2
+ ,12
+ ,60735
+ ,315394
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,209
+ ,19
+ ,64016
+ ,316647
+ ,384
+ ,17
+ ,54683
+ ,309836
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,365
+ ,32
+ ,87192
+ ,346611
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,49
+ ,14
+ ,64107
+ ,322031
+ ,3
+ ,8
+ ,60761
+ ,315656
+ ,133
+ ,4
+ ,65990
+ ,339445
+ ,32
+ ,0
+ ,59988
+ ,314964
+ ,368
+ ,20
+ ,61167
+ ,297141
+ ,1
+ ,5
+ ,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
+ ,1
+ ,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
+ ,4
+ ,60379
+ ,313729
+ ,1
+ ,1
+ ,60727
+ ,315388
+ ,314
+ ,4
+ ,60720
+ ,315371
+ ,844
+ ,20
+ ,60925
+ ,296139
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,26
+ ,1
+ ,60896
+ ,313880
+ ,125
+ ,10
+ ,59734
+ ,317698
+ ,304
+ ,12
+ ,62969
+ ,295580
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,621
+ ,13
+ ,60720
+ ,308256
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,119
+ ,3
+ ,59118
+ ,303677
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1595
+ ,10
+ ,60720
+ ,319369
+ ,312
+ ,3
+ ,58598
+ ,318690
+ ,60
+ ,7
+ ,61124
+ ,314049
+ ,587
+ ,10
+ ,59595
+ ,325699
+ ,135
+ ,1
+ ,62065
+ ,314210
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,514
+ ,15
+ ,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
+ ,28
+ ,61600
+ ,308336
+ ,180
+ ,9
+ ,59635
+ ,316386
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,218
+ ,7
+ ,60720
+ ,315553
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,448
+ ,7
+ ,59781
+ ,323361
+ ,227
+ ,7
+ ,76644
+ ,336639
+ ,174
+ ,3
+ ,64820
+ ,307424
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,121
+ ,11
+ ,56178
+ ,295370
+ ,607
+ ,7
+ ,60436
+ ,322340
+ ,2212
+ ,10
+ ,60720
+ ,319864
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,530
+ ,18
+ ,73433
+ ,317291
+ ,571
+ ,14
+ ,41477
+ ,280398
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,78
+ ,12
+ ,62700
+ ,317330
+ ,2489
+ ,29
+ ,67804
+ ,238125
+ ,131
+ ,3
+ ,59661
+ ,327071
+ ,923
+ ,6
+ ,58620
+ ,309038
+ ,72
+ ,3
+ ,60398
+ ,314210
+ ,572
+ ,8
+ ,58580
+ ,307930
+ ,397
+ ,10
+ ,62710
+ ,322327
+ ,450
+ ,6
+ ,59325
+ ,292136
+ ,622
+ ,8
+ ,60950
+ ,263276
+ ,694
+ ,6
+ ,68060
+ ,367655
+ ,3425
+ ,9
+ ,83620
+ ,283910
+ ,562
+ ,8
+ ,58456
+ ,283587
+ ,4917
+ ,26
+ ,52811
+ ,243650
+ ,1442
+ ,239
+ ,121173
+ ,438493
+ ,529
+ ,7
+ ,63870
+ ,296261
+ ,2126
+ ,41
+ ,21001
+ ,230621
+ ,1061
+ ,3
+ ,70415
+ ,304252
+ ,776
+ ,8
+ ,64230
+ ,333505
+ ,611
+ ,6
+ ,59190
+ ,296919
+ ,1526
+ ,21
+ ,69351
+ ,278990
+ ,592
+ ,7
+ ,64270
+ ,276898
+ ,1182
+ ,11
+ ,70694
+ ,327007
+ ,621
+ ,11
+ ,68005
+ ,317046
+ ,989
+ ,12
+ ,58930
+ ,304555
+ ,438
+ ,9
+ ,58320
+ ,298096
+ ,726
+ ,3
+ ,69980
+ ,231861
+ ,1303
+ ,57
+ ,69863
+ ,309422
+ ,7419
+ ,21
+ ,63255
+ ,286963
+ ,1164
+ ,15
+ ,57320
+ ,269753
+ ,3310
+ ,32
+ ,75230
+ ,448243
+ ,1920
+ ,11
+ ,79420
+ ,165404
+ ,965
+ ,2
+ ,73490
+ ,204325
+ ,3256
+ ,23
+ ,35250
+ ,407159
+ ,1135
+ ,20
+ ,62285
+ ,290476
+ ,1270
+ ,24
+ ,69206
+ ,275311
+ ,661
+ ,1
+ ,65920
+ ,246541
+ ,1013
+ ,1
+ ,69770
+ ,253468
+ ,2844
+ ,74
+ ,72683
+ ,240897
+ ,11528
+ ,68
+ ,-14545
+ ,-83265
+ ,6526
+ ,20
+ ,55830
+ ,-42143
+ ,2264
+ ,20
+ ,55174
+ ,272713
+ ,5109
+ ,82
+ ,67038
+ ,215362
+ ,3999
+ ,21
+ ,51252
+ ,42754
+ ,35624
+ ,244
+ ,157278
+ ,306275
+ ,9252
+ ,32
+ ,79510
+ ,253537
+ ,15236
+ ,86
+ ,77440
+ ,372631
+ ,18073
+ ,69
+ ,27284
+ ,-7170)
+ ,dim=c(4
+ ,431)
+ ,dimnames=list(c('Kosten'
+ ,'Oders'
+ ,'Divid.'
+ ,'Wealth')
+ ,1:431))
> y <- array(NA,dim=c(4,431),dimnames=list(c('Kosten','Oders','Divid.','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 Kosten Oders Divid.
1 6282154 162556 807 213118
2 4321023 29790 444 81767
3 4111912 87550 412 153198
4 223193 84738 428 -26007
5 1491348 54660 315 126942
6 1629616 42634 168 157214
7 1398893 40949 263 129352
8 1926517 45187 267 234817
9 983660 37704 228 60448
10 1443586 16275 129 47818
11 1073089 25830 104 245546
12 984885 12679 122 48020
13 1405225 18014 393 -1710
14 227132 43556 190 32648
15 929118 24811 280 95350
16 1071292 6575 63 151352
17 638830 7123 102 288170
18 856956 21950 265 114337
19 992426 37597 234 37884
20 444477 17821 277 122844
21 857217 12988 73 82340
22 711969 22330 67 79801
23 702380 13326 103 165548
24 358589 16189 290 116384
25 297978 7146 83 134028
26 585715 15824 56 63838
27 657954 27664 236 74996
28 209458 11920 73 31080
29 786690 8568 34 32168
30 439798 14416 139 49857
31 688779 3369 26 87161
32 574339 11819 70 106113
33 741409 6984 40 80570
34 597793 4519 42 102129
35 644190 2220 12 301670
36 377934 18562 211 102313
37 640273 10327 74 88577
38 697458 5336 80 112477
39 550608 2365 83 191778
40 207393 4069 131 79804
41 301607 8636 203 128294
42 345783 13718 56 96448
43 501749 4525 89 93811
44 379983 6869 88 117520
45 387475 4628 39 69159
46 377305 3689 25 101792
47 370837 4891 49 210568
48 430866 7489 149 136996
49 469107 4901 58 121920
50 194493 2284 41 76403
51 530670 3160 90 108094
52 518365 4150 136 134759
53 491303 7285 97 188873
54 527021 1134 63 146216
55 233773 4658 114 156608
56 405972 2384 77 61348
57 652925 3748 6 50350
58 446211 5371 47 87720
59 341340 1285 51 99489
60 387699 9327 85 87419
61 493408 5565 43 94355
62 146494 1528 32 60326
63 414462 3122 25 94670
64 364304 7561 77 82425
65 355178 2675 54 59017
66 357760 13253 251 90829
67 261216 880 15 80791
68 397144 2053 44 100423
69 374943 1424 73 131116
70 424898 4036 85 100269
71 202055 3045 49 27330
72 378525 5119 38 39039
73 310768 1431 35 106885
74 325738 554 9 79285
75 394510 1975 34 118881
76 247060 1765 20 77623
77 368078 1012 29 114768
78 236761 810 11 74015
79 312378 1280 52 69465
80 339836 666 13 117869
81 347385 1380 29 60982
82 426280 4677 66 90131
83 352850 876 33 138971
84 301881 814 15 39625
85 377516 514 15 102725
86 357312 5692 68 64239
87 458343 3642 100 90262
88 354228 540 13 103960
89 308636 2099 45 106611
90 386212 567 14 103345
91 393343 2001 36 95551
92 378509 2949 40 82903
93 452469 2253 68 63593
94 364839 6533 29 126910
95 358649 1889 43 37527
96 376641 3055 30 60247
97 429112 272 9 112995
98 330546 1414 22 70184
99 403560 2564 19 130140
100 317892 1383 9 73221
101 307528 1261 31 76114
102 235133 975 19 90534
103 299243 3366 55 108479
104 314073 576 8 113761
105 368186 1686 28 68696
106 269661 746 29 71561
107 125390 3192 48 59831
108 510834 2045 16 97890
109 321896 5702 47 101481
110 249898 1932 20 72954
111 408881 936 22 67939
112 158492 3437 33 48022
113 292154 5131 44 86111
114 289513 2397 13 74020
115 378049 1389 6 57530
116 343466 1503 35 56364
117 332743 402 8 84990
118 442882 2239 17 88590
119 214215 2234 11 77200
120 315688 837 21 61262
121 375195 10579 92 110309
122 334280 875 12 67000
123 355864 1585 112 93099
124 480382 1659 25 107577
125 353058 2647 17 62920
126 217193 3294 23 75832
127 315380 0 0 60720
128 314533 94 10 60793
129 318056 422 23 57935
130 315380 0 0 60720
131 314353 34 7 60630
132 369448 1558 25 55637
133 315380 0 0 60720
134 312846 43 20 60887
135 312075 645 4 60720
136 315009 316 4 60505
137 318903 115 10 60945
138 314887 5 1 60720
139 314913 897 4 60720
140 315380 0 0 60720
141 325506 389 8 58990
142 315380 0 0 60720
143 298568 1002 11 56750
144 315834 36 4 60894
145 329784 460 15 63346
146 312878 309 9 56535
147 315380 0 0 60720
148 314987 9 7 60835
149 325249 271 2 60720
150 315877 14 0 61016
151 291650 520 7 58650
152 305959 1766 46 60438
153 315380 0 5 60720
154 297765 458 7 58625
155 315245 20 2 60938
156 315380 0 0 60720
157 315380 0 0 60720
158 315236 98 2 61490
159 336425 405 5 60845
160 315380 0 0 60720
161 315380 0 0 60720
162 315380 0 0 60720
163 315380 0 0 60720
164 306268 483 7 60830
165 302187 454 24 63261
166 314882 47 1 60720
167 315380 0 0 60720
168 382712 757 18 45689
169 341570 4655 55 60720
170 315380 0 0 60720
171 315380 0 0 60720
172 312412 36 3 61564
173 315380 0 0 60720
174 309596 203 9 61938
175 315380 0 0 60720
176 315547 126 8 60951
177 313267 400 113 60720
178 316176 71 0 60745
179 315380 0 0 60720
180 315380 0 0 60720
181 359335 972 19 71642
182 330068 531 11 71641
183 314289 2461 25 55792
184 297413 378 16 71873
185 314806 23 5 62555
186 333210 638 11 60370
187 352108 2300 23 64873
188 313332 149 6 62041
189 291787 226 5 65745
190 315380 0 0 60720
191 318745 275 7 59500
192 315380 0 0 60720
193 315366 141 7 61630
194 315380 0 0 60720
195 315688 28 3 60890
196 315380 0 0 60720
197 409642 4980 89 113521
198 315380 0 0 60720
199 315380 0 0 60720
200 269587 472 19 80045
201 315380 0 0 60720
202 315380 0 0 60720
203 315380 0 0 60720
204 300962 203 12 50804
205 325479 496 12 87390
206 316155 10 5 61656
207 318574 63 2 65688
208 315380 0 0 60720
209 343613 1136 26 48522
210 306948 265 3 60720
211 315380 0 0 60720
212 315380 0 0 60720
213 330059 267 11 57640
214 288985 474 10 61977
215 304485 534 5 62620
216 315380 0 2 60720
217 315688 15 6 60831
218 317736 397 7 60646
219 315380 0 2 60720
220 322331 1866 28 56225
221 296656 288 3 60510
222 315380 0 0 60720
223 315354 3 1 60698
224 312161 468 20 60720
225 315576 20 1 60805
226 314922 278 22 61404
227 314551 61 9 60720
228 315380 0 0 60720
229 312339 192 2 65276
230 315380 0 0 60720
231 298700 317 7 63915
232 321376 738 9 60720
233 315380 0 0 60720
234 303230 368 13 61686
235 315380 0 0 60720
236 315487 2 0 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 7 59818
245 324385 2332 2 72680
246 315380 0 0 60720
247 315380 0 0 60720
248 308989 131 15 61808
249 315380 0 0 60720
250 315380 0 0 60720
251 296702 206 9 53110
252 315380 0 0 60720
253 307322 167 1 64245
254 304376 622 38 73007
255 253588 2328 57 82732
256 315380 0 0 60720
257 309560 365 7 54820
258 298466 364 26 47705
259 315380 0 0 60720
260 315380 0 0 60720
261 315380 0 0 60720
262 315380 0 0 60720
263 343929 226 13 72835
264 331955 307 10 58856
265 315380 0 0 60720
266 315380 0 0 60720
267 315380 0 0 60720
268 381180 188 9 77655
269 315380 0 0 60720
270 331420 138 26 69817
271 315380 0 0 60720
272 315380 0 0 60720
273 315380 0 0 60720
274 310201 125 19 60798
275 315380 0 0 60720
276 320016 282 12 62452
277 320398 335 23 64175
278 315380 0 0 60720
279 291841 1324 29 67440
280 310670 176 8 68136
281 315380 0 0 60720
282 315380 0 0 60720
283 313491 249 26 56726
284 315380 0 0 60720
285 331323 333 9 70811
286 315380 0 0 60720
287 319210 601 5 60720
288 318098 30 3 62045
289 315380 0 0 60720
290 292754 249 13 54323
291 315380 0 0 60720
292 325176 165 12 62841
293 365959 453 19 81125
294 315380 0 0 60720
295 302409 53 10 59506
296 340968 382 9 59365
297 315380 0 0 60720
298 315380 0 0 60720
299 315380 0 0 60720
300 315380 0 9 60720
301 313164 30 4 60798
302 301164 290 1 58790
303 315380 0 1 60720
304 315380 0 0 60720
305 344425 366 14 61808
306 315394 2 12 60735
307 315380 0 0 60720
308 316647 209 19 64016
309 309836 384 17 54683
310 315380 0 0 60720
311 315380 0 0 60720
312 346611 365 32 87192
313 315380 0 0 60720
314 322031 49 14 64107
315 315656 3 8 60761
316 339445 133 4 65990
317 314964 32 0 59988
318 297141 368 20 61167
319 315372 1 5 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 1 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 4 60379
335 315388 1 1 60727
336 315371 314 4 60720
337 296139 844 20 60925
338 315380 0 0 60720
339 313880 26 1 60896
340 317698 125 10 59734
341 295580 304 12 62969
342 315380 0 0 60720
343 315380 0 0 60720
344 315380 0 0 60720
345 308256 621 13 60720
346 315380 0 0 60720
347 303677 119 3 59118
348 315380 0 0 60720
349 315380 0 0 60720
350 319369 1595 10 60720
351 318690 312 3 58598
352 314049 60 7 61124
353 325699 587 10 59595
354 314210 135 1 62065
355 315380 0 0 60720
356 315380 0 0 60720
357 322378 514 15 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 28 61600
365 316386 180 9 59635
366 315380 0 0 60720
367 315380 0 0 60720
368 315380 0 0 60720
369 315380 0 0 60720
370 315553 218 7 60720
371 315380 0 0 60720
372 323361 448 7 59781
373 336639 227 7 76644
374 307424 174 3 64820
375 315380 0 0 60720
376 315380 0 0 60720
377 295370 121 11 56178
378 322340 607 7 60436
379 319864 2212 10 60720
380 315380 0 0 60720
381 315380 0 0 60720
382 317291 530 18 73433
383 280398 571 14 41477
384 315380 0 0 60720
385 317330 78 12 62700
386 238125 2489 29 67804
387 327071 131 3 59661
388 309038 923 6 58620
389 314210 72 3 60398
390 307930 572 8 58580
391 322327 397 10 62710
392 292136 450 6 59325
393 263276 622 8 60950
394 367655 694 6 68060
395 283910 3425 9 83620
396 283587 562 8 58456
397 243650 4917 26 52811
398 438493 1442 239 121173
399 296261 529 7 63870
400 230621 2126 41 21001
401 304252 1061 3 70415
402 333505 776 8 64230
403 296919 611 6 59190
404 278990 1526 21 69351
405 276898 592 7 64270
406 327007 1182 11 70694
407 317046 621 11 68005
408 304555 989 12 58930
409 298096 438 9 58320
410 231861 726 3 69980
411 309422 1303 57 69863
412 286963 7419 21 63255
413 269753 1164 15 57320
414 448243 3310 32 75230
415 165404 1920 11 79420
416 204325 965 2 73490
417 407159 3256 23 35250
418 290476 1135 20 62285
419 275311 1270 24 69206
420 246541 661 1 65920
421 253468 1013 1 69770
422 240897 2844 74 72683
423 -83265 11528 68 -14545
424 -42143 6526 20 55830
425 272713 2264 20 55174
426 215362 5109 82 67038
427 42754 3999 21 51252
428 306275 35624 244 157278
429 253537 9252 32 79510
430 372631 15236 86 77440
431 -7170 18073 69 27284
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) Kosten Oders Divid.
1.650e+05 2.202e+01 1.009e+03 1.568e+00
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-2198539 -21927 41819 55153 2923814
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.650e+05 3.063e+04 5.386 1.19e-07 ***
Kosten 2.202e+01 2.155e+00 10.217 < 2e-16 ***
Oders 1.009e+03 3.497e+02 2.886 0.004096 **
Divid. 1.569e+00 4.115e-01 3.812 0.000158 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 253500 on 427 degrees of freedom
Multiple R-squared: 0.6556, Adjusted R-squared: 0.6532
F-statistic: 271 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 2.430592e-32 1.215296e-32
[2,] 1.0000000 8.739682e-53 4.369841e-53
[3,] 1.0000000 2.606180e-53 1.303090e-53
[4,] 1.0000000 3.344876e-75 1.672438e-75
[5,] 1.0000000 3.789859e-78 1.894929e-78
[6,] 1.0000000 8.074852e-84 4.037426e-84
[7,] 1.0000000 3.228560e-105 1.614280e-105
[8,] 1.0000000 2.615482e-107 1.307741e-107
[9,] 1.0000000 8.428079e-116 4.214039e-116
[10,] 1.0000000 2.058138e-123 1.029069e-123
[11,] 1.0000000 1.769463e-133 8.847316e-134
[12,] 1.0000000 1.349026e-139 6.745128e-140
[13,] 1.0000000 1.674390e-144 8.371948e-145
[14,] 1.0000000 9.979506e-153 4.989753e-153
[15,] 1.0000000 4.714147e-161 2.357073e-161
[16,] 1.0000000 2.522661e-164 1.261330e-164
[17,] 1.0000000 7.254403e-165 3.627202e-165
[18,] 1.0000000 4.523952e-171 2.261976e-171
[19,] 1.0000000 1.148221e-170 5.741105e-171
[20,] 1.0000000 7.161661e-173 3.580830e-173
[21,] 1.0000000 6.188206e-176 3.094103e-176
[22,] 1.0000000 3.089147e-175 1.544573e-175
[23,] 1.0000000 4.706165e-189 2.353083e-189
[24,] 1.0000000 4.328376e-189 2.164188e-189
[25,] 1.0000000 6.442047e-195 3.221023e-195
[26,] 1.0000000 2.904653e-196 1.452326e-196
[27,] 1.0000000 1.152967e-205 5.764835e-206
[28,] 1.0000000 4.251535e-208 2.125767e-208
[29,] 1.0000000 2.166124e-207 1.083062e-207
[30,] 1.0000000 1.983912e-208 9.919558e-209
[31,] 1.0000000 4.716584e-214 2.358292e-214
[32,] 1.0000000 6.244748e-220 3.122374e-220
[33,] 1.0000000 4.578503e-219 2.289252e-219
[34,] 1.0000000 1.143927e-219 5.719637e-220
[35,] 1.0000000 3.044899e-221 1.522449e-221
[36,] 1.0000000 1.875420e-220 9.377099e-221
[37,] 1.0000000 8.053626e-221 4.026813e-221
[38,] 1.0000000 9.407400e-220 4.703700e-220
[39,] 1.0000000 2.088861e-219 1.044430e-219
[40,] 1.0000000 1.808484e-218 9.042421e-219
[41,] 1.0000000 4.585386e-219 2.292693e-219
[42,] 1.0000000 2.174066e-218 1.087033e-218
[43,] 1.0000000 8.208906e-218 4.104453e-218
[44,] 1.0000000 4.343402e-218 2.171701e-218
[45,] 1.0000000 1.574703e-218 7.873517e-219
[46,] 1.0000000 4.484256e-218 2.242128e-218
[47,] 1.0000000 2.730826e-217 1.365413e-217
[48,] 1.0000000 8.423766e-217 4.211883e-217
[49,] 1.0000000 5.051767e-220 2.525883e-220
[50,] 1.0000000 1.098346e-219 5.491728e-220
[51,] 1.0000000 3.518083e-233 1.759042e-233
[52,] 1.0000000 1.416428e-233 7.082138e-234
[53,] 1.0000000 1.301640e-232 6.508202e-233
[54,] 1.0000000 3.366758e-232 1.683379e-232
[55,] 1.0000000 6.717500e-234 3.358750e-234
[56,] 1.0000000 1.972486e-235 9.862431e-236
[57,] 1.0000000 5.452362e-235 2.726181e-235
[58,] 1.0000000 3.004395e-234 1.502198e-234
[59,] 1.0000000 1.671496e-233 8.357482e-234
[60,] 1.0000000 8.480176e-234 4.240088e-234
[61,] 1.0000000 2.273922e-233 1.136961e-233
[62,] 1.0000000 2.142692e-232 1.071346e-232
[63,] 1.0000000 1.214715e-231 6.073577e-232
[64,] 1.0000000 6.707206e-231 3.353603e-231
[65,] 1.0000000 4.723453e-230 2.361727e-230
[66,] 1.0000000 3.817530e-231 1.908765e-231
[67,] 1.0000000 1.310705e-230 6.553526e-231
[68,] 1.0000000 1.124158e-229 5.620790e-230
[69,] 1.0000000 1.349869e-228 6.749343e-229
[70,] 1.0000000 3.164417e-228 1.582209e-228
[71,] 1.0000000 2.854783e-227 1.427392e-227
[72,] 1.0000000 3.837427e-227 1.918714e-227
[73,] 1.0000000 4.036207e-226 2.018104e-226
[74,] 1.0000000 1.525577e-225 7.627887e-226
[75,] 1.0000000 9.075956e-225 4.537978e-225
[76,] 1.0000000 1.497523e-224 7.487613e-225
[77,] 1.0000000 1.459550e-224 7.297752e-225
[78,] 1.0000000 8.179804e-224 4.089902e-224
[79,] 1.0000000 7.945820e-223 3.972910e-223
[80,] 1.0000000 2.123714e-222 1.061857e-222
[81,] 1.0000000 1.694008e-222 8.470042e-223
[82,] 1.0000000 1.267032e-221 6.335159e-222
[83,] 1.0000000 2.738350e-221 1.369175e-221
[84,] 1.0000000 2.664020e-220 1.332010e-220
[85,] 1.0000000 2.421778e-219 1.210889e-219
[86,] 1.0000000 1.577767e-218 7.888834e-219
[87,] 1.0000000 5.308984e-220 2.654492e-220
[88,] 1.0000000 5.402230e-219 2.701115e-219
[89,] 1.0000000 1.923829e-219 9.619147e-220
[90,] 1.0000000 1.707649e-219 8.538247e-220
[91,] 1.0000000 1.533171e-218 7.665856e-219
[92,] 1.0000000 1.506117e-217 7.530586e-218
[93,] 1.0000000 1.317935e-216 6.589677e-217
[94,] 1.0000000 1.295777e-215 6.478886e-216
[95,] 1.0000000 1.170218e-214 5.851090e-215
[96,] 1.0000000 9.191774e-216 4.595887e-216
[97,] 1.0000000 8.432706e-216 4.216353e-216
[98,] 1.0000000 1.403658e-216 7.018291e-217
[99,] 1.0000000 6.191573e-216 3.095787e-216
[100,] 1.0000000 2.187586e-215 1.093793e-215
[101,] 1.0000000 1.346190e-217 6.730948e-218
[102,] 1.0000000 5.649664e-219 2.824832e-219
[103,] 1.0000000 4.299186e-218 2.149593e-218
[104,] 1.0000000 8.620852e-218 4.310426e-218
[105,] 1.0000000 6.134868e-218 3.067434e-218
[106,] 1.0000000 3.811163e-218 1.905582e-218
[107,] 1.0000000 2.946839e-217 1.473420e-217
[108,] 1.0000000 2.387761e-216 1.193881e-216
[109,] 1.0000000 1.584638e-216 7.923190e-217
[110,] 1.0000000 6.365781e-216 3.182891e-216
[111,] 1.0000000 4.585815e-215 2.292908e-215
[112,] 1.0000000 3.686476e-215 1.843238e-215
[113,] 1.0000000 3.978613e-216 1.989306e-216
[114,] 1.0000000 4.088253e-215 2.044126e-215
[115,] 1.0000000 1.975573e-214 9.877866e-215
[116,] 1.0000000 1.832932e-213 9.164659e-214
[117,] 1.0000000 1.589629e-212 7.948143e-213
[118,] 1.0000000 2.198678e-212 1.099339e-212
[119,] 1.0000000 6.318306e-212 3.159153e-212
[120,] 1.0000000 2.063023e-212 1.031511e-212
[121,] 1.0000000 1.926414e-211 9.632071e-212
[122,] 1.0000000 1.902403e-210 9.512017e-211
[123,] 1.0000000 1.818865e-209 9.094323e-210
[124,] 1.0000000 1.710070e-208 8.550349e-209
[125,] 1.0000000 1.668946e-207 8.344731e-208
[126,] 1.0000000 1.435411e-207 7.177056e-208
[127,] 1.0000000 1.367575e-206 6.837874e-207
[128,] 1.0000000 1.417370e-205 7.086848e-206
[129,] 1.0000000 1.401973e-204 7.009863e-205
[130,] 1.0000000 1.358313e-203 6.791565e-204
[131,] 1.0000000 1.322004e-202 6.610022e-203
[132,] 1.0000000 1.269751e-201 6.348753e-202
[133,] 1.0000000 1.214274e-200 6.071372e-201
[134,] 1.0000000 1.158771e-199 5.793855e-200
[135,] 1.0000000 9.589873e-199 4.794937e-199
[136,] 1.0000000 9.151115e-198 4.575558e-198
[137,] 1.0000000 9.369098e-197 4.684549e-197
[138,] 1.0000000 9.036756e-196 4.518378e-196
[139,] 1.0000000 8.271089e-195 4.135544e-195
[140,] 1.0000000 7.584350e-194 3.792175e-194
[141,] 1.0000000 7.191716e-193 3.595858e-193
[142,] 1.0000000 6.975854e-192 3.487927e-192
[143,] 1.0000000 6.006059e-191 3.003029e-191
[144,] 1.0000000 5.677643e-190 2.838821e-190
[145,] 1.0000000 5.527930e-189 2.763965e-189
[146,] 1.0000000 5.706216e-188 2.853108e-188
[147,] 1.0000000 5.419894e-187 2.709947e-187
[148,] 1.0000000 5.317814e-186 2.658907e-186
[149,] 1.0000000 4.999803e-185 2.499902e-185
[150,] 1.0000000 4.652301e-184 2.326150e-184
[151,] 1.0000000 4.320326e-183 2.160163e-183
[152,] 1.0000000 4.057165e-182 2.028583e-182
[153,] 1.0000000 2.837675e-181 1.418837e-181
[154,] 1.0000000 2.619488e-180 1.309744e-180
[155,] 1.0000000 2.412398e-179 1.206199e-179
[156,] 1.0000000 2.216283e-178 1.108142e-178
[157,] 1.0000000 2.031012e-177 1.015506e-177
[158,] 1.0000000 1.939728e-176 9.698638e-177
[159,] 1.0000000 1.809404e-175 9.047021e-176
[160,] 1.0000000 1.649790e-174 8.248948e-175
[161,] 1.0000000 1.493476e-173 7.467381e-174
[162,] 1.0000000 1.850795e-174 9.253976e-175
[163,] 1.0000000 3.303986e-174 1.651993e-174
[164,] 1.0000000 3.027321e-173 1.513660e-173
[165,] 1.0000000 2.765752e-172 1.382876e-172
[166,] 1.0000000 2.574255e-171 1.287127e-171
[167,] 1.0000000 2.337443e-170 1.168721e-170
[168,] 1.0000000 2.194040e-169 1.097020e-169
[169,] 1.0000000 1.979190e-168 9.895948e-169
[170,] 1.0000000 1.805319e-167 9.026594e-168
[171,] 1.0000000 1.695989e-166 8.479946e-167
[172,] 1.0000000 1.505326e-165 7.526629e-166
[173,] 1.0000000 1.338383e-164 6.691913e-165
[174,] 1.0000000 1.186093e-163 5.930463e-164
[175,] 1.0000000 7.911250e-163 3.955625e-163
[176,] 1.0000000 7.250046e-162 3.625023e-162
[177,] 1.0000000 4.225910e-161 2.112955e-161
[178,] 1.0000000 2.568503e-160 1.284252e-160
[179,] 1.0000000 2.293050e-159 1.146525e-159
[180,] 1.0000000 1.530856e-158 7.654282e-159
[181,] 1.0000000 5.726082e-158 2.863041e-158
[182,] 1.0000000 5.092835e-157 2.546418e-157
[183,] 1.0000000 3.600176e-156 1.800088e-156
[184,] 1.0000000 3.101590e-155 1.550795e-155
[185,] 1.0000000 2.539664e-154 1.269832e-154
[186,] 1.0000000 2.172516e-153 1.086258e-153
[187,] 1.0000000 1.883443e-152 9.417215e-153
[188,] 1.0000000 1.599390e-151 7.996952e-152
[189,] 1.0000000 1.358950e-150 6.794751e-151
[190,] 1.0000000 1.145603e-149 5.728014e-150
[191,] 1.0000000 8.787550e-149 4.393775e-149
[192,] 1.0000000 7.351060e-148 3.675530e-148
[193,] 1.0000000 6.126847e-147 3.063423e-147
[194,] 1.0000000 6.344816e-147 3.172408e-147
[195,] 1.0000000 5.301083e-146 2.650541e-146
[196,] 1.0000000 4.412611e-145 2.206306e-145
[197,] 1.0000000 3.659338e-144 1.829669e-144
[198,] 1.0000000 2.808778e-143 1.404389e-143
[199,] 1.0000000 1.235292e-142 6.176459e-143
[200,] 1.0000000 1.027787e-141 5.138934e-142
[201,] 1.0000000 8.550471e-141 4.275235e-141
[202,] 1.0000000 6.981267e-140 3.490634e-140
[203,] 1.0000000 1.005138e-139 5.025692e-140
[204,] 1.0000000 8.403893e-139 4.201947e-139
[205,] 1.0000000 6.831375e-138 3.415688e-138
[206,] 1.0000000 5.531530e-137 2.765765e-137
[207,] 1.0000000 3.352656e-136 1.676328e-136
[208,] 1.0000000 2.520675e-135 1.260337e-135
[209,] 1.0000000 2.070219e-134 1.035110e-134
[210,] 1.0000000 1.655285e-133 8.276427e-134
[211,] 1.0000000 1.322122e-132 6.610610e-133
[212,] 1.0000000 1.021302e-131 5.106509e-132
[213,] 1.0000000 8.066985e-131 4.033492e-131
[214,] 1.0000000 3.866636e-130 1.933318e-130
[215,] 1.0000000 3.046118e-129 1.523059e-129
[216,] 1.0000000 2.377357e-128 1.188679e-128
[217,] 1.0000000 1.849042e-127 9.245209e-128
[218,] 1.0000000 1.451266e-126 7.256329e-127
[219,] 1.0000000 1.119107e-125 5.595536e-126
[220,] 1.0000000 8.723133e-125 4.361566e-125
[221,] 1.0000000 6.723391e-124 3.361695e-124
[222,] 1.0000000 5.117062e-123 2.558531e-123
[223,] 1.0000000 3.917908e-122 1.958954e-122
[224,] 1.0000000 2.957740e-121 1.478870e-121
[225,] 1.0000000 2.147152e-120 1.073576e-120
[226,] 1.0000000 1.481676e-119 7.408380e-120
[227,] 1.0000000 1.105263e-118 5.526315e-119
[228,] 1.0000000 8.419570e-118 4.209785e-118
[229,] 1.0000000 6.227574e-117 3.113787e-117
[230,] 1.0000000 4.585707e-116 2.292853e-116
[231,] 1.0000000 3.363278e-115 1.681639e-115
[232,] 1.0000000 2.445823e-114 1.222912e-114
[233,] 1.0000000 1.778540e-113 8.892699e-114
[234,] 1.0000000 1.287785e-112 6.438925e-113
[235,] 1.0000000 9.284482e-112 4.642241e-112
[236,] 1.0000000 6.759561e-111 3.379780e-111
[237,] 1.0000000 4.831331e-110 2.415665e-110
[238,] 1.0000000 3.405378e-109 1.702689e-109
[239,] 1.0000000 2.471831e-108 1.235915e-108
[240,] 1.0000000 1.743174e-107 8.715870e-108
[241,] 1.0000000 1.223942e-106 6.119712e-107
[242,] 1.0000000 8.778635e-106 4.389318e-106
[243,] 1.0000000 6.109333e-105 3.054667e-105
[244,] 1.0000000 4.232923e-104 2.116461e-104
[245,] 1.0000000 2.900174e-103 1.450087e-103
[246,] 1.0000000 1.991501e-102 9.957507e-103
[247,] 1.0000000 1.367456e-101 6.837281e-102
[248,] 1.0000000 7.697318e-101 3.848659e-101
[249,] 1.0000000 6.133853e-101 3.066926e-101
[250,] 1.0000000 4.198894e-100 2.099447e-100
[251,] 1.0000000 2.626027e-99 1.313014e-99
[252,] 1.0000000 1.494175e-98 7.470875e-99
[253,] 1.0000000 1.010003e-97 5.050017e-98
[254,] 1.0000000 6.796205e-97 3.398102e-97
[255,] 1.0000000 4.552238e-96 2.276119e-96
[256,] 1.0000000 3.035225e-95 1.517612e-95
[257,] 1.0000000 2.040084e-94 1.020042e-94
[258,] 1.0000000 1.070511e-93 5.352554e-94
[259,] 1.0000000 7.045406e-93 3.522703e-93
[260,] 1.0000000 4.615288e-92 2.307644e-92
[261,] 1.0000000 3.009272e-91 1.504636e-91
[262,] 1.0000000 1.604508e-90 8.022538e-91
[263,] 1.0000000 1.037089e-89 5.185444e-90
[264,] 1.0000000 6.846043e-89 3.423022e-89
[265,] 1.0000000 4.383118e-88 2.191559e-88
[266,] 1.0000000 2.792913e-87 1.396457e-87
[267,] 1.0000000 1.771145e-86 8.855723e-87
[268,] 1.0000000 1.142786e-85 5.713931e-86
[269,] 1.0000000 7.177221e-85 3.588611e-85
[270,] 1.0000000 4.469982e-84 2.234991e-84
[271,] 1.0000000 2.824435e-83 1.412218e-83
[272,] 1.0000000 1.747989e-82 8.739945e-83
[273,] 1.0000000 9.987421e-82 4.993710e-82
[274,] 1.0000000 5.987048e-81 2.993524e-81
[275,] 1.0000000 3.657469e-80 1.828735e-80
[276,] 1.0000000 2.223306e-79 1.111653e-79
[277,] 1.0000000 1.274009e-78 6.370044e-79
[278,] 1.0000000 7.669229e-78 3.834615e-78
[279,] 1.0000000 4.693385e-77 2.346693e-77
[280,] 1.0000000 2.797238e-76 1.398619e-76
[281,] 1.0000000 1.582610e-75 7.913049e-76
[282,] 1.0000000 9.371309e-75 4.685654e-75
[283,] 1.0000000 5.501093e-74 2.750546e-74
[284,] 1.0000000 3.270733e-73 1.635366e-73
[285,] 1.0000000 1.900239e-72 9.501195e-73
[286,] 1.0000000 1.074963e-71 5.374813e-72
[287,] 1.0000000 6.230254e-71 3.115127e-71
[288,] 1.0000000 3.563613e-70 1.781807e-70
[289,] 1.0000000 2.076082e-69 1.038041e-69
[290,] 1.0000000 8.466843e-69 4.233422e-69
[291,] 1.0000000 4.775010e-68 2.387505e-68
[292,] 1.0000000 2.678632e-67 1.339316e-67
[293,] 1.0000000 1.494603e-66 7.473014e-67
[294,] 1.0000000 8.325508e-66 4.162754e-66
[295,] 1.0000000 4.640352e-65 2.320176e-65
[296,] 1.0000000 2.602008e-64 1.301004e-64
[297,] 1.0000000 1.421324e-63 7.106620e-64
[298,] 1.0000000 7.717969e-63 3.858985e-63
[299,] 1.0000000 3.200789e-62 1.600394e-62
[300,] 1.0000000 1.728956e-61 8.644781e-62
[301,] 1.0000000 9.245426e-61 4.622713e-61
[302,] 1.0000000 5.032402e-60 2.516201e-60
[303,] 1.0000000 2.446646e-59 1.223323e-59
[304,] 1.0000000 1.286819e-58 6.434096e-59
[305,] 1.0000000 6.729447e-58 3.364724e-58
[306,] 1.0000000 3.171111e-57 1.585556e-57
[307,] 1.0000000 1.642811e-56 8.214056e-57
[308,] 1.0000000 8.563712e-56 4.281856e-56
[309,] 1.0000000 4.393795e-55 2.196897e-55
[310,] 1.0000000 2.102368e-54 1.051184e-54
[311,] 1.0000000 1.055659e-53 5.278297e-54
[312,] 1.0000000 5.449743e-53 2.724871e-53
[313,] 1.0000000 2.731083e-52 1.365541e-52
[314,] 1.0000000 1.357996e-51 6.789981e-52
[315,] 1.0000000 6.711486e-51 3.355743e-51
[316,] 1.0000000 3.296681e-50 1.648340e-50
[317,] 1.0000000 1.609365e-49 8.046823e-50
[318,] 1.0000000 7.807862e-49 3.903931e-49
[319,] 1.0000000 3.764353e-48 1.882176e-48
[320,] 1.0000000 1.822163e-47 9.110813e-48
[321,] 1.0000000 8.674707e-47 4.337353e-47
[322,] 1.0000000 4.103382e-46 2.051691e-46
[323,] 1.0000000 1.928533e-45 9.642665e-46
[324,] 1.0000000 9.005099e-45 4.502550e-45
[325,] 1.0000000 4.177378e-44 2.088689e-44
[326,] 1.0000000 1.925085e-43 9.625426e-44
[327,] 1.0000000 8.812592e-43 4.406296e-43
[328,] 1.0000000 4.013298e-42 2.006649e-42
[329,] 1.0000000 1.813349e-41 9.066746e-42
[330,] 1.0000000 8.058085e-41 4.029042e-41
[331,] 1.0000000 3.710338e-40 1.855169e-40
[332,] 1.0000000 1.641984e-39 8.209921e-40
[333,] 1.0000000 7.267902e-39 3.633951e-39
[334,] 1.0000000 3.095162e-38 1.547581e-38
[335,] 1.0000000 1.374656e-37 6.873279e-38
[336,] 1.0000000 5.919262e-37 2.959631e-37
[337,] 1.0000000 2.530407e-36 1.265204e-36
[338,] 1.0000000 1.073822e-35 5.369108e-36
[339,] 1.0000000 4.577995e-35 2.288997e-35
[340,] 1.0000000 1.914210e-34 9.571050e-35
[341,] 1.0000000 8.121518e-34 4.060759e-34
[342,] 1.0000000 3.345165e-33 1.672583e-33
[343,] 1.0000000 1.367273e-32 6.836365e-33
[344,] 1.0000000 5.011547e-32 2.505774e-32
[345,] 1.0000000 1.897232e-31 9.486159e-32
[346,] 1.0000000 7.657861e-31 3.828930e-31
[347,] 1.0000000 2.728176e-30 1.364088e-30
[348,] 1.0000000 1.088639e-29 5.443194e-30
[349,] 1.0000000 4.246226e-29 2.123113e-29
[350,] 1.0000000 1.642496e-28 8.212478e-29
[351,] 1.0000000 6.318022e-28 3.159011e-28
[352,] 1.0000000 2.409193e-27 1.204596e-27
[353,] 1.0000000 9.107831e-27 4.553916e-27
[354,] 1.0000000 3.413206e-26 1.706603e-26
[355,] 1.0000000 1.270191e-25 6.350954e-26
[356,] 1.0000000 4.676202e-25 2.338101e-25
[357,] 1.0000000 1.705948e-24 8.529739e-25
[358,] 1.0000000 6.081439e-24 3.040719e-24
[359,] 1.0000000 2.134975e-23 1.067487e-23
[360,] 1.0000000 7.577712e-23 3.788856e-23
[361,] 1.0000000 2.663818e-22 1.331909e-22
[362,] 1.0000000 9.273145e-22 4.636572e-22
[363,] 1.0000000 3.196226e-21 1.598113e-21
[364,] 1.0000000 1.085630e-20 5.428148e-21
[365,] 1.0000000 3.666904e-20 1.833452e-20
[366,] 1.0000000 1.139622e-19 5.698109e-20
[367,] 1.0000000 3.983320e-19 1.991660e-19
[368,] 1.0000000 1.374135e-18 6.870676e-19
[369,] 1.0000000 4.456475e-18 2.228238e-18
[370,] 1.0000000 1.429197e-17 7.145986e-18
[371,] 1.0000000 4.676550e-17 2.338275e-17
[372,] 1.0000000 1.379755e-16 6.898773e-17
[373,] 1.0000000 3.754474e-16 1.877237e-16
[374,] 1.0000000 1.150483e-15 5.752416e-16
[375,] 1.0000000 3.481903e-15 1.740951e-15
[376,] 1.0000000 1.120447e-14 5.602236e-15
[377,] 1.0000000 2.962514e-14 1.481257e-14
[378,] 1.0000000 8.658496e-14 4.329248e-14
[379,] 1.0000000 2.549114e-13 1.274557e-13
[380,] 1.0000000 6.812795e-13 3.406397e-13
[381,] 1.0000000 1.781479e-12 8.907395e-13
[382,] 1.0000000 4.844002e-12 2.422001e-12
[383,] 1.0000000 1.336663e-11 6.683313e-12
[384,] 1.0000000 3.604832e-11 1.802416e-11
[385,] 1.0000000 9.508511e-11 4.754255e-11
[386,] 1.0000000 2.667857e-10 1.333928e-10
[387,] 1.0000000 7.603123e-10 3.801562e-10
[388,] 1.0000000 1.526677e-09 7.633383e-10
[389,] 1.0000000 4.142153e-09 2.071077e-09
[390,] 1.0000000 1.118990e-08 5.594948e-09
[391,] 1.0000000 2.889582e-08 1.444791e-08
[392,] 1.0000000 5.852297e-08 2.926148e-08
[393,] 0.9999999 1.525094e-07 7.625468e-08
[394,] 0.9999998 3.216808e-07 1.608404e-07
[395,] 0.9999996 8.186697e-07 4.093349e-07
[396,] 0.9999991 1.759357e-06 8.796784e-07
[397,] 0.9999980 4.096960e-06 2.048480e-06
[398,] 0.9999949 1.011551e-05 5.057754e-06
[399,] 0.9999878 2.439833e-05 1.219916e-05
[400,] 0.9999732 5.355622e-05 2.677811e-05
[401,] 0.9999412 1.176231e-04 5.881155e-05
[402,] 0.9998805 2.390808e-04 1.195404e-04
[403,] 0.9997563 4.874979e-04 2.437490e-04
[404,] 0.9994846 1.030896e-03 5.154479e-04
[405,] 0.9989024 2.195194e-03 1.097597e-03
[406,] 0.9982032 3.593522e-03 1.796761e-03
[407,] 0.9964342 7.131593e-03 3.565796e-03
[408,] 0.9968653 6.269431e-03 3.134716e-03
[409,] 0.9949696 1.006073e-02 5.030367e-03
[410,] 0.9905244 1.895117e-02 9.475586e-03
[411,] 0.9977198 4.560433e-03 2.280217e-03
[412,] 0.9952421 9.515868e-03 4.757934e-03
[413,] 0.9890591 2.188182e-02 1.094091e-02
[414,] 0.9754772 4.904553e-02 2.452276e-02
[415,] 0.9482408 1.035184e-01 5.175922e-02
[416,] 0.8942305 2.115390e-01 1.057695e-01
[417,] 0.8019684 3.960632e-01 1.980316e-01
[418,] 0.8245668 3.508663e-01 1.754332e-01
> postscript(file="/var/www/html/freestat/rcomp/tmp/1uzlk1291211442.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/2m8k51291211442.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/3m8k51291211442.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/4m8k51291211442.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/5fzj81291211442.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
1389636.7423 2923814.3860 1363346.2455 -2198539.3196 -394043.3532
6 7 8 9 10
109869.7630 -135939.7551 128928.3808 -336334.3876 715092.5993
11 12 13 14 15
-150657.5105 342307.7668 449675.2072 -1139741.2853 -214255.2259
16 17 18 19 20
460574.8057 -237906.6182 -238071.5057 -295877.9417 -585102.6254
21 22 23 24 25
203463.1027 -137416.1514 -119601.7096 -638049.8630 -318324.0983
26 27 28 29 30
-84297.5193 -471896.4012 -340384.4471 348299.6853 -261059.4968
31 32 33 34 35
286668.0415 -87937.6349 255919.3586 130738.8627 -54941.9881
36 37 38 39 40
-569143.8636 34310.7675 157833.6426 -51015.8179 -304565.6193
41 42 43 44 45
-459620.0553 -329011.9718 172.3323 -209375.2579 -27239.0962
46 47 48 49 50
-53789.6599 -281553.3224 -264256.0037 -53547.0125 -181997.0873
51 52 53 54 55
35733.0747 -86617.6159 -228211.1136 44146.0101 -394458.1521
56 57 58 59 60
14558.8138 320392.3041 -22047.0390 -59458.7651 -205533.6443
61 62 63 64 65
14509.4419 -179052.6824 7020.6734 -174141.9735 -15771.8893
66 67 68 69 70
-494796.9414 -65004.9696 -14962.5874 -100725.8195 -72004.9465
71 72 73 74 75
-122294.6437 1252.3272 -88696.6011 15111.9538 -34736.9165
76 77 78 79 80
-98722.3509 -28470.4246 -73253.8756 -42228.7857 -37810.1071
81 82 83 84 85
27095.5975 -49657.3672 -82702.7231 41679.9422 24950.3027
86 87 88 89 90
-102379.3908 -29329.7313 1171.4484 -115198.1764 32516.3322
91 92 93 94 95
-1902.9294 -21807.0740 69502.3235 -172301.4464 49811.5996
96 97 98 99 100
19619.0604 71821.8316 2140.2336 -41173.8667 -1473.8500
101 102 103 104 105
-35893.9515 -112498.0489 -165507.6930 -50102.0373 30070.0947
106 107 108 109 110
-53263.1917 -252162.8870 131136.6578 -175232.6525 -92237.8474
111 112 113 114 115
94519.8432 -190793.4409 -165268.7848 -57466.9920 86189.5153
116 117 118 119 120
21656.0303 17524.5950 72490.9475 -132145.5451 14988.2296
121 122 123 124 125
-288568.0383 32827.5084 -103083.9599 84905.4238 13946.7576
126 127 128 129 130
-162469.9198 55152.5245 42028.6252 29692.3360 55152.5245
131 132 133 134 135
46453.0982 57660.3793 55152.5245 31224.0439 33610.3153
136 137 138 139 140
44124.6501 45697.8939 53540.1528 30900.3954 55152.5245
141 142 143 144 145
51353.5430 55152.5245 11405.4575 50503.8779 40171.2046
146 147 148 149 150
43328.0424 55152.5245 47315.9545 57036.7210 54877.0451
151 152 153 154 155
16156.0728 -39133.1674 50106.0545 23675.2486 52216.7029
156 157 158 159 160
55152.5245 55152.5245 49624.7042 62038.6981 55152.5245
161 162 163 164 165
55152.5245 55152.5245 55152.5245 28169.4235 3755.9701
166 167 168 169 170
52610.4994 55152.5245 111226.8553 -76651.0561 55152.5245
171 172 173 174 175
55152.5245 47040.3092 55152.5245 33905.3464 55152.5245
176 177 178 179 180
44108.9001 -69816.9214 54346.2088 55152.5245 55152.5245
181 182 183 184 185
41401.1884 29918.9688 -17621.7771 -4778.0021 46147.5877
186 187 188 189 190
48383.3242 11517.2049 41696.5179 13656.0689 55152.5245
191 192 193 194 195
47311.6999 55152.5245 43541.9820 55152.5245 51549.5701
196 197 198 199 200
55152.5245 -132865.9295 55152.5245 55152.5245 -50518.7341
201 202 203 204 205
55152.5245 55152.5245 55152.5245 39706.6063 389.6725
206 207 208 209 210
49192.8281 47148.8882 55152.5245 51266.1839 37858.5204
211 212 213 214 215
55152.5245 55152.5245 57681.9685 6257.6673 24474.6896
216 217 218 219 220
53133.9365 48900.4290 41819.3564 53133.9365 -187.5434
221 222 223 224 225
27389.5375 55152.5245 54085.6898 21444.3644 53765.6010
226 227 228 229 230
25296.9108 43896.9295 55152.5245 38720.0837 55152.5245
231 232 233 234 235
19417.3320 35817.3986 55152.5245 20264.8536 55152.5245
236 237 238 239 240
55179.4190 55152.5245 48924.8316 55152.5245 55152.5245
241 242 243 244 245
55152.5245 46558.3712 55152.5245 49230.8366 -7960.0214
246 247 248 249 250
55152.5245 55152.5245 29031.6010 55152.5245 55152.5245
251 252 253 254 255
34791.5919 55152.5245 36879.8403 -27169.8902 -149946.2072
256 257 258 259 260
55152.5245 43485.6542 24396.6177 55152.5245 55152.5245
261 262 263 264 265
55152.5245 55152.5245 46603.3938 57799.4030 55152.5245
266 267 268 269 270
55152.5245 55152.5245 81168.2234 55152.5245 27644.5278
271 272 273 274 275
55152.5245 55152.5245 55152.5245 27922.6546 55152.5245
276 277 278 279 280
38752.0513 24162.5503 55152.5245 -37344.5929 26861.7964
281 282 283 284 285
55152.5245 55152.5245 27804.4032 55152.5245 38853.4518
286 287 288 289 290
55152.5245 40704.7058 52103.9773 55152.5245 23957.2150
291 292 293 294 295
55152.5245 45877.7435 44577.6335 55152.5245 32825.8604
296 297 298 299 300
65372.1884 55152.5245 55152.5245 55152.5245 46068.8784
301 302 303 304 305
48116.5426 36569.8314 54143.2305 55152.5245 60303.2396
306 307 308 309 310
42987.4385 55152.5245 27472.0702 33465.2963 55152.5245
311 312 313 314 315
55152.5245 4530.3833 55152.5245 41282.3000 47223.8192
316 317 318 319 320
63986.5389 55180.1335 7924.8219 50077.6074 55152.5245
321 322 323 324 325
55152.5245 55152.5245 55152.5245 55152.5245 55152.5245
326 327 328 329 330
50777.7514 55152.5245 55152.5245 55152.5245 55152.5245
331 332 333 334 335
55152.5245 55152.5245 55152.5245 47885.6972 54118.2358
336 337 338 339 340
44193.4641 -3177.0170 55152.5245 51794.7786 46172.1334
341 342 343 344 345
13020.8195 55152.5245 55152.5245 55152.5245 21236.0425
346 347 348 349 350
55152.5245 40314.4511 55152.5245 55152.5245 13933.7738
351 352 353 354 355
51894.0439 44801.8786 44219.9617 47891.5630 55152.5245
356 357 358 359 360
55152.5245 7368.8801 55152.5245 55152.5245 55152.5245
361 362 363 364 365
51108.1960 55152.5245 55152.5245 -20345.3694 44813.8486
366 367 368 369 370
55152.5245 55152.5245 55152.5245 55152.5245 43461.0754
371 372 373 374 375
55152.5245 47678.2739 39372.9095 33907.2837 55152.5245
376 377 378 379 380
55152.5245 28500.3161 42129.4647 845.1762 55152.5245
381 382 383 384 385
55152.5245 7288.2607 23651.2435 55152.5245 40168.2486
386 387 388 389 390
-117279.6312 62592.5952 25728.1516 49874.5640 30391.7616
391 392 393 394 395
40145.1899 18133.7508 -19080.2469 74580.5292 -96721.9437
396 397 398 399 400
6463.4051 -138664.7201 -189519.7904 12381.6147 -55495.2651
401 402 403 404 405
2431.9986 42613.8359 19583.9874 -49565.6943 -8995.7460
406 407 408 409 410
14011.1662 20618.4595 13250.1409 22906.3494 -61901.5139
411 412 413 414 415
-51361.8573 -161769.0794 -25907.2560 60088.3931 -177525.6236
416 417 418 419 420
-99195.2033 91983.5806 -17379.6380 -50409.1687 -37404.0006
421 422 423 424 425
-44265.0154 -175393.8602 -547870.2678 -458560.1402 -48844.9394
426 427 428 429 430
-250014.5348 -311858.7540 -1135948.9209 -272147.0098 -336049.2741
431
-682483.1396
> postscript(file="/var/www/html/freestat/rcomp/tmp/6fzj81291211442.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 1389636.7423 NA
1 2923814.3860 1389636.7423
2 1363346.2455 2923814.3860
3 -2198539.3196 1363346.2455
4 -394043.3532 -2198539.3196
5 109869.7630 -394043.3532
6 -135939.7551 109869.7630
7 128928.3808 -135939.7551
8 -336334.3876 128928.3808
9 715092.5993 -336334.3876
10 -150657.5105 715092.5993
11 342307.7668 -150657.5105
12 449675.2072 342307.7668
13 -1139741.2853 449675.2072
14 -214255.2259 -1139741.2853
15 460574.8057 -214255.2259
16 -237906.6182 460574.8057
17 -238071.5057 -237906.6182
18 -295877.9417 -238071.5057
19 -585102.6254 -295877.9417
20 203463.1027 -585102.6254
21 -137416.1514 203463.1027
22 -119601.7096 -137416.1514
23 -638049.8630 -119601.7096
24 -318324.0983 -638049.8630
25 -84297.5193 -318324.0983
26 -471896.4012 -84297.5193
27 -340384.4471 -471896.4012
28 348299.6853 -340384.4471
29 -261059.4968 348299.6853
30 286668.0415 -261059.4968
31 -87937.6349 286668.0415
32 255919.3586 -87937.6349
33 130738.8627 255919.3586
34 -54941.9881 130738.8627
35 -569143.8636 -54941.9881
36 34310.7675 -569143.8636
37 157833.6426 34310.7675
38 -51015.8179 157833.6426
39 -304565.6193 -51015.8179
40 -459620.0553 -304565.6193
41 -329011.9718 -459620.0553
42 172.3323 -329011.9718
43 -209375.2579 172.3323
44 -27239.0962 -209375.2579
45 -53789.6599 -27239.0962
46 -281553.3224 -53789.6599
47 -264256.0037 -281553.3224
48 -53547.0125 -264256.0037
49 -181997.0873 -53547.0125
50 35733.0747 -181997.0873
51 -86617.6159 35733.0747
52 -228211.1136 -86617.6159
53 44146.0101 -228211.1136
54 -394458.1521 44146.0101
55 14558.8138 -394458.1521
56 320392.3041 14558.8138
57 -22047.0390 320392.3041
58 -59458.7651 -22047.0390
59 -205533.6443 -59458.7651
60 14509.4419 -205533.6443
61 -179052.6824 14509.4419
62 7020.6734 -179052.6824
63 -174141.9735 7020.6734
64 -15771.8893 -174141.9735
65 -494796.9414 -15771.8893
66 -65004.9696 -494796.9414
67 -14962.5874 -65004.9696
68 -100725.8195 -14962.5874
69 -72004.9465 -100725.8195
70 -122294.6437 -72004.9465
71 1252.3272 -122294.6437
72 -88696.6011 1252.3272
73 15111.9538 -88696.6011
74 -34736.9165 15111.9538
75 -98722.3509 -34736.9165
76 -28470.4246 -98722.3509
77 -73253.8756 -28470.4246
78 -42228.7857 -73253.8756
79 -37810.1071 -42228.7857
80 27095.5975 -37810.1071
81 -49657.3672 27095.5975
82 -82702.7231 -49657.3672
83 41679.9422 -82702.7231
84 24950.3027 41679.9422
85 -102379.3908 24950.3027
86 -29329.7313 -102379.3908
87 1171.4484 -29329.7313
88 -115198.1764 1171.4484
89 32516.3322 -115198.1764
90 -1902.9294 32516.3322
91 -21807.0740 -1902.9294
92 69502.3235 -21807.0740
93 -172301.4464 69502.3235
94 49811.5996 -172301.4464
95 19619.0604 49811.5996
96 71821.8316 19619.0604
97 2140.2336 71821.8316
98 -41173.8667 2140.2336
99 -1473.8500 -41173.8667
100 -35893.9515 -1473.8500
101 -112498.0489 -35893.9515
102 -165507.6930 -112498.0489
103 -50102.0373 -165507.6930
104 30070.0947 -50102.0373
105 -53263.1917 30070.0947
106 -252162.8870 -53263.1917
107 131136.6578 -252162.8870
108 -175232.6525 131136.6578
109 -92237.8474 -175232.6525
110 94519.8432 -92237.8474
111 -190793.4409 94519.8432
112 -165268.7848 -190793.4409
113 -57466.9920 -165268.7848
114 86189.5153 -57466.9920
115 21656.0303 86189.5153
116 17524.5950 21656.0303
117 72490.9475 17524.5950
118 -132145.5451 72490.9475
119 14988.2296 -132145.5451
120 -288568.0383 14988.2296
121 32827.5084 -288568.0383
122 -103083.9599 32827.5084
123 84905.4238 -103083.9599
124 13946.7576 84905.4238
125 -162469.9198 13946.7576
126 55152.5245 -162469.9198
127 42028.6252 55152.5245
128 29692.3360 42028.6252
129 55152.5245 29692.3360
130 46453.0982 55152.5245
131 57660.3793 46453.0982
132 55152.5245 57660.3793
133 31224.0439 55152.5245
134 33610.3153 31224.0439
135 44124.6501 33610.3153
136 45697.8939 44124.6501
137 53540.1528 45697.8939
138 30900.3954 53540.1528
139 55152.5245 30900.3954
140 51353.5430 55152.5245
141 55152.5245 51353.5430
142 11405.4575 55152.5245
143 50503.8779 11405.4575
144 40171.2046 50503.8779
145 43328.0424 40171.2046
146 55152.5245 43328.0424
147 47315.9545 55152.5245
148 57036.7210 47315.9545
149 54877.0451 57036.7210
150 16156.0728 54877.0451
151 -39133.1674 16156.0728
152 50106.0545 -39133.1674
153 23675.2486 50106.0545
154 52216.7029 23675.2486
155 55152.5245 52216.7029
156 55152.5245 55152.5245
157 49624.7042 55152.5245
158 62038.6981 49624.7042
159 55152.5245 62038.6981
160 55152.5245 55152.5245
161 55152.5245 55152.5245
162 55152.5245 55152.5245
163 28169.4235 55152.5245
164 3755.9701 28169.4235
165 52610.4994 3755.9701
166 55152.5245 52610.4994
167 111226.8553 55152.5245
168 -76651.0561 111226.8553
169 55152.5245 -76651.0561
170 55152.5245 55152.5245
171 47040.3092 55152.5245
172 55152.5245 47040.3092
173 33905.3464 55152.5245
174 55152.5245 33905.3464
175 44108.9001 55152.5245
176 -69816.9214 44108.9001
177 54346.2088 -69816.9214
178 55152.5245 54346.2088
179 55152.5245 55152.5245
180 41401.1884 55152.5245
181 29918.9688 41401.1884
182 -17621.7771 29918.9688
183 -4778.0021 -17621.7771
184 46147.5877 -4778.0021
185 48383.3242 46147.5877
186 11517.2049 48383.3242
187 41696.5179 11517.2049
188 13656.0689 41696.5179
189 55152.5245 13656.0689
190 47311.6999 55152.5245
191 55152.5245 47311.6999
192 43541.9820 55152.5245
193 55152.5245 43541.9820
194 51549.5701 55152.5245
195 55152.5245 51549.5701
196 -132865.9295 55152.5245
197 55152.5245 -132865.9295
198 55152.5245 55152.5245
199 -50518.7341 55152.5245
200 55152.5245 -50518.7341
201 55152.5245 55152.5245
202 55152.5245 55152.5245
203 39706.6063 55152.5245
204 389.6725 39706.6063
205 49192.8281 389.6725
206 47148.8882 49192.8281
207 55152.5245 47148.8882
208 51266.1839 55152.5245
209 37858.5204 51266.1839
210 55152.5245 37858.5204
211 55152.5245 55152.5245
212 57681.9685 55152.5245
213 6257.6673 57681.9685
214 24474.6896 6257.6673
215 53133.9365 24474.6896
216 48900.4290 53133.9365
217 41819.3564 48900.4290
218 53133.9365 41819.3564
219 -187.5434 53133.9365
220 27389.5375 -187.5434
221 55152.5245 27389.5375
222 54085.6898 55152.5245
223 21444.3644 54085.6898
224 53765.6010 21444.3644
225 25296.9108 53765.6010
226 43896.9295 25296.9108
227 55152.5245 43896.9295
228 38720.0837 55152.5245
229 55152.5245 38720.0837
230 19417.3320 55152.5245
231 35817.3986 19417.3320
232 55152.5245 35817.3986
233 20264.8536 55152.5245
234 55152.5245 20264.8536
235 55179.4190 55152.5245
236 55152.5245 55179.4190
237 48924.8316 55152.5245
238 55152.5245 48924.8316
239 55152.5245 55152.5245
240 55152.5245 55152.5245
241 46558.3712 55152.5245
242 55152.5245 46558.3712
243 49230.8366 55152.5245
244 -7960.0214 49230.8366
245 55152.5245 -7960.0214
246 55152.5245 55152.5245
247 29031.6010 55152.5245
248 55152.5245 29031.6010
249 55152.5245 55152.5245
250 34791.5919 55152.5245
251 55152.5245 34791.5919
252 36879.8403 55152.5245
253 -27169.8902 36879.8403
254 -149946.2072 -27169.8902
255 55152.5245 -149946.2072
256 43485.6542 55152.5245
257 24396.6177 43485.6542
258 55152.5245 24396.6177
259 55152.5245 55152.5245
260 55152.5245 55152.5245
261 55152.5245 55152.5245
262 46603.3938 55152.5245
263 57799.4030 46603.3938
264 55152.5245 57799.4030
265 55152.5245 55152.5245
266 55152.5245 55152.5245
267 81168.2234 55152.5245
268 55152.5245 81168.2234
269 27644.5278 55152.5245
270 55152.5245 27644.5278
271 55152.5245 55152.5245
272 55152.5245 55152.5245
273 27922.6546 55152.5245
274 55152.5245 27922.6546
275 38752.0513 55152.5245
276 24162.5503 38752.0513
277 55152.5245 24162.5503
278 -37344.5929 55152.5245
279 26861.7964 -37344.5929
280 55152.5245 26861.7964
281 55152.5245 55152.5245
282 27804.4032 55152.5245
283 55152.5245 27804.4032
284 38853.4518 55152.5245
285 55152.5245 38853.4518
286 40704.7058 55152.5245
287 52103.9773 40704.7058
288 55152.5245 52103.9773
289 23957.2150 55152.5245
290 55152.5245 23957.2150
291 45877.7435 55152.5245
292 44577.6335 45877.7435
293 55152.5245 44577.6335
294 32825.8604 55152.5245
295 65372.1884 32825.8604
296 55152.5245 65372.1884
297 55152.5245 55152.5245
298 55152.5245 55152.5245
299 46068.8784 55152.5245
300 48116.5426 46068.8784
301 36569.8314 48116.5426
302 54143.2305 36569.8314
303 55152.5245 54143.2305
304 60303.2396 55152.5245
305 42987.4385 60303.2396
306 55152.5245 42987.4385
307 27472.0702 55152.5245
308 33465.2963 27472.0702
309 55152.5245 33465.2963
310 55152.5245 55152.5245
311 4530.3833 55152.5245
312 55152.5245 4530.3833
313 41282.3000 55152.5245
314 47223.8192 41282.3000
315 63986.5389 47223.8192
316 55180.1335 63986.5389
317 7924.8219 55180.1335
318 50077.6074 7924.8219
319 55152.5245 50077.6074
320 55152.5245 55152.5245
321 55152.5245 55152.5245
322 55152.5245 55152.5245
323 55152.5245 55152.5245
324 55152.5245 55152.5245
325 50777.7514 55152.5245
326 55152.5245 50777.7514
327 55152.5245 55152.5245
328 55152.5245 55152.5245
329 55152.5245 55152.5245
330 55152.5245 55152.5245
331 55152.5245 55152.5245
332 55152.5245 55152.5245
333 47885.6972 55152.5245
334 54118.2358 47885.6972
335 44193.4641 54118.2358
336 -3177.0170 44193.4641
337 55152.5245 -3177.0170
338 51794.7786 55152.5245
339 46172.1334 51794.7786
340 13020.8195 46172.1334
341 55152.5245 13020.8195
342 55152.5245 55152.5245
343 55152.5245 55152.5245
344 21236.0425 55152.5245
345 55152.5245 21236.0425
346 40314.4511 55152.5245
347 55152.5245 40314.4511
348 55152.5245 55152.5245
349 13933.7738 55152.5245
350 51894.0439 13933.7738
351 44801.8786 51894.0439
352 44219.9617 44801.8786
353 47891.5630 44219.9617
354 55152.5245 47891.5630
355 55152.5245 55152.5245
356 7368.8801 55152.5245
357 55152.5245 7368.8801
358 55152.5245 55152.5245
359 55152.5245 55152.5245
360 51108.1960 55152.5245
361 55152.5245 51108.1960
362 55152.5245 55152.5245
363 -20345.3694 55152.5245
364 44813.8486 -20345.3694
365 55152.5245 44813.8486
366 55152.5245 55152.5245
367 55152.5245 55152.5245
368 55152.5245 55152.5245
369 43461.0754 55152.5245
370 55152.5245 43461.0754
371 47678.2739 55152.5245
372 39372.9095 47678.2739
373 33907.2837 39372.9095
374 55152.5245 33907.2837
375 55152.5245 55152.5245
376 28500.3161 55152.5245
377 42129.4647 28500.3161
378 845.1762 42129.4647
379 55152.5245 845.1762
380 55152.5245 55152.5245
381 7288.2607 55152.5245
382 23651.2435 7288.2607
383 55152.5245 23651.2435
384 40168.2486 55152.5245
385 -117279.6312 40168.2486
386 62592.5952 -117279.6312
387 25728.1516 62592.5952
388 49874.5640 25728.1516
389 30391.7616 49874.5640
390 40145.1899 30391.7616
391 18133.7508 40145.1899
392 -19080.2469 18133.7508
393 74580.5292 -19080.2469
394 -96721.9437 74580.5292
395 6463.4051 -96721.9437
396 -138664.7201 6463.4051
397 -189519.7904 -138664.7201
398 12381.6147 -189519.7904
399 -55495.2651 12381.6147
400 2431.9986 -55495.2651
401 42613.8359 2431.9986
402 19583.9874 42613.8359
403 -49565.6943 19583.9874
404 -8995.7460 -49565.6943
405 14011.1662 -8995.7460
406 20618.4595 14011.1662
407 13250.1409 20618.4595
408 22906.3494 13250.1409
409 -61901.5139 22906.3494
410 -51361.8573 -61901.5139
411 -161769.0794 -51361.8573
412 -25907.2560 -161769.0794
413 60088.3931 -25907.2560
414 -177525.6236 60088.3931
415 -99195.2033 -177525.6236
416 91983.5806 -99195.2033
417 -17379.6380 91983.5806
418 -50409.1687 -17379.6380
419 -37404.0006 -50409.1687
420 -44265.0154 -37404.0006
421 -175393.8602 -44265.0154
422 -547870.2678 -175393.8602
423 -458560.1402 -547870.2678
424 -48844.9394 -458560.1402
425 -250014.5348 -48844.9394
426 -311858.7540 -250014.5348
427 -1135948.9209 -311858.7540
428 -272147.0098 -1135948.9209
429 -336049.2741 -272147.0098
430 -682483.1396 -336049.2741
431 NA -682483.1396
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] 2923814.3860 1389636.7423
[2,] 1363346.2455 2923814.3860
[3,] -2198539.3196 1363346.2455
[4,] -394043.3532 -2198539.3196
[5,] 109869.7630 -394043.3532
[6,] -135939.7551 109869.7630
[7,] 128928.3808 -135939.7551
[8,] -336334.3876 128928.3808
[9,] 715092.5993 -336334.3876
[10,] -150657.5105 715092.5993
[11,] 342307.7668 -150657.5105
[12,] 449675.2072 342307.7668
[13,] -1139741.2853 449675.2072
[14,] -214255.2259 -1139741.2853
[15,] 460574.8057 -214255.2259
[16,] -237906.6182 460574.8057
[17,] -238071.5057 -237906.6182
[18,] -295877.9417 -238071.5057
[19,] -585102.6254 -295877.9417
[20,] 203463.1027 -585102.6254
[21,] -137416.1514 203463.1027
[22,] -119601.7096 -137416.1514
[23,] -638049.8630 -119601.7096
[24,] -318324.0983 -638049.8630
[25,] -84297.5193 -318324.0983
[26,] -471896.4012 -84297.5193
[27,] -340384.4471 -471896.4012
[28,] 348299.6853 -340384.4471
[29,] -261059.4968 348299.6853
[30,] 286668.0415 -261059.4968
[31,] -87937.6349 286668.0415
[32,] 255919.3586 -87937.6349
[33,] 130738.8627 255919.3586
[34,] -54941.9881 130738.8627
[35,] -569143.8636 -54941.9881
[36,] 34310.7675 -569143.8636
[37,] 157833.6426 34310.7675
[38,] -51015.8179 157833.6426
[39,] -304565.6193 -51015.8179
[40,] -459620.0553 -304565.6193
[41,] -329011.9718 -459620.0553
[42,] 172.3323 -329011.9718
[43,] -209375.2579 172.3323
[44,] -27239.0962 -209375.2579
[45,] -53789.6599 -27239.0962
[46,] -281553.3224 -53789.6599
[47,] -264256.0037 -281553.3224
[48,] -53547.0125 -264256.0037
[49,] -181997.0873 -53547.0125
[50,] 35733.0747 -181997.0873
[51,] -86617.6159 35733.0747
[52,] -228211.1136 -86617.6159
[53,] 44146.0101 -228211.1136
[54,] -394458.1521 44146.0101
[55,] 14558.8138 -394458.1521
[56,] 320392.3041 14558.8138
[57,] -22047.0390 320392.3041
[58,] -59458.7651 -22047.0390
[59,] -205533.6443 -59458.7651
[60,] 14509.4419 -205533.6443
[61,] -179052.6824 14509.4419
[62,] 7020.6734 -179052.6824
[63,] -174141.9735 7020.6734
[64,] -15771.8893 -174141.9735
[65,] -494796.9414 -15771.8893
[66,] -65004.9696 -494796.9414
[67,] -14962.5874 -65004.9696
[68,] -100725.8195 -14962.5874
[69,] -72004.9465 -100725.8195
[70,] -122294.6437 -72004.9465
[71,] 1252.3272 -122294.6437
[72,] -88696.6011 1252.3272
[73,] 15111.9538 -88696.6011
[74,] -34736.9165 15111.9538
[75,] -98722.3509 -34736.9165
[76,] -28470.4246 -98722.3509
[77,] -73253.8756 -28470.4246
[78,] -42228.7857 -73253.8756
[79,] -37810.1071 -42228.7857
[80,] 27095.5975 -37810.1071
[81,] -49657.3672 27095.5975
[82,] -82702.7231 -49657.3672
[83,] 41679.9422 -82702.7231
[84,] 24950.3027 41679.9422
[85,] -102379.3908 24950.3027
[86,] -29329.7313 -102379.3908
[87,] 1171.4484 -29329.7313
[88,] -115198.1764 1171.4484
[89,] 32516.3322 -115198.1764
[90,] -1902.9294 32516.3322
[91,] -21807.0740 -1902.9294
[92,] 69502.3235 -21807.0740
[93,] -172301.4464 69502.3235
[94,] 49811.5996 -172301.4464
[95,] 19619.0604 49811.5996
[96,] 71821.8316 19619.0604
[97,] 2140.2336 71821.8316
[98,] -41173.8667 2140.2336
[99,] -1473.8500 -41173.8667
[100,] -35893.9515 -1473.8500
[101,] -112498.0489 -35893.9515
[102,] -165507.6930 -112498.0489
[103,] -50102.0373 -165507.6930
[104,] 30070.0947 -50102.0373
[105,] -53263.1917 30070.0947
[106,] -252162.8870 -53263.1917
[107,] 131136.6578 -252162.8870
[108,] -175232.6525 131136.6578
[109,] -92237.8474 -175232.6525
[110,] 94519.8432 -92237.8474
[111,] -190793.4409 94519.8432
[112,] -165268.7848 -190793.4409
[113,] -57466.9920 -165268.7848
[114,] 86189.5153 -57466.9920
[115,] 21656.0303 86189.5153
[116,] 17524.5950 21656.0303
[117,] 72490.9475 17524.5950
[118,] -132145.5451 72490.9475
[119,] 14988.2296 -132145.5451
[120,] -288568.0383 14988.2296
[121,] 32827.5084 -288568.0383
[122,] -103083.9599 32827.5084
[123,] 84905.4238 -103083.9599
[124,] 13946.7576 84905.4238
[125,] -162469.9198 13946.7576
[126,] 55152.5245 -162469.9198
[127,] 42028.6252 55152.5245
[128,] 29692.3360 42028.6252
[129,] 55152.5245 29692.3360
[130,] 46453.0982 55152.5245
[131,] 57660.3793 46453.0982
[132,] 55152.5245 57660.3793
[133,] 31224.0439 55152.5245
[134,] 33610.3153 31224.0439
[135,] 44124.6501 33610.3153
[136,] 45697.8939 44124.6501
[137,] 53540.1528 45697.8939
[138,] 30900.3954 53540.1528
[139,] 55152.5245 30900.3954
[140,] 51353.5430 55152.5245
[141,] 55152.5245 51353.5430
[142,] 11405.4575 55152.5245
[143,] 50503.8779 11405.4575
[144,] 40171.2046 50503.8779
[145,] 43328.0424 40171.2046
[146,] 55152.5245 43328.0424
[147,] 47315.9545 55152.5245
[148,] 57036.7210 47315.9545
[149,] 54877.0451 57036.7210
[150,] 16156.0728 54877.0451
[151,] -39133.1674 16156.0728
[152,] 50106.0545 -39133.1674
[153,] 23675.2486 50106.0545
[154,] 52216.7029 23675.2486
[155,] 55152.5245 52216.7029
[156,] 55152.5245 55152.5245
[157,] 49624.7042 55152.5245
[158,] 62038.6981 49624.7042
[159,] 55152.5245 62038.6981
[160,] 55152.5245 55152.5245
[161,] 55152.5245 55152.5245
[162,] 55152.5245 55152.5245
[163,] 28169.4235 55152.5245
[164,] 3755.9701 28169.4235
[165,] 52610.4994 3755.9701
[166,] 55152.5245 52610.4994
[167,] 111226.8553 55152.5245
[168,] -76651.0561 111226.8553
[169,] 55152.5245 -76651.0561
[170,] 55152.5245 55152.5245
[171,] 47040.3092 55152.5245
[172,] 55152.5245 47040.3092
[173,] 33905.3464 55152.5245
[174,] 55152.5245 33905.3464
[175,] 44108.9001 55152.5245
[176,] -69816.9214 44108.9001
[177,] 54346.2088 -69816.9214
[178,] 55152.5245 54346.2088
[179,] 55152.5245 55152.5245
[180,] 41401.1884 55152.5245
[181,] 29918.9688 41401.1884
[182,] -17621.7771 29918.9688
[183,] -4778.0021 -17621.7771
[184,] 46147.5877 -4778.0021
[185,] 48383.3242 46147.5877
[186,] 11517.2049 48383.3242
[187,] 41696.5179 11517.2049
[188,] 13656.0689 41696.5179
[189,] 55152.5245 13656.0689
[190,] 47311.6999 55152.5245
[191,] 55152.5245 47311.6999
[192,] 43541.9820 55152.5245
[193,] 55152.5245 43541.9820
[194,] 51549.5701 55152.5245
[195,] 55152.5245 51549.5701
[196,] -132865.9295 55152.5245
[197,] 55152.5245 -132865.9295
[198,] 55152.5245 55152.5245
[199,] -50518.7341 55152.5245
[200,] 55152.5245 -50518.7341
[201,] 55152.5245 55152.5245
[202,] 55152.5245 55152.5245
[203,] 39706.6063 55152.5245
[204,] 389.6725 39706.6063
[205,] 49192.8281 389.6725
[206,] 47148.8882 49192.8281
[207,] 55152.5245 47148.8882
[208,] 51266.1839 55152.5245
[209,] 37858.5204 51266.1839
[210,] 55152.5245 37858.5204
[211,] 55152.5245 55152.5245
[212,] 57681.9685 55152.5245
[213,] 6257.6673 57681.9685
[214,] 24474.6896 6257.6673
[215,] 53133.9365 24474.6896
[216,] 48900.4290 53133.9365
[217,] 41819.3564 48900.4290
[218,] 53133.9365 41819.3564
[219,] -187.5434 53133.9365
[220,] 27389.5375 -187.5434
[221,] 55152.5245 27389.5375
[222,] 54085.6898 55152.5245
[223,] 21444.3644 54085.6898
[224,] 53765.6010 21444.3644
[225,] 25296.9108 53765.6010
[226,] 43896.9295 25296.9108
[227,] 55152.5245 43896.9295
[228,] 38720.0837 55152.5245
[229,] 55152.5245 38720.0837
[230,] 19417.3320 55152.5245
[231,] 35817.3986 19417.3320
[232,] 55152.5245 35817.3986
[233,] 20264.8536 55152.5245
[234,] 55152.5245 20264.8536
[235,] 55179.4190 55152.5245
[236,] 55152.5245 55179.4190
[237,] 48924.8316 55152.5245
[238,] 55152.5245 48924.8316
[239,] 55152.5245 55152.5245
[240,] 55152.5245 55152.5245
[241,] 46558.3712 55152.5245
[242,] 55152.5245 46558.3712
[243,] 49230.8366 55152.5245
[244,] -7960.0214 49230.8366
[245,] 55152.5245 -7960.0214
[246,] 55152.5245 55152.5245
[247,] 29031.6010 55152.5245
[248,] 55152.5245 29031.6010
[249,] 55152.5245 55152.5245
[250,] 34791.5919 55152.5245
[251,] 55152.5245 34791.5919
[252,] 36879.8403 55152.5245
[253,] -27169.8902 36879.8403
[254,] -149946.2072 -27169.8902
[255,] 55152.5245 -149946.2072
[256,] 43485.6542 55152.5245
[257,] 24396.6177 43485.6542
[258,] 55152.5245 24396.6177
[259,] 55152.5245 55152.5245
[260,] 55152.5245 55152.5245
[261,] 55152.5245 55152.5245
[262,] 46603.3938 55152.5245
[263,] 57799.4030 46603.3938
[264,] 55152.5245 57799.4030
[265,] 55152.5245 55152.5245
[266,] 55152.5245 55152.5245
[267,] 81168.2234 55152.5245
[268,] 55152.5245 81168.2234
[269,] 27644.5278 55152.5245
[270,] 55152.5245 27644.5278
[271,] 55152.5245 55152.5245
[272,] 55152.5245 55152.5245
[273,] 27922.6546 55152.5245
[274,] 55152.5245 27922.6546
[275,] 38752.0513 55152.5245
[276,] 24162.5503 38752.0513
[277,] 55152.5245 24162.5503
[278,] -37344.5929 55152.5245
[279,] 26861.7964 -37344.5929
[280,] 55152.5245 26861.7964
[281,] 55152.5245 55152.5245
[282,] 27804.4032 55152.5245
[283,] 55152.5245 27804.4032
[284,] 38853.4518 55152.5245
[285,] 55152.5245 38853.4518
[286,] 40704.7058 55152.5245
[287,] 52103.9773 40704.7058
[288,] 55152.5245 52103.9773
[289,] 23957.2150 55152.5245
[290,] 55152.5245 23957.2150
[291,] 45877.7435 55152.5245
[292,] 44577.6335 45877.7435
[293,] 55152.5245 44577.6335
[294,] 32825.8604 55152.5245
[295,] 65372.1884 32825.8604
[296,] 55152.5245 65372.1884
[297,] 55152.5245 55152.5245
[298,] 55152.5245 55152.5245
[299,] 46068.8784 55152.5245
[300,] 48116.5426 46068.8784
[301,] 36569.8314 48116.5426
[302,] 54143.2305 36569.8314
[303,] 55152.5245 54143.2305
[304,] 60303.2396 55152.5245
[305,] 42987.4385 60303.2396
[306,] 55152.5245 42987.4385
[307,] 27472.0702 55152.5245
[308,] 33465.2963 27472.0702
[309,] 55152.5245 33465.2963
[310,] 55152.5245 55152.5245
[311,] 4530.3833 55152.5245
[312,] 55152.5245 4530.3833
[313,] 41282.3000 55152.5245
[314,] 47223.8192 41282.3000
[315,] 63986.5389 47223.8192
[316,] 55180.1335 63986.5389
[317,] 7924.8219 55180.1335
[318,] 50077.6074 7924.8219
[319,] 55152.5245 50077.6074
[320,] 55152.5245 55152.5245
[321,] 55152.5245 55152.5245
[322,] 55152.5245 55152.5245
[323,] 55152.5245 55152.5245
[324,] 55152.5245 55152.5245
[325,] 50777.7514 55152.5245
[326,] 55152.5245 50777.7514
[327,] 55152.5245 55152.5245
[328,] 55152.5245 55152.5245
[329,] 55152.5245 55152.5245
[330,] 55152.5245 55152.5245
[331,] 55152.5245 55152.5245
[332,] 55152.5245 55152.5245
[333,] 47885.6972 55152.5245
[334,] 54118.2358 47885.6972
[335,] 44193.4641 54118.2358
[336,] -3177.0170 44193.4641
[337,] 55152.5245 -3177.0170
[338,] 51794.7786 55152.5245
[339,] 46172.1334 51794.7786
[340,] 13020.8195 46172.1334
[341,] 55152.5245 13020.8195
[342,] 55152.5245 55152.5245
[343,] 55152.5245 55152.5245
[344,] 21236.0425 55152.5245
[345,] 55152.5245 21236.0425
[346,] 40314.4511 55152.5245
[347,] 55152.5245 40314.4511
[348,] 55152.5245 55152.5245
[349,] 13933.7738 55152.5245
[350,] 51894.0439 13933.7738
[351,] 44801.8786 51894.0439
[352,] 44219.9617 44801.8786
[353,] 47891.5630 44219.9617
[354,] 55152.5245 47891.5630
[355,] 55152.5245 55152.5245
[356,] 7368.8801 55152.5245
[357,] 55152.5245 7368.8801
[358,] 55152.5245 55152.5245
[359,] 55152.5245 55152.5245
[360,] 51108.1960 55152.5245
[361,] 55152.5245 51108.1960
[362,] 55152.5245 55152.5245
[363,] -20345.3694 55152.5245
[364,] 44813.8486 -20345.3694
[365,] 55152.5245 44813.8486
[366,] 55152.5245 55152.5245
[367,] 55152.5245 55152.5245
[368,] 55152.5245 55152.5245
[369,] 43461.0754 55152.5245
[370,] 55152.5245 43461.0754
[371,] 47678.2739 55152.5245
[372,] 39372.9095 47678.2739
[373,] 33907.2837 39372.9095
[374,] 55152.5245 33907.2837
[375,] 55152.5245 55152.5245
[376,] 28500.3161 55152.5245
[377,] 42129.4647 28500.3161
[378,] 845.1762 42129.4647
[379,] 55152.5245 845.1762
[380,] 55152.5245 55152.5245
[381,] 7288.2607 55152.5245
[382,] 23651.2435 7288.2607
[383,] 55152.5245 23651.2435
[384,] 40168.2486 55152.5245
[385,] -117279.6312 40168.2486
[386,] 62592.5952 -117279.6312
[387,] 25728.1516 62592.5952
[388,] 49874.5640 25728.1516
[389,] 30391.7616 49874.5640
[390,] 40145.1899 30391.7616
[391,] 18133.7508 40145.1899
[392,] -19080.2469 18133.7508
[393,] 74580.5292 -19080.2469
[394,] -96721.9437 74580.5292
[395,] 6463.4051 -96721.9437
[396,] -138664.7201 6463.4051
[397,] -189519.7904 -138664.7201
[398,] 12381.6147 -189519.7904
[399,] -55495.2651 12381.6147
[400,] 2431.9986 -55495.2651
[401,] 42613.8359 2431.9986
[402,] 19583.9874 42613.8359
[403,] -49565.6943 19583.9874
[404,] -8995.7460 -49565.6943
[405,] 14011.1662 -8995.7460
[406,] 20618.4595 14011.1662
[407,] 13250.1409 20618.4595
[408,] 22906.3494 13250.1409
[409,] -61901.5139 22906.3494
[410,] -51361.8573 -61901.5139
[411,] -161769.0794 -51361.8573
[412,] -25907.2560 -161769.0794
[413,] 60088.3931 -25907.2560
[414,] -177525.6236 60088.3931
[415,] -99195.2033 -177525.6236
[416,] 91983.5806 -99195.2033
[417,] -17379.6380 91983.5806
[418,] -50409.1687 -17379.6380
[419,] -37404.0006 -50409.1687
[420,] -44265.0154 -37404.0006
[421,] -175393.8602 -44265.0154
[422,] -547870.2678 -175393.8602
[423,] -458560.1402 -547870.2678
[424,] -48844.9394 -458560.1402
[425,] -250014.5348 -48844.9394
[426,] -311858.7540 -250014.5348
[427,] -1135948.9209 -311858.7540
[428,] -272147.0098 -1135948.9209
[429,] -336049.2741 -272147.0098
[430,] -682483.1396 -336049.2741
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 2923814.3860 1389636.7423
2 1363346.2455 2923814.3860
3 -2198539.3196 1363346.2455
4 -394043.3532 -2198539.3196
5 109869.7630 -394043.3532
6 -135939.7551 109869.7630
7 128928.3808 -135939.7551
8 -336334.3876 128928.3808
9 715092.5993 -336334.3876
10 -150657.5105 715092.5993
11 342307.7668 -150657.5105
12 449675.2072 342307.7668
13 -1139741.2853 449675.2072
14 -214255.2259 -1139741.2853
15 460574.8057 -214255.2259
16 -237906.6182 460574.8057
17 -238071.5057 -237906.6182
18 -295877.9417 -238071.5057
19 -585102.6254 -295877.9417
20 203463.1027 -585102.6254
21 -137416.1514 203463.1027
22 -119601.7096 -137416.1514
23 -638049.8630 -119601.7096
24 -318324.0983 -638049.8630
25 -84297.5193 -318324.0983
26 -471896.4012 -84297.5193
27 -340384.4471 -471896.4012
28 348299.6853 -340384.4471
29 -261059.4968 348299.6853
30 286668.0415 -261059.4968
31 -87937.6349 286668.0415
32 255919.3586 -87937.6349
33 130738.8627 255919.3586
34 -54941.9881 130738.8627
35 -569143.8636 -54941.9881
36 34310.7675 -569143.8636
37 157833.6426 34310.7675
38 -51015.8179 157833.6426
39 -304565.6193 -51015.8179
40 -459620.0553 -304565.6193
41 -329011.9718 -459620.0553
42 172.3323 -329011.9718
43 -209375.2579 172.3323
44 -27239.0962 -209375.2579
45 -53789.6599 -27239.0962
46 -281553.3224 -53789.6599
47 -264256.0037 -281553.3224
48 -53547.0125 -264256.0037
49 -181997.0873 -53547.0125
50 35733.0747 -181997.0873
51 -86617.6159 35733.0747
52 -228211.1136 -86617.6159
53 44146.0101 -228211.1136
54 -394458.1521 44146.0101
55 14558.8138 -394458.1521
56 320392.3041 14558.8138
57 -22047.0390 320392.3041
58 -59458.7651 -22047.0390
59 -205533.6443 -59458.7651
60 14509.4419 -205533.6443
61 -179052.6824 14509.4419
62 7020.6734 -179052.6824
63 -174141.9735 7020.6734
64 -15771.8893 -174141.9735
65 -494796.9414 -15771.8893
66 -65004.9696 -494796.9414
67 -14962.5874 -65004.9696
68 -100725.8195 -14962.5874
69 -72004.9465 -100725.8195
70 -122294.6437 -72004.9465
71 1252.3272 -122294.6437
72 -88696.6011 1252.3272
73 15111.9538 -88696.6011
74 -34736.9165 15111.9538
75 -98722.3509 -34736.9165
76 -28470.4246 -98722.3509
77 -73253.8756 -28470.4246
78 -42228.7857 -73253.8756
79 -37810.1071 -42228.7857
80 27095.5975 -37810.1071
81 -49657.3672 27095.5975
82 -82702.7231 -49657.3672
83 41679.9422 -82702.7231
84 24950.3027 41679.9422
85 -102379.3908 24950.3027
86 -29329.7313 -102379.3908
87 1171.4484 -29329.7313
88 -115198.1764 1171.4484
89 32516.3322 -115198.1764
90 -1902.9294 32516.3322
91 -21807.0740 -1902.9294
92 69502.3235 -21807.0740
93 -172301.4464 69502.3235
94 49811.5996 -172301.4464
95 19619.0604 49811.5996
96 71821.8316 19619.0604
97 2140.2336 71821.8316
98 -41173.8667 2140.2336
99 -1473.8500 -41173.8667
100 -35893.9515 -1473.8500
101 -112498.0489 -35893.9515
102 -165507.6930 -112498.0489
103 -50102.0373 -165507.6930
104 30070.0947 -50102.0373
105 -53263.1917 30070.0947
106 -252162.8870 -53263.1917
107 131136.6578 -252162.8870
108 -175232.6525 131136.6578
109 -92237.8474 -175232.6525
110 94519.8432 -92237.8474
111 -190793.4409 94519.8432
112 -165268.7848 -190793.4409
113 -57466.9920 -165268.7848
114 86189.5153 -57466.9920
115 21656.0303 86189.5153
116 17524.5950 21656.0303
117 72490.9475 17524.5950
118 -132145.5451 72490.9475
119 14988.2296 -132145.5451
120 -288568.0383 14988.2296
121 32827.5084 -288568.0383
122 -103083.9599 32827.5084
123 84905.4238 -103083.9599
124 13946.7576 84905.4238
125 -162469.9198 13946.7576
126 55152.5245 -162469.9198
127 42028.6252 55152.5245
128 29692.3360 42028.6252
129 55152.5245 29692.3360
130 46453.0982 55152.5245
131 57660.3793 46453.0982
132 55152.5245 57660.3793
133 31224.0439 55152.5245
134 33610.3153 31224.0439
135 44124.6501 33610.3153
136 45697.8939 44124.6501
137 53540.1528 45697.8939
138 30900.3954 53540.1528
139 55152.5245 30900.3954
140 51353.5430 55152.5245
141 55152.5245 51353.5430
142 11405.4575 55152.5245
143 50503.8779 11405.4575
144 40171.2046 50503.8779
145 43328.0424 40171.2046
146 55152.5245 43328.0424
147 47315.9545 55152.5245
148 57036.7210 47315.9545
149 54877.0451 57036.7210
150 16156.0728 54877.0451
151 -39133.1674 16156.0728
152 50106.0545 -39133.1674
153 23675.2486 50106.0545
154 52216.7029 23675.2486
155 55152.5245 52216.7029
156 55152.5245 55152.5245
157 49624.7042 55152.5245
158 62038.6981 49624.7042
159 55152.5245 62038.6981
160 55152.5245 55152.5245
161 55152.5245 55152.5245
162 55152.5245 55152.5245
163 28169.4235 55152.5245
164 3755.9701 28169.4235
165 52610.4994 3755.9701
166 55152.5245 52610.4994
167 111226.8553 55152.5245
168 -76651.0561 111226.8553
169 55152.5245 -76651.0561
170 55152.5245 55152.5245
171 47040.3092 55152.5245
172 55152.5245 47040.3092
173 33905.3464 55152.5245
174 55152.5245 33905.3464
175 44108.9001 55152.5245
176 -69816.9214 44108.9001
177 54346.2088 -69816.9214
178 55152.5245 54346.2088
179 55152.5245 55152.5245
180 41401.1884 55152.5245
181 29918.9688 41401.1884
182 -17621.7771 29918.9688
183 -4778.0021 -17621.7771
184 46147.5877 -4778.0021
185 48383.3242 46147.5877
186 11517.2049 48383.3242
187 41696.5179 11517.2049
188 13656.0689 41696.5179
189 55152.5245 13656.0689
190 47311.6999 55152.5245
191 55152.5245 47311.6999
192 43541.9820 55152.5245
193 55152.5245 43541.9820
194 51549.5701 55152.5245
195 55152.5245 51549.5701
196 -132865.9295 55152.5245
197 55152.5245 -132865.9295
198 55152.5245 55152.5245
199 -50518.7341 55152.5245
200 55152.5245 -50518.7341
201 55152.5245 55152.5245
202 55152.5245 55152.5245
203 39706.6063 55152.5245
204 389.6725 39706.6063
205 49192.8281 389.6725
206 47148.8882 49192.8281
207 55152.5245 47148.8882
208 51266.1839 55152.5245
209 37858.5204 51266.1839
210 55152.5245 37858.5204
211 55152.5245 55152.5245
212 57681.9685 55152.5245
213 6257.6673 57681.9685
214 24474.6896 6257.6673
215 53133.9365 24474.6896
216 48900.4290 53133.9365
217 41819.3564 48900.4290
218 53133.9365 41819.3564
219 -187.5434 53133.9365
220 27389.5375 -187.5434
221 55152.5245 27389.5375
222 54085.6898 55152.5245
223 21444.3644 54085.6898
224 53765.6010 21444.3644
225 25296.9108 53765.6010
226 43896.9295 25296.9108
227 55152.5245 43896.9295
228 38720.0837 55152.5245
229 55152.5245 38720.0837
230 19417.3320 55152.5245
231 35817.3986 19417.3320
232 55152.5245 35817.3986
233 20264.8536 55152.5245
234 55152.5245 20264.8536
235 55179.4190 55152.5245
236 55152.5245 55179.4190
237 48924.8316 55152.5245
238 55152.5245 48924.8316
239 55152.5245 55152.5245
240 55152.5245 55152.5245
241 46558.3712 55152.5245
242 55152.5245 46558.3712
243 49230.8366 55152.5245
244 -7960.0214 49230.8366
245 55152.5245 -7960.0214
246 55152.5245 55152.5245
247 29031.6010 55152.5245
248 55152.5245 29031.6010
249 55152.5245 55152.5245
250 34791.5919 55152.5245
251 55152.5245 34791.5919
252 36879.8403 55152.5245
253 -27169.8902 36879.8403
254 -149946.2072 -27169.8902
255 55152.5245 -149946.2072
256 43485.6542 55152.5245
257 24396.6177 43485.6542
258 55152.5245 24396.6177
259 55152.5245 55152.5245
260 55152.5245 55152.5245
261 55152.5245 55152.5245
262 46603.3938 55152.5245
263 57799.4030 46603.3938
264 55152.5245 57799.4030
265 55152.5245 55152.5245
266 55152.5245 55152.5245
267 81168.2234 55152.5245
268 55152.5245 81168.2234
269 27644.5278 55152.5245
270 55152.5245 27644.5278
271 55152.5245 55152.5245
272 55152.5245 55152.5245
273 27922.6546 55152.5245
274 55152.5245 27922.6546
275 38752.0513 55152.5245
276 24162.5503 38752.0513
277 55152.5245 24162.5503
278 -37344.5929 55152.5245
279 26861.7964 -37344.5929
280 55152.5245 26861.7964
281 55152.5245 55152.5245
282 27804.4032 55152.5245
283 55152.5245 27804.4032
284 38853.4518 55152.5245
285 55152.5245 38853.4518
286 40704.7058 55152.5245
287 52103.9773 40704.7058
288 55152.5245 52103.9773
289 23957.2150 55152.5245
290 55152.5245 23957.2150
291 45877.7435 55152.5245
292 44577.6335 45877.7435
293 55152.5245 44577.6335
294 32825.8604 55152.5245
295 65372.1884 32825.8604
296 55152.5245 65372.1884
297 55152.5245 55152.5245
298 55152.5245 55152.5245
299 46068.8784 55152.5245
300 48116.5426 46068.8784
301 36569.8314 48116.5426
302 54143.2305 36569.8314
303 55152.5245 54143.2305
304 60303.2396 55152.5245
305 42987.4385 60303.2396
306 55152.5245 42987.4385
307 27472.0702 55152.5245
308 33465.2963 27472.0702
309 55152.5245 33465.2963
310 55152.5245 55152.5245
311 4530.3833 55152.5245
312 55152.5245 4530.3833
313 41282.3000 55152.5245
314 47223.8192 41282.3000
315 63986.5389 47223.8192
316 55180.1335 63986.5389
317 7924.8219 55180.1335
318 50077.6074 7924.8219
319 55152.5245 50077.6074
320 55152.5245 55152.5245
321 55152.5245 55152.5245
322 55152.5245 55152.5245
323 55152.5245 55152.5245
324 55152.5245 55152.5245
325 50777.7514 55152.5245
326 55152.5245 50777.7514
327 55152.5245 55152.5245
328 55152.5245 55152.5245
329 55152.5245 55152.5245
330 55152.5245 55152.5245
331 55152.5245 55152.5245
332 55152.5245 55152.5245
333 47885.6972 55152.5245
334 54118.2358 47885.6972
335 44193.4641 54118.2358
336 -3177.0170 44193.4641
337 55152.5245 -3177.0170
338 51794.7786 55152.5245
339 46172.1334 51794.7786
340 13020.8195 46172.1334
341 55152.5245 13020.8195
342 55152.5245 55152.5245
343 55152.5245 55152.5245
344 21236.0425 55152.5245
345 55152.5245 21236.0425
346 40314.4511 55152.5245
347 55152.5245 40314.4511
348 55152.5245 55152.5245
349 13933.7738 55152.5245
350 51894.0439 13933.7738
351 44801.8786 51894.0439
352 44219.9617 44801.8786
353 47891.5630 44219.9617
354 55152.5245 47891.5630
355 55152.5245 55152.5245
356 7368.8801 55152.5245
357 55152.5245 7368.8801
358 55152.5245 55152.5245
359 55152.5245 55152.5245
360 51108.1960 55152.5245
361 55152.5245 51108.1960
362 55152.5245 55152.5245
363 -20345.3694 55152.5245
364 44813.8486 -20345.3694
365 55152.5245 44813.8486
366 55152.5245 55152.5245
367 55152.5245 55152.5245
368 55152.5245 55152.5245
369 43461.0754 55152.5245
370 55152.5245 43461.0754
371 47678.2739 55152.5245
372 39372.9095 47678.2739
373 33907.2837 39372.9095
374 55152.5245 33907.2837
375 55152.5245 55152.5245
376 28500.3161 55152.5245
377 42129.4647 28500.3161
378 845.1762 42129.4647
379 55152.5245 845.1762
380 55152.5245 55152.5245
381 7288.2607 55152.5245
382 23651.2435 7288.2607
383 55152.5245 23651.2435
384 40168.2486 55152.5245
385 -117279.6312 40168.2486
386 62592.5952 -117279.6312
387 25728.1516 62592.5952
388 49874.5640 25728.1516
389 30391.7616 49874.5640
390 40145.1899 30391.7616
391 18133.7508 40145.1899
392 -19080.2469 18133.7508
393 74580.5292 -19080.2469
394 -96721.9437 74580.5292
395 6463.4051 -96721.9437
396 -138664.7201 6463.4051
397 -189519.7904 -138664.7201
398 12381.6147 -189519.7904
399 -55495.2651 12381.6147
400 2431.9986 -55495.2651
401 42613.8359 2431.9986
402 19583.9874 42613.8359
403 -49565.6943 19583.9874
404 -8995.7460 -49565.6943
405 14011.1662 -8995.7460
406 20618.4595 14011.1662
407 13250.1409 20618.4595
408 22906.3494 13250.1409
409 -61901.5139 22906.3494
410 -51361.8573 -61901.5139
411 -161769.0794 -51361.8573
412 -25907.2560 -161769.0794
413 60088.3931 -25907.2560
414 -177525.6236 60088.3931
415 -99195.2033 -177525.6236
416 91983.5806 -99195.2033
417 -17379.6380 91983.5806
418 -50409.1687 -17379.6380
419 -37404.0006 -50409.1687
420 -44265.0154 -37404.0006
421 -175393.8602 -44265.0154
422 -547870.2678 -175393.8602
423 -458560.1402 -547870.2678
424 -48844.9394 -458560.1402
425 -250014.5348 -48844.9394
426 -311858.7540 -250014.5348
427 -1135948.9209 -311858.7540
428 -272147.0098 -1135948.9209
429 -336049.2741 -272147.0098
430 -682483.1396 -336049.2741
> 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/7q81b1291211442.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/8q81b1291211442.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/9j0ie1291211442.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/10j0ie1291211442.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/114izj1291211442.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/12p1x71291211442.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/13wkcj1291211442.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/14pbt41291211442.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/15auas1291211442.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/1663811291211442.tab")
+ }
>
> try(system("convert tmp/1uzlk1291211442.ps tmp/1uzlk1291211442.png",intern=TRUE))
character(0)
> try(system("convert tmp/2m8k51291211442.ps tmp/2m8k51291211442.png",intern=TRUE))
character(0)
> try(system("convert tmp/3m8k51291211442.ps tmp/3m8k51291211442.png",intern=TRUE))
character(0)
> try(system("convert tmp/4m8k51291211442.ps tmp/4m8k51291211442.png",intern=TRUE))
character(0)
> try(system("convert tmp/5fzj81291211442.ps tmp/5fzj81291211442.png",intern=TRUE))
character(0)
> try(system("convert tmp/6fzj81291211442.ps tmp/6fzj81291211442.png",intern=TRUE))
character(0)
> try(system("convert tmp/7q81b1291211442.ps tmp/7q81b1291211442.png",intern=TRUE))
character(0)
> try(system("convert tmp/8q81b1291211442.ps tmp/8q81b1291211442.png",intern=TRUE))
character(0)
> try(system("convert tmp/9j0ie1291211442.ps tmp/9j0ie1291211442.png",intern=TRUE))
character(0)
> try(system("convert tmp/10j0ie1291211442.ps tmp/10j0ie1291211442.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
13.372 3.174 13.841