R version 2.12.0 (2010-10-15)
Copyright (C) 2010 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i486-pc-linux-gnu (32-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> x <- array(list(162556
+ ,807
+ ,131182692
+ ,6282154
+ ,29790
+ ,444
+ ,13226760
+ ,4321023
+ ,87550
+ ,412
+ ,36070600
+ ,4111912
+ ,84738
+ ,428
+ ,36267864
+ ,223193
+ ,54660
+ ,315
+ ,17217900
+ ,1491348
+ ,42634
+ ,168
+ ,7162512
+ ,1629616
+ ,40949
+ ,263
+ ,10769587
+ ,1398893
+ ,45187
+ ,267
+ ,12064929
+ ,1926517
+ ,37704
+ ,228
+ ,8596512
+ ,983660
+ ,16275
+ ,129
+ ,2099475
+ ,1443586
+ ,25830
+ ,104
+ ,2686320
+ ,1073089
+ ,12679
+ ,122
+ ,1546838
+ ,984885
+ ,18014
+ ,393
+ ,7079502
+ ,1405225
+ ,43556
+ ,190
+ ,8275640
+ ,227132
+ ,24811
+ ,280
+ ,6947080
+ ,929118
+ ,6575
+ ,63
+ ,414225
+ ,1071292
+ ,7123
+ ,102
+ ,726546
+ ,638830
+ ,21950
+ ,265
+ ,5816750
+ ,856956
+ ,37597
+ ,234
+ ,8797698
+ ,992426
+ ,17821
+ ,277
+ ,4936417
+ ,444477
+ ,12988
+ ,73
+ ,948124
+ ,857217
+ ,22330
+ ,67
+ ,1496110
+ ,711969
+ ,13326
+ ,103
+ ,1372578
+ ,702380
+ ,16189
+ ,290
+ ,4694810
+ ,358589
+ ,7146
+ ,83
+ ,593118
+ ,297978
+ ,15824
+ ,56
+ ,886144
+ ,585715
+ ,27664
+ ,236
+ ,6528704
+ ,657954
+ ,11920
+ ,73
+ ,870160
+ ,209458
+ ,8568
+ ,34
+ ,291312
+ ,786690
+ ,14416
+ ,139
+ ,2003824
+ ,439798
+ ,3369
+ ,26
+ ,87594
+ ,688779
+ ,11819
+ ,70
+ ,827330
+ ,574339
+ ,6984
+ ,40
+ ,279360
+ ,741409
+ ,4519
+ ,42
+ ,189798
+ ,597793
+ ,2220
+ ,12
+ ,26640
+ ,644190
+ ,18562
+ ,211
+ ,3916582
+ ,377934
+ ,10327
+ ,74
+ ,764198
+ ,640273
+ ,5336
+ ,80
+ ,426880
+ ,697458
+ ,2365
+ ,83
+ ,196295
+ ,550608
+ ,4069
+ ,131
+ ,533039
+ ,207393
+ ,8636
+ ,203
+ ,1753108
+ ,301607
+ ,13718
+ ,56
+ ,768208
+ ,345783
+ ,4525
+ ,89
+ ,402725
+ ,501749
+ ,6869
+ ,88
+ ,604472
+ ,379983
+ ,4628
+ ,39
+ ,180492
+ ,387475
+ ,3689
+ ,25
+ ,92225
+ ,377305
+ ,4891
+ ,49
+ ,239659
+ ,370837
+ ,7489
+ ,149
+ ,1115861
+ ,430866
+ ,4901
+ ,58
+ ,284258
+ ,469107
+ ,2284
+ ,41
+ ,93644
+ ,194493
+ ,3160
+ ,90
+ ,284400
+ ,530670
+ ,4150
+ ,136
+ ,564400
+ ,518365
+ ,7285
+ ,97
+ ,706645
+ ,491303
+ ,1134
+ ,63
+ ,71442
+ ,527021
+ ,4658
+ ,114
+ ,531012
+ ,233773
+ ,2384
+ ,77
+ ,183568
+ ,405972
+ ,3748
+ ,6
+ ,22488
+ ,652925
+ ,5371
+ ,47
+ ,252437
+ ,446211
+ ,1285
+ ,51
+ ,65535
+ ,341340
+ ,9327
+ ,85
+ ,792795
+ ,387699
+ ,5565
+ ,43
+ ,239295
+ ,493408
+ ,1528
+ ,32
+ ,48896
+ ,146494
+ ,3122
+ ,25
+ ,78050
+ ,414462
+ ,7561
+ ,77
+ ,582197
+ ,364304
+ ,2675
+ ,54
+ ,144450
+ ,355178
+ ,13253
+ ,251
+ ,3326503
+ ,357760
+ ,880
+ ,15
+ ,13200
+ ,261216
+ ,2053
+ ,44
+ ,90332
+ ,397144
+ ,1424
+ ,73
+ ,103952
+ ,374943
+ ,4036
+ ,85
+ ,343060
+ ,424898
+ ,3045
+ ,49
+ ,149205
+ ,202055
+ ,5119
+ ,38
+ ,194522
+ ,378525
+ ,1431
+ ,35
+ ,50085
+ ,310768
+ ,554
+ ,9
+ ,4986
+ ,325738
+ ,1975
+ ,34
+ ,67150
+ ,394510
+ ,1765
+ ,20
+ ,35300
+ ,247060
+ ,1012
+ ,29
+ ,29348
+ ,368078
+ ,810
+ ,11
+ ,8910
+ ,236761
+ ,1280
+ ,52
+ ,66560
+ ,312378
+ ,666
+ ,13
+ ,8658
+ ,339836
+ ,1380
+ ,29
+ ,40020
+ ,347385
+ ,4677
+ ,66
+ ,308682
+ ,426280
+ ,876
+ ,33
+ ,28908
+ ,352850
+ ,814
+ ,15
+ ,12210
+ ,301881
+ ,514
+ ,15
+ ,7710
+ ,377516
+ ,5692
+ ,68
+ ,387056
+ ,357312
+ ,3642
+ ,100
+ ,364200
+ ,458343
+ ,540
+ ,13
+ ,7020
+ ,354228
+ ,2099
+ ,45
+ ,94455
+ ,308636
+ ,567
+ ,14
+ ,7938
+ ,386212
+ ,2001
+ ,36
+ ,72036
+ ,393343
+ ,2949
+ ,40
+ ,117960
+ ,378509
+ ,2253
+ ,68
+ ,153204
+ ,452469
+ ,6533
+ ,29
+ ,189457
+ ,364839
+ ,1889
+ ,43
+ ,81227
+ ,358649
+ ,3055
+ ,30
+ ,91650
+ ,376641
+ ,272
+ ,9
+ ,2448
+ ,429112
+ ,1414
+ ,22
+ ,31108
+ ,330546
+ ,2564
+ ,19
+ ,48716
+ ,403560
+ ,1383
+ ,9
+ ,12447
+ ,317892
+ ,1261
+ ,31
+ ,39091
+ ,307528
+ ,975
+ ,19
+ ,18525
+ ,235133
+ ,3366
+ ,55
+ ,185130
+ ,299243
+ ,576
+ ,8
+ ,4608
+ ,314073
+ ,1686
+ ,28
+ ,47208
+ ,368186
+ ,746
+ ,29
+ ,21634
+ ,269661
+ ,3192
+ ,48
+ ,153216
+ ,125390
+ ,2045
+ ,16
+ ,32720
+ ,510834
+ ,5702
+ ,47
+ ,267994
+ ,321896
+ ,1932
+ ,20
+ ,38640
+ ,249898
+ ,936
+ ,22
+ ,20592
+ ,408881
+ ,3437
+ ,33
+ ,113421
+ ,158492
+ ,5131
+ ,44
+ ,225764
+ ,292154
+ ,2397
+ ,13
+ ,31161
+ ,289513
+ ,1389
+ ,6
+ ,8334
+ ,378049
+ ,1503
+ ,35
+ ,52605
+ ,343466
+ ,402
+ ,8
+ ,3216
+ ,332743
+ ,2239
+ ,17
+ ,38063
+ ,442882
+ ,2234
+ ,11
+ ,24574
+ ,214215
+ ,837
+ ,21
+ ,17577
+ ,315688
+ ,10579
+ ,92
+ ,973268
+ ,375195
+ ,875
+ ,12
+ ,10500
+ ,334280
+ ,1585
+ ,112
+ ,177520
+ ,355864
+ ,1659
+ ,25
+ ,41475
+ ,480382
+ ,2647
+ ,17
+ ,44999
+ ,353058
+ ,3294
+ ,23
+ ,75762
+ ,217193
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,94
+ ,10
+ ,940
+ ,314533
+ ,422
+ ,23
+ ,9706
+ ,318056
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,34
+ ,7
+ ,238
+ ,314353
+ ,1558
+ ,25
+ ,38950
+ ,369448
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,43
+ ,20
+ ,860
+ ,312846
+ ,645
+ ,4
+ ,2580
+ ,312075
+ ,316
+ ,4
+ ,1264
+ ,315009
+ ,115
+ ,10
+ ,1150
+ ,318903
+ ,5
+ ,1
+ ,5
+ ,314887
+ ,897
+ ,4
+ ,3588
+ ,314913
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,389
+ ,8
+ ,3112
+ ,325506
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,1002
+ ,11
+ ,11022
+ ,298568
+ ,36
+ ,4
+ ,144
+ ,315834
+ ,460
+ ,15
+ ,6900
+ ,329784
+ ,309
+ ,9
+ ,2781
+ ,312878
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,9
+ ,7
+ ,63
+ ,314987
+ ,271
+ ,2
+ ,542
+ ,325249
+ ,14
+ ,0
+ ,0
+ ,315877
+ ,520
+ ,7
+ ,3640
+ ,291650
+ ,1766
+ ,46
+ ,81236
+ ,305959
+ ,0
+ ,5
+ ,0
+ ,315380
+ ,458
+ ,7
+ ,3206
+ ,297765
+ ,20
+ ,2
+ ,40
+ ,315245
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,98
+ ,2
+ ,196
+ ,315236
+ ,405
+ ,5
+ ,2025
+ ,336425
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,483
+ ,7
+ ,3381
+ ,306268
+ ,454
+ ,24
+ ,10896
+ ,302187
+ ,47
+ ,1
+ ,47
+ ,314882
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,757
+ ,18
+ ,13626
+ ,382712
+ ,4655
+ ,55
+ ,256025
+ ,341570
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,36
+ ,3
+ ,108
+ ,312412
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,203
+ ,9
+ ,1827
+ ,309596
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,126
+ ,8
+ ,1008
+ ,315547
+ ,400
+ ,113
+ ,45200
+ ,313267
+ ,71
+ ,0
+ ,0
+ ,316176
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,972
+ ,19
+ ,18468
+ ,359335
+ ,531
+ ,11
+ ,5841
+ ,330068
+ ,2461
+ ,25
+ ,61525
+ ,314289
+ ,378
+ ,16
+ ,6048
+ ,297413
+ ,23
+ ,5
+ ,115
+ ,314806
+ ,638
+ ,11
+ ,7018
+ ,333210
+ ,2300
+ ,23
+ ,52900
+ ,352108
+ ,149
+ ,6
+ ,894
+ ,313332
+ ,226
+ ,5
+ ,1130
+ ,291787
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,275
+ ,7
+ ,1925
+ ,318745
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,141
+ ,7
+ ,987
+ ,315366
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,28
+ ,3
+ ,84
+ ,315688
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,4980
+ ,89
+ ,443220
+ ,409642
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,472
+ ,19
+ ,8968
+ ,269587
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,203
+ ,12
+ ,2436
+ ,300962
+ ,496
+ ,12
+ ,5952
+ ,325479
+ ,10
+ ,5
+ ,50
+ ,316155
+ ,63
+ ,2
+ ,126
+ ,318574
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,1136
+ ,26
+ ,29536
+ ,343613
+ ,265
+ ,3
+ ,795
+ ,306948
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,267
+ ,11
+ ,2937
+ ,330059
+ ,474
+ ,10
+ ,4740
+ ,288985
+ ,534
+ ,5
+ ,2670
+ ,304485
+ ,0
+ ,2
+ ,0
+ ,315380
+ ,15
+ ,6
+ ,90
+ ,315688
+ ,397
+ ,7
+ ,2779
+ ,317736
+ ,0
+ ,2
+ ,0
+ ,315380
+ ,1866
+ ,28
+ ,52248
+ ,322331
+ ,288
+ ,3
+ ,864
+ ,296656
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,3
+ ,1
+ ,3
+ ,315354
+ ,468
+ ,20
+ ,9360
+ ,312161
+ ,20
+ ,1
+ ,20
+ ,315576
+ ,278
+ ,22
+ ,6116
+ ,314922
+ ,61
+ ,9
+ ,549
+ ,314551
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,192
+ ,2
+ ,384
+ ,312339
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,317
+ ,7
+ ,2219
+ ,298700
+ ,738
+ ,9
+ ,6642
+ ,321376
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,368
+ ,13
+ ,4784
+ ,303230
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,2
+ ,0
+ ,0
+ ,315487
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,53
+ ,6
+ ,318
+ ,315793
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,94
+ ,3
+ ,282
+ ,312887
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,24
+ ,7
+ ,168
+ ,315637
+ ,2332
+ ,2
+ ,4664
+ ,324385
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,131
+ ,15
+ ,1965
+ ,308989
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,206
+ ,9
+ ,1854
+ ,296702
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,167
+ ,1
+ ,167
+ ,307322
+ ,622
+ ,38
+ ,23636
+ ,304376
+ ,2328
+ ,57
+ ,132696
+ ,253588
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,365
+ ,7
+ ,2555
+ ,309560
+ ,364
+ ,26
+ ,9464
+ ,298466
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,226
+ ,13
+ ,2938
+ ,343929
+ ,307
+ ,10
+ ,3070
+ ,331955
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,188
+ ,9
+ ,1692
+ ,381180
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,138
+ ,26
+ ,3588
+ ,331420
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,125
+ ,19
+ ,2375
+ ,310201
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,282
+ ,12
+ ,3384
+ ,320016
+ ,335
+ ,23
+ ,7705
+ ,320398
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,1324
+ ,29
+ ,38396
+ ,291841
+ ,176
+ ,8
+ ,1408
+ ,310670
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,249
+ ,26
+ ,6474
+ ,313491
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,333
+ ,9
+ ,2997
+ ,331323
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,601
+ ,5
+ ,3005
+ ,319210
+ ,30
+ ,3
+ ,90
+ ,318098
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,249
+ ,13
+ ,3237
+ ,292754
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,165
+ ,12
+ ,1980
+ ,325176
+ ,453
+ ,19
+ ,8607
+ ,365959
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,53
+ ,10
+ ,530
+ ,302409
+ ,382
+ ,9
+ ,3438
+ ,340968
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,9
+ ,0
+ ,315380
+ ,30
+ ,4
+ ,120
+ ,313164
+ ,290
+ ,1
+ ,290
+ ,301164
+ ,0
+ ,1
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,366
+ ,14
+ ,5124
+ ,344425
+ ,2
+ ,12
+ ,24
+ ,315394
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,209
+ ,19
+ ,3971
+ ,316647
+ ,384
+ ,17
+ ,6528
+ ,309836
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,365
+ ,32
+ ,11680
+ ,346611
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,49
+ ,14
+ ,686
+ ,322031
+ ,3
+ ,8
+ ,24
+ ,315656
+ ,133
+ ,4
+ ,532
+ ,339445
+ ,32
+ ,0
+ ,0
+ ,314964
+ ,368
+ ,20
+ ,7360
+ ,297141
+ ,1
+ ,5
+ ,5
+ ,315372
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,22
+ ,1
+ ,22
+ ,312502
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,96
+ ,4
+ ,384
+ ,313729
+ ,1
+ ,1
+ ,1
+ ,315388
+ ,314
+ ,4
+ ,1256
+ ,315371
+ ,844
+ ,20
+ ,16880
+ ,296139
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,26
+ ,1
+ ,26
+ ,313880
+ ,125
+ ,10
+ ,1250
+ ,317698
+ ,304
+ ,12
+ ,3648
+ ,295580
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,621
+ ,13
+ ,8073
+ ,308256
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,119
+ ,3
+ ,357
+ ,303677
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,1595
+ ,10
+ ,15950
+ ,319369
+ ,312
+ ,3
+ ,936
+ ,318690
+ ,60
+ ,7
+ ,420
+ ,314049
+ ,587
+ ,10
+ ,5870
+ ,325699
+ ,135
+ ,1
+ ,135
+ ,314210
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,514
+ ,15
+ ,7710
+ ,322378
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,1
+ ,4
+ ,4
+ ,315398
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,1763
+ ,28
+ ,49364
+ ,308336
+ ,180
+ ,9
+ ,1620
+ ,316386
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,218
+ ,7
+ ,1526
+ ,315553
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,448
+ ,7
+ ,3136
+ ,323361
+ ,227
+ ,7
+ ,1589
+ ,336639
+ ,174
+ ,3
+ ,522
+ ,307424
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,121
+ ,11
+ ,1331
+ ,295370
+ ,607
+ ,7
+ ,4249
+ ,322340
+ ,2212
+ ,10
+ ,22120
+ ,319864
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,530
+ ,18
+ ,9540
+ ,317291
+ ,571
+ ,14
+ ,7994
+ ,280398
+ ,0
+ ,0
+ ,0
+ ,315380
+ ,78
+ ,12
+ ,936
+ ,317330
+ ,2489
+ ,29
+ ,72181
+ ,238125
+ ,131
+ ,3
+ ,393
+ ,327071
+ ,923
+ ,6
+ ,5538
+ ,309038
+ ,72
+ ,3
+ ,216
+ ,314210
+ ,572
+ ,8
+ ,4576
+ ,307930
+ ,397
+ ,10
+ ,3970
+ ,322327
+ ,450
+ ,6
+ ,2700
+ ,292136
+ ,622
+ ,8
+ ,4976
+ ,263276
+ ,694
+ ,6
+ ,4164
+ ,367655
+ ,3425
+ ,9
+ ,30825
+ ,283910
+ ,562
+ ,8
+ ,4496
+ ,283587
+ ,4917
+ ,26
+ ,127842
+ ,243650
+ ,1442
+ ,239
+ ,344638
+ ,438493
+ ,529
+ ,7
+ ,3703
+ ,296261
+ ,2126
+ ,41
+ ,87166
+ ,230621
+ ,1061
+ ,3
+ ,3183
+ ,304252
+ ,776
+ ,8
+ ,6208
+ ,333505
+ ,611
+ ,6
+ ,3666
+ ,296919
+ ,1526
+ ,21
+ ,32046
+ ,278990
+ ,592
+ ,7
+ ,4144
+ ,276898
+ ,1182
+ ,11
+ ,13002
+ ,327007
+ ,621
+ ,11
+ ,6831
+ ,317046
+ ,989
+ ,12
+ ,11868
+ ,304555
+ ,438
+ ,9
+ ,3942
+ ,298096
+ ,726
+ ,3
+ ,2178
+ ,231861
+ ,1303
+ ,57
+ ,74271
+ ,309422
+ ,7419
+ ,21
+ ,155799
+ ,286963
+ ,1164
+ ,15
+ ,17460
+ ,269753
+ ,3310
+ ,32
+ ,105920
+ ,448243
+ ,1920
+ ,11
+ ,21120
+ ,165404
+ ,965
+ ,2
+ ,1930
+ ,204325
+ ,3256
+ ,23
+ ,74888
+ ,407159
+ ,1135
+ ,20
+ ,22700
+ ,290476
+ ,1270
+ ,24
+ ,30480
+ ,275311
+ ,661
+ ,1
+ ,661
+ ,246541
+ ,1013
+ ,1
+ ,1013
+ ,253468
+ ,2844
+ ,74
+ ,210456
+ ,240897
+ ,11528
+ ,68
+ ,783904
+ ,-83265
+ ,6526
+ ,20
+ ,130520
+ ,-42143
+ ,2264
+ ,20
+ ,45280
+ ,272713
+ ,5109
+ ,82
+ ,418938
+ ,215362
+ ,3999
+ ,21
+ ,83979
+ ,42754
+ ,35624
+ ,244
+ ,8692256
+ ,306275
+ ,9252
+ ,32
+ ,296064
+ ,253537
+ ,15236
+ ,86
+ ,1310296
+ ,372631
+ ,18073
+ ,69
+ ,1247037
+ ,-7170)
+ ,dim=c(4
+ ,431)
+ ,dimnames=list(c('Costs'
+ ,'Orders'
+ ,'C_O'
+ ,'Wealth
')
+ ,1:431))
> y <- array(NA,dim=c(4,431),dimnames=list(c('Costs','Orders','C_O','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
> n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test
> par1 <- as.numeric(par1)
> x <- t(y)
> k <- length(x[1,])
> n <- length(x[,1])
> x1 <- cbind(x[,par1], x[,1:k!=par1])
> mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1])
> colnames(x1) <- mycolnames #colnames(x)[par1]
> x <- x1
> if (par3 == 'First Differences'){
+ x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep='')))
+ for (i in 1:n-1) {
+ for (j in 1:k) {
+ x2[i,j] <- x[i+1,j] - x[i,j]
+ }
+ }
+ x <- x2
+ }
> if (par2 == 'Include Monthly Dummies'){
+ x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep ='')))
+ for (i in 1:11){
+ x2[seq(i,n,12),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> if (par2 == 'Include Quarterly Dummies'){
+ x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep ='')))
+ for (i in 1:3){
+ x2[seq(i,n,4),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> k <- length(x[1,])
> if (par3 == 'Linear Trend'){
+ x <- cbind(x, c(1:n))
+ colnames(x)[k+1] <- 't'
+ }
> x
Wealth\r Costs Orders C_O
1 6282154 162556 807 131182692
2 4321023 29790 444 13226760
3 4111912 87550 412 36070600
4 223193 84738 428 36267864
5 1491348 54660 315 17217900
6 1629616 42634 168 7162512
7 1398893 40949 263 10769587
8 1926517 45187 267 12064929
9 983660 37704 228 8596512
10 1443586 16275 129 2099475
11 1073089 25830 104 2686320
12 984885 12679 122 1546838
13 1405225 18014 393 7079502
14 227132 43556 190 8275640
15 929118 24811 280 6947080
16 1071292 6575 63 414225
17 638830 7123 102 726546
18 856956 21950 265 5816750
19 992426 37597 234 8797698
20 444477 17821 277 4936417
21 857217 12988 73 948124
22 711969 22330 67 1496110
23 702380 13326 103 1372578
24 358589 16189 290 4694810
25 297978 7146 83 593118
26 585715 15824 56 886144
27 657954 27664 236 6528704
28 209458 11920 73 870160
29 786690 8568 34 291312
30 439798 14416 139 2003824
31 688779 3369 26 87594
32 574339 11819 70 827330
33 741409 6984 40 279360
34 597793 4519 42 189798
35 644190 2220 12 26640
36 377934 18562 211 3916582
37 640273 10327 74 764198
38 697458 5336 80 426880
39 550608 2365 83 196295
40 207393 4069 131 533039
41 301607 8636 203 1753108
42 345783 13718 56 768208
43 501749 4525 89 402725
44 379983 6869 88 604472
45 387475 4628 39 180492
46 377305 3689 25 92225
47 370837 4891 49 239659
48 430866 7489 149 1115861
49 469107 4901 58 284258
50 194493 2284 41 93644
51 530670 3160 90 284400
52 518365 4150 136 564400
53 491303 7285 97 706645
54 527021 1134 63 71442
55 233773 4658 114 531012
56 405972 2384 77 183568
57 652925 3748 6 22488
58 446211 5371 47 252437
59 341340 1285 51 65535
60 387699 9327 85 792795
61 493408 5565 43 239295
62 146494 1528 32 48896
63 414462 3122 25 78050
64 364304 7561 77 582197
65 355178 2675 54 144450
66 357760 13253 251 3326503
67 261216 880 15 13200
68 397144 2053 44 90332
69 374943 1424 73 103952
70 424898 4036 85 343060
71 202055 3045 49 149205
72 378525 5119 38 194522
73 310768 1431 35 50085
74 325738 554 9 4986
75 394510 1975 34 67150
76 247060 1765 20 35300
77 368078 1012 29 29348
78 236761 810 11 8910
79 312378 1280 52 66560
80 339836 666 13 8658
81 347385 1380 29 40020
82 426280 4677 66 308682
83 352850 876 33 28908
84 301881 814 15 12210
85 377516 514 15 7710
86 357312 5692 68 387056
87 458343 3642 100 364200
88 354228 540 13 7020
89 308636 2099 45 94455
90 386212 567 14 7938
91 393343 2001 36 72036
92 378509 2949 40 117960
93 452469 2253 68 153204
94 364839 6533 29 189457
95 358649 1889 43 81227
96 376641 3055 30 91650
97 429112 272 9 2448
98 330546 1414 22 31108
99 403560 2564 19 48716
100 317892 1383 9 12447
101 307528 1261 31 39091
102 235133 975 19 18525
103 299243 3366 55 185130
104 314073 576 8 4608
105 368186 1686 28 47208
106 269661 746 29 21634
107 125390 3192 48 153216
108 510834 2045 16 32720
109 321896 5702 47 267994
110 249898 1932 20 38640
111 408881 936 22 20592
112 158492 3437 33 113421
113 292154 5131 44 225764
114 289513 2397 13 31161
115 378049 1389 6 8334
116 343466 1503 35 52605
117 332743 402 8 3216
118 442882 2239 17 38063
119 214215 2234 11 24574
120 315688 837 21 17577
121 375195 10579 92 973268
122 334280 875 12 10500
123 355864 1585 112 177520
124 480382 1659 25 41475
125 353058 2647 17 44999
126 217193 3294 23 75762
127 315380 0 0 0
128 314533 94 10 940
129 318056 422 23 9706
130 315380 0 0 0
131 314353 34 7 238
132 369448 1558 25 38950
133 315380 0 0 0
134 312846 43 20 860
135 312075 645 4 2580
136 315009 316 4 1264
137 318903 115 10 1150
138 314887 5 1 5
139 314913 897 4 3588
140 315380 0 0 0
141 325506 389 8 3112
142 315380 0 0 0
143 298568 1002 11 11022
144 315834 36 4 144
145 329784 460 15 6900
146 312878 309 9 2781
147 315380 0 0 0
148 314987 9 7 63
149 325249 271 2 542
150 315877 14 0 0
151 291650 520 7 3640
152 305959 1766 46 81236
153 315380 0 5 0
154 297765 458 7 3206
155 315245 20 2 40
156 315380 0 0 0
157 315380 0 0 0
158 315236 98 2 196
159 336425 405 5 2025
160 315380 0 0 0
161 315380 0 0 0
162 315380 0 0 0
163 315380 0 0 0
164 306268 483 7 3381
165 302187 454 24 10896
166 314882 47 1 47
167 315380 0 0 0
168 382712 757 18 13626
169 341570 4655 55 256025
170 315380 0 0 0
171 315380 0 0 0
172 312412 36 3 108
173 315380 0 0 0
174 309596 203 9 1827
175 315380 0 0 0
176 315547 126 8 1008
177 313267 400 113 45200
178 316176 71 0 0
179 315380 0 0 0
180 315380 0 0 0
181 359335 972 19 18468
182 330068 531 11 5841
183 314289 2461 25 61525
184 297413 378 16 6048
185 314806 23 5 115
186 333210 638 11 7018
187 352108 2300 23 52900
188 313332 149 6 894
189 291787 226 5 1130
190 315380 0 0 0
191 318745 275 7 1925
192 315380 0 0 0
193 315366 141 7 987
194 315380 0 0 0
195 315688 28 3 84
196 315380 0 0 0
197 409642 4980 89 443220
198 315380 0 0 0
199 315380 0 0 0
200 269587 472 19 8968
201 315380 0 0 0
202 315380 0 0 0
203 315380 0 0 0
204 300962 203 12 2436
205 325479 496 12 5952
206 316155 10 5 50
207 318574 63 2 126
208 315380 0 0 0
209 343613 1136 26 29536
210 306948 265 3 795
211 315380 0 0 0
212 315380 0 0 0
213 330059 267 11 2937
214 288985 474 10 4740
215 304485 534 5 2670
216 315380 0 2 0
217 315688 15 6 90
218 317736 397 7 2779
219 315380 0 2 0
220 322331 1866 28 52248
221 296656 288 3 864
222 315380 0 0 0
223 315354 3 1 3
224 312161 468 20 9360
225 315576 20 1 20
226 314922 278 22 6116
227 314551 61 9 549
228 315380 0 0 0
229 312339 192 2 384
230 315380 0 0 0
231 298700 317 7 2219
232 321376 738 9 6642
233 315380 0 0 0
234 303230 368 13 4784
235 315380 0 0 0
236 315487 2 0 0
237 315380 0 0 0
238 315793 53 6 318
239 315380 0 0 0
240 315380 0 0 0
241 315380 0 0 0
242 312887 94 3 282
243 315380 0 0 0
244 315637 24 7 168
245 324385 2332 2 4664
246 315380 0 0 0
247 315380 0 0 0
248 308989 131 15 1965
249 315380 0 0 0
250 315380 0 0 0
251 296702 206 9 1854
252 315380 0 0 0
253 307322 167 1 167
254 304376 622 38 23636
255 253588 2328 57 132696
256 315380 0 0 0
257 309560 365 7 2555
258 298466 364 26 9464
259 315380 0 0 0
260 315380 0 0 0
261 315380 0 0 0
262 315380 0 0 0
263 343929 226 13 2938
264 331955 307 10 3070
265 315380 0 0 0
266 315380 0 0 0
267 315380 0 0 0
268 381180 188 9 1692
269 315380 0 0 0
270 331420 138 26 3588
271 315380 0 0 0
272 315380 0 0 0
273 315380 0 0 0
274 310201 125 19 2375
275 315380 0 0 0
276 320016 282 12 3384
277 320398 335 23 7705
278 315380 0 0 0
279 291841 1324 29 38396
280 310670 176 8 1408
281 315380 0 0 0
282 315380 0 0 0
283 313491 249 26 6474
284 315380 0 0 0
285 331323 333 9 2997
286 315380 0 0 0
287 319210 601 5 3005
288 318098 30 3 90
289 315380 0 0 0
290 292754 249 13 3237
291 315380 0 0 0
292 325176 165 12 1980
293 365959 453 19 8607
294 315380 0 0 0
295 302409 53 10 530
296 340968 382 9 3438
297 315380 0 0 0
298 315380 0 0 0
299 315380 0 0 0
300 315380 0 9 0
301 313164 30 4 120
302 301164 290 1 290
303 315380 0 1 0
304 315380 0 0 0
305 344425 366 14 5124
306 315394 2 12 24
307 315380 0 0 0
308 316647 209 19 3971
309 309836 384 17 6528
310 315380 0 0 0
311 315380 0 0 0
312 346611 365 32 11680
313 315380 0 0 0
314 322031 49 14 686
315 315656 3 8 24
316 339445 133 4 532
317 314964 32 0 0
318 297141 368 20 7360
319 315372 1 5 5
320 315380 0 0 0
321 315380 0 0 0
322 315380 0 0 0
323 315380 0 0 0
324 315380 0 0 0
325 315380 0 0 0
326 312502 22 1 22
327 315380 0 0 0
328 315380 0 0 0
329 315380 0 0 0
330 315380 0 0 0
331 315380 0 0 0
332 315380 0 0 0
333 315380 0 0 0
334 313729 96 4 384
335 315388 1 1 1
336 315371 314 4 1256
337 296139 844 20 16880
338 315380 0 0 0
339 313880 26 1 26
340 317698 125 10 1250
341 295580 304 12 3648
342 315380 0 0 0
343 315380 0 0 0
344 315380 0 0 0
345 308256 621 13 8073
346 315380 0 0 0
347 303677 119 3 357
348 315380 0 0 0
349 315380 0 0 0
350 319369 1595 10 15950
351 318690 312 3 936
352 314049 60 7 420
353 325699 587 10 5870
354 314210 135 1 135
355 315380 0 0 0
356 315380 0 0 0
357 322378 514 15 7710
358 315380 0 0 0
359 315380 0 0 0
360 315380 0 0 0
361 315398 1 4 4
362 315380 0 0 0
363 315380 0 0 0
364 308336 1763 28 49364
365 316386 180 9 1620
366 315380 0 0 0
367 315380 0 0 0
368 315380 0 0 0
369 315380 0 0 0
370 315553 218 7 1526
371 315380 0 0 0
372 323361 448 7 3136
373 336639 227 7 1589
374 307424 174 3 522
375 315380 0 0 0
376 315380 0 0 0
377 295370 121 11 1331
378 322340 607 7 4249
379 319864 2212 10 22120
380 315380 0 0 0
381 315380 0 0 0
382 317291 530 18 9540
383 280398 571 14 7994
384 315380 0 0 0
385 317330 78 12 936
386 238125 2489 29 72181
387 327071 131 3 393
388 309038 923 6 5538
389 314210 72 3 216
390 307930 572 8 4576
391 322327 397 10 3970
392 292136 450 6 2700
393 263276 622 8 4976
394 367655 694 6 4164
395 283910 3425 9 30825
396 283587 562 8 4496
397 243650 4917 26 127842
398 438493 1442 239 344638
399 296261 529 7 3703
400 230621 2126 41 87166
401 304252 1061 3 3183
402 333505 776 8 6208
403 296919 611 6 3666
404 278990 1526 21 32046
405 276898 592 7 4144
406 327007 1182 11 13002
407 317046 621 11 6831
408 304555 989 12 11868
409 298096 438 9 3942
410 231861 726 3 2178
411 309422 1303 57 74271
412 286963 7419 21 155799
413 269753 1164 15 17460
414 448243 3310 32 105920
415 165404 1920 11 21120
416 204325 965 2 1930
417 407159 3256 23 74888
418 290476 1135 20 22700
419 275311 1270 24 30480
420 246541 661 1 661
421 253468 1013 1 1013
422 240897 2844 74 210456
423 -83265 11528 68 783904
424 -42143 6526 20 130520
425 272713 2264 20 45280
426 215362 5109 82 418938
427 42754 3999 21 83979
428 306275 35624 244 8692256
429 253537 9252 32 296064
430 372631 15236 86 1310296
431 -7170 18073 69 1247037
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) Costs Orders C_O
2.943e+05 1.814e+00 1.799e+03 3.285e-02
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-2186068 -19503 14165 21036 2739587
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 2.943e+05 1.302e+04 22.612 < 2e-16 ***
Costs 1.814e+00 2.986e+00 0.608 0.544
Orders 1.799e+03 3.233e+02 5.563 4.69e-08 ***
C_O 3.285e-02 3.629e-03 9.051 < 2e-16 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 236100 on 427 degrees of freedom
Multiple R-squared: 0.7012, Adjusted R-squared: 0.6991
F-statistic: 334.1 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.986130e-79 1.493065e-79
[2,] 1.0000000 4.462198e-85 2.231099e-85
[3,] 1.0000000 2.291864e-87 1.145932e-87
[4,] 1.0000000 4.578807e-95 2.289404e-95
[5,] 1.0000000 6.767462e-98 3.383731e-98
[6,] 1.0000000 6.080998e-102 3.040499e-102
[7,] 1.0000000 1.316877e-121 6.584386e-122
[8,] 1.0000000 3.789018e-128 1.894509e-128
[9,] 1.0000000 7.189170e-133 3.594585e-133
[10,] 1.0000000 3.310738e-141 1.655369e-141
[11,] 1.0000000 2.984591e-142 1.492296e-142
[12,] 1.0000000 1.793196e-146 8.965980e-147
[13,] 1.0000000 7.622464e-151 3.811232e-151
[14,] 1.0000000 9.851845e-156 4.925922e-156
[15,] 1.0000000 1.269400e-159 6.347002e-160
[16,] 1.0000000 1.821285e-160 9.106426e-161
[17,] 1.0000000 4.134923e-162 2.067461e-162
[18,] 1.0000000 1.190540e-166 5.952701e-167
[19,] 1.0000000 1.702515e-166 8.512577e-167
[20,] 1.0000000 9.288106e-167 4.644053e-167
[21,] 1.0000000 2.134658e-169 1.067329e-169
[22,] 1.0000000 9.227937e-170 4.613968e-170
[23,] 1.0000000 1.960419e-175 9.802097e-176
[24,] 1.0000000 3.971279e-175 1.985640e-175
[25,] 1.0000000 9.371045e-179 4.685523e-179
[26,] 1.0000000 8.593132e-180 4.296566e-180
[27,] 1.0000000 7.196158e-186 3.598079e-186
[28,] 1.0000000 6.571173e-188 3.285587e-188
[29,] 1.0000000 7.834955e-192 3.917477e-192
[30,] 1.0000000 1.207013e-192 6.035065e-193
[31,] 1.0000000 8.506888e-197 4.253444e-197
[32,] 1.0000000 2.868729e-202 1.434365e-202
[33,] 1.0000000 3.404503e-203 1.702251e-203
[34,] 1.0000000 3.058060e-205 1.529030e-205
[35,] 1.0000000 1.198440e-206 5.992201e-207
[36,] 1.0000000 2.798552e-206 1.399276e-206
[37,] 1.0000000 1.772512e-206 8.862559e-207
[38,] 1.0000000 9.970932e-206 4.985466e-206
[39,] 1.0000000 4.437559e-205 2.218779e-205
[40,] 1.0000000 2.269889e-204 1.134944e-204
[41,] 1.0000000 1.417048e-203 7.085238e-204
[42,] 1.0000000 5.109524e-203 2.554762e-203
[43,] 1.0000000 2.432195e-203 1.216098e-203
[44,] 1.0000000 1.435655e-203 7.178273e-204
[45,] 1.0000000 2.865201e-204 1.432600e-204
[46,] 1.0000000 2.185618e-204 1.092809e-204
[47,] 1.0000000 3.435965e-205 1.717982e-205
[48,] 1.0000000 3.153286e-206 1.576643e-206
[49,] 1.0000000 1.200474e-206 6.002372e-207
[50,] 1.0000000 7.319867e-206 3.659934e-206
[51,] 1.0000000 3.918681e-215 1.959340e-215
[52,] 1.0000000 5.929913e-216 2.964956e-216
[53,] 1.0000000 6.400927e-215 3.200464e-215
[54,] 1.0000000 6.892235e-215 3.446117e-215
[55,] 1.0000000 1.343014e-217 6.715068e-218
[56,] 1.0000000 3.188109e-219 1.594055e-219
[57,] 1.0000000 1.927854e-219 9.639271e-220
[58,] 1.0000000 4.648559e-219 2.324280e-219
[59,] 1.0000000 4.066386e-218 2.033193e-218
[60,] 1.0000000 6.231747e-219 3.115873e-219
[61,] 1.0000000 4.011241e-218 2.005621e-218
[62,] 1.0000000 1.507455e-217 7.537273e-218
[63,] 1.0000000 1.523518e-216 7.617588e-217
[64,] 1.0000000 2.787387e-216 1.393694e-216
[65,] 1.0000000 2.051241e-216 1.025620e-216
[66,] 1.0000000 2.163337e-216 1.081668e-216
[67,] 1.0000000 2.378781e-215 1.189390e-215
[68,] 1.0000000 2.680152e-214 1.340076e-214
[69,] 1.0000000 7.538850e-214 3.769425e-214
[70,] 1.0000000 4.115888e-213 2.057944e-213
[71,] 1.0000000 3.026770e-212 1.513385e-212
[72,] 1.0000000 1.016613e-211 5.083064e-212
[73,] 1.0000000 9.380857e-211 4.690428e-211
[74,] 1.0000000 9.218530e-210 4.609265e-210
[75,] 1.0000000 8.397379e-209 4.198689e-209
[76,] 1.0000000 3.488472e-209 1.744236e-209
[77,] 1.0000000 3.307912e-208 1.653956e-208
[78,] 1.0000000 3.492358e-207 1.746179e-207
[79,] 1.0000000 1.790205e-206 8.951024e-207
[80,] 1.0000000 6.220207e-206 3.110103e-206
[81,] 1.0000000 3.026237e-206 1.513118e-206
[82,] 1.0000000 2.356200e-205 1.178100e-205
[83,] 1.0000000 2.369083e-204 1.184542e-204
[84,] 1.0000000 8.479372e-204 4.239686e-204
[85,] 1.0000000 1.883330e-203 9.416648e-204
[86,] 1.0000000 4.612910e-203 2.306455e-203
[87,] 1.0000000 1.305544e-203 6.527718e-204
[88,] 1.0000000 2.313255e-204 1.156627e-204
[89,] 1.0000000 1.442250e-203 7.211252e-204
[90,] 1.0000000 2.050897e-203 1.025448e-203
[91,] 1.0000000 1.031310e-203 5.156550e-204
[92,] 1.0000000 9.240375e-203 4.620188e-203
[93,] 1.0000000 2.988342e-203 1.494171e-203
[94,] 1.0000000 2.816072e-202 1.408036e-202
[95,] 1.0000000 2.879394e-201 1.439697e-201
[96,] 1.0000000 7.668690e-201 3.834345e-201
[97,] 1.0000000 7.106003e-200 3.553001e-200
[98,] 1.0000000 7.317150e-199 3.658575e-199
[99,] 1.0000000 2.705310e-198 1.352655e-198
[100,] 1.0000000 1.496459e-197 7.482293e-198
[101,] 1.0000000 1.007305e-199 5.036526e-200
[102,] 1.0000000 2.703604e-204 1.351802e-204
[103,] 1.0000000 7.188054e-204 3.594027e-204
[104,] 1.0000000 4.479580e-203 2.239790e-203
[105,] 1.0000000 3.491234e-203 1.745617e-203
[106,] 1.0000000 5.852013e-204 2.926006e-204
[107,] 1.0000000 4.222707e-203 2.111353e-203
[108,] 1.0000000 4.384562e-202 2.192281e-202
[109,] 1.0000000 5.952982e-202 2.976491e-202
[110,] 1.0000000 4.603021e-201 2.301511e-201
[111,] 1.0000000 4.312343e-200 2.156171e-200
[112,] 1.0000000 3.031507e-202 1.515753e-202
[113,] 1.0000000 8.597236e-202 4.298618e-202
[114,] 1.0000000 9.174526e-201 4.587263e-201
[115,] 1.0000000 5.885963e-203 2.942981e-203
[116,] 1.0000000 4.863894e-202 2.431947e-202
[117,] 1.0000000 4.330530e-201 2.165265e-201
[118,] 1.0000000 9.145711e-205 4.572855e-205
[119,] 1.0000000 1.136468e-204 5.682340e-205
[120,] 1.0000000 5.551988e-204 2.775994e-204
[121,] 1.0000000 6.115386e-203 3.057693e-203
[122,] 1.0000000 6.711203e-202 3.355602e-202
[123,] 1.0000000 7.350718e-201 3.675359e-201
[124,] 1.0000000 7.961281e-200 3.980640e-200
[125,] 1.0000000 8.582814e-199 4.291407e-199
[126,] 1.0000000 1.695524e-198 8.477620e-199
[127,] 1.0000000 1.817175e-197 9.085876e-198
[128,] 1.0000000 1.870359e-196 9.351793e-197
[129,] 1.0000000 1.946685e-195 9.733425e-196
[130,] 1.0000000 2.045662e-194 1.022831e-194
[131,] 1.0000000 2.154833e-193 1.077416e-193
[132,] 1.0000000 2.250069e-192 1.125034e-192
[133,] 1.0000000 2.148018e-191 1.074009e-191
[134,] 1.0000000 2.217806e-190 1.108903e-190
[135,] 1.0000000 2.132703e-189 1.066351e-189
[136,] 1.0000000 2.179410e-188 1.089705e-188
[137,] 1.0000000 2.222206e-187 1.111103e-187
[138,] 1.0000000 2.252437e-186 1.126219e-186
[139,] 1.0000000 2.083101e-185 1.041550e-185
[140,] 1.0000000 2.094921e-184 1.047461e-184
[141,] 1.0000000 2.085630e-183 1.042815e-183
[142,] 1.0000000 2.067551e-182 1.033776e-182
[143,] 1.0000000 1.894099e-181 9.470493e-182
[144,] 1.0000000 1.856067e-180 9.280335e-181
[145,] 1.0000000 1.678252e-179 8.391258e-180
[146,] 1.0000000 1.567856e-178 7.839280e-179
[147,] 1.0000000 1.516768e-177 7.583839e-178
[148,] 1.0000000 1.408460e-176 7.042298e-177
[149,] 1.0000000 1.348274e-175 6.741372e-176
[150,] 1.0000000 1.282397e-174 6.411987e-175
[151,] 1.0000000 1.213640e-173 6.068200e-174
[152,] 1.0000000 1.143735e-172 5.718674e-173
[153,] 1.0000000 8.505684e-172 4.252842e-172
[154,] 1.0000000 7.936772e-171 3.968386e-171
[155,] 1.0000000 7.369179e-170 3.684589e-170
[156,] 1.0000000 6.808308e-169 3.404154e-169
[157,] 1.0000000 6.259049e-168 3.129524e-168
[158,] 1.0000000 5.740002e-167 2.870001e-167
[159,] 1.0000000 4.835615e-166 2.417807e-166
[160,] 1.0000000 4.383781e-165 2.191890e-165
[161,] 1.0000000 3.951705e-164 1.975853e-164
[162,] 1.0000000 7.005321e-164 3.502661e-164
[163,] 1.0000000 9.194036e-164 4.597018e-164
[164,] 1.0000000 8.281807e-163 4.140903e-163
[165,] 1.0000000 7.423851e-162 3.711925e-162
[166,] 1.0000000 6.621165e-161 3.310582e-161
[167,] 1.0000000 5.877762e-160 2.938881e-160
[168,] 1.0000000 5.171270e-159 2.585635e-159
[169,] 1.0000000 4.546255e-158 2.273128e-158
[170,] 1.0000000 3.989373e-157 1.994686e-157
[171,] 1.0000000 1.258968e-156 6.294840e-157
[172,] 1.0000000 1.091832e-155 5.459162e-156
[173,] 1.0000000 9.451799e-155 4.725899e-155
[174,] 1.0000000 8.143100e-154 4.071550e-154
[175,] 1.0000000 2.971764e-153 1.485882e-153
[176,] 1.0000000 2.244949e-152 1.122475e-152
[177,] 1.0000000 1.328886e-151 6.644429e-152
[178,] 1.0000000 1.035106e-150 5.175531e-151
[179,] 1.0000000 8.772781e-150 4.386390e-150
[180,] 1.0000000 6.052441e-149 3.026221e-149
[181,] 1.0000000 1.160190e-148 5.800949e-149
[182,] 1.0000000 9.772404e-148 4.886202e-148
[183,] 1.0000000 7.314336e-147 3.657168e-147
[184,] 1.0000000 6.087254e-146 3.043627e-146
[185,] 1.0000000 4.986315e-145 2.493158e-145
[186,] 1.0000000 4.110760e-144 2.055380e-144
[187,] 1.0000000 3.382190e-143 1.691095e-143
[188,] 1.0000000 2.762007e-142 1.381003e-142
[189,] 1.0000000 2.248398e-141 1.124199e-141
[190,] 1.0000000 1.818864e-140 9.094320e-141
[191,] 1.0000000 1.353912e-141 6.769562e-142
[192,] 1.0000000 1.108295e-140 5.541474e-141
[193,] 1.0000000 9.029825e-140 4.514912e-140
[194,] 1.0000000 4.233923e-139 2.116961e-139
[195,] 1.0000000 3.425685e-138 1.712842e-138
[196,] 1.0000000 2.758770e-137 1.379385e-137
[197,] 1.0000000 2.211303e-136 1.105651e-136
[198,] 1.0000000 1.666528e-135 8.332639e-136
[199,] 1.0000000 1.232781e-134 6.163904e-135
[200,] 1.0000000 9.765283e-134 4.882641e-134
[201,] 1.0000000 7.644545e-133 3.822273e-133
[202,] 1.0000000 5.988747e-132 2.994374e-132
[203,] 1.0000000 2.961893e-131 1.480946e-131
[204,] 1.0000000 2.303442e-130 1.151721e-130
[205,] 1.0000000 1.783031e-129 8.915157e-130
[206,] 1.0000000 1.373758e-128 6.868788e-129
[207,] 1.0000000 9.781531e-128 4.890765e-128
[208,] 1.0000000 6.737143e-127 3.368571e-127
[209,] 1.0000000 5.134625e-126 2.567312e-126
[210,] 1.0000000 3.889340e-125 1.944670e-125
[211,] 1.0000000 2.934482e-124 1.467241e-124
[212,] 1.0000000 2.158640e-123 1.079320e-123
[213,] 1.0000000 1.612373e-122 8.061863e-123
[214,] 1.0000000 8.934440e-122 4.467220e-122
[215,] 1.0000000 6.358999e-121 3.179500e-121
[216,] 1.0000000 4.686473e-120 2.343236e-120
[217,] 1.0000000 3.439496e-119 1.719748e-119
[218,] 1.0000000 2.510906e-118 1.255453e-118
[219,] 1.0000000 1.824943e-117 9.124716e-118
[220,] 1.0000000 1.317299e-116 6.586496e-117
[221,] 1.0000000 9.492152e-116 4.746076e-116
[222,] 1.0000000 6.800895e-115 3.400447e-115
[223,] 1.0000000 4.854533e-114 2.427266e-114
[224,] 1.0000000 3.445589e-113 1.722794e-113
[225,] 1.0000000 2.356110e-112 1.178055e-112
[226,] 1.0000000 1.518527e-111 7.592635e-112
[227,] 1.0000000 1.063306e-110 5.316531e-111
[228,] 1.0000000 7.279884e-110 3.639942e-110
[229,] 1.0000000 5.050107e-109 2.525053e-109
[230,] 1.0000000 3.486456e-108 1.743228e-108
[231,] 1.0000000 2.395845e-107 1.197922e-107
[232,] 1.0000000 1.641330e-106 8.206648e-107
[233,] 1.0000000 1.117254e-105 5.586269e-106
[234,] 1.0000000 7.569187e-105 3.784593e-105
[235,] 1.0000000 5.103704e-104 2.551852e-104
[236,] 1.0000000 3.429373e-103 1.714686e-103
[237,] 1.0000000 2.290427e-102 1.145214e-102
[238,] 1.0000000 1.524643e-101 7.623213e-102
[239,] 1.0000000 4.330007e-101 2.165004e-101
[240,] 1.0000000 2.864148e-100 1.432074e-100
[241,] 1.0000000 1.885491e-99 9.427455e-100
[242,] 1.0000000 1.215567e-98 6.077833e-99
[243,] 1.0000000 7.926572e-98 3.963286e-98
[244,] 1.0000000 5.144047e-97 2.572024e-97
[245,] 1.0000000 3.121452e-96 1.560726e-96
[246,] 1.0000000 2.007089e-95 1.003545e-95
[247,] 1.0000000 1.282687e-94 6.413433e-95
[248,] 1.0000000 7.701568e-94 3.850784e-94
[249,] 1.0000000 2.862040e-93 1.431020e-93
[250,] 1.0000000 1.811293e-92 9.056463e-93
[251,] 1.0000000 1.142812e-91 5.714061e-92
[252,] 1.0000000 6.506173e-91 3.253087e-91
[253,] 1.0000000 4.060129e-90 2.030065e-90
[254,] 1.0000000 2.521336e-89 1.260668e-89
[255,] 1.0000000 1.558094e-88 7.790470e-89
[256,] 1.0000000 9.581291e-88 4.790645e-88
[257,] 1.0000000 4.937649e-87 2.468825e-87
[258,] 1.0000000 2.771753e-86 1.385876e-86
[259,] 1.0000000 1.681873e-85 8.409367e-86
[260,] 1.0000000 1.015500e-84 5.077498e-85
[261,] 1.0000000 6.101105e-84 3.050553e-84
[262,] 1.0000000 1.416014e-83 7.080069e-84
[263,] 1.0000000 8.473021e-83 4.236510e-83
[264,] 1.0000000 4.970925e-82 2.485463e-82
[265,] 1.0000000 2.944936e-81 1.472468e-81
[266,] 1.0000000 1.735914e-80 8.679570e-81
[267,] 1.0000000 1.018095e-79 5.090475e-80
[268,] 1.0000000 5.841676e-79 2.920838e-79
[269,] 1.0000000 3.392031e-78 1.696016e-78
[270,] 1.0000000 1.943950e-77 9.719751e-78
[271,] 1.0000000 1.115525e-76 5.577624e-77
[272,] 1.0000000 6.379004e-76 3.189502e-76
[273,] 1.0000000 3.536915e-75 1.768457e-75
[274,] 1.0000000 2.003613e-74 1.001807e-74
[275,] 1.0000000 1.128393e-73 5.641966e-74
[276,] 1.0000000 6.322021e-73 3.161010e-73
[277,] 1.0000000 3.487873e-72 1.743937e-72
[278,] 1.0000000 1.934068e-71 9.670340e-72
[279,] 1.0000000 9.845030e-71 4.922515e-71
[280,] 1.0000000 5.404826e-70 2.702413e-70
[281,] 1.0000000 2.807663e-69 1.403831e-69
[282,] 1.0000000 1.522402e-68 7.612008e-69
[283,] 1.0000000 8.228484e-68 4.114242e-68
[284,] 1.0000000 4.086438e-67 2.043219e-67
[285,] 1.0000000 2.186839e-66 1.093419e-66
[286,] 1.0000000 1.145877e-65 5.729387e-66
[287,] 1.0000000 3.669879e-65 1.834940e-65
[288,] 1.0000000 1.939583e-64 9.697914e-65
[289,] 1.0000000 9.887100e-64 4.943550e-64
[290,] 1.0000000 4.322249e-63 2.161124e-63
[291,] 1.0000000 2.250807e-62 1.125404e-62
[292,] 1.0000000 1.165672e-61 5.828359e-62
[293,] 1.0000000 6.003624e-61 3.001812e-61
[294,] 1.0000000 3.081624e-60 1.540812e-60
[295,] 1.0000000 1.572979e-59 7.864895e-60
[296,] 1.0000000 7.948673e-59 3.974337e-59
[297,] 1.0000000 4.006110e-58 2.003055e-58
[298,] 1.0000000 2.006532e-57 1.003266e-57
[299,] 1.0000000 8.376554e-57 4.188277e-57
[300,] 1.0000000 4.160463e-56 2.080231e-56
[301,] 1.0000000 2.051138e-55 1.025569e-55
[302,] 1.0000000 1.007133e-54 5.035663e-55
[303,] 1.0000000 4.911339e-54 2.455669e-54
[304,] 1.0000000 2.379701e-53 1.189851e-53
[305,] 1.0000000 1.146312e-52 5.731560e-53
[306,] 1.0000000 4.921522e-52 2.460761e-52
[307,] 1.0000000 2.344870e-51 1.172435e-51
[308,] 1.0000000 1.109473e-50 5.547366e-51
[309,] 1.0000000 5.237164e-50 2.618582e-50
[310,] 1.0000000 2.171328e-49 1.085664e-49
[311,] 1.0000000 1.010581e-48 5.052903e-49
[312,] 1.0000000 4.507939e-48 2.253970e-48
[313,] 1.0000000 2.079157e-47 1.039578e-47
[314,] 1.0000000 9.508772e-47 4.754386e-47
[315,] 1.0000000 4.321863e-46 2.160931e-46
[316,] 1.0000000 1.952136e-45 9.760681e-46
[317,] 1.0000000 8.762435e-45 4.381218e-45
[318,] 1.0000000 3.908387e-44 1.954193e-44
[319,] 1.0000000 1.732249e-43 8.661247e-44
[320,] 1.0000000 7.646822e-43 3.823411e-43
[321,] 1.0000000 3.345998e-42 1.672999e-42
[322,] 1.0000000 1.454637e-41 7.273184e-42
[323,] 1.0000000 6.282726e-41 3.141363e-41
[324,] 1.0000000 2.695792e-40 1.347896e-40
[325,] 1.0000000 1.149076e-39 5.745382e-40
[326,] 1.0000000 4.865357e-39 2.432678e-39
[327,] 1.0000000 2.046267e-38 1.023134e-38
[328,] 1.0000000 8.563538e-38 4.281769e-38
[329,] 1.0000000 3.555081e-37 1.777541e-37
[330,] 1.0000000 1.453380e-36 7.266900e-37
[331,] 1.0000000 5.915038e-36 2.957519e-36
[332,] 1.0000000 2.404253e-35 1.202127e-35
[333,] 1.0000000 9.718749e-35 4.859374e-35
[334,] 1.0000000 3.895587e-34 1.947793e-34
[335,] 1.0000000 1.518998e-33 7.594991e-34
[336,] 1.0000000 6.003939e-33 3.001970e-33
[337,] 1.0000000 2.355889e-32 1.177945e-32
[338,] 1.0000000 9.176669e-32 4.588335e-32
[339,] 1.0000000 3.549111e-31 1.774555e-31
[340,] 1.0000000 1.362048e-30 6.810239e-31
[341,] 1.0000000 5.189148e-30 2.594574e-30
[342,] 1.0000000 1.961855e-29 9.809276e-30
[343,] 1.0000000 7.360371e-29 3.680185e-29
[344,] 1.0000000 2.380863e-28 1.190431e-28
[345,] 1.0000000 8.659573e-28 4.329787e-28
[346,] 1.0000000 3.188870e-27 1.594435e-27
[347,] 1.0000000 1.103430e-26 5.517151e-27
[348,] 1.0000000 3.981060e-26 1.990530e-26
[349,] 1.0000000 1.426829e-25 7.134143e-26
[350,] 1.0000000 5.071716e-25 2.535858e-25
[351,] 1.0000000 1.748621e-24 8.743104e-25
[352,] 1.0000000 6.112894e-24 3.056447e-24
[353,] 1.0000000 2.118770e-23 1.059385e-23
[354,] 1.0000000 7.280529e-23 3.640264e-23
[355,] 1.0000000 2.485938e-22 1.242969e-22
[356,] 1.0000000 8.393155e-22 4.196578e-22
[357,] 1.0000000 2.808421e-21 1.404210e-21
[358,] 1.0000000 8.991045e-21 4.495522e-21
[359,] 1.0000000 2.954851e-20 1.477426e-20
[360,] 1.0000000 9.622641e-20 4.811321e-20
[361,] 1.0000000 3.104179e-19 1.552089e-19
[362,] 1.0000000 9.918263e-19 4.959132e-19
[363,] 1.0000000 3.138352e-18 1.569176e-18
[364,] 1.0000000 9.817852e-18 4.908926e-18
[365,] 1.0000000 3.045494e-17 1.522747e-17
[366,] 1.0000000 9.051956e-17 4.525978e-17
[367,] 1.0000000 2.570163e-16 1.285082e-16
[368,] 1.0000000 7.790569e-16 3.895284e-16
[369,] 1.0000000 2.319963e-15 1.159982e-15
[370,] 1.0000000 6.834008e-15 3.417004e-15
[371,] 1.0000000 1.995787e-14 9.978936e-15
[372,] 1.0000000 5.534588e-14 2.767294e-14
[373,] 1.0000000 1.310007e-13 6.550035e-14
[374,] 1.0000000 3.703638e-13 1.851819e-13
[375,] 1.0000000 1.034688e-12 5.173440e-13
[376,] 1.0000000 2.842002e-12 1.421001e-12
[377,] 1.0000000 7.744822e-12 3.872411e-12
[378,] 1.0000000 2.089218e-11 1.044609e-11
[379,] 1.0000000 5.591317e-11 2.795658e-11
[380,] 1.0000000 1.404408e-10 7.022041e-11
[381,] 1.0000000 3.550831e-10 1.775415e-10
[382,] 1.0000000 8.950069e-10 4.475034e-10
[383,] 1.0000000 2.271859e-09 1.135929e-09
[384,] 1.0000000 5.678978e-09 2.839489e-09
[385,] 1.0000000 1.370072e-08 6.850361e-09
[386,] 1.0000000 3.383562e-08 1.691781e-08
[387,] 1.0000000 8.033980e-08 4.016990e-08
[388,] 0.9999999 1.448434e-07 7.242171e-08
[389,] 0.9999998 3.084942e-07 1.542471e-07
[390,] 0.9999996 7.276590e-07 3.638295e-07
[391,] 0.9999992 1.628309e-06 8.141547e-07
[392,] 0.9999982 3.514232e-06 1.757116e-06
[393,] 0.9999961 7.870591e-06 3.935296e-06
[394,] 0.9999919 1.612224e-05 8.061119e-06
[395,] 0.9999830 3.390226e-05 1.695113e-05
[396,] 0.9999671 6.589578e-05 3.294789e-05
[397,] 0.9999313 1.373156e-04 6.865778e-05
[398,] 0.9998585 2.830867e-04 1.415433e-04
[399,] 0.9997132 5.736757e-04 2.868378e-04
[400,] 0.9994810 1.037927e-03 5.189637e-04
[401,] 0.9990341 1.931895e-03 9.659476e-04
[402,] 0.9982226 3.554817e-03 1.777408e-03
[403,] 0.9967418 6.516351e-03 3.258176e-03
[404,] 0.9941257 1.174856e-02 5.874278e-03
[405,] 0.9896106 2.077874e-02 1.038937e-02
[406,] 0.9867937 2.641261e-02 1.320630e-02
[407,] 0.9772497 4.550061e-02 2.275031e-02
[408,] 0.9829130 3.417392e-02 1.708696e-02
[409,] 0.9715303 5.693933e-02 2.846967e-02
[410,] 0.9520033 9.599342e-02 4.799671e-02
[411,] 0.9563626 8.727480e-02 4.363740e-02
[412,] 0.9290470 1.419060e-01 7.095298e-02
[413,] 0.8864334 2.271332e-01 1.135666e-01
[414,] 0.8241096 3.517809e-01 1.758904e-01
[415,] 0.7574574 4.850853e-01 2.425426e-01
[416,] 0.6392711 7.214578e-01 3.607289e-01
[417,] 0.6939126 6.121748e-01 3.060874e-01
[418,] 0.6362986 7.274027e-01 3.637014e-01
> postscript(file="/var/www/rcomp/tmp/1xk311292955502.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/rcomp/tmp/28u241292955502.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/rcomp/tmp/38u241292955502.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/rcomp/tmp/48u241292955502.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/rcomp/tmp/5jljo1292955502.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
-6.791538e+04 2.739587e+06 1.732806e+06 -2.186068e+06 -3.431192e+04
6 7 8 9 10
7.204795e+05 2.034611e+05 6.736496e+05 -7.155026e+04 8.187372e+05
11 12 13 14 15
4.565894e+05 3.973043e+05 1.388178e+05 -7.598112e+05 -1.420388e+05
16 17 18 19 20
6.381048e+05 1.242467e+05 -1.449007e+05 -7.999024e+04 -5.425513e+05
21 22 23 24 25
3.768709e+05 2.074633e+05 1.535217e+05 -6.409222e+05 -1.780921e+05
26 27 28 29 30
1.328343e+05 -3.255006e+05 -2.663893e+05 4.060816e+05 -1.965210e+05
31 32 33 34 35
3.386834e+05 1.054774e+05 3.532760e+05 2.134773e+05 3.233606e+05
36 37 38 39 40
-4.582349e+05 1.689982e+05 2.355283e+05 9.624783e+04 -3.474496e+05
41 42 43 44 45
-4.310941e+05 -9.940253e+04 2.589757e+04 -1.049501e+05 8.662826e+03
46 47 48 49 50
2.827525e+04 -2.838113e+04 -1.817006e+05 5.221901e+04 -1.808105e+05
51 52 53 54 55
5.938359e+04 -4.664730e+04 -1.392787e+04 1.149660e+05 -2.914968e+05
56 57 58 59 60
-3.721347e+04 3.402509e+05 4.929928e+04 -4.921277e+04 -1.024846e+05
61 62 63 64 65
1.037701e+05 -2.097813e+05 6.692660e+04 -1.013692e+05 -4.588457e+04
66 67 68 69 70
-5.213330e+05 -6.213626e+04 1.697286e+04 -5.669127e+04 -4.091216e+04
71 72 73 74 75
-1.908425e+05 1.596308e+02 -5.076590e+04 1.403811e+04 3.322705e+04
76 77 78 79 80
-8.761648e+04 1.877659e+04 -7.912928e+04 -7.999789e+04 2.061823e+04
81 82 83 84 85
-2.934650e+03 -5.392035e+03 -3.384267e+03 -2.131900e+04 5.500812e+04
86 87 88 89 90
-8.237319e+04 -3.442450e+04 3.529264e+04 -7.355256e+04 6.539898e+04
91 92 93 94 95
2.825534e+04 2.998713e+03 2.670530e+04 2.612294e+02 -1.912698e+04
96 97 98 99 100
1.978783e+04 1.180071e+05 -6.952977e+03 6.879170e+04 4.442967e+03
101 102 103 104 105
-4.614226e+04 -9.576060e+04 -1.062081e+05 4.144130e+03 1.887356e+04
106 107 108 109 110
-7.890440e+04 -2.661074e+05 1.829283e+05 -7.612731e+04 -8.519118e+04
111 112 113 114 115
7.259471e+04 -2.051650e+05 -9.805053e+04 -3.358455e+04 7.011974e+04
116 117 118 119 120
-1.828131e+04 2.317555e+04 1.126503e+05 -1.047734e+05 -1.852112e+04
121 122 123 124 125
-1.357783e+05 1.642105e+04 -1.486211e+05 1.367024e+05 2.185824e+04
126 127 128 129 130
-1.269823e+05 2.103565e+04 2.002088e+03 -1.873865e+04 2.103565e+04
131 132 133 134 135
7.349548e+03 2.603461e+04 2.103565e+04 -1.757489e+04 9.281619e+03
136 137 138 139 140
1.285575e+04 6.327089e+03 1.873490e+04 1.162930e+04 2.103565e+04
141 142 143 144 145
1.596555e+04 2.103565e+04 -1.774000e+04 1.422555e+04 7.400699e+03
146 147 148 149 150
1.695051e+03 2.103565e+04 8.034654e+03 2.679814e+04 2.150725e+04
151 152 153 154 155
-1.634696e+04 -7.698966e+04 1.204308e+04 -1.010521e+04 1.726602e+04
156 157 158 159 160
2.103565e+04 2.103565e+04 1.711038e+04 3.228677e+04 2.103565e+04
161 162 163 164 165
2.103565e+04 2.103565e+04 2.103565e+04 -1.653319e+03 -3.650331e+04
166 167 168 169 170
1.865232e+04 2.103565e+04 5.417336e+04 -6.854867e+04 2.103565e+04
171 172 173 174 175
2.103565e+04 1.260324e+04 2.103565e+04 -1.363295e+03 2.103565e+04
176 177 178 179 180
6.552824e+03 -1.865200e+05 2.170283e+04 2.103565e+04 2.103565e+04
181 182 183 184 185
2.844871e+04 1.478473e+04 -3.150429e+04 -2.659205e+04 1.142357e+04
186 187 188 189 190
1.769393e+04 1.048717e+04 7.896868e+03 -1.199707e+04 2.103565e+04
191 192 193 194 195
1.124888e+04 2.103565e+04 8.143813e+03 2.103565e+04 1.589455e+04
196 197 198 199 200
2.103565e+04 -6.836521e+04 2.103565e+04 2.103565e+04 -6.008006e+04
201 202 203 204 205
2.103565e+04 2.103565e+04 2.103565e+04 -1.541284e+04 8.457068e+03
206 207 208 209 210
1.279829e+04 2.051418e+04 2.103565e+04 -5.240170e+02 6.701203e+03
211 212 213 214 215
2.103565e+04 2.103565e+04 1.535010e+04 -2.436018e+04 9.153589e+01
216 217 218 219 220
1.743862e+04 1.052239e+04 9.990486e+03 1.743862e+04 -2.747358e+04
221 222 223 224 225
-3.634792e+03 2.103565e+04 1.920559e+04 -1.931020e+04 1.939619e+04
226 227 228 229 230
-1.969494e+04 3.891318e+03 2.103565e+04 1.403666e+04 2.103565e+04
231 232 233 234 235
-8.881974e+03 9.287883e+03 2.103565e+04 -1.531985e+04 2.103565e+04
236 237 238 239 240
2.113902e+04 2.103565e+04 1.055096e+04 2.103565e+04 2.103565e+04
241 242 243 244 245
2.103565e+04 1.296730e+04 2.103565e+04 8.653990e+03 2.205947e+04
246 247 248 249 250
2.103565e+04 2.103565e+04 -1.263528e+04 2.103565e+04 2.103565e+04
251 252 253 254 255
-1.426362e+04 2.103565e+04 1.087066e+04 -6.021682e+04 -1.518544e+05
256 257 258 259 260
2.103565e+04 1.879902e+03 -4.361101e+04 2.103565e+04 2.103565e+04
261 262 263 264 265
2.103565e+04 2.103565e+04 2.569742e+04 1.896767e+04 2.103565e+04
266 267 268 269 270
2.103565e+04 2.103565e+04 7.025235e+04 2.103565e+04 -1.005395e+04
271 272 273 274 275
2.103565e+04 2.103565e+04 2.103565e+04 -1.861992e+04 2.103565e+04
276 277 278 279 280
3.466686e+03 -1.617307e+04 2.103565e+04 -5.832370e+04 1.571969e+03
281 282 283 284 285
2.103565e+04 2.103565e+04 -2.827914e+04 2.103565e+04 2.008941e+04
286 287 288 289 290
2.103565e+04 1.468397e+04 1.830072e+04 2.103565e+04 -2.552913e+04
291 292 293 294 295
2.103565e+04 8.885081e+03 3.633827e+04 2.103565e+04 -1.003406e+04
296 297 298 299 300
2.963102e+04 2.103565e+04 2.103565e+04 2.103565e+04 4.849025e+03
301 302 303 304 305
1.156722e+04 4.485463e+03 1.923713e+04 2.103565e+04 2.406910e+04
306 307 308 309 310
-5.369338e+02 2.103565e+04 -1.237875e+04 -1.599422e+04 2.103565e+04
311 312 313 314 315
2.103565e+04 -6.331702e+03 2.103565e+04 2.396020e+03 6.917307e+03
316 317 318 319 320
3.764782e+04 2.056159e+04 -3.408307e+04 1.203310e+04 2.103565e+04
321 322 323 324 325
2.103565e+04 2.103565e+04 2.103565e+04 2.103565e+04 2.103565e+04
326 327 328 329 330
1.631850e+04 2.103565e+04 2.103565e+04 2.103565e+04 2.103565e+04
331 332 333 334 335
2.103565e+04 2.103565e+04 2.103565e+04 1.200381e+04 1.924329e+04
336 337 338 339 340
1.322164e+04 -3.626141e+04 2.103565e+04 1.768911e+04 5.100661e+03
341 342 343 344 345
-2.101790e+04 2.103565e+04 2.103565e+04 2.103565e+04 -1.086091e+04
346 347 348 349 350
2.103565e+04 3.709479e+03 2.103565e+04 2.103565e+04 3.621747e+03
351 352 353 354 355
1.835330e+04 6.992397e+03 1.211169e+04 1.781777e+04 2.103565e+04
356 357 358 359 360
2.103565e+04 -1.298819e+02 2.103565e+04 2.103565e+04 2.103565e+04
361 362 363 364 365
1.385765e+04 2.103565e+04 2.103565e+04 -4.118696e+04 5.475234e+03
366 367 368 369 370
2.103565e+04 2.103565e+04 2.103565e+04 2.103565e+04 8.173406e+03
371 372 373 374 375
2.103565e+04 1.551123e+04 2.924101e+04 7.351272e+03 2.103565e+04
376 377 378 379 380
2.103565e+04 -1.902126e+04 1.416519e+04 2.794641e+03 2.103565e+04
381 382 383 384 385
2.103565e+04 -1.070157e+04 -4.042411e+04 2.103565e+04 1.231220e+03
386 387 388 389 390
-1.152632e+05 2.708052e+04 2.046047e+03 1.433238e+04 -1.990561e+03
391 392 393 394 395
9.146820e+03 -1.390456e+04 -4.674842e+04 6.112366e+04 -3.384755e+04
396 397 398 399 400
-2.631279e+04 -1.105763e+05 -2.996338e+05 -1.175435e+04 -1.441830e+05
401 402 403 404 405
2.482578e+03 2.316071e+04 -9.445398e+03 -5.694448e+04 -3.124614e+04
406 407 408 409 410
1.030738e+04 1.566919e+03 -1.355572e+04 -1.335913e+04 -6.926762e+04
411 412 413 414 415
-9.224149e+04 -6.372845e+04 -5.425446e+04 8.686138e+04 -1.529012e+05
416 417 418 419 420
-9.543057e+04 6.308139e+04 -4.264355e+04 -6.550312e+04 -5.082283e+04
421 422 423 424 425
-4.454602e+04 -1.986108e+05 -5.465750e+05 -3.885853e+05 -6.319665e+04
426 427 428 429 430
-2.494920e+05 -2.993732e+05 -7.770830e+05 -1.248715e+05 -1.470718e+05
431
-4.993672e+05
> postscript(file="/var/www/rcomp/tmp/6jljo1292955502.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 -6.791538e+04 NA
1 2.739587e+06 -6.791538e+04
2 1.732806e+06 2.739587e+06
3 -2.186068e+06 1.732806e+06
4 -3.431192e+04 -2.186068e+06
5 7.204795e+05 -3.431192e+04
6 2.034611e+05 7.204795e+05
7 6.736496e+05 2.034611e+05
8 -7.155026e+04 6.736496e+05
9 8.187372e+05 -7.155026e+04
10 4.565894e+05 8.187372e+05
11 3.973043e+05 4.565894e+05
12 1.388178e+05 3.973043e+05
13 -7.598112e+05 1.388178e+05
14 -1.420388e+05 -7.598112e+05
15 6.381048e+05 -1.420388e+05
16 1.242467e+05 6.381048e+05
17 -1.449007e+05 1.242467e+05
18 -7.999024e+04 -1.449007e+05
19 -5.425513e+05 -7.999024e+04
20 3.768709e+05 -5.425513e+05
21 2.074633e+05 3.768709e+05
22 1.535217e+05 2.074633e+05
23 -6.409222e+05 1.535217e+05
24 -1.780921e+05 -6.409222e+05
25 1.328343e+05 -1.780921e+05
26 -3.255006e+05 1.328343e+05
27 -2.663893e+05 -3.255006e+05
28 4.060816e+05 -2.663893e+05
29 -1.965210e+05 4.060816e+05
30 3.386834e+05 -1.965210e+05
31 1.054774e+05 3.386834e+05
32 3.532760e+05 1.054774e+05
33 2.134773e+05 3.532760e+05
34 3.233606e+05 2.134773e+05
35 -4.582349e+05 3.233606e+05
36 1.689982e+05 -4.582349e+05
37 2.355283e+05 1.689982e+05
38 9.624783e+04 2.355283e+05
39 -3.474496e+05 9.624783e+04
40 -4.310941e+05 -3.474496e+05
41 -9.940253e+04 -4.310941e+05
42 2.589757e+04 -9.940253e+04
43 -1.049501e+05 2.589757e+04
44 8.662826e+03 -1.049501e+05
45 2.827525e+04 8.662826e+03
46 -2.838113e+04 2.827525e+04
47 -1.817006e+05 -2.838113e+04
48 5.221901e+04 -1.817006e+05
49 -1.808105e+05 5.221901e+04
50 5.938359e+04 -1.808105e+05
51 -4.664730e+04 5.938359e+04
52 -1.392787e+04 -4.664730e+04
53 1.149660e+05 -1.392787e+04
54 -2.914968e+05 1.149660e+05
55 -3.721347e+04 -2.914968e+05
56 3.402509e+05 -3.721347e+04
57 4.929928e+04 3.402509e+05
58 -4.921277e+04 4.929928e+04
59 -1.024846e+05 -4.921277e+04
60 1.037701e+05 -1.024846e+05
61 -2.097813e+05 1.037701e+05
62 6.692660e+04 -2.097813e+05
63 -1.013692e+05 6.692660e+04
64 -4.588457e+04 -1.013692e+05
65 -5.213330e+05 -4.588457e+04
66 -6.213626e+04 -5.213330e+05
67 1.697286e+04 -6.213626e+04
68 -5.669127e+04 1.697286e+04
69 -4.091216e+04 -5.669127e+04
70 -1.908425e+05 -4.091216e+04
71 1.596308e+02 -1.908425e+05
72 -5.076590e+04 1.596308e+02
73 1.403811e+04 -5.076590e+04
74 3.322705e+04 1.403811e+04
75 -8.761648e+04 3.322705e+04
76 1.877659e+04 -8.761648e+04
77 -7.912928e+04 1.877659e+04
78 -7.999789e+04 -7.912928e+04
79 2.061823e+04 -7.999789e+04
80 -2.934650e+03 2.061823e+04
81 -5.392035e+03 -2.934650e+03
82 -3.384267e+03 -5.392035e+03
83 -2.131900e+04 -3.384267e+03
84 5.500812e+04 -2.131900e+04
85 -8.237319e+04 5.500812e+04
86 -3.442450e+04 -8.237319e+04
87 3.529264e+04 -3.442450e+04
88 -7.355256e+04 3.529264e+04
89 6.539898e+04 -7.355256e+04
90 2.825534e+04 6.539898e+04
91 2.998713e+03 2.825534e+04
92 2.670530e+04 2.998713e+03
93 2.612294e+02 2.670530e+04
94 -1.912698e+04 2.612294e+02
95 1.978783e+04 -1.912698e+04
96 1.180071e+05 1.978783e+04
97 -6.952977e+03 1.180071e+05
98 6.879170e+04 -6.952977e+03
99 4.442967e+03 6.879170e+04
100 -4.614226e+04 4.442967e+03
101 -9.576060e+04 -4.614226e+04
102 -1.062081e+05 -9.576060e+04
103 4.144130e+03 -1.062081e+05
104 1.887356e+04 4.144130e+03
105 -7.890440e+04 1.887356e+04
106 -2.661074e+05 -7.890440e+04
107 1.829283e+05 -2.661074e+05
108 -7.612731e+04 1.829283e+05
109 -8.519118e+04 -7.612731e+04
110 7.259471e+04 -8.519118e+04
111 -2.051650e+05 7.259471e+04
112 -9.805053e+04 -2.051650e+05
113 -3.358455e+04 -9.805053e+04
114 7.011974e+04 -3.358455e+04
115 -1.828131e+04 7.011974e+04
116 2.317555e+04 -1.828131e+04
117 1.126503e+05 2.317555e+04
118 -1.047734e+05 1.126503e+05
119 -1.852112e+04 -1.047734e+05
120 -1.357783e+05 -1.852112e+04
121 1.642105e+04 -1.357783e+05
122 -1.486211e+05 1.642105e+04
123 1.367024e+05 -1.486211e+05
124 2.185824e+04 1.367024e+05
125 -1.269823e+05 2.185824e+04
126 2.103565e+04 -1.269823e+05
127 2.002088e+03 2.103565e+04
128 -1.873865e+04 2.002088e+03
129 2.103565e+04 -1.873865e+04
130 7.349548e+03 2.103565e+04
131 2.603461e+04 7.349548e+03
132 2.103565e+04 2.603461e+04
133 -1.757489e+04 2.103565e+04
134 9.281619e+03 -1.757489e+04
135 1.285575e+04 9.281619e+03
136 6.327089e+03 1.285575e+04
137 1.873490e+04 6.327089e+03
138 1.162930e+04 1.873490e+04
139 2.103565e+04 1.162930e+04
140 1.596555e+04 2.103565e+04
141 2.103565e+04 1.596555e+04
142 -1.774000e+04 2.103565e+04
143 1.422555e+04 -1.774000e+04
144 7.400699e+03 1.422555e+04
145 1.695051e+03 7.400699e+03
146 2.103565e+04 1.695051e+03
147 8.034654e+03 2.103565e+04
148 2.679814e+04 8.034654e+03
149 2.150725e+04 2.679814e+04
150 -1.634696e+04 2.150725e+04
151 -7.698966e+04 -1.634696e+04
152 1.204308e+04 -7.698966e+04
153 -1.010521e+04 1.204308e+04
154 1.726602e+04 -1.010521e+04
155 2.103565e+04 1.726602e+04
156 2.103565e+04 2.103565e+04
157 1.711038e+04 2.103565e+04
158 3.228677e+04 1.711038e+04
159 2.103565e+04 3.228677e+04
160 2.103565e+04 2.103565e+04
161 2.103565e+04 2.103565e+04
162 2.103565e+04 2.103565e+04
163 -1.653319e+03 2.103565e+04
164 -3.650331e+04 -1.653319e+03
165 1.865232e+04 -3.650331e+04
166 2.103565e+04 1.865232e+04
167 5.417336e+04 2.103565e+04
168 -6.854867e+04 5.417336e+04
169 2.103565e+04 -6.854867e+04
170 2.103565e+04 2.103565e+04
171 1.260324e+04 2.103565e+04
172 2.103565e+04 1.260324e+04
173 -1.363295e+03 2.103565e+04
174 2.103565e+04 -1.363295e+03
175 6.552824e+03 2.103565e+04
176 -1.865200e+05 6.552824e+03
177 2.170283e+04 -1.865200e+05
178 2.103565e+04 2.170283e+04
179 2.103565e+04 2.103565e+04
180 2.844871e+04 2.103565e+04
181 1.478473e+04 2.844871e+04
182 -3.150429e+04 1.478473e+04
183 -2.659205e+04 -3.150429e+04
184 1.142357e+04 -2.659205e+04
185 1.769393e+04 1.142357e+04
186 1.048717e+04 1.769393e+04
187 7.896868e+03 1.048717e+04
188 -1.199707e+04 7.896868e+03
189 2.103565e+04 -1.199707e+04
190 1.124888e+04 2.103565e+04
191 2.103565e+04 1.124888e+04
192 8.143813e+03 2.103565e+04
193 2.103565e+04 8.143813e+03
194 1.589455e+04 2.103565e+04
195 2.103565e+04 1.589455e+04
196 -6.836521e+04 2.103565e+04
197 2.103565e+04 -6.836521e+04
198 2.103565e+04 2.103565e+04
199 -6.008006e+04 2.103565e+04
200 2.103565e+04 -6.008006e+04
201 2.103565e+04 2.103565e+04
202 2.103565e+04 2.103565e+04
203 -1.541284e+04 2.103565e+04
204 8.457068e+03 -1.541284e+04
205 1.279829e+04 8.457068e+03
206 2.051418e+04 1.279829e+04
207 2.103565e+04 2.051418e+04
208 -5.240170e+02 2.103565e+04
209 6.701203e+03 -5.240170e+02
210 2.103565e+04 6.701203e+03
211 2.103565e+04 2.103565e+04
212 1.535010e+04 2.103565e+04
213 -2.436018e+04 1.535010e+04
214 9.153589e+01 -2.436018e+04
215 1.743862e+04 9.153589e+01
216 1.052239e+04 1.743862e+04
217 9.990486e+03 1.052239e+04
218 1.743862e+04 9.990486e+03
219 -2.747358e+04 1.743862e+04
220 -3.634792e+03 -2.747358e+04
221 2.103565e+04 -3.634792e+03
222 1.920559e+04 2.103565e+04
223 -1.931020e+04 1.920559e+04
224 1.939619e+04 -1.931020e+04
225 -1.969494e+04 1.939619e+04
226 3.891318e+03 -1.969494e+04
227 2.103565e+04 3.891318e+03
228 1.403666e+04 2.103565e+04
229 2.103565e+04 1.403666e+04
230 -8.881974e+03 2.103565e+04
231 9.287883e+03 -8.881974e+03
232 2.103565e+04 9.287883e+03
233 -1.531985e+04 2.103565e+04
234 2.103565e+04 -1.531985e+04
235 2.113902e+04 2.103565e+04
236 2.103565e+04 2.113902e+04
237 1.055096e+04 2.103565e+04
238 2.103565e+04 1.055096e+04
239 2.103565e+04 2.103565e+04
240 2.103565e+04 2.103565e+04
241 1.296730e+04 2.103565e+04
242 2.103565e+04 1.296730e+04
243 8.653990e+03 2.103565e+04
244 2.205947e+04 8.653990e+03
245 2.103565e+04 2.205947e+04
246 2.103565e+04 2.103565e+04
247 -1.263528e+04 2.103565e+04
248 2.103565e+04 -1.263528e+04
249 2.103565e+04 2.103565e+04
250 -1.426362e+04 2.103565e+04
251 2.103565e+04 -1.426362e+04
252 1.087066e+04 2.103565e+04
253 -6.021682e+04 1.087066e+04
254 -1.518544e+05 -6.021682e+04
255 2.103565e+04 -1.518544e+05
256 1.879902e+03 2.103565e+04
257 -4.361101e+04 1.879902e+03
258 2.103565e+04 -4.361101e+04
259 2.103565e+04 2.103565e+04
260 2.103565e+04 2.103565e+04
261 2.103565e+04 2.103565e+04
262 2.569742e+04 2.103565e+04
263 1.896767e+04 2.569742e+04
264 2.103565e+04 1.896767e+04
265 2.103565e+04 2.103565e+04
266 2.103565e+04 2.103565e+04
267 7.025235e+04 2.103565e+04
268 2.103565e+04 7.025235e+04
269 -1.005395e+04 2.103565e+04
270 2.103565e+04 -1.005395e+04
271 2.103565e+04 2.103565e+04
272 2.103565e+04 2.103565e+04
273 -1.861992e+04 2.103565e+04
274 2.103565e+04 -1.861992e+04
275 3.466686e+03 2.103565e+04
276 -1.617307e+04 3.466686e+03
277 2.103565e+04 -1.617307e+04
278 -5.832370e+04 2.103565e+04
279 1.571969e+03 -5.832370e+04
280 2.103565e+04 1.571969e+03
281 2.103565e+04 2.103565e+04
282 -2.827914e+04 2.103565e+04
283 2.103565e+04 -2.827914e+04
284 2.008941e+04 2.103565e+04
285 2.103565e+04 2.008941e+04
286 1.468397e+04 2.103565e+04
287 1.830072e+04 1.468397e+04
288 2.103565e+04 1.830072e+04
289 -2.552913e+04 2.103565e+04
290 2.103565e+04 -2.552913e+04
291 8.885081e+03 2.103565e+04
292 3.633827e+04 8.885081e+03
293 2.103565e+04 3.633827e+04
294 -1.003406e+04 2.103565e+04
295 2.963102e+04 -1.003406e+04
296 2.103565e+04 2.963102e+04
297 2.103565e+04 2.103565e+04
298 2.103565e+04 2.103565e+04
299 4.849025e+03 2.103565e+04
300 1.156722e+04 4.849025e+03
301 4.485463e+03 1.156722e+04
302 1.923713e+04 4.485463e+03
303 2.103565e+04 1.923713e+04
304 2.406910e+04 2.103565e+04
305 -5.369338e+02 2.406910e+04
306 2.103565e+04 -5.369338e+02
307 -1.237875e+04 2.103565e+04
308 -1.599422e+04 -1.237875e+04
309 2.103565e+04 -1.599422e+04
310 2.103565e+04 2.103565e+04
311 -6.331702e+03 2.103565e+04
312 2.103565e+04 -6.331702e+03
313 2.396020e+03 2.103565e+04
314 6.917307e+03 2.396020e+03
315 3.764782e+04 6.917307e+03
316 2.056159e+04 3.764782e+04
317 -3.408307e+04 2.056159e+04
318 1.203310e+04 -3.408307e+04
319 2.103565e+04 1.203310e+04
320 2.103565e+04 2.103565e+04
321 2.103565e+04 2.103565e+04
322 2.103565e+04 2.103565e+04
323 2.103565e+04 2.103565e+04
324 2.103565e+04 2.103565e+04
325 1.631850e+04 2.103565e+04
326 2.103565e+04 1.631850e+04
327 2.103565e+04 2.103565e+04
328 2.103565e+04 2.103565e+04
329 2.103565e+04 2.103565e+04
330 2.103565e+04 2.103565e+04
331 2.103565e+04 2.103565e+04
332 2.103565e+04 2.103565e+04
333 1.200381e+04 2.103565e+04
334 1.924329e+04 1.200381e+04
335 1.322164e+04 1.924329e+04
336 -3.626141e+04 1.322164e+04
337 2.103565e+04 -3.626141e+04
338 1.768911e+04 2.103565e+04
339 5.100661e+03 1.768911e+04
340 -2.101790e+04 5.100661e+03
341 2.103565e+04 -2.101790e+04
342 2.103565e+04 2.103565e+04
343 2.103565e+04 2.103565e+04
344 -1.086091e+04 2.103565e+04
345 2.103565e+04 -1.086091e+04
346 3.709479e+03 2.103565e+04
347 2.103565e+04 3.709479e+03
348 2.103565e+04 2.103565e+04
349 3.621747e+03 2.103565e+04
350 1.835330e+04 3.621747e+03
351 6.992397e+03 1.835330e+04
352 1.211169e+04 6.992397e+03
353 1.781777e+04 1.211169e+04
354 2.103565e+04 1.781777e+04
355 2.103565e+04 2.103565e+04
356 -1.298819e+02 2.103565e+04
357 2.103565e+04 -1.298819e+02
358 2.103565e+04 2.103565e+04
359 2.103565e+04 2.103565e+04
360 1.385765e+04 2.103565e+04
361 2.103565e+04 1.385765e+04
362 2.103565e+04 2.103565e+04
363 -4.118696e+04 2.103565e+04
364 5.475234e+03 -4.118696e+04
365 2.103565e+04 5.475234e+03
366 2.103565e+04 2.103565e+04
367 2.103565e+04 2.103565e+04
368 2.103565e+04 2.103565e+04
369 8.173406e+03 2.103565e+04
370 2.103565e+04 8.173406e+03
371 1.551123e+04 2.103565e+04
372 2.924101e+04 1.551123e+04
373 7.351272e+03 2.924101e+04
374 2.103565e+04 7.351272e+03
375 2.103565e+04 2.103565e+04
376 -1.902126e+04 2.103565e+04
377 1.416519e+04 -1.902126e+04
378 2.794641e+03 1.416519e+04
379 2.103565e+04 2.794641e+03
380 2.103565e+04 2.103565e+04
381 -1.070157e+04 2.103565e+04
382 -4.042411e+04 -1.070157e+04
383 2.103565e+04 -4.042411e+04
384 1.231220e+03 2.103565e+04
385 -1.152632e+05 1.231220e+03
386 2.708052e+04 -1.152632e+05
387 2.046047e+03 2.708052e+04
388 1.433238e+04 2.046047e+03
389 -1.990561e+03 1.433238e+04
390 9.146820e+03 -1.990561e+03
391 -1.390456e+04 9.146820e+03
392 -4.674842e+04 -1.390456e+04
393 6.112366e+04 -4.674842e+04
394 -3.384755e+04 6.112366e+04
395 -2.631279e+04 -3.384755e+04
396 -1.105763e+05 -2.631279e+04
397 -2.996338e+05 -1.105763e+05
398 -1.175435e+04 -2.996338e+05
399 -1.441830e+05 -1.175435e+04
400 2.482578e+03 -1.441830e+05
401 2.316071e+04 2.482578e+03
402 -9.445398e+03 2.316071e+04
403 -5.694448e+04 -9.445398e+03
404 -3.124614e+04 -5.694448e+04
405 1.030738e+04 -3.124614e+04
406 1.566919e+03 1.030738e+04
407 -1.355572e+04 1.566919e+03
408 -1.335913e+04 -1.355572e+04
409 -6.926762e+04 -1.335913e+04
410 -9.224149e+04 -6.926762e+04
411 -6.372845e+04 -9.224149e+04
412 -5.425446e+04 -6.372845e+04
413 8.686138e+04 -5.425446e+04
414 -1.529012e+05 8.686138e+04
415 -9.543057e+04 -1.529012e+05
416 6.308139e+04 -9.543057e+04
417 -4.264355e+04 6.308139e+04
418 -6.550312e+04 -4.264355e+04
419 -5.082283e+04 -6.550312e+04
420 -4.454602e+04 -5.082283e+04
421 -1.986108e+05 -4.454602e+04
422 -5.465750e+05 -1.986108e+05
423 -3.885853e+05 -5.465750e+05
424 -6.319665e+04 -3.885853e+05
425 -2.494920e+05 -6.319665e+04
426 -2.993732e+05 -2.494920e+05
427 -7.770830e+05 -2.993732e+05
428 -1.248715e+05 -7.770830e+05
429 -1.470718e+05 -1.248715e+05
430 -4.993672e+05 -1.470718e+05
431 NA -4.993672e+05
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] 2.739587e+06 -6.791538e+04
[2,] 1.732806e+06 2.739587e+06
[3,] -2.186068e+06 1.732806e+06
[4,] -3.431192e+04 -2.186068e+06
[5,] 7.204795e+05 -3.431192e+04
[6,] 2.034611e+05 7.204795e+05
[7,] 6.736496e+05 2.034611e+05
[8,] -7.155026e+04 6.736496e+05
[9,] 8.187372e+05 -7.155026e+04
[10,] 4.565894e+05 8.187372e+05
[11,] 3.973043e+05 4.565894e+05
[12,] 1.388178e+05 3.973043e+05
[13,] -7.598112e+05 1.388178e+05
[14,] -1.420388e+05 -7.598112e+05
[15,] 6.381048e+05 -1.420388e+05
[16,] 1.242467e+05 6.381048e+05
[17,] -1.449007e+05 1.242467e+05
[18,] -7.999024e+04 -1.449007e+05
[19,] -5.425513e+05 -7.999024e+04
[20,] 3.768709e+05 -5.425513e+05
[21,] 2.074633e+05 3.768709e+05
[22,] 1.535217e+05 2.074633e+05
[23,] -6.409222e+05 1.535217e+05
[24,] -1.780921e+05 -6.409222e+05
[25,] 1.328343e+05 -1.780921e+05
[26,] -3.255006e+05 1.328343e+05
[27,] -2.663893e+05 -3.255006e+05
[28,] 4.060816e+05 -2.663893e+05
[29,] -1.965210e+05 4.060816e+05
[30,] 3.386834e+05 -1.965210e+05
[31,] 1.054774e+05 3.386834e+05
[32,] 3.532760e+05 1.054774e+05
[33,] 2.134773e+05 3.532760e+05
[34,] 3.233606e+05 2.134773e+05
[35,] -4.582349e+05 3.233606e+05
[36,] 1.689982e+05 -4.582349e+05
[37,] 2.355283e+05 1.689982e+05
[38,] 9.624783e+04 2.355283e+05
[39,] -3.474496e+05 9.624783e+04
[40,] -4.310941e+05 -3.474496e+05
[41,] -9.940253e+04 -4.310941e+05
[42,] 2.589757e+04 -9.940253e+04
[43,] -1.049501e+05 2.589757e+04
[44,] 8.662826e+03 -1.049501e+05
[45,] 2.827525e+04 8.662826e+03
[46,] -2.838113e+04 2.827525e+04
[47,] -1.817006e+05 -2.838113e+04
[48,] 5.221901e+04 -1.817006e+05
[49,] -1.808105e+05 5.221901e+04
[50,] 5.938359e+04 -1.808105e+05
[51,] -4.664730e+04 5.938359e+04
[52,] -1.392787e+04 -4.664730e+04
[53,] 1.149660e+05 -1.392787e+04
[54,] -2.914968e+05 1.149660e+05
[55,] -3.721347e+04 -2.914968e+05
[56,] 3.402509e+05 -3.721347e+04
[57,] 4.929928e+04 3.402509e+05
[58,] -4.921277e+04 4.929928e+04
[59,] -1.024846e+05 -4.921277e+04
[60,] 1.037701e+05 -1.024846e+05
[61,] -2.097813e+05 1.037701e+05
[62,] 6.692660e+04 -2.097813e+05
[63,] -1.013692e+05 6.692660e+04
[64,] -4.588457e+04 -1.013692e+05
[65,] -5.213330e+05 -4.588457e+04
[66,] -6.213626e+04 -5.213330e+05
[67,] 1.697286e+04 -6.213626e+04
[68,] -5.669127e+04 1.697286e+04
[69,] -4.091216e+04 -5.669127e+04
[70,] -1.908425e+05 -4.091216e+04
[71,] 1.596308e+02 -1.908425e+05
[72,] -5.076590e+04 1.596308e+02
[73,] 1.403811e+04 -5.076590e+04
[74,] 3.322705e+04 1.403811e+04
[75,] -8.761648e+04 3.322705e+04
[76,] 1.877659e+04 -8.761648e+04
[77,] -7.912928e+04 1.877659e+04
[78,] -7.999789e+04 -7.912928e+04
[79,] 2.061823e+04 -7.999789e+04
[80,] -2.934650e+03 2.061823e+04
[81,] -5.392035e+03 -2.934650e+03
[82,] -3.384267e+03 -5.392035e+03
[83,] -2.131900e+04 -3.384267e+03
[84,] 5.500812e+04 -2.131900e+04
[85,] -8.237319e+04 5.500812e+04
[86,] -3.442450e+04 -8.237319e+04
[87,] 3.529264e+04 -3.442450e+04
[88,] -7.355256e+04 3.529264e+04
[89,] 6.539898e+04 -7.355256e+04
[90,] 2.825534e+04 6.539898e+04
[91,] 2.998713e+03 2.825534e+04
[92,] 2.670530e+04 2.998713e+03
[93,] 2.612294e+02 2.670530e+04
[94,] -1.912698e+04 2.612294e+02
[95,] 1.978783e+04 -1.912698e+04
[96,] 1.180071e+05 1.978783e+04
[97,] -6.952977e+03 1.180071e+05
[98,] 6.879170e+04 -6.952977e+03
[99,] 4.442967e+03 6.879170e+04
[100,] -4.614226e+04 4.442967e+03
[101,] -9.576060e+04 -4.614226e+04
[102,] -1.062081e+05 -9.576060e+04
[103,] 4.144130e+03 -1.062081e+05
[104,] 1.887356e+04 4.144130e+03
[105,] -7.890440e+04 1.887356e+04
[106,] -2.661074e+05 -7.890440e+04
[107,] 1.829283e+05 -2.661074e+05
[108,] -7.612731e+04 1.829283e+05
[109,] -8.519118e+04 -7.612731e+04
[110,] 7.259471e+04 -8.519118e+04
[111,] -2.051650e+05 7.259471e+04
[112,] -9.805053e+04 -2.051650e+05
[113,] -3.358455e+04 -9.805053e+04
[114,] 7.011974e+04 -3.358455e+04
[115,] -1.828131e+04 7.011974e+04
[116,] 2.317555e+04 -1.828131e+04
[117,] 1.126503e+05 2.317555e+04
[118,] -1.047734e+05 1.126503e+05
[119,] -1.852112e+04 -1.047734e+05
[120,] -1.357783e+05 -1.852112e+04
[121,] 1.642105e+04 -1.357783e+05
[122,] -1.486211e+05 1.642105e+04
[123,] 1.367024e+05 -1.486211e+05
[124,] 2.185824e+04 1.367024e+05
[125,] -1.269823e+05 2.185824e+04
[126,] 2.103565e+04 -1.269823e+05
[127,] 2.002088e+03 2.103565e+04
[128,] -1.873865e+04 2.002088e+03
[129,] 2.103565e+04 -1.873865e+04
[130,] 7.349548e+03 2.103565e+04
[131,] 2.603461e+04 7.349548e+03
[132,] 2.103565e+04 2.603461e+04
[133,] -1.757489e+04 2.103565e+04
[134,] 9.281619e+03 -1.757489e+04
[135,] 1.285575e+04 9.281619e+03
[136,] 6.327089e+03 1.285575e+04
[137,] 1.873490e+04 6.327089e+03
[138,] 1.162930e+04 1.873490e+04
[139,] 2.103565e+04 1.162930e+04
[140,] 1.596555e+04 2.103565e+04
[141,] 2.103565e+04 1.596555e+04
[142,] -1.774000e+04 2.103565e+04
[143,] 1.422555e+04 -1.774000e+04
[144,] 7.400699e+03 1.422555e+04
[145,] 1.695051e+03 7.400699e+03
[146,] 2.103565e+04 1.695051e+03
[147,] 8.034654e+03 2.103565e+04
[148,] 2.679814e+04 8.034654e+03
[149,] 2.150725e+04 2.679814e+04
[150,] -1.634696e+04 2.150725e+04
[151,] -7.698966e+04 -1.634696e+04
[152,] 1.204308e+04 -7.698966e+04
[153,] -1.010521e+04 1.204308e+04
[154,] 1.726602e+04 -1.010521e+04
[155,] 2.103565e+04 1.726602e+04
[156,] 2.103565e+04 2.103565e+04
[157,] 1.711038e+04 2.103565e+04
[158,] 3.228677e+04 1.711038e+04
[159,] 2.103565e+04 3.228677e+04
[160,] 2.103565e+04 2.103565e+04
[161,] 2.103565e+04 2.103565e+04
[162,] 2.103565e+04 2.103565e+04
[163,] -1.653319e+03 2.103565e+04
[164,] -3.650331e+04 -1.653319e+03
[165,] 1.865232e+04 -3.650331e+04
[166,] 2.103565e+04 1.865232e+04
[167,] 5.417336e+04 2.103565e+04
[168,] -6.854867e+04 5.417336e+04
[169,] 2.103565e+04 -6.854867e+04
[170,] 2.103565e+04 2.103565e+04
[171,] 1.260324e+04 2.103565e+04
[172,] 2.103565e+04 1.260324e+04
[173,] -1.363295e+03 2.103565e+04
[174,] 2.103565e+04 -1.363295e+03
[175,] 6.552824e+03 2.103565e+04
[176,] -1.865200e+05 6.552824e+03
[177,] 2.170283e+04 -1.865200e+05
[178,] 2.103565e+04 2.170283e+04
[179,] 2.103565e+04 2.103565e+04
[180,] 2.844871e+04 2.103565e+04
[181,] 1.478473e+04 2.844871e+04
[182,] -3.150429e+04 1.478473e+04
[183,] -2.659205e+04 -3.150429e+04
[184,] 1.142357e+04 -2.659205e+04
[185,] 1.769393e+04 1.142357e+04
[186,] 1.048717e+04 1.769393e+04
[187,] 7.896868e+03 1.048717e+04
[188,] -1.199707e+04 7.896868e+03
[189,] 2.103565e+04 -1.199707e+04
[190,] 1.124888e+04 2.103565e+04
[191,] 2.103565e+04 1.124888e+04
[192,] 8.143813e+03 2.103565e+04
[193,] 2.103565e+04 8.143813e+03
[194,] 1.589455e+04 2.103565e+04
[195,] 2.103565e+04 1.589455e+04
[196,] -6.836521e+04 2.103565e+04
[197,] 2.103565e+04 -6.836521e+04
[198,] 2.103565e+04 2.103565e+04
[199,] -6.008006e+04 2.103565e+04
[200,] 2.103565e+04 -6.008006e+04
[201,] 2.103565e+04 2.103565e+04
[202,] 2.103565e+04 2.103565e+04
[203,] -1.541284e+04 2.103565e+04
[204,] 8.457068e+03 -1.541284e+04
[205,] 1.279829e+04 8.457068e+03
[206,] 2.051418e+04 1.279829e+04
[207,] 2.103565e+04 2.051418e+04
[208,] -5.240170e+02 2.103565e+04
[209,] 6.701203e+03 -5.240170e+02
[210,] 2.103565e+04 6.701203e+03
[211,] 2.103565e+04 2.103565e+04
[212,] 1.535010e+04 2.103565e+04
[213,] -2.436018e+04 1.535010e+04
[214,] 9.153589e+01 -2.436018e+04
[215,] 1.743862e+04 9.153589e+01
[216,] 1.052239e+04 1.743862e+04
[217,] 9.990486e+03 1.052239e+04
[218,] 1.743862e+04 9.990486e+03
[219,] -2.747358e+04 1.743862e+04
[220,] -3.634792e+03 -2.747358e+04
[221,] 2.103565e+04 -3.634792e+03
[222,] 1.920559e+04 2.103565e+04
[223,] -1.931020e+04 1.920559e+04
[224,] 1.939619e+04 -1.931020e+04
[225,] -1.969494e+04 1.939619e+04
[226,] 3.891318e+03 -1.969494e+04
[227,] 2.103565e+04 3.891318e+03
[228,] 1.403666e+04 2.103565e+04
[229,] 2.103565e+04 1.403666e+04
[230,] -8.881974e+03 2.103565e+04
[231,] 9.287883e+03 -8.881974e+03
[232,] 2.103565e+04 9.287883e+03
[233,] -1.531985e+04 2.103565e+04
[234,] 2.103565e+04 -1.531985e+04
[235,] 2.113902e+04 2.103565e+04
[236,] 2.103565e+04 2.113902e+04
[237,] 1.055096e+04 2.103565e+04
[238,] 2.103565e+04 1.055096e+04
[239,] 2.103565e+04 2.103565e+04
[240,] 2.103565e+04 2.103565e+04
[241,] 1.296730e+04 2.103565e+04
[242,] 2.103565e+04 1.296730e+04
[243,] 8.653990e+03 2.103565e+04
[244,] 2.205947e+04 8.653990e+03
[245,] 2.103565e+04 2.205947e+04
[246,] 2.103565e+04 2.103565e+04
[247,] -1.263528e+04 2.103565e+04
[248,] 2.103565e+04 -1.263528e+04
[249,] 2.103565e+04 2.103565e+04
[250,] -1.426362e+04 2.103565e+04
[251,] 2.103565e+04 -1.426362e+04
[252,] 1.087066e+04 2.103565e+04
[253,] -6.021682e+04 1.087066e+04
[254,] -1.518544e+05 -6.021682e+04
[255,] 2.103565e+04 -1.518544e+05
[256,] 1.879902e+03 2.103565e+04
[257,] -4.361101e+04 1.879902e+03
[258,] 2.103565e+04 -4.361101e+04
[259,] 2.103565e+04 2.103565e+04
[260,] 2.103565e+04 2.103565e+04
[261,] 2.103565e+04 2.103565e+04
[262,] 2.569742e+04 2.103565e+04
[263,] 1.896767e+04 2.569742e+04
[264,] 2.103565e+04 1.896767e+04
[265,] 2.103565e+04 2.103565e+04
[266,] 2.103565e+04 2.103565e+04
[267,] 7.025235e+04 2.103565e+04
[268,] 2.103565e+04 7.025235e+04
[269,] -1.005395e+04 2.103565e+04
[270,] 2.103565e+04 -1.005395e+04
[271,] 2.103565e+04 2.103565e+04
[272,] 2.103565e+04 2.103565e+04
[273,] -1.861992e+04 2.103565e+04
[274,] 2.103565e+04 -1.861992e+04
[275,] 3.466686e+03 2.103565e+04
[276,] -1.617307e+04 3.466686e+03
[277,] 2.103565e+04 -1.617307e+04
[278,] -5.832370e+04 2.103565e+04
[279,] 1.571969e+03 -5.832370e+04
[280,] 2.103565e+04 1.571969e+03
[281,] 2.103565e+04 2.103565e+04
[282,] -2.827914e+04 2.103565e+04
[283,] 2.103565e+04 -2.827914e+04
[284,] 2.008941e+04 2.103565e+04
[285,] 2.103565e+04 2.008941e+04
[286,] 1.468397e+04 2.103565e+04
[287,] 1.830072e+04 1.468397e+04
[288,] 2.103565e+04 1.830072e+04
[289,] -2.552913e+04 2.103565e+04
[290,] 2.103565e+04 -2.552913e+04
[291,] 8.885081e+03 2.103565e+04
[292,] 3.633827e+04 8.885081e+03
[293,] 2.103565e+04 3.633827e+04
[294,] -1.003406e+04 2.103565e+04
[295,] 2.963102e+04 -1.003406e+04
[296,] 2.103565e+04 2.963102e+04
[297,] 2.103565e+04 2.103565e+04
[298,] 2.103565e+04 2.103565e+04
[299,] 4.849025e+03 2.103565e+04
[300,] 1.156722e+04 4.849025e+03
[301,] 4.485463e+03 1.156722e+04
[302,] 1.923713e+04 4.485463e+03
[303,] 2.103565e+04 1.923713e+04
[304,] 2.406910e+04 2.103565e+04
[305,] -5.369338e+02 2.406910e+04
[306,] 2.103565e+04 -5.369338e+02
[307,] -1.237875e+04 2.103565e+04
[308,] -1.599422e+04 -1.237875e+04
[309,] 2.103565e+04 -1.599422e+04
[310,] 2.103565e+04 2.103565e+04
[311,] -6.331702e+03 2.103565e+04
[312,] 2.103565e+04 -6.331702e+03
[313,] 2.396020e+03 2.103565e+04
[314,] 6.917307e+03 2.396020e+03
[315,] 3.764782e+04 6.917307e+03
[316,] 2.056159e+04 3.764782e+04
[317,] -3.408307e+04 2.056159e+04
[318,] 1.203310e+04 -3.408307e+04
[319,] 2.103565e+04 1.203310e+04
[320,] 2.103565e+04 2.103565e+04
[321,] 2.103565e+04 2.103565e+04
[322,] 2.103565e+04 2.103565e+04
[323,] 2.103565e+04 2.103565e+04
[324,] 2.103565e+04 2.103565e+04
[325,] 1.631850e+04 2.103565e+04
[326,] 2.103565e+04 1.631850e+04
[327,] 2.103565e+04 2.103565e+04
[328,] 2.103565e+04 2.103565e+04
[329,] 2.103565e+04 2.103565e+04
[330,] 2.103565e+04 2.103565e+04
[331,] 2.103565e+04 2.103565e+04
[332,] 2.103565e+04 2.103565e+04
[333,] 1.200381e+04 2.103565e+04
[334,] 1.924329e+04 1.200381e+04
[335,] 1.322164e+04 1.924329e+04
[336,] -3.626141e+04 1.322164e+04
[337,] 2.103565e+04 -3.626141e+04
[338,] 1.768911e+04 2.103565e+04
[339,] 5.100661e+03 1.768911e+04
[340,] -2.101790e+04 5.100661e+03
[341,] 2.103565e+04 -2.101790e+04
[342,] 2.103565e+04 2.103565e+04
[343,] 2.103565e+04 2.103565e+04
[344,] -1.086091e+04 2.103565e+04
[345,] 2.103565e+04 -1.086091e+04
[346,] 3.709479e+03 2.103565e+04
[347,] 2.103565e+04 3.709479e+03
[348,] 2.103565e+04 2.103565e+04
[349,] 3.621747e+03 2.103565e+04
[350,] 1.835330e+04 3.621747e+03
[351,] 6.992397e+03 1.835330e+04
[352,] 1.211169e+04 6.992397e+03
[353,] 1.781777e+04 1.211169e+04
[354,] 2.103565e+04 1.781777e+04
[355,] 2.103565e+04 2.103565e+04
[356,] -1.298819e+02 2.103565e+04
[357,] 2.103565e+04 -1.298819e+02
[358,] 2.103565e+04 2.103565e+04
[359,] 2.103565e+04 2.103565e+04
[360,] 1.385765e+04 2.103565e+04
[361,] 2.103565e+04 1.385765e+04
[362,] 2.103565e+04 2.103565e+04
[363,] -4.118696e+04 2.103565e+04
[364,] 5.475234e+03 -4.118696e+04
[365,] 2.103565e+04 5.475234e+03
[366,] 2.103565e+04 2.103565e+04
[367,] 2.103565e+04 2.103565e+04
[368,] 2.103565e+04 2.103565e+04
[369,] 8.173406e+03 2.103565e+04
[370,] 2.103565e+04 8.173406e+03
[371,] 1.551123e+04 2.103565e+04
[372,] 2.924101e+04 1.551123e+04
[373,] 7.351272e+03 2.924101e+04
[374,] 2.103565e+04 7.351272e+03
[375,] 2.103565e+04 2.103565e+04
[376,] -1.902126e+04 2.103565e+04
[377,] 1.416519e+04 -1.902126e+04
[378,] 2.794641e+03 1.416519e+04
[379,] 2.103565e+04 2.794641e+03
[380,] 2.103565e+04 2.103565e+04
[381,] -1.070157e+04 2.103565e+04
[382,] -4.042411e+04 -1.070157e+04
[383,] 2.103565e+04 -4.042411e+04
[384,] 1.231220e+03 2.103565e+04
[385,] -1.152632e+05 1.231220e+03
[386,] 2.708052e+04 -1.152632e+05
[387,] 2.046047e+03 2.708052e+04
[388,] 1.433238e+04 2.046047e+03
[389,] -1.990561e+03 1.433238e+04
[390,] 9.146820e+03 -1.990561e+03
[391,] -1.390456e+04 9.146820e+03
[392,] -4.674842e+04 -1.390456e+04
[393,] 6.112366e+04 -4.674842e+04
[394,] -3.384755e+04 6.112366e+04
[395,] -2.631279e+04 -3.384755e+04
[396,] -1.105763e+05 -2.631279e+04
[397,] -2.996338e+05 -1.105763e+05
[398,] -1.175435e+04 -2.996338e+05
[399,] -1.441830e+05 -1.175435e+04
[400,] 2.482578e+03 -1.441830e+05
[401,] 2.316071e+04 2.482578e+03
[402,] -9.445398e+03 2.316071e+04
[403,] -5.694448e+04 -9.445398e+03
[404,] -3.124614e+04 -5.694448e+04
[405,] 1.030738e+04 -3.124614e+04
[406,] 1.566919e+03 1.030738e+04
[407,] -1.355572e+04 1.566919e+03
[408,] -1.335913e+04 -1.355572e+04
[409,] -6.926762e+04 -1.335913e+04
[410,] -9.224149e+04 -6.926762e+04
[411,] -6.372845e+04 -9.224149e+04
[412,] -5.425446e+04 -6.372845e+04
[413,] 8.686138e+04 -5.425446e+04
[414,] -1.529012e+05 8.686138e+04
[415,] -9.543057e+04 -1.529012e+05
[416,] 6.308139e+04 -9.543057e+04
[417,] -4.264355e+04 6.308139e+04
[418,] -6.550312e+04 -4.264355e+04
[419,] -5.082283e+04 -6.550312e+04
[420,] -4.454602e+04 -5.082283e+04
[421,] -1.986108e+05 -4.454602e+04
[422,] -5.465750e+05 -1.986108e+05
[423,] -3.885853e+05 -5.465750e+05
[424,] -6.319665e+04 -3.885853e+05
[425,] -2.494920e+05 -6.319665e+04
[426,] -2.993732e+05 -2.494920e+05
[427,] -7.770830e+05 -2.993732e+05
[428,] -1.248715e+05 -7.770830e+05
[429,] -1.470718e+05 -1.248715e+05
[430,] -4.993672e+05 -1.470718e+05
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 2.739587e+06 -6.791538e+04
2 1.732806e+06 2.739587e+06
3 -2.186068e+06 1.732806e+06
4 -3.431192e+04 -2.186068e+06
5 7.204795e+05 -3.431192e+04
6 2.034611e+05 7.204795e+05
7 6.736496e+05 2.034611e+05
8 -7.155026e+04 6.736496e+05
9 8.187372e+05 -7.155026e+04
10 4.565894e+05 8.187372e+05
11 3.973043e+05 4.565894e+05
12 1.388178e+05 3.973043e+05
13 -7.598112e+05 1.388178e+05
14 -1.420388e+05 -7.598112e+05
15 6.381048e+05 -1.420388e+05
16 1.242467e+05 6.381048e+05
17 -1.449007e+05 1.242467e+05
18 -7.999024e+04 -1.449007e+05
19 -5.425513e+05 -7.999024e+04
20 3.768709e+05 -5.425513e+05
21 2.074633e+05 3.768709e+05
22 1.535217e+05 2.074633e+05
23 -6.409222e+05 1.535217e+05
24 -1.780921e+05 -6.409222e+05
25 1.328343e+05 -1.780921e+05
26 -3.255006e+05 1.328343e+05
27 -2.663893e+05 -3.255006e+05
28 4.060816e+05 -2.663893e+05
29 -1.965210e+05 4.060816e+05
30 3.386834e+05 -1.965210e+05
31 1.054774e+05 3.386834e+05
32 3.532760e+05 1.054774e+05
33 2.134773e+05 3.532760e+05
34 3.233606e+05 2.134773e+05
35 -4.582349e+05 3.233606e+05
36 1.689982e+05 -4.582349e+05
37 2.355283e+05 1.689982e+05
38 9.624783e+04 2.355283e+05
39 -3.474496e+05 9.624783e+04
40 -4.310941e+05 -3.474496e+05
41 -9.940253e+04 -4.310941e+05
42 2.589757e+04 -9.940253e+04
43 -1.049501e+05 2.589757e+04
44 8.662826e+03 -1.049501e+05
45 2.827525e+04 8.662826e+03
46 -2.838113e+04 2.827525e+04
47 -1.817006e+05 -2.838113e+04
48 5.221901e+04 -1.817006e+05
49 -1.808105e+05 5.221901e+04
50 5.938359e+04 -1.808105e+05
51 -4.664730e+04 5.938359e+04
52 -1.392787e+04 -4.664730e+04
53 1.149660e+05 -1.392787e+04
54 -2.914968e+05 1.149660e+05
55 -3.721347e+04 -2.914968e+05
56 3.402509e+05 -3.721347e+04
57 4.929928e+04 3.402509e+05
58 -4.921277e+04 4.929928e+04
59 -1.024846e+05 -4.921277e+04
60 1.037701e+05 -1.024846e+05
61 -2.097813e+05 1.037701e+05
62 6.692660e+04 -2.097813e+05
63 -1.013692e+05 6.692660e+04
64 -4.588457e+04 -1.013692e+05
65 -5.213330e+05 -4.588457e+04
66 -6.213626e+04 -5.213330e+05
67 1.697286e+04 -6.213626e+04
68 -5.669127e+04 1.697286e+04
69 -4.091216e+04 -5.669127e+04
70 -1.908425e+05 -4.091216e+04
71 1.596308e+02 -1.908425e+05
72 -5.076590e+04 1.596308e+02
73 1.403811e+04 -5.076590e+04
74 3.322705e+04 1.403811e+04
75 -8.761648e+04 3.322705e+04
76 1.877659e+04 -8.761648e+04
77 -7.912928e+04 1.877659e+04
78 -7.999789e+04 -7.912928e+04
79 2.061823e+04 -7.999789e+04
80 -2.934650e+03 2.061823e+04
81 -5.392035e+03 -2.934650e+03
82 -3.384267e+03 -5.392035e+03
83 -2.131900e+04 -3.384267e+03
84 5.500812e+04 -2.131900e+04
85 -8.237319e+04 5.500812e+04
86 -3.442450e+04 -8.237319e+04
87 3.529264e+04 -3.442450e+04
88 -7.355256e+04 3.529264e+04
89 6.539898e+04 -7.355256e+04
90 2.825534e+04 6.539898e+04
91 2.998713e+03 2.825534e+04
92 2.670530e+04 2.998713e+03
93 2.612294e+02 2.670530e+04
94 -1.912698e+04 2.612294e+02
95 1.978783e+04 -1.912698e+04
96 1.180071e+05 1.978783e+04
97 -6.952977e+03 1.180071e+05
98 6.879170e+04 -6.952977e+03
99 4.442967e+03 6.879170e+04
100 -4.614226e+04 4.442967e+03
101 -9.576060e+04 -4.614226e+04
102 -1.062081e+05 -9.576060e+04
103 4.144130e+03 -1.062081e+05
104 1.887356e+04 4.144130e+03
105 -7.890440e+04 1.887356e+04
106 -2.661074e+05 -7.890440e+04
107 1.829283e+05 -2.661074e+05
108 -7.612731e+04 1.829283e+05
109 -8.519118e+04 -7.612731e+04
110 7.259471e+04 -8.519118e+04
111 -2.051650e+05 7.259471e+04
112 -9.805053e+04 -2.051650e+05
113 -3.358455e+04 -9.805053e+04
114 7.011974e+04 -3.358455e+04
115 -1.828131e+04 7.011974e+04
116 2.317555e+04 -1.828131e+04
117 1.126503e+05 2.317555e+04
118 -1.047734e+05 1.126503e+05
119 -1.852112e+04 -1.047734e+05
120 -1.357783e+05 -1.852112e+04
121 1.642105e+04 -1.357783e+05
122 -1.486211e+05 1.642105e+04
123 1.367024e+05 -1.486211e+05
124 2.185824e+04 1.367024e+05
125 -1.269823e+05 2.185824e+04
126 2.103565e+04 -1.269823e+05
127 2.002088e+03 2.103565e+04
128 -1.873865e+04 2.002088e+03
129 2.103565e+04 -1.873865e+04
130 7.349548e+03 2.103565e+04
131 2.603461e+04 7.349548e+03
132 2.103565e+04 2.603461e+04
133 -1.757489e+04 2.103565e+04
134 9.281619e+03 -1.757489e+04
135 1.285575e+04 9.281619e+03
136 6.327089e+03 1.285575e+04
137 1.873490e+04 6.327089e+03
138 1.162930e+04 1.873490e+04
139 2.103565e+04 1.162930e+04
140 1.596555e+04 2.103565e+04
141 2.103565e+04 1.596555e+04
142 -1.774000e+04 2.103565e+04
143 1.422555e+04 -1.774000e+04
144 7.400699e+03 1.422555e+04
145 1.695051e+03 7.400699e+03
146 2.103565e+04 1.695051e+03
147 8.034654e+03 2.103565e+04
148 2.679814e+04 8.034654e+03
149 2.150725e+04 2.679814e+04
150 -1.634696e+04 2.150725e+04
151 -7.698966e+04 -1.634696e+04
152 1.204308e+04 -7.698966e+04
153 -1.010521e+04 1.204308e+04
154 1.726602e+04 -1.010521e+04
155 2.103565e+04 1.726602e+04
156 2.103565e+04 2.103565e+04
157 1.711038e+04 2.103565e+04
158 3.228677e+04 1.711038e+04
159 2.103565e+04 3.228677e+04
160 2.103565e+04 2.103565e+04
161 2.103565e+04 2.103565e+04
162 2.103565e+04 2.103565e+04
163 -1.653319e+03 2.103565e+04
164 -3.650331e+04 -1.653319e+03
165 1.865232e+04 -3.650331e+04
166 2.103565e+04 1.865232e+04
167 5.417336e+04 2.103565e+04
168 -6.854867e+04 5.417336e+04
169 2.103565e+04 -6.854867e+04
170 2.103565e+04 2.103565e+04
171 1.260324e+04 2.103565e+04
172 2.103565e+04 1.260324e+04
173 -1.363295e+03 2.103565e+04
174 2.103565e+04 -1.363295e+03
175 6.552824e+03 2.103565e+04
176 -1.865200e+05 6.552824e+03
177 2.170283e+04 -1.865200e+05
178 2.103565e+04 2.170283e+04
179 2.103565e+04 2.103565e+04
180 2.844871e+04 2.103565e+04
181 1.478473e+04 2.844871e+04
182 -3.150429e+04 1.478473e+04
183 -2.659205e+04 -3.150429e+04
184 1.142357e+04 -2.659205e+04
185 1.769393e+04 1.142357e+04
186 1.048717e+04 1.769393e+04
187 7.896868e+03 1.048717e+04
188 -1.199707e+04 7.896868e+03
189 2.103565e+04 -1.199707e+04
190 1.124888e+04 2.103565e+04
191 2.103565e+04 1.124888e+04
192 8.143813e+03 2.103565e+04
193 2.103565e+04 8.143813e+03
194 1.589455e+04 2.103565e+04
195 2.103565e+04 1.589455e+04
196 -6.836521e+04 2.103565e+04
197 2.103565e+04 -6.836521e+04
198 2.103565e+04 2.103565e+04
199 -6.008006e+04 2.103565e+04
200 2.103565e+04 -6.008006e+04
201 2.103565e+04 2.103565e+04
202 2.103565e+04 2.103565e+04
203 -1.541284e+04 2.103565e+04
204 8.457068e+03 -1.541284e+04
205 1.279829e+04 8.457068e+03
206 2.051418e+04 1.279829e+04
207 2.103565e+04 2.051418e+04
208 -5.240170e+02 2.103565e+04
209 6.701203e+03 -5.240170e+02
210 2.103565e+04 6.701203e+03
211 2.103565e+04 2.103565e+04
212 1.535010e+04 2.103565e+04
213 -2.436018e+04 1.535010e+04
214 9.153589e+01 -2.436018e+04
215 1.743862e+04 9.153589e+01
216 1.052239e+04 1.743862e+04
217 9.990486e+03 1.052239e+04
218 1.743862e+04 9.990486e+03
219 -2.747358e+04 1.743862e+04
220 -3.634792e+03 -2.747358e+04
221 2.103565e+04 -3.634792e+03
222 1.920559e+04 2.103565e+04
223 -1.931020e+04 1.920559e+04
224 1.939619e+04 -1.931020e+04
225 -1.969494e+04 1.939619e+04
226 3.891318e+03 -1.969494e+04
227 2.103565e+04 3.891318e+03
228 1.403666e+04 2.103565e+04
229 2.103565e+04 1.403666e+04
230 -8.881974e+03 2.103565e+04
231 9.287883e+03 -8.881974e+03
232 2.103565e+04 9.287883e+03
233 -1.531985e+04 2.103565e+04
234 2.103565e+04 -1.531985e+04
235 2.113902e+04 2.103565e+04
236 2.103565e+04 2.113902e+04
237 1.055096e+04 2.103565e+04
238 2.103565e+04 1.055096e+04
239 2.103565e+04 2.103565e+04
240 2.103565e+04 2.103565e+04
241 1.296730e+04 2.103565e+04
242 2.103565e+04 1.296730e+04
243 8.653990e+03 2.103565e+04
244 2.205947e+04 8.653990e+03
245 2.103565e+04 2.205947e+04
246 2.103565e+04 2.103565e+04
247 -1.263528e+04 2.103565e+04
248 2.103565e+04 -1.263528e+04
249 2.103565e+04 2.103565e+04
250 -1.426362e+04 2.103565e+04
251 2.103565e+04 -1.426362e+04
252 1.087066e+04 2.103565e+04
253 -6.021682e+04 1.087066e+04
254 -1.518544e+05 -6.021682e+04
255 2.103565e+04 -1.518544e+05
256 1.879902e+03 2.103565e+04
257 -4.361101e+04 1.879902e+03
258 2.103565e+04 -4.361101e+04
259 2.103565e+04 2.103565e+04
260 2.103565e+04 2.103565e+04
261 2.103565e+04 2.103565e+04
262 2.569742e+04 2.103565e+04
263 1.896767e+04 2.569742e+04
264 2.103565e+04 1.896767e+04
265 2.103565e+04 2.103565e+04
266 2.103565e+04 2.103565e+04
267 7.025235e+04 2.103565e+04
268 2.103565e+04 7.025235e+04
269 -1.005395e+04 2.103565e+04
270 2.103565e+04 -1.005395e+04
271 2.103565e+04 2.103565e+04
272 2.103565e+04 2.103565e+04
273 -1.861992e+04 2.103565e+04
274 2.103565e+04 -1.861992e+04
275 3.466686e+03 2.103565e+04
276 -1.617307e+04 3.466686e+03
277 2.103565e+04 -1.617307e+04
278 -5.832370e+04 2.103565e+04
279 1.571969e+03 -5.832370e+04
280 2.103565e+04 1.571969e+03
281 2.103565e+04 2.103565e+04
282 -2.827914e+04 2.103565e+04
283 2.103565e+04 -2.827914e+04
284 2.008941e+04 2.103565e+04
285 2.103565e+04 2.008941e+04
286 1.468397e+04 2.103565e+04
287 1.830072e+04 1.468397e+04
288 2.103565e+04 1.830072e+04
289 -2.552913e+04 2.103565e+04
290 2.103565e+04 -2.552913e+04
291 8.885081e+03 2.103565e+04
292 3.633827e+04 8.885081e+03
293 2.103565e+04 3.633827e+04
294 -1.003406e+04 2.103565e+04
295 2.963102e+04 -1.003406e+04
296 2.103565e+04 2.963102e+04
297 2.103565e+04 2.103565e+04
298 2.103565e+04 2.103565e+04
299 4.849025e+03 2.103565e+04
300 1.156722e+04 4.849025e+03
301 4.485463e+03 1.156722e+04
302 1.923713e+04 4.485463e+03
303 2.103565e+04 1.923713e+04
304 2.406910e+04 2.103565e+04
305 -5.369338e+02 2.406910e+04
306 2.103565e+04 -5.369338e+02
307 -1.237875e+04 2.103565e+04
308 -1.599422e+04 -1.237875e+04
309 2.103565e+04 -1.599422e+04
310 2.103565e+04 2.103565e+04
311 -6.331702e+03 2.103565e+04
312 2.103565e+04 -6.331702e+03
313 2.396020e+03 2.103565e+04
314 6.917307e+03 2.396020e+03
315 3.764782e+04 6.917307e+03
316 2.056159e+04 3.764782e+04
317 -3.408307e+04 2.056159e+04
318 1.203310e+04 -3.408307e+04
319 2.103565e+04 1.203310e+04
320 2.103565e+04 2.103565e+04
321 2.103565e+04 2.103565e+04
322 2.103565e+04 2.103565e+04
323 2.103565e+04 2.103565e+04
324 2.103565e+04 2.103565e+04
325 1.631850e+04 2.103565e+04
326 2.103565e+04 1.631850e+04
327 2.103565e+04 2.103565e+04
328 2.103565e+04 2.103565e+04
329 2.103565e+04 2.103565e+04
330 2.103565e+04 2.103565e+04
331 2.103565e+04 2.103565e+04
332 2.103565e+04 2.103565e+04
333 1.200381e+04 2.103565e+04
334 1.924329e+04 1.200381e+04
335 1.322164e+04 1.924329e+04
336 -3.626141e+04 1.322164e+04
337 2.103565e+04 -3.626141e+04
338 1.768911e+04 2.103565e+04
339 5.100661e+03 1.768911e+04
340 -2.101790e+04 5.100661e+03
341 2.103565e+04 -2.101790e+04
342 2.103565e+04 2.103565e+04
343 2.103565e+04 2.103565e+04
344 -1.086091e+04 2.103565e+04
345 2.103565e+04 -1.086091e+04
346 3.709479e+03 2.103565e+04
347 2.103565e+04 3.709479e+03
348 2.103565e+04 2.103565e+04
349 3.621747e+03 2.103565e+04
350 1.835330e+04 3.621747e+03
351 6.992397e+03 1.835330e+04
352 1.211169e+04 6.992397e+03
353 1.781777e+04 1.211169e+04
354 2.103565e+04 1.781777e+04
355 2.103565e+04 2.103565e+04
356 -1.298819e+02 2.103565e+04
357 2.103565e+04 -1.298819e+02
358 2.103565e+04 2.103565e+04
359 2.103565e+04 2.103565e+04
360 1.385765e+04 2.103565e+04
361 2.103565e+04 1.385765e+04
362 2.103565e+04 2.103565e+04
363 -4.118696e+04 2.103565e+04
364 5.475234e+03 -4.118696e+04
365 2.103565e+04 5.475234e+03
366 2.103565e+04 2.103565e+04
367 2.103565e+04 2.103565e+04
368 2.103565e+04 2.103565e+04
369 8.173406e+03 2.103565e+04
370 2.103565e+04 8.173406e+03
371 1.551123e+04 2.103565e+04
372 2.924101e+04 1.551123e+04
373 7.351272e+03 2.924101e+04
374 2.103565e+04 7.351272e+03
375 2.103565e+04 2.103565e+04
376 -1.902126e+04 2.103565e+04
377 1.416519e+04 -1.902126e+04
378 2.794641e+03 1.416519e+04
379 2.103565e+04 2.794641e+03
380 2.103565e+04 2.103565e+04
381 -1.070157e+04 2.103565e+04
382 -4.042411e+04 -1.070157e+04
383 2.103565e+04 -4.042411e+04
384 1.231220e+03 2.103565e+04
385 -1.152632e+05 1.231220e+03
386 2.708052e+04 -1.152632e+05
387 2.046047e+03 2.708052e+04
388 1.433238e+04 2.046047e+03
389 -1.990561e+03 1.433238e+04
390 9.146820e+03 -1.990561e+03
391 -1.390456e+04 9.146820e+03
392 -4.674842e+04 -1.390456e+04
393 6.112366e+04 -4.674842e+04
394 -3.384755e+04 6.112366e+04
395 -2.631279e+04 -3.384755e+04
396 -1.105763e+05 -2.631279e+04
397 -2.996338e+05 -1.105763e+05
398 -1.175435e+04 -2.996338e+05
399 -1.441830e+05 -1.175435e+04
400 2.482578e+03 -1.441830e+05
401 2.316071e+04 2.482578e+03
402 -9.445398e+03 2.316071e+04
403 -5.694448e+04 -9.445398e+03
404 -3.124614e+04 -5.694448e+04
405 1.030738e+04 -3.124614e+04
406 1.566919e+03 1.030738e+04
407 -1.355572e+04 1.566919e+03
408 -1.335913e+04 -1.355572e+04
409 -6.926762e+04 -1.335913e+04
410 -9.224149e+04 -6.926762e+04
411 -6.372845e+04 -9.224149e+04
412 -5.425446e+04 -6.372845e+04
413 8.686138e+04 -5.425446e+04
414 -1.529012e+05 8.686138e+04
415 -9.543057e+04 -1.529012e+05
416 6.308139e+04 -9.543057e+04
417 -4.264355e+04 6.308139e+04
418 -6.550312e+04 -4.264355e+04
419 -5.082283e+04 -6.550312e+04
420 -4.454602e+04 -5.082283e+04
421 -1.986108e+05 -4.454602e+04
422 -5.465750e+05 -1.986108e+05
423 -3.885853e+05 -5.465750e+05
424 -6.319665e+04 -3.885853e+05
425 -2.494920e+05 -6.319665e+04
426 -2.993732e+05 -2.494920e+05
427 -7.770830e+05 -2.993732e+05
428 -1.248715e+05 -7.770830e+05
429 -1.470718e+05 -1.248715e+05
430 -4.993672e+05 -1.470718e+05
> 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/rcomp/tmp/7bu1r1292955502.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/rcomp/tmp/8bu1r1292955502.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/rcomp/tmp/94l0c1292955502.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/rcomp/tmp/104l0c1292955502.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/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/www/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/rcomp/tmp/11iejv1292955503.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/rcomp/tmp/123ezj1292955503.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/rcomp/tmp/13h6xs1292955503.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/rcomp/tmp/143pey1292955503.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/rcomp/tmp/1567c31292955503.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/rcomp/tmp/162hau1292955503.tab")
+ }
>
> try(system("convert tmp/1xk311292955502.ps tmp/1xk311292955502.png",intern=TRUE))
character(0)
> try(system("convert tmp/28u241292955502.ps tmp/28u241292955502.png",intern=TRUE))
character(0)
> try(system("convert tmp/38u241292955502.ps tmp/38u241292955502.png",intern=TRUE))
character(0)
> try(system("convert tmp/48u241292955502.ps tmp/48u241292955502.png",intern=TRUE))
character(0)
> try(system("convert tmp/5jljo1292955502.ps tmp/5jljo1292955502.png",intern=TRUE))
character(0)
> try(system("convert tmp/6jljo1292955502.ps tmp/6jljo1292955502.png",intern=TRUE))
character(0)
> try(system("convert tmp/7bu1r1292955502.ps tmp/7bu1r1292955502.png",intern=TRUE))
character(0)
> try(system("convert tmp/8bu1r1292955502.ps tmp/8bu1r1292955502.png",intern=TRUE))
character(0)
> try(system("convert tmp/94l0c1292955502.ps tmp/94l0c1292955502.png",intern=TRUE))
character(0)
> try(system("convert tmp/104l0c1292955502.ps tmp/104l0c1292955502.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
11.530 1.710 13.283