R version 2.9.0 (2009-04-17)
Copyright (C) 2009 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> x <- array(list(162556
+ ,807
+ ,213118
+ ,6282154
+ ,56289
+ ,29790
+ ,444
+ ,81767
+ ,4321023
+ ,28328
+ ,87550
+ ,412
+ ,153198
+ ,4111912
+ ,17936
+ ,84738
+ ,428
+ ,-26007
+ ,223193
+ ,4145
+ ,54660
+ ,315
+ ,126942
+ ,1491348
+ ,3040
+ ,42634
+ ,168
+ ,157214
+ ,1629616
+ ,2964
+ ,40949
+ ,263
+ ,129352
+ ,1398893
+ ,2865
+ ,45187
+ ,267
+ ,234817
+ ,1926517
+ ,2854
+ ,37704
+ ,228
+ ,60448
+ ,983660
+ ,2167
+ ,16275
+ ,129
+ ,47818
+ ,1443586
+ ,1974
+ ,25830
+ ,104
+ ,245546
+ ,1073089
+ ,1910
+ ,12679
+ ,122
+ ,48020
+ ,984885
+ ,1871
+ ,18014
+ ,393
+ ,-1710
+ ,1405225
+ ,943
+ ,43556
+ ,190
+ ,32648
+ ,227132
+ ,929
+ ,24811
+ ,280
+ ,95350
+ ,929118
+ ,822
+ ,6575
+ ,63
+ ,151352
+ ,1071292
+ ,819
+ ,7123
+ ,102
+ ,288170
+ ,638830
+ ,769
+ ,21950
+ ,265
+ ,114337
+ ,856956
+ ,745
+ ,37597
+ ,234
+ ,37884
+ ,992426
+ ,652
+ ,17821
+ ,277
+ ,122844
+ ,444477
+ ,643
+ ,12988
+ ,73
+ ,82340
+ ,857217
+ ,601
+ ,22330
+ ,67
+ ,79801
+ ,711969
+ ,446
+ ,13326
+ ,103
+ ,165548
+ ,702380
+ ,436
+ ,16189
+ ,290
+ ,116384
+ ,358589
+ ,379
+ ,7146
+ ,83
+ ,134028
+ ,297978
+ ,305
+ ,15824
+ ,56
+ ,63838
+ ,585715
+ ,284
+ ,27664
+ ,236
+ ,74996
+ ,657954
+ ,247
+ ,11920
+ ,73
+ ,31080
+ ,209458
+ ,238
+ ,8568
+ ,34
+ ,32168
+ ,786690
+ ,223
+ ,14416
+ ,139
+ ,49857
+ ,439798
+ ,220
+ ,3369
+ ,26
+ ,87161
+ ,688779
+ ,217
+ ,11819
+ ,70
+ ,106113
+ ,574339
+ ,199
+ ,6984
+ ,40
+ ,80570
+ ,741409
+ ,195
+ ,4519
+ ,42
+ ,102129
+ ,597793
+ ,163
+ ,2220
+ ,12
+ ,301670
+ ,644190
+ ,154
+ ,18562
+ ,211
+ ,102313
+ ,377934
+ ,143
+ ,10327
+ ,74
+ ,88577
+ ,640273
+ ,135
+ ,5336
+ ,80
+ ,112477
+ ,697458
+ ,132
+ ,2365
+ ,83
+ ,191778
+ ,550608
+ ,120
+ ,4069
+ ,131
+ ,79804
+ ,207393
+ ,119
+ ,8636
+ ,203
+ ,128294
+ ,301607
+ ,108
+ ,13718
+ ,56
+ ,96448
+ ,345783
+ ,104
+ ,4525
+ ,89
+ ,93811
+ ,501749
+ ,101
+ ,6869
+ ,88
+ ,117520
+ ,379983
+ ,90
+ ,4628
+ ,39
+ ,69159
+ ,387475
+ ,85
+ ,3689
+ ,25
+ ,101792
+ ,377305
+ ,74
+ ,4891
+ ,49
+ ,210568
+ ,370837
+ ,73
+ ,7489
+ ,149
+ ,136996
+ ,430866
+ ,70
+ ,4901
+ ,58
+ ,121920
+ ,469107
+ ,67
+ ,2284
+ ,41
+ ,76403
+ ,194493
+ ,66
+ ,3160
+ ,90
+ ,108094
+ ,530670
+ ,66
+ ,4150
+ ,136
+ ,134759
+ ,518365
+ ,66
+ ,7285
+ ,97
+ ,188873
+ ,491303
+ ,59
+ ,1134
+ ,63
+ ,146216
+ ,527021
+ ,58
+ ,4658
+ ,114
+ ,156608
+ ,233773
+ ,58
+ ,2384
+ ,77
+ ,61348
+ ,405972
+ ,54
+ ,3748
+ ,6
+ ,50350
+ ,652925
+ ,53
+ ,5371
+ ,47
+ ,87720
+ ,446211
+ ,49
+ ,1285
+ ,51
+ ,99489
+ ,341340
+ ,49
+ ,9327
+ ,85
+ ,87419
+ ,387699
+ ,44
+ ,5565
+ ,43
+ ,94355
+ ,493408
+ ,39
+ ,1528
+ ,32
+ ,60326
+ ,146494
+ ,38
+ ,3122
+ ,25
+ ,94670
+ ,414462
+ ,38
+ ,7561
+ ,77
+ ,82425
+ ,364304
+ ,37
+ ,2675
+ ,54
+ ,59017
+ ,355178
+ ,36
+ ,13253
+ ,251
+ ,90829
+ ,357760
+ ,35
+ ,880
+ ,15
+ ,80791
+ ,261216
+ ,34
+ ,2053
+ ,44
+ ,100423
+ ,397144
+ ,33
+ ,1424
+ ,73
+ ,131116
+ ,374943
+ ,32
+ ,4036
+ ,85
+ ,100269
+ ,424898
+ ,32
+ ,3045
+ ,49
+ ,27330
+ ,202055
+ ,31
+ ,5119
+ ,38
+ ,39039
+ ,378525
+ ,31
+ ,1431
+ ,35
+ ,106885
+ ,310768
+ ,30
+ ,554
+ ,9
+ ,79285
+ ,325738
+ ,30
+ ,1975
+ ,34
+ ,118881
+ ,394510
+ ,30
+ ,1765
+ ,20
+ ,77623
+ ,247060
+ ,28
+ ,1012
+ ,29
+ ,114768
+ ,368078
+ ,26
+ ,810
+ ,11
+ ,74015
+ ,236761
+ ,26
+ ,1280
+ ,52
+ ,69465
+ ,312378
+ ,26
+ ,666
+ ,13
+ ,117869
+ ,339836
+ ,25
+ ,1380
+ ,29
+ ,60982
+ ,347385
+ ,25
+ ,4677
+ ,66
+ ,90131
+ ,426280
+ ,24
+ ,876
+ ,33
+ ,138971
+ ,352850
+ ,24
+ ,814
+ ,15
+ ,39625
+ ,301881
+ ,23
+ ,514
+ ,15
+ ,102725
+ ,377516
+ ,23
+ ,5692
+ ,68
+ ,64239
+ ,357312
+ ,22
+ ,3642
+ ,100
+ ,90262
+ ,458343
+ ,22
+ ,540
+ ,13
+ ,103960
+ ,354228
+ ,22
+ ,2099
+ ,45
+ ,106611
+ ,308636
+ ,20
+ ,567
+ ,14
+ ,103345
+ ,386212
+ ,20
+ ,2001
+ ,36
+ ,95551
+ ,393343
+ ,19
+ ,2949
+ ,40
+ ,82903
+ ,378509
+ ,18
+ ,2253
+ ,68
+ ,63593
+ ,452469
+ ,17
+ ,6533
+ ,29
+ ,126910
+ ,364839
+ ,16
+ ,1889
+ ,43
+ ,37527
+ ,358649
+ ,14
+ ,3055
+ ,30
+ ,60247
+ ,376641
+ ,14
+ ,272
+ ,9
+ ,112995
+ ,429112
+ ,14
+ ,1414
+ ,22
+ ,70184
+ ,330546
+ ,13
+ ,2564
+ ,19
+ ,130140
+ ,403560
+ ,13
+ ,1383
+ ,9
+ ,73221
+ ,317892
+ ,13
+ ,1261
+ ,31
+ ,76114
+ ,307528
+ ,13
+ ,975
+ ,19
+ ,90534
+ ,235133
+ ,13
+ ,3366
+ ,55
+ ,108479
+ ,299243
+ ,12
+ ,576
+ ,8
+ ,113761
+ ,314073
+ ,12
+ ,1686
+ ,28
+ ,68696
+ ,368186
+ ,11
+ ,746
+ ,29
+ ,71561
+ ,269661
+ ,10
+ ,3192
+ ,48
+ ,59831
+ ,125390
+ ,10
+ ,2045
+ ,16
+ ,97890
+ ,510834
+ ,10
+ ,5702
+ ,47
+ ,101481
+ ,321896
+ ,10
+ ,1932
+ ,20
+ ,72954
+ ,249898
+ ,10
+ ,936
+ ,22
+ ,67939
+ ,408881
+ ,9
+ ,3437
+ ,33
+ ,48022
+ ,158492
+ ,9
+ ,5131
+ ,44
+ ,86111
+ ,292154
+ ,9
+ ,2397
+ ,13
+ ,74020
+ ,289513
+ ,6
+ ,1389
+ ,6
+ ,57530
+ ,378049
+ ,6
+ ,1503
+ ,35
+ ,56364
+ ,343466
+ ,6
+ ,402
+ ,8
+ ,84990
+ ,332743
+ ,5
+ ,2239
+ ,17
+ ,88590
+ ,442882
+ ,5
+ ,2234
+ ,11
+ ,77200
+ ,214215
+ ,5
+ ,837
+ ,21
+ ,61262
+ ,315688
+ ,4
+ ,10579
+ ,92
+ ,110309
+ ,375195
+ ,4
+ ,875
+ ,12
+ ,67000
+ ,334280
+ ,3
+ ,1585
+ ,112
+ ,93099
+ ,355864
+ ,2
+ ,1659
+ ,25
+ ,107577
+ ,480382
+ ,2
+ ,2647
+ ,17
+ ,62920
+ ,353058
+ ,1
+ ,3294
+ ,23
+ ,75832
+ ,217193
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,94
+ ,10
+ ,60793
+ ,314533
+ ,0
+ ,422
+ ,23
+ ,57935
+ ,318056
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,34
+ ,7
+ ,60630
+ ,314353
+ ,0
+ ,1558
+ ,25
+ ,55637
+ ,369448
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,43
+ ,20
+ ,60887
+ ,312846
+ ,0
+ ,645
+ ,4
+ ,60720
+ ,312075
+ ,0
+ ,316
+ ,4
+ ,60505
+ ,315009
+ ,0
+ ,115
+ ,10
+ ,60945
+ ,318903
+ ,0
+ ,5
+ ,1
+ ,60720
+ ,314887
+ ,0
+ ,897
+ ,4
+ ,60720
+ ,314913
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,389
+ ,8
+ ,58990
+ ,325506
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1002
+ ,11
+ ,56750
+ ,298568
+ ,0
+ ,36
+ ,4
+ ,60894
+ ,315834
+ ,0
+ ,460
+ ,15
+ ,63346
+ ,329784
+ ,0
+ ,309
+ ,9
+ ,56535
+ ,312878
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,9
+ ,7
+ ,60835
+ ,314987
+ ,0
+ ,271
+ ,2
+ ,60720
+ ,325249
+ ,0
+ ,14
+ ,0
+ ,61016
+ ,315877
+ ,0
+ ,520
+ ,7
+ ,58650
+ ,291650
+ ,0
+ ,1766
+ ,46
+ ,60438
+ ,305959
+ ,0
+ ,0
+ ,5
+ ,60720
+ ,315380
+ ,0
+ ,458
+ ,7
+ ,58625
+ ,297765
+ ,0
+ ,20
+ ,2
+ ,60938
+ ,315245
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,98
+ ,2
+ ,61490
+ ,315236
+ ,0
+ ,405
+ ,5
+ ,60845
+ ,336425
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,483
+ ,7
+ ,60830
+ ,306268
+ ,0
+ ,454
+ ,24
+ ,63261
+ ,302187
+ ,0
+ ,47
+ ,1
+ ,60720
+ ,314882
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,757
+ ,18
+ ,45689
+ ,382712
+ ,0
+ ,4655
+ ,55
+ ,60720
+ ,341570
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,36
+ ,3
+ ,61564
+ ,312412
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,203
+ ,9
+ ,61938
+ ,309596
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,126
+ ,8
+ ,60951
+ ,315547
+ ,0
+ ,400
+ ,113
+ ,60720
+ ,313267
+ ,0
+ ,71
+ ,0
+ ,60745
+ ,316176
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,972
+ ,19
+ ,71642
+ ,359335
+ ,0
+ ,531
+ ,11
+ ,71641
+ ,330068
+ ,0
+ ,2461
+ ,25
+ ,55792
+ ,314289
+ ,0
+ ,378
+ ,16
+ ,71873
+ ,297413
+ ,0
+ ,23
+ ,5
+ ,62555
+ ,314806
+ ,0
+ ,638
+ ,11
+ ,60370
+ ,333210
+ ,0
+ ,2300
+ ,23
+ ,64873
+ ,352108
+ ,0
+ ,149
+ ,6
+ ,62041
+ ,313332
+ ,0
+ ,226
+ ,5
+ ,65745
+ ,291787
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,275
+ ,7
+ ,59500
+ ,318745
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,141
+ ,7
+ ,61630
+ ,315366
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,28
+ ,3
+ ,60890
+ ,315688
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,4980
+ ,89
+ ,113521
+ ,409642
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,472
+ ,19
+ ,80045
+ ,269587
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,203
+ ,12
+ ,50804
+ ,300962
+ ,0
+ ,496
+ ,12
+ ,87390
+ ,325479
+ ,0
+ ,10
+ ,5
+ ,61656
+ ,316155
+ ,0
+ ,63
+ ,2
+ ,65688
+ ,318574
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1136
+ ,26
+ ,48522
+ ,343613
+ ,0
+ ,265
+ ,3
+ ,60720
+ ,306948
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,267
+ ,11
+ ,57640
+ ,330059
+ ,0
+ ,474
+ ,10
+ ,61977
+ ,288985
+ ,0
+ ,534
+ ,5
+ ,62620
+ ,304485
+ ,0
+ ,0
+ ,2
+ ,60720
+ ,315380
+ ,0
+ ,15
+ ,6
+ ,60831
+ ,315688
+ ,0
+ ,397
+ ,7
+ ,60646
+ ,317736
+ ,0
+ ,0
+ ,2
+ ,60720
+ ,315380
+ ,0
+ ,1866
+ ,28
+ ,56225
+ ,322331
+ ,0
+ ,288
+ ,3
+ ,60510
+ ,296656
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,3
+ ,1
+ ,60698
+ ,315354
+ ,0
+ ,468
+ ,20
+ ,60720
+ ,312161
+ ,0
+ ,20
+ ,1
+ ,60805
+ ,315576
+ ,0
+ ,278
+ ,22
+ ,61404
+ ,314922
+ ,0
+ ,61
+ ,9
+ ,60720
+ ,314551
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,192
+ ,2
+ ,65276
+ ,312339
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,317
+ ,7
+ ,63915
+ ,298700
+ ,0
+ ,738
+ ,9
+ ,60720
+ ,321376
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,368
+ ,13
+ ,61686
+ ,303230
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,2
+ ,0
+ ,60743
+ ,315487
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,53
+ ,6
+ ,60349
+ ,315793
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,94
+ ,3
+ ,61360
+ ,312887
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,24
+ ,7
+ ,59818
+ ,315637
+ ,0
+ ,2332
+ ,2
+ ,72680
+ ,324385
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,131
+ ,15
+ ,61808
+ ,308989
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,206
+ ,9
+ ,53110
+ ,296702
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,167
+ ,1
+ ,64245
+ ,307322
+ ,0
+ ,622
+ ,38
+ ,73007
+ ,304376
+ ,0
+ ,2328
+ ,57
+ ,82732
+ ,253588
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,365
+ ,7
+ ,54820
+ ,309560
+ ,0
+ ,364
+ ,26
+ ,47705
+ ,298466
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,226
+ ,13
+ ,72835
+ ,343929
+ ,0
+ ,307
+ ,10
+ ,58856
+ ,331955
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,188
+ ,9
+ ,77655
+ ,381180
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,138
+ ,26
+ ,69817
+ ,331420
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,125
+ ,19
+ ,60798
+ ,310201
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,282
+ ,12
+ ,62452
+ ,320016
+ ,0
+ ,335
+ ,23
+ ,64175
+ ,320398
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1324
+ ,29
+ ,67440
+ ,291841
+ ,0
+ ,176
+ ,8
+ ,68136
+ ,310670
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,249
+ ,26
+ ,56726
+ ,313491
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,333
+ ,9
+ ,70811
+ ,331323
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,601
+ ,5
+ ,60720
+ ,319210
+ ,0
+ ,30
+ ,3
+ ,62045
+ ,318098
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,249
+ ,13
+ ,54323
+ ,292754
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,165
+ ,12
+ ,62841
+ ,325176
+ ,0
+ ,453
+ ,19
+ ,81125
+ ,365959
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,53
+ ,10
+ ,59506
+ ,302409
+ ,0
+ ,382
+ ,9
+ ,59365
+ ,340968
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,9
+ ,60720
+ ,315380
+ ,0
+ ,30
+ ,4
+ ,60798
+ ,313164
+ ,0
+ ,290
+ ,1
+ ,58790
+ ,301164
+ ,0
+ ,0
+ ,1
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,366
+ ,14
+ ,61808
+ ,344425
+ ,0
+ ,2
+ ,12
+ ,60735
+ ,315394
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,209
+ ,19
+ ,64016
+ ,316647
+ ,0
+ ,384
+ ,17
+ ,54683
+ ,309836
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,365
+ ,32
+ ,87192
+ ,346611
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,49
+ ,14
+ ,64107
+ ,322031
+ ,0
+ ,3
+ ,8
+ ,60761
+ ,315656
+ ,0
+ ,133
+ ,4
+ ,65990
+ ,339445
+ ,0
+ ,32
+ ,0
+ ,59988
+ ,314964
+ ,0
+ ,368
+ ,20
+ ,61167
+ ,297141
+ ,0
+ ,1
+ ,5
+ ,60719
+ ,315372
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,22
+ ,1
+ ,60722
+ ,312502
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,96
+ ,4
+ ,60379
+ ,313729
+ ,0
+ ,1
+ ,1
+ ,60727
+ ,315388
+ ,0
+ ,314
+ ,4
+ ,60720
+ ,315371
+ ,0
+ ,844
+ ,20
+ ,60925
+ ,296139
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,26
+ ,1
+ ,60896
+ ,313880
+ ,0
+ ,125
+ ,10
+ ,59734
+ ,317698
+ ,0
+ ,304
+ ,12
+ ,62969
+ ,295580
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,621
+ ,13
+ ,60720
+ ,308256
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,119
+ ,3
+ ,59118
+ ,303677
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1595
+ ,10
+ ,60720
+ ,319369
+ ,0
+ ,312
+ ,3
+ ,58598
+ ,318690
+ ,0
+ ,60
+ ,7
+ ,61124
+ ,314049
+ ,0
+ ,587
+ ,10
+ ,59595
+ ,325699
+ ,0
+ ,135
+ ,1
+ ,62065
+ ,314210
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,514
+ ,15
+ ,78780
+ ,322378
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,4
+ ,60722
+ ,315398
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1763
+ ,28
+ ,61600
+ ,308336
+ ,0
+ ,180
+ ,9
+ ,59635
+ ,316386
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,218
+ ,7
+ ,60720
+ ,315553
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,448
+ ,7
+ ,59781
+ ,323361
+ ,0
+ ,227
+ ,7
+ ,76644
+ ,336639
+ ,0
+ ,174
+ ,3
+ ,64820
+ ,307424
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,121
+ ,11
+ ,56178
+ ,295370
+ ,0
+ ,607
+ ,7
+ ,60436
+ ,322340
+ ,0
+ ,2212
+ ,10
+ ,60720
+ ,319864
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,530
+ ,18
+ ,73433
+ ,317291
+ ,0
+ ,571
+ ,14
+ ,41477
+ ,280398
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,78
+ ,12
+ ,62700
+ ,317330
+ ,0
+ ,2489
+ ,29
+ ,67804
+ ,238125
+ ,0
+ ,131
+ ,3
+ ,59661
+ ,327071
+ ,0
+ ,923
+ ,6
+ ,58620
+ ,309038
+ ,0
+ ,72
+ ,3
+ ,60398
+ ,314210
+ ,0
+ ,572
+ ,8
+ ,58580
+ ,307930
+ ,-1
+ ,397
+ ,10
+ ,62710
+ ,322327
+ ,-1
+ ,450
+ ,6
+ ,59325
+ ,292136
+ ,-1
+ ,622
+ ,8
+ ,60950
+ ,263276
+ ,-2
+ ,694
+ ,6
+ ,68060
+ ,367655
+ ,-2
+ ,3425
+ ,9
+ ,83620
+ ,283910
+ ,-3
+ ,562
+ ,8
+ ,58456
+ ,283587
+ ,-3
+ ,4917
+ ,26
+ ,52811
+ ,243650
+ ,-3
+ ,1442
+ ,239
+ ,121173
+ ,438493
+ ,-3
+ ,529
+ ,7
+ ,63870
+ ,296261
+ ,-3
+ ,2126
+ ,41
+ ,21001
+ ,230621
+ ,-3
+ ,1061
+ ,3
+ ,70415
+ ,304252
+ ,-4
+ ,776
+ ,8
+ ,64230
+ ,333505
+ ,-4
+ ,611
+ ,6
+ ,59190
+ ,296919
+ ,-4
+ ,1526
+ ,21
+ ,69351
+ ,278990
+ ,-4
+ ,592
+ ,7
+ ,64270
+ ,276898
+ ,-4
+ ,1182
+ ,11
+ ,70694
+ ,327007
+ ,-5
+ ,621
+ ,11
+ ,68005
+ ,317046
+ ,-5
+ ,989
+ ,12
+ ,58930
+ ,304555
+ ,-5
+ ,438
+ ,9
+ ,58320
+ ,298096
+ ,-5
+ ,726
+ ,3
+ ,69980
+ ,231861
+ ,-5
+ ,1303
+ ,57
+ ,69863
+ ,309422
+ ,-5
+ ,7419
+ ,21
+ ,63255
+ ,286963
+ ,-5
+ ,1164
+ ,15
+ ,57320
+ ,269753
+ ,-6
+ ,3310
+ ,32
+ ,75230
+ ,448243
+ ,-6
+ ,1920
+ ,11
+ ,79420
+ ,165404
+ ,-6
+ ,965
+ ,2
+ ,73490
+ ,204325
+ ,-8
+ ,3256
+ ,23
+ ,35250
+ ,407159
+ ,-8
+ ,1135
+ ,20
+ ,62285
+ ,290476
+ ,-9
+ ,1270
+ ,24
+ ,69206
+ ,275311
+ ,-10
+ ,661
+ ,1
+ ,65920
+ ,246541
+ ,-12
+ ,1013
+ ,1
+ ,69770
+ ,253468
+ ,-14
+ ,2844
+ ,74
+ ,72683
+ ,240897
+ ,-15
+ ,11528
+ ,68
+ ,-14545
+ ,-83265
+ ,-25
+ ,6526
+ ,20
+ ,55830
+ ,-42143
+ ,-28
+ ,2264
+ ,20
+ ,55174
+ ,272713
+ ,-31
+ ,5109
+ ,82
+ ,67038
+ ,215362
+ ,-31
+ ,3999
+ ,21
+ ,51252
+ ,42754
+ ,-58
+ ,35624
+ ,244
+ ,157278
+ ,306275
+ ,-70
+ ,9252
+ ,32
+ ,79510
+ ,253537
+ ,-129
+ ,15236
+ ,86
+ ,77440
+ ,372631
+ ,-158
+ ,18073
+ ,69
+ ,27284
+ ,-7170
+ ,-366)
+ ,dim=c(5
+ ,431)
+ ,dimnames=list(c('kosten'
+ ,'orders'
+ ,'dividenden'
+ ,'rijkdom'
+ ,'CC-score')
+ ,1:431))
> y <- array(NA,dim=c(5,431),dimnames=list(c('kosten','orders','dividenden','rijkdom','CC-score'),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 = '5'
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Prof. Dr. P. Wessa
> #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #Technical description: Write here your technical program description (don't use hard returns!)
> library(lattice)
> library(lmtest)
Loading required package: zoo
Attaching package: 'zoo'
The following object(s) are masked from package:base :
as.Date.numeric
> n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test
> par1 <- as.numeric(par1)
> x <- t(y)
> k <- length(x[1,])
> n <- length(x[,1])
> x1 <- cbind(x[,par1], x[,1:k!=par1])
> mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1])
> colnames(x1) <- mycolnames #colnames(x)[par1]
> x <- x1
> if (par3 == 'First Differences'){
+ x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep='')))
+ for (i in 1:n-1) {
+ for (j in 1:k) {
+ x2[i,j] <- x[i+1,j] - x[i,j]
+ }
+ }
+ x <- x2
+ }
> if (par2 == 'Include Monthly Dummies'){
+ x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep ='')))
+ for (i in 1:11){
+ x2[seq(i,n,12),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> if (par2 == 'Include Quarterly Dummies'){
+ x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep ='')))
+ for (i in 1:3){
+ x2[seq(i,n,4),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> k <- length(x[1,])
> if (par3 == 'Linear Trend'){
+ x <- cbind(x, c(1:n))
+ colnames(x)[k+1] <- 't'
+ }
> x
CC-score kosten orders dividenden rijkdom
1 56289 162556 807 213118 6282154
2 28328 29790 444 81767 4321023
3 17936 87550 412 153198 4111912
4 4145 84738 428 -26007 223193
5 3040 54660 315 126942 1491348
6 2964 42634 168 157214 1629616
7 2865 40949 263 129352 1398893
8 2854 45187 267 234817 1926517
9 2167 37704 228 60448 983660
10 1974 16275 129 47818 1443586
11 1910 25830 104 245546 1073089
12 1871 12679 122 48020 984885
13 943 18014 393 -1710 1405225
14 929 43556 190 32648 227132
15 822 24811 280 95350 929118
16 819 6575 63 151352 1071292
17 769 7123 102 288170 638830
18 745 21950 265 114337 856956
19 652 37597 234 37884 992426
20 643 17821 277 122844 444477
21 601 12988 73 82340 857217
22 446 22330 67 79801 711969
23 436 13326 103 165548 702380
24 379 16189 290 116384 358589
25 305 7146 83 134028 297978
26 284 15824 56 63838 585715
27 247 27664 236 74996 657954
28 238 11920 73 31080 209458
29 223 8568 34 32168 786690
30 220 14416 139 49857 439798
31 217 3369 26 87161 688779
32 199 11819 70 106113 574339
33 195 6984 40 80570 741409
34 163 4519 42 102129 597793
35 154 2220 12 301670 644190
36 143 18562 211 102313 377934
37 135 10327 74 88577 640273
38 132 5336 80 112477 697458
39 120 2365 83 191778 550608
40 119 4069 131 79804 207393
41 108 8636 203 128294 301607
42 104 13718 56 96448 345783
43 101 4525 89 93811 501749
44 90 6869 88 117520 379983
45 85 4628 39 69159 387475
46 74 3689 25 101792 377305
47 73 4891 49 210568 370837
48 70 7489 149 136996 430866
49 67 4901 58 121920 469107
50 66 2284 41 76403 194493
51 66 3160 90 108094 530670
52 66 4150 136 134759 518365
53 59 7285 97 188873 491303
54 58 1134 63 146216 527021
55 58 4658 114 156608 233773
56 54 2384 77 61348 405972
57 53 3748 6 50350 652925
58 49 5371 47 87720 446211
59 49 1285 51 99489 341340
60 44 9327 85 87419 387699
61 39 5565 43 94355 493408
62 38 1528 32 60326 146494
63 38 3122 25 94670 414462
64 37 7561 77 82425 364304
65 36 2675 54 59017 355178
66 35 13253 251 90829 357760
67 34 880 15 80791 261216
68 33 2053 44 100423 397144
69 32 1424 73 131116 374943
70 32 4036 85 100269 424898
71 31 3045 49 27330 202055
72 31 5119 38 39039 378525
73 30 1431 35 106885 310768
74 30 554 9 79285 325738
75 30 1975 34 118881 394510
76 28 1765 20 77623 247060
77 26 1012 29 114768 368078
78 26 810 11 74015 236761
79 26 1280 52 69465 312378
80 25 666 13 117869 339836
81 25 1380 29 60982 347385
82 24 4677 66 90131 426280
83 24 876 33 138971 352850
84 23 814 15 39625 301881
85 23 514 15 102725 377516
86 22 5692 68 64239 357312
87 22 3642 100 90262 458343
88 22 540 13 103960 354228
89 20 2099 45 106611 308636
90 20 567 14 103345 386212
91 19 2001 36 95551 393343
92 18 2949 40 82903 378509
93 17 2253 68 63593 452469
94 16 6533 29 126910 364839
95 14 1889 43 37527 358649
96 14 3055 30 60247 376641
97 14 272 9 112995 429112
98 13 1414 22 70184 330546
99 13 2564 19 130140 403560
100 13 1383 9 73221 317892
101 13 1261 31 76114 307528
102 13 975 19 90534 235133
103 12 3366 55 108479 299243
104 12 576 8 113761 314073
105 11 1686 28 68696 368186
106 10 746 29 71561 269661
107 10 3192 48 59831 125390
108 10 2045 16 97890 510834
109 10 5702 47 101481 321896
110 10 1932 20 72954 249898
111 9 936 22 67939 408881
112 9 3437 33 48022 158492
113 9 5131 44 86111 292154
114 6 2397 13 74020 289513
115 6 1389 6 57530 378049
116 6 1503 35 56364 343466
117 5 402 8 84990 332743
118 5 2239 17 88590 442882
119 5 2234 11 77200 214215
120 4 837 21 61262 315688
121 4 10579 92 110309 375195
122 3 875 12 67000 334280
123 2 1585 112 93099 355864
124 2 1659 25 107577 480382
125 1 2647 17 62920 353058
126 1 3294 23 75832 217193
127 0 0 0 60720 315380
128 0 94 10 60793 314533
129 0 422 23 57935 318056
130 0 0 0 60720 315380
131 0 34 7 60630 314353
132 0 1558 25 55637 369448
133 0 0 0 60720 315380
134 0 43 20 60887 312846
135 0 645 4 60720 312075
136 0 316 4 60505 315009
137 0 115 10 60945 318903
138 0 5 1 60720 314887
139 0 897 4 60720 314913
140 0 0 0 60720 315380
141 0 389 8 58990 325506
142 0 0 0 60720 315380
143 0 1002 11 56750 298568
144 0 36 4 60894 315834
145 0 460 15 63346 329784
146 0 309 9 56535 312878
147 0 0 0 60720 315380
148 0 9 7 60835 314987
149 0 271 2 60720 325249
150 0 14 0 61016 315877
151 0 520 7 58650 291650
152 0 1766 46 60438 305959
153 0 0 5 60720 315380
154 0 458 7 58625 297765
155 0 20 2 60938 315245
156 0 0 0 60720 315380
157 0 0 0 60720 315380
158 0 98 2 61490 315236
159 0 405 5 60845 336425
160 0 0 0 60720 315380
161 0 0 0 60720 315380
162 0 0 0 60720 315380
163 0 0 0 60720 315380
164 0 483 7 60830 306268
165 0 454 24 63261 302187
166 0 47 1 60720 314882
167 0 0 0 60720 315380
168 0 757 18 45689 382712
169 0 4655 55 60720 341570
170 0 0 0 60720 315380
171 0 0 0 60720 315380
172 0 36 3 61564 312412
173 0 0 0 60720 315380
174 0 203 9 61938 309596
175 0 0 0 60720 315380
176 0 126 8 60951 315547
177 0 400 113 60720 313267
178 0 71 0 60745 316176
179 0 0 0 60720 315380
180 0 0 0 60720 315380
181 0 972 19 71642 359335
182 0 531 11 71641 330068
183 0 2461 25 55792 314289
184 0 378 16 71873 297413
185 0 23 5 62555 314806
186 0 638 11 60370 333210
187 0 2300 23 64873 352108
188 0 149 6 62041 313332
189 0 226 5 65745 291787
190 0 0 0 60720 315380
191 0 275 7 59500 318745
192 0 0 0 60720 315380
193 0 141 7 61630 315366
194 0 0 0 60720 315380
195 0 28 3 60890 315688
196 0 0 0 60720 315380
197 0 4980 89 113521 409642
198 0 0 0 60720 315380
199 0 0 0 60720 315380
200 0 472 19 80045 269587
201 0 0 0 60720 315380
202 0 0 0 60720 315380
203 0 0 0 60720 315380
204 0 203 12 50804 300962
205 0 496 12 87390 325479
206 0 10 5 61656 316155
207 0 63 2 65688 318574
208 0 0 0 60720 315380
209 0 1136 26 48522 343613
210 0 265 3 60720 306948
211 0 0 60720 315380 0
212 267 0 60720 315380 0
213 474 11 57640 330059 0
214 534 10 61977 288985 0
215 0 5 62620 304485 0
216 15 2 60720 315380 0
217 397 6 60831 315688 0
218 0 7 60646 317736 0
219 1866 2 60720 315380 0
220 288 28 56225 322331 0
221 0 3 60510 296656 0
222 3 0 60720 315380 0
223 468 1 60698 315354 0
224 20 20 60720 312161 0
225 278 1 60805 315576 0
226 61 22 61404 314922 0
227 0 9 60720 314551 0
228 192 0 60720 315380 0
229 0 2 65276 312339 0
230 317 0 60720 315380 0
231 738 7 63915 298700 0
232 0 9 60720 321376 0
233 368 0 60720 315380 0
234 0 13 61686 303230 0
235 2 0 60720 315380 0
236 0 0 60743 315487 0
237 53 0 60720 315380 0
238 0 6 60349 315793 0
239 0 0 60720 315380 0
240 0 0 60720 315380 0
241 94 0 60720 315380 0
242 0 3 61360 312887 0
243 24 0 60720 315380 0
244 2332 7 59818 315637 0
245 0 2 72680 324385 0
246 0 0 60720 315380 0
247 131 0 60720 315380 0
248 0 15 61808 308989 0
249 0 0 60720 315380 0
250 206 0 60720 315380 0
251 0 9 53110 296702 0
252 167 0 60720 315380 0
253 622 1 64245 307322 0
254 2328 38 73007 304376 0
255 0 57 82732 253588 0
256 365 0 60720 315380 0
257 364 7 54820 309560 0
258 0 26 47705 298466 0
259 0 0 60720 315380 0
260 0 0 60720 315380 0
261 0 0 60720 315380 0
262 226 0 60720 315380 0
263 307 13 72835 343929 0
264 0 10 58856 331955 0
265 0 0 60720 315380 0
266 0 0 60720 315380 0
267 188 0 60720 315380 0
268 0 9 77655 381180 0
269 138 0 60720 315380 0
270 0 26 69817 331420 0
271 0 0 60720 315380 0
272 0 0 60720 315380 0
273 125 0 60720 315380 0
274 0 19 60798 310201 0
275 282 0 60720 315380 0
276 335 12 62452 320016 0
277 0 23 64175 320398 0
278 1324 0 60720 315380 0
279 176 29 67440 291841 0
280 0 8 68136 310670 0
281 0 0 60720 315380 0
282 249 0 60720 315380 0
283 0 26 56726 313491 0
284 333 0 60720 315380 0
285 0 9 70811 331323 0
286 601 0 60720 315380 0
287 30 5 60720 319210 0
288 0 3 62045 318098 0
289 249 0 60720 315380 0
290 0 13 54323 292754 0
291 165 0 60720 315380 0
292 453 12 62841 325176 0
293 0 19 81125 365959 0
294 53 0 60720 315380 0
295 382 10 59506 302409 0
296 0 9 59365 340968 0
297 0 0 60720 315380 0
298 0 0 60720 315380 0
299 0 0 60720 315380 0
300 30 9 60720 315380 0
301 290 4 60798 313164 0
302 0 1 58790 301164 0
303 0 1 60720 315380 0
304 366 0 60720 315380 0
305 2 14 61808 344425 0
306 0 12 60735 315394 0
307 209 0 60720 315380 0
308 384 19 64016 316647 0
309 0 17 54683 309836 0
310 0 0 60720 315380 0
311 365 0 60720 315380 0
312 0 32 87192 346611 0
313 49 0 60720 315380 0
314 3 14 64107 322031 0
315 133 8 60761 315656 0
316 32 4 65990 339445 0
317 368 0 59988 314964 0
318 1 20 61167 297141 0
319 0 5 60719 315372 0
320 0 0 60720 315380 0
321 0 0 60720 315380 0
322 0 0 60720 315380 0
323 0 0 60720 315380 0
324 0 0 60720 315380 0
325 22 0 60720 315380 0
326 0 1 60722 312502 0
327 0 0 60720 315380 0
328 0 0 60720 315380 0
329 0 0 60720 315380 0
330 0 0 60720 315380 0
331 0 0 60720 315380 0
332 0 0 60720 315380 0
333 96 0 60720 315380 0
334 1 4 60379 313729 0
335 314 1 60727 315388 0
336 844 4 60720 315371 0
337 0 20 60925 296139 0
338 26 0 60720 315380 0
339 125 1 60896 313880 0
340 304 10 59734 317698 0
341 0 12 62969 295580 0
342 0 0 60720 315380 0
343 0 0 60720 315380 0
344 621 0 60720 315380 0
345 0 13 60720 308256 0
346 119 0 60720 315380 0
347 0 3 59118 303677 0
348 0 0 60720 315380 0
349 1595 0 60720 315380 0
350 312 10 60720 319369 0
351 60 3 58598 318690 0
352 587 7 61124 314049 0
353 135 10 59595 325699 0
354 0 1 62065 314210 0
355 0 0 60720 315380 0
356 514 0 60720 315380 0
357 0 15 78780 322378 0
358 0 0 60720 315380 0
359 0 0 60720 315380 0
360 1 0 60720 315380 0
361 0 4 60722 315398 0
362 0 0 60720 315380 0
363 1763 0 60720 315380 0
364 180 28 61600 308336 0
365 0 9 59635 316386 0
366 0 0 60720 315380 0
367 0 0 60720 315380 0
368 0 0 60720 315380 0
369 218 0 60720 315380 0
370 0 7 60720 315553 0
371 448 0 60720 315380 0
372 227 7 59781 323361 0
373 174 7 76644 336639 0
374 0 3 64820 307424 0
375 0 0 60720 315380 0
376 121 0 60720 315380 0
377 607 11 56178 295370 0
378 2212 7 60436 322340 0
379 0 10 60720 319864 0
380 0 0 60720 315380 0
381 530 0 60720 315380 0
382 571 18 73433 317291 0
383 0 14 41477 280398 0
384 78 0 60720 315380 0
385 2489 12 62700 317330 0
386 131 29 67804 238125 0
387 923 3 59661 327071 0
388 72 6 58620 309038 0
389 572 3 60398 314210 0
390 397 8 58580 307930 -1
391 450 10 62710 322327 -1
392 622 6 59325 292136 -1
393 694 8 60950 263276 -2
394 3425 6 68060 367655 -2
395 562 9 83620 283910 -3
396 4917 8 58456 283587 -3
397 1442 26 52811 243650 -3
398 529 239 121173 438493 -3
399 2126 7 63870 296261 -3
400 1061 41 21001 230621 -3
401 776 3 70415 304252 -4
402 611 8 64230 333505 -4
403 1526 6 59190 296919 -4
404 592 21 69351 278990 -4
405 1182 7 64270 276898 -4
406 621 11 70694 327007 -5
407 989 11 68005 317046 -5
408 438 12 58930 304555 -5
409 726 9 58320 298096 -5
410 1303 3 69980 231861 -5
411 7419 57 69863 309422 -5
412 1164 21 63255 286963 -5
413 3310 15 57320 269753 -6
414 1920 32 75230 448243 -6
415 965 11 79420 165404 -6
416 3256 2 73490 204325 -8
417 1135 23 35250 407159 -8
418 1270 20 62285 290476 -9
419 661 24 69206 275311 -10
420 1013 1 65920 246541 -12
421 2844 1 69770 253468 -14
422 11528 74 72683 240897 -15
423 6526 68 -14545 -83265 -25
424 2264 20 55830 -42143 -28
425 5109 20 55174 272713 -31
426 3999 82 67038 215362 -31
427 35624 21 51252 42754 -58
428 9252 244 157278 306275 -70
429 15236 32 79510 253537 -129
430 18073 86 77440 372631 -158
431 162556 69 27284 -7170 -366
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) kosten orders dividenden rijkdom
2.225e+03 6.075e-02 2.215e-01 -4.925e-02 5.793e-03
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-14401.8 -1061.7 -144.3 182.6 153933.4
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 2.225e+03 9.701e+02 2.294 0.022287 *
kosten 6.075e-02 5.822e-02 1.043 0.297355
orders 2.215e-01 3.277e-02 6.758 4.62e-11 ***
dividenden -4.925e-02 8.248e-03 -5.971 4.96e-09 ***
rijkdom 5.793e-03 1.557e-03 3.721 0.000225 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 7992 on 426 degrees of freedom
Multiple R-squared: 0.1639, Adjusted R-squared: 0.156
F-statistic: 20.88 on 4 and 426 DF, p-value: 9.929e-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,] 2.125167e-01 4.250334e-01 0.7874833
[2,] 1.298748e-01 2.597497e-01 0.8701252
[3,] 8.070540e-02 1.614108e-01 0.9192946
[4,] 1.132189e-01 2.264379e-01 0.8867811
[5,] 8.523411e-02 1.704682e-01 0.9147659
[6,] 7.627887e-02 1.525577e-01 0.9237211
[7,] 5.666892e-02 1.133378e-01 0.9433311
[8,] 3.196655e-02 6.393311e-02 0.9680334
[9,] 2.288814e-02 4.577629e-02 0.9771119
[10,] 2.338291e-02 4.676583e-02 0.9766171
[11,] 1.289508e-02 2.579017e-02 0.9871049
[12,] 7.005531e-03 1.401106e-02 0.9929945
[13,] 4.044243e-03 8.088486e-03 0.9959558
[14,] 2.606657e-03 5.213314e-03 0.9973933
[15,] 1.593433e-03 3.186867e-03 0.9984066
[16,] 9.500028e-04 1.900006e-03 0.9990500
[17,] 5.234383e-04 1.046877e-03 0.9994766
[18,] 4.890164e-04 9.780328e-04 0.9995110
[19,] 3.064900e-04 6.129800e-04 0.9996935
[20,] 1.488652e-04 2.977304e-04 0.9998511
[21,] 1.353729e-04 2.707459e-04 0.9998646
[22,] 6.953701e-05 1.390740e-04 0.9999305
[23,] 3.868887e-05 7.737774e-05 0.9999613
[24,] 2.109271e-05 4.218543e-05 0.9999789
[25,] 1.098999e-05 2.197998e-05 0.9999890
[26,] 5.324534e-06 1.064907e-05 0.9999947
[27,] 2.818646e-06 5.637291e-06 0.9999972
[28,] 1.829723e-06 3.659447e-06 0.9999982
[29,] 9.085901e-07 1.817180e-06 0.9999991
[30,] 4.203005e-07 8.406010e-07 0.9999996
[31,] 1.982056e-07 3.964112e-07 0.9999998
[32,] 1.082174e-07 2.164347e-07 0.9999999
[33,] 7.587511e-08 1.517502e-07 0.9999999
[34,] 3.899582e-08 7.799165e-08 1.0000000
[35,] 2.102492e-08 4.204985e-08 1.0000000
[36,] 1.008029e-08 2.016058e-08 1.0000000
[37,] 5.139158e-09 1.027832e-08 1.0000000
[38,] 2.736776e-09 5.473553e-09 1.0000000
[39,] 1.463929e-09 2.927857e-09 1.0000000
[40,] 7.785330e-10 1.557066e-09 1.0000000
[41,] 3.451712e-10 6.903424e-10 1.0000000
[42,] 1.583951e-10 3.167902e-10 1.0000000
[43,] 1.046328e-10 2.092655e-10 1.0000000
[44,] 4.589607e-11 9.179213e-11 1.0000000
[45,] 2.000633e-11 4.001267e-11 1.0000000
[46,] 9.019674e-12 1.803935e-11 1.0000000
[47,] 4.253908e-12 8.507817e-12 1.0000000
[48,] 2.162127e-12 4.324254e-12 1.0000000
[49,] 9.709263e-13 1.941853e-12 1.0000000
[50,] 4.053908e-13 8.107816e-13 1.0000000
[51,] 1.716025e-13 3.432049e-13 1.0000000
[52,] 8.134551e-14 1.626910e-13 1.0000000
[53,] 3.324879e-14 6.649758e-14 1.0000000
[54,] 1.348461e-14 2.696923e-14 1.0000000
[55,] 9.065586e-15 1.813117e-14 1.0000000
[56,] 3.821229e-15 7.642458e-15 1.0000000
[57,] 1.549156e-15 3.098312e-15 1.0000000
[58,] 6.562047e-16 1.312409e-15 1.0000000
[59,] 2.451638e-16 4.903276e-16 1.0000000
[60,] 1.197098e-16 2.394195e-16 1.0000000
[61,] 4.874148e-17 9.748295e-17 1.0000000
[62,] 2.054608e-17 4.109215e-17 1.0000000
[63,] 7.892697e-18 1.578539e-17 1.0000000
[64,] 4.283150e-18 8.566301e-18 1.0000000
[65,] 1.631164e-18 3.262328e-18 1.0000000
[66,] 6.865030e-19 1.373006e-18 1.0000000
[67,] 2.808978e-19 5.617956e-19 1.0000000
[68,] 1.110217e-19 2.220435e-19 1.0000000
[69,] 4.915851e-20 9.831702e-20 1.0000000
[70,] 1.940275e-20 3.880549e-20 1.0000000
[71,] 8.615550e-21 1.723110e-20 1.0000000
[72,] 3.329965e-21 6.659930e-21 1.0000000
[73,] 1.307677e-21 2.615353e-21 1.0000000
[74,] 4.759898e-22 9.519796e-22 1.0000000
[75,] 1.649766e-22 3.299531e-22 1.0000000
[76,] 6.549897e-23 1.309979e-22 1.0000000
[77,] 2.470359e-23 4.940717e-23 1.0000000
[78,] 9.032674e-24 1.806535e-23 1.0000000
[79,] 3.072709e-24 6.145417e-24 1.0000000
[80,] 1.063313e-24 2.126627e-24 1.0000000
[81,] 3.839309e-25 7.678617e-25 1.0000000
[82,] 1.391529e-25 2.783058e-25 1.0000000
[83,] 4.944555e-26 9.889109e-26 1.0000000
[84,] 1.688219e-26 3.376439e-26 1.0000000
[85,] 5.561028e-27 1.112206e-26 1.0000000
[86,] 1.816844e-27 3.633688e-27 1.0000000
[87,] 6.045227e-28 1.209045e-27 1.0000000
[88,] 1.940506e-28 3.881013e-28 1.0000000
[89,] 6.092513e-29 1.218503e-28 1.0000000
[90,] 2.245242e-29 4.490483e-29 1.0000000
[91,] 7.293385e-30 1.458677e-29 1.0000000
[92,] 2.642190e-30 5.284380e-30 1.0000000
[93,] 8.558370e-31 1.711674e-30 1.0000000
[94,] 2.805789e-31 5.611579e-31 1.0000000
[95,] 1.028237e-31 2.056474e-31 1.0000000
[96,] 3.398433e-32 6.796865e-32 1.0000000
[97,] 1.139655e-32 2.279310e-32 1.0000000
[98,] 3.453843e-33 6.907686e-33 1.0000000
[99,] 1.156283e-33 2.312565e-33 1.0000000
[100,] 6.609495e-34 1.321899e-33 1.0000000
[101,] 2.794897e-34 5.589794e-34 1.0000000
[102,] 8.499258e-35 1.699852e-34 1.0000000
[103,] 2.905298e-35 5.810596e-35 1.0000000
[104,] 8.901899e-36 1.780380e-35 1.0000000
[105,] 4.999616e-36 9.999232e-36 1.0000000
[106,] 1.575010e-36 3.150020e-36 1.0000000
[107,] 4.876073e-37 9.752145e-37 1.0000000
[108,] 1.393284e-37 2.786569e-37 1.0000000
[109,] 4.002322e-38 8.004643e-38 1.0000000
[110,] 1.191275e-38 2.382551e-38 1.0000000
[111,] 4.008940e-39 8.017881e-39 1.0000000
[112,] 1.508357e-39 3.016715e-39 1.0000000
[113,] 4.346884e-40 8.693768e-40 1.0000000
[114,] 1.210674e-40 2.421347e-40 1.0000000
[115,] 3.394383e-41 6.788767e-41 1.0000000
[116,] 1.024693e-41 2.049386e-41 1.0000000
[117,] 4.924799e-42 9.849597e-42 1.0000000
[118,] 1.335838e-42 2.671676e-42 1.0000000
[119,] 5.401411e-43 1.080282e-42 1.0000000
[120,] 1.496811e-43 2.993621e-43 1.0000000
[121,] 4.138401e-44 8.276802e-44 1.0000000
[122,] 1.135942e-44 2.271884e-44 1.0000000
[123,] 3.087154e-45 6.174308e-45 1.0000000
[124,] 8.365772e-46 1.673154e-45 1.0000000
[125,] 2.190263e-46 4.380527e-46 1.0000000
[126,] 5.843829e-47 1.168766e-46 1.0000000
[127,] 1.567106e-47 3.134211e-47 1.0000000
[128,] 4.141103e-48 8.282206e-48 1.0000000
[129,] 1.084649e-48 2.169299e-48 1.0000000
[130,] 2.828215e-49 5.656429e-49 1.0000000
[131,] 7.324404e-50 1.464881e-49 1.0000000
[132,] 1.884864e-50 3.769729e-50 1.0000000
[133,] 4.820098e-51 9.640196e-51 1.0000000
[134,] 1.217029e-51 2.434059e-51 1.0000000
[135,] 3.076106e-52 6.152212e-52 1.0000000
[136,] 8.021139e-53 1.604228e-52 1.0000000
[137,] 2.007301e-53 4.014603e-53 1.0000000
[138,] 4.979964e-54 9.959927e-54 1.0000000
[139,] 1.234700e-54 2.469401e-54 1.0000000
[140,] 3.033337e-55 6.066673e-55 1.0000000
[141,] 7.437712e-56 1.487542e-55 1.0000000
[142,] 1.795917e-56 3.591833e-56 1.0000000
[143,] 4.338172e-57 8.676345e-57 1.0000000
[144,] 1.085420e-57 2.170841e-57 1.0000000
[145,] 2.751911e-58 5.503821e-58 1.0000000
[146,] 6.558370e-59 1.311674e-58 1.0000000
[147,] 1.595341e-59 3.190682e-59 1.0000000
[148,] 3.755854e-60 7.511709e-60 1.0000000
[149,] 8.784980e-61 1.756996e-60 1.0000000
[150,] 2.043755e-61 4.087510e-61 1.0000000
[151,] 4.733893e-62 9.467786e-62 1.0000000
[152,] 1.082159e-62 2.164319e-62 1.0000000
[153,] 2.477840e-63 4.955681e-63 1.0000000
[154,] 5.643535e-64 1.128707e-63 1.0000000
[155,] 1.278599e-64 2.557198e-64 1.0000000
[156,] 2.881588e-65 5.763176e-65 1.0000000
[157,] 6.547938e-66 1.309588e-65 1.0000000
[158,] 1.506076e-66 3.012153e-66 1.0000000
[159,] 3.343327e-67 6.686653e-67 1.0000000
[160,] 7.379812e-68 1.475962e-67 1.0000000
[161,] 1.634993e-68 3.269986e-68 1.0000000
[162,] 4.199971e-69 8.399942e-69 1.0000000
[163,] 9.149807e-70 1.829961e-69 1.0000000
[164,] 1.983188e-70 3.966375e-70 1.0000000
[165,] 4.291285e-71 8.582570e-71 1.0000000
[166,] 9.207111e-72 1.841422e-71 1.0000000
[167,] 1.982459e-72 3.964918e-72 1.0000000
[168,] 4.210498e-73 8.420996e-73 1.0000000
[169,] 8.923368e-74 1.784674e-73 1.0000000
[170,] 2.064496e-74 4.128992e-74 1.0000000
[171,] 4.317101e-75 8.634201e-75 1.0000000
[172,] 8.992734e-76 1.798547e-75 1.0000000
[173,] 1.863978e-76 3.727956e-76 1.0000000
[174,] 3.923746e-77 7.847491e-77 1.0000000
[175,] 8.146689e-78 1.629338e-77 1.0000000
[176,] 1.962605e-78 3.925210e-78 1.0000000
[177,] 4.138857e-79 8.277715e-79 1.0000000
[178,] 8.424276e-80 1.684855e-79 1.0000000
[179,] 1.683556e-80 3.367113e-80 1.0000000
[180,] 3.436523e-81 6.873046e-81 1.0000000
[181,] 6.884869e-82 1.376974e-81 1.0000000
[182,] 1.415035e-82 2.830070e-82 1.0000000
[183,] 2.797711e-83 5.595422e-83 1.0000000
[184,] 5.493742e-84 1.098748e-83 1.0000000
[185,] 1.075712e-84 2.151424e-84 1.0000000
[186,] 2.101082e-85 4.202165e-85 1.0000000
[187,] 4.075116e-86 8.150232e-86 1.0000000
[188,] 7.875212e-87 1.575042e-86 1.0000000
[189,] 1.513430e-87 3.026860e-87 1.0000000
[190,] 3.504022e-88 7.008044e-88 1.0000000
[191,] 6.756597e-89 1.351319e-88 1.0000000
[192,] 1.301737e-89 2.603473e-89 1.0000000
[193,] 2.933585e-90 5.867171e-90 1.0000000
[194,] 5.619687e-91 1.123937e-90 1.0000000
[195,] 1.080354e-91 2.160709e-91 1.0000000
[196,] 2.096977e-92 4.193955e-92 1.0000000
[197,] 3.941754e-93 7.883507e-93 1.0000000
[198,] 7.454274e-94 1.490855e-93 1.0000000
[199,] 1.473384e-94 2.946768e-94 1.0000000
[200,] 3.152911e-95 6.305822e-95 1.0000000
[201,] 1.143386e-95 2.286772e-95 1.0000000
[202,] 1.572090e-94 3.144180e-94 1.0000000
[203,] 1.719520e-31 3.439040e-31 1.0000000
[204,] 3.179882e-31 6.359764e-31 1.0000000
[205,] 1.359689e-31 2.719378e-31 1.0000000
[206,] 6.022542e-32 1.204508e-31 1.0000000
[207,] 2.553133e-32 5.106265e-32 1.0000000
[208,] 1.101506e-32 2.203012e-32 1.0000000
[209,] 4.613961e-33 9.227921e-33 1.0000000
[210,] 1.933284e-33 3.866567e-33 1.0000000
[211,] 8.020659e-34 1.604132e-33 1.0000000
[212,] 4.720117e-34 9.440234e-34 1.0000000
[213,] 1.973781e-34 3.947562e-34 1.0000000
[214,] 8.337541e-35 1.667508e-34 1.0000000
[215,] 3.437695e-35 6.875391e-35 1.0000000
[216,] 1.399912e-35 2.799823e-35 1.0000000
[217,] 5.727277e-36 1.145455e-35 1.0000000
[218,] 2.300448e-36 4.600896e-36 1.0000000
[219,] 9.328884e-37 1.865777e-36 1.0000000
[220,] 3.745340e-37 7.490679e-37 1.0000000
[221,] 1.483174e-37 2.966349e-37 1.0000000
[222,] 6.147085e-38 1.229417e-37 1.0000000
[223,] 2.414678e-38 4.829357e-38 1.0000000
[224,] 9.593374e-39 1.918675e-38 1.0000000
[225,] 3.743117e-39 7.486234e-39 1.0000000
[226,] 1.452639e-39 2.905278e-39 1.0000000
[227,] 5.709466e-40 1.141893e-39 1.0000000
[228,] 2.199261e-40 4.398522e-40 1.0000000
[229,] 8.429253e-41 1.685851e-40 1.0000000
[230,] 3.206270e-41 6.412539e-41 1.0000000
[231,] 1.214409e-41 2.428818e-41 1.0000000
[232,] 4.584536e-42 9.169072e-42 1.0000000
[233,] 1.722015e-42 3.444031e-42 1.0000000
[234,] 6.414213e-43 1.282843e-42 1.0000000
[235,] 2.394097e-43 4.788194e-43 1.0000000
[236,] 8.849297e-44 1.769859e-43 1.0000000
[237,] 7.057370e-44 1.411474e-43 1.0000000
[238,] 3.093903e-44 6.187806e-44 1.0000000
[239,] 1.135660e-44 2.271320e-44 1.0000000
[240,] 4.130207e-45 8.260415e-45 1.0000000
[241,] 1.517328e-45 3.034656e-45 1.0000000
[242,] 5.486182e-46 1.097236e-45 1.0000000
[243,] 1.966769e-46 3.933537e-46 1.0000000
[244,] 7.167588e-47 1.433518e-46 1.0000000
[245,] 2.543980e-47 5.087961e-47 1.0000000
[246,] 9.127289e-48 1.825458e-47 1.0000000
[247,] 5.595379e-48 1.119076e-47 1.0000000
[248,] 4.996868e-48 9.993736e-48 1.0000000
[249,] 1.768589e-48 3.537179e-48 1.0000000
[250,] 6.498884e-49 1.299777e-48 1.0000000
[251,] 2.443329e-49 4.886658e-49 1.0000000
[252,] 8.486887e-50 1.697377e-49 1.0000000
[253,] 2.932990e-50 5.865980e-50 1.0000000
[254,] 1.008479e-50 2.016958e-50 1.0000000
[255,] 3.441039e-51 6.882079e-51 1.0000000
[256,] 1.330236e-51 2.660472e-51 1.0000000
[257,] 4.502382e-52 9.004763e-52 1.0000000
[258,] 1.517786e-52 3.035572e-52 1.0000000
[259,] 5.090485e-53 1.018097e-52 1.0000000
[260,] 1.694250e-53 3.388499e-53 1.0000000
[261,] 7.244327e-54 1.448865e-53 1.0000000
[262,] 2.388934e-54 4.777867e-54 1.0000000
[263,] 8.901083e-55 1.780217e-54 1.0000000
[264,] 2.912523e-55 5.825045e-55 1.0000000
[265,] 9.481096e-56 1.896219e-55 1.0000000
[266,] 3.065366e-56 6.130732e-56 1.0000000
[267,] 9.926118e-57 1.985224e-56 1.0000000
[268,] 3.195349e-57 6.390697e-57 1.0000000
[269,] 1.036771e-57 2.073542e-57 1.0000000
[270,] 3.404980e-58 6.809960e-58 1.0000000
[271,] 1.408534e-58 2.817067e-58 1.0000000
[272,] 4.872769e-59 9.745538e-59 1.0000000
[273,] 1.667587e-59 3.335174e-59 1.0000000
[274,] 5.193284e-60 1.038657e-59 1.0000000
[275,] 1.613981e-60 3.227962e-60 1.0000000
[276,] 4.991982e-61 9.983964e-61 1.0000000
[277,] 1.544566e-61 3.089132e-61 1.0000000
[278,] 5.353445e-62 1.070689e-61 1.0000000
[279,] 1.705631e-62 3.411263e-62 1.0000000
[280,] 5.149150e-63 1.029830e-62 1.0000000
[281,] 1.554649e-63 3.109299e-63 1.0000000
[282,] 4.661483e-64 9.322965e-64 1.0000000
[283,] 1.400195e-64 2.800391e-64 1.0000000
[284,] 4.136964e-65 8.273927e-65 1.0000000
[285,] 1.267296e-65 2.534592e-65 1.0000000
[286,] 5.972430e-66 1.194486e-65 1.0000000
[287,] 1.738281e-66 3.476561e-66 1.0000000
[288,] 5.118820e-67 1.023764e-66 1.0000000
[289,] 1.487909e-67 2.975818e-67 1.0000000
[290,] 4.264271e-68 8.528542e-68 1.0000000
[291,] 1.215266e-68 2.430531e-68 1.0000000
[292,] 3.443857e-69 6.887715e-69 1.0000000
[293,] 9.720773e-70 1.944155e-69 1.0000000
[294,] 2.748689e-70 5.497378e-70 1.0000000
[295,] 7.661794e-71 1.532359e-70 1.0000000
[296,] 2.122818e-71 4.245636e-71 1.0000000
[297,] 5.946448e-72 1.189290e-71 1.0000000
[298,] 1.658230e-72 3.316459e-72 1.0000000
[299,] 4.532656e-73 9.065311e-73 1.0000000
[300,] 1.231090e-73 2.462179e-73 1.0000000
[301,] 3.467139e-74 6.934278e-74 1.0000000
[302,] 9.405537e-75 1.881107e-74 1.0000000
[303,] 2.503292e-75 5.006585e-75 1.0000000
[304,] 6.737848e-76 1.347570e-75 1.0000000
[305,] 4.635612e-76 9.271223e-76 1.0000000
[306,] 1.215673e-76 2.431347e-76 1.0000000
[307,] 3.279483e-77 6.558966e-77 1.0000000
[308,] 8.548824e-78 1.709765e-77 1.0000000
[309,] 2.303353e-78 4.606706e-78 1.0000000
[310,] 6.052766e-79 1.210553e-78 1.0000000
[311,] 1.563370e-79 3.126740e-79 1.0000000
[312,] 3.962455e-80 7.924909e-80 1.0000000
[313,] 9.964648e-81 1.992930e-80 1.0000000
[314,] 2.490101e-81 4.980202e-81 1.0000000
[315,] 6.183183e-82 1.236637e-81 1.0000000
[316,] 1.525562e-82 3.051124e-82 1.0000000
[317,] 3.739836e-83 7.479671e-83 1.0000000
[318,] 9.106880e-84 1.821376e-83 1.0000000
[319,] 2.204781e-84 4.409562e-84 1.0000000
[320,] 5.300201e-85 1.060040e-84 1.0000000
[321,] 1.265751e-85 2.531503e-85 1.0000000
[322,] 3.002711e-86 6.005422e-86 1.0000000
[323,] 7.075667e-87 1.415133e-86 1.0000000
[324,] 1.656109e-87 3.312218e-87 1.0000000
[325,] 3.849965e-88 7.699931e-88 1.0000000
[326,] 8.897520e-89 1.779504e-88 1.0000000
[327,] 2.040650e-89 4.081300e-89 1.0000000
[328,] 4.739729e-90 9.479459e-90 1.0000000
[329,] 1.257923e-90 2.515846e-90 1.0000000
[330,] 2.867742e-91 5.735484e-91 1.0000000
[331,] 6.400568e-92 1.280114e-91 1.0000000
[332,] 1.422070e-92 2.844140e-92 1.0000000
[333,] 3.234830e-93 6.469660e-93 1.0000000
[334,] 7.302499e-94 1.460500e-93 1.0000000
[335,] 1.584240e-94 3.168479e-94 1.0000000
[336,] 3.411563e-95 6.823127e-95 1.0000000
[337,] 7.942871e-96 1.588574e-95 1.0000000
[338,] 1.698175e-96 3.396351e-96 1.0000000
[339,] 3.586179e-97 7.172358e-97 1.0000000
[340,] 7.468370e-98 1.493674e-97 1.0000000
[341,] 1.549369e-98 3.098739e-98 1.0000000
[342,] 5.799935e-99 1.159987e-98 1.0000000
[343,] 1.233566e-99 2.467132e-99 1.0000000
[344,] 2.511657e-100 5.023314e-100 1.0000000
[345,] 5.508980e-101 1.101796e-100 1.0000000
[346,] 1.125702e-101 2.251404e-101 1.0000000
[347,] 2.260961e-102 4.521922e-102 1.0000000
[348,] 4.456139e-103 8.912278e-103 1.0000000
[349,] 9.264020e-104 1.852804e-103 1.0000000
[350,] 3.225471e-104 6.450942e-104 1.0000000
[351,] 6.220965e-105 1.244193e-104 1.0000000
[352,] 1.189691e-105 2.379383e-105 1.0000000
[353,] 2.255730e-106 4.511461e-106 1.0000000
[354,] 4.253237e-107 8.506474e-107 1.0000000
[355,] 7.925552e-108 1.585110e-107 1.0000000
[356,] 3.224558e-108 6.449116e-108 1.0000000
[357,] 6.221687e-109 1.244337e-108 1.0000000
[358,] 1.136705e-109 2.273409e-109 1.0000000
[359,] 2.054750e-110 4.109500e-110 1.0000000
[360,] 3.680542e-111 7.361085e-111 1.0000000
[361,] 6.532314e-112 1.306463e-111 1.0000000
[362,] 1.166145e-112 2.332291e-112 1.0000000
[363,] 2.047116e-113 4.094233e-113 1.0000000
[364,] 3.749704e-114 7.499408e-114 1.0000000
[365,] 6.666230e-115 1.333246e-114 1.0000000
[366,] 1.813893e-115 3.627787e-115 1.0000000
[367,] 3.268060e-116 6.536119e-116 1.0000000
[368,] 5.469190e-117 1.093838e-116 1.0000000
[369,] 9.139755e-118 1.827951e-117 1.0000000
[370,] 1.678753e-118 3.357506e-118 1.0000000
[371,] 1.189169e-118 2.378339e-118 1.0000000
[372,] 1.974428e-119 3.948856e-119 1.0000000
[373,] 3.189461e-120 6.378921e-120 1.0000000
[374,] 5.640105e-121 1.128021e-120 1.0000000
[375,] 1.364488e-121 2.728977e-121 1.0000000
[376,] 2.611075e-122 5.222151e-122 1.0000000
[377,] 4.105721e-123 8.211441e-123 1.0000000
[378,] 3.859383e-123 7.718767e-123 1.0000000
[379,] 7.379655e-124 1.475931e-123 1.0000000
[380,] 1.655556e-124 3.311112e-124 1.0000000
[381,] 2.507606e-125 5.015212e-125 1.0000000
[382,] 4.419909e-126 8.839819e-126 1.0000000
[383,] 6.785805e-127 1.357161e-126 1.0000000
[384,] 1.138074e-127 2.276149e-127 1.0000000
[385,] 1.762340e-128 3.524680e-128 1.0000000
[386,] 2.444690e-129 4.889380e-129 1.0000000
[387,] 1.626040e-128 3.252080e-128 1.0000000
[388,] 5.083736e-129 1.016747e-128 1.0000000
[389,] 3.453745e-127 6.907490e-127 1.0000000
[390,] 6.878716e-128 1.375743e-127 1.0000000
[391,] 4.061757e-123 8.123514e-123 1.0000000
[392,] 1.649220e-123 3.298439e-123 1.0000000
[393,] 1.298665e-123 2.597331e-123 1.0000000
[394,] 2.095980e-124 4.191960e-124 1.0000000
[395,] 3.717499e-125 7.434998e-125 1.0000000
[396,] 8.612823e-126 1.722565e-125 1.0000000
[397,] 1.454992e-126 2.909985e-126 1.0000000
[398,] 2.328971e-127 4.657943e-127 1.0000000
[399,] 4.811506e-128 9.623012e-128 1.0000000
[400,] 1.057440e-128 2.114881e-128 1.0000000
[401,] 1.462247e-129 2.924493e-129 1.0000000
[402,] 2.153549e-130 4.307097e-130 1.0000000
[403,] 2.733865e-131 5.467730e-131 1.0000000
[404,] 1.267842e-124 2.535683e-124 1.0000000
[405,] 4.409326e-125 8.818653e-125 1.0000000
[406,] 1.514407e-124 3.028815e-124 1.0000000
[407,] 7.863393e-121 1.572679e-120 1.0000000
[408,] 1.156639e-121 2.313278e-121 1.0000000
[409,] 8.431701e-122 1.686340e-121 1.0000000
[410,] 5.847202e-118 1.169440e-117 1.0000000
[411,] 2.992952e-116 5.985904e-116 1.0000000
[412,] 9.868378e-115 1.973676e-114 1.0000000
[413,] 6.394888e-114 1.278978e-113 1.0000000
[414,] 6.690517e-110 1.338103e-109 1.0000000
[415,] 1.131871e-84 2.263741e-84 1.0000000
[416,] 5.888222e-51 1.177644e-50 1.0000000
> postscript(file="/var/www/html/rcomp/tmp/1menx1291297725.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> plot(x[,1], type='l', main='Actuals and Interpolation', ylab='value of Actuals and Interpolation (dots)', xlab='time or index')
> points(x[,1]-mysum$resid)
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/2menx1291297725.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> plot(mysum$resid, type='b', pch=19, main='Residuals', ylab='value of Residuals', xlab='time or index')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/3f5401291297725.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> hist(mysum$resid, main='Residual Histogram', xlab='values of Residuals')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/4f5401291297725.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> densityplot(~mysum$resid,col='black',main='Residual Density Plot', xlab='values of Residuals')
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/5f5401291297725.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
18112.030820 3189.149011 -5975.240229 -5896.762010 -4963.274386
6 7 8 9 10
-3586.218801 -3639.496594 -1771.217901 -5120.670902 -7276.418018
11 12 13 14 15
3969.352621 -4492.098869 -10688.605698 -3692.210291 -3658.988881
16 17 18 19 20
-571.565104 8580.303693 -2205.676265 -7792.585933 749.071535
21 22 23 24 25
-3340.109525 -3344.907924 1462.809853 760.744286 2502.076015
26 27 28 29 30
-3164.005247 -3829.109271 -2410.225134 -5503.404364 -3004.098847
31 32 33 34 35
-1916.099554 -860.802328 -2790.326038 -779.212044 8916.934765
36 37 38 39 40
-407.005211 -2080.687204 -935.992092 3988.236079 346.561017
41 42 43 44 45
1884.552701 -220.002535 -705.251195 1014.684702 -1268.559636
46 47 48 49 50
446.715966 5762.193305 1607.948969 818.275166 329.161773
51 52 53 54 55
-121.639573 1192.597994 3825.734858 1898.113007 3883.404235
56 57 58 59 60
-1663.508514 -3703.968891 -777.589508 656.930989 -707.185995
61 62 63 64 65
-745.130968 -164.659171 -120.864169 -715.556170 -1514.636577
66 67 68 69 70
-650.065442 217.786117 318.539775 1989.611071 19.621278
71 72 73 74 75
-2214.551765 -2783.762310 1173.976909 -213.043438 1246.830179
76 77 78 79 80
82.895705 1252.997086 22.872089 -676.919991 1592.900422
81 82 83 84 85
-1299.502738 -530.426761 2538.618480 -2052.259268 635.543149
86 87 88 89 90
-1470.195210 -656.375725 829.144472 1120.037178 609.702605
91 92 93 94 95
91.538895 -504.930538 -1849.356020 1524.364018 -2564.966168
96 97 98 99 100
-1618.168056 849.475477 -761.256688 1699.459918 -533.611160
101 102 103 104 105
-328.546992 821.087835 1179.268373 1533.402505 -1072.452550
106 107 108 109 110
-304.685718 -199.410741 -481.158937 561.232612 -191.644522
111 112 113 114 115
-1300.599370 -985.347261 10.897599 -399.340499 -1661.614585
116 117 118 119 120
-1532.042011 11.808421 -562.540244 202.842125 -1088.325075
121 122 123 124 125
375.001446 -914.745069 179.341066 185.809164 -1335.232702
126 127 128 129 130
47.161008 -1061.735897 -1061.158843 -1245.133403 -1061.735897
131 132 133 134 135
-1063.834583 -1725.494028 -1061.735897 -1045.872342 -1082.659339
136 137 138 139 140
-1090.258677 -1080.264849 -1059.405046 -1114.409716 -1061.735897
141 142 143 144 145
-1231.006356 -1061.735897 -1223.175503 -1058.869179 -1047.115252
146 147 148 149 150
-1274.122439 -1061.735897 -1055.892234 -1135.815656 -1050.887300
151 152 153 154 155
-1059.353040 -1138.519584 -1062.843185 -1092.243375 -1051.874989
156 157 158 159 160
-1061.735897 -1061.735897 -1029.374768 -1203.209376 -1061.735897
161 162 163 164 165
-1061.735897 -1061.735897 -1061.735897 -1034.423174 -893.054439
166 167 168 169 170
-1061.927612 -1061.735897 -2242.074409 -1508.435528 -1061.735897
171 172 173 174 175
-1061.735897 -1005.825011 -1061.735897 -982.565463 -1061.735897
176 177 178 179 180
-1060.752605 -1098.819809 -1069.429333 -1061.735897 -1061.735897
181 182 183 184 185
-841.713596 -643.649635 -1453.169244 -434.857855 -970.539261
186 187 188 189 190
-1223.462113 -1214.789853 -995.191188 -692.406099 -1061.735897
191 192 193 194 195
-1159.573192 -1061.735897 -1026.952296 -1061.735897 -1057.512909
196 197 198 199 200
-1061.735897 670.444619 -1061.735897 -1061.735897 122.450100
201 202 203 204 205
-1061.735897 -1061.735897 -1061.735897 -1481.573495 160.496827
206 207 208 209 210
-1021.841361 -839.829927 -1061.735897 -1900.832563 -1029.650703
211 212 213 214 215
-139.270215 127.729785 1739.108604 -1184.234238 -1096.934730
216 217 218 219 220
-124.391717 247.952847 -7.271879 1726.608283 1484.825340
221 222 223 224 225
-1015.125087 -136.270215 332.260570 -279.024713 129.498368
226 227 228 229 230
-253.640739 -180.646183 52.729785 -1298.125172 177.729785
231 232 233 234 235
-930.761739 155.492980 228.729785 -952.389572 -137.270215
236 237 238 239 240
-139.093865 -86.270215 -37.133233 -139.270215 -139.270215
241 242 243 244 245
-45.270215 -403.968451 -115.270215 2404.716804 -2344.517538
246 247 248 249 250
-139.270215 -8.270215 -695.891470 -139.270215 66.729785
251 252 253 254 255
625.561964 27.729785 -694.834735 -1076.587638 -8060.784634
256 257 258 259 260
225.729785 1244.262432 1908.386400 -139.270215 -139.270215
261 262 263 264 265
-139.270215 86.729785 -1109.947060 1089.257151 -139.270215
266 267 268 269 270
-139.270215 48.729785 -649.474595 -1.270215 -1365.460633
271 272 273 274 275
-139.270215 -139.270215 -14.270215 -412.769924 142.729785
276 277 278 279 280
39.764643 -658.661062 1184.729785 -2612.549884 -2014.058505
281 282 283 284 285
-139.270215 109.729785 650.616351 193.729785 -1589.334055
286 287 288 289 290
461.729785 79.057971 -299.019094 109.729785 162.247354
291 292 293 294 295
25.729785 325.753632 -2168.191760 -86.270215 -127.865024
296 297 298 299 300
1420.496699 -139.270215 -139.270215 -139.270215 -109.816972
301 302 303 304 305
24.072522 -412.072342 -139.330966 226.729785 1051.433365
306 307 308 309 310
-142.631572 69.729785 -423.947404 923.587931 -139.270215
311 312 313 314 315
225.729785 -4465.476016 -90.270215 -559.628047 -2.242662
316 317 318 319 320
-89.365450 370.348242 -1136.768660 -139.746521 -139.270215
321 322 323 324 325
-139.270215 -139.270215 -139.270215 -139.270215 -117.270215
326 327 328 329 330
-281.518718 -139.270215 -139.270215 -139.270215 -139.270215
331 332 333 334 335
-139.270215 -139.270215 -43.270215 -144.309852 173.512840
336 337 338 339 340
704.043521 -1133.525590 -113.270215 -127.184237 496.643627
341 342 343 344 345
-1613.230251 -139.270215 -139.270215 481.729785 -490.925235
346 347 348 349 350
-20.270215 -361.063762 -139.270215 1455.729785 368.585152
351 352 353 354 355
553.501829 292.282379 752.484755 -494.815251 -139.270215
356 357 358 359 360
374.729785 -3795.045538 -139.270215 -139.270215 -138.270215
361 362 363 364 365
-139.069613 -139.270215 1623.729785 -502.779063 150.011155
366 367 368 369 370
-139.270215 -139.270215 -139.270215 78.729785 -131.175020
371 372 373 374 375
308.729785 688.326689 -2445.155383 -1439.270720 -139.270215
376 377 378 379 380
-18.270215 487.406100 2477.986548 80.964476 -139.270215
381 382 383 384 385
390.729785 -2290.634823 2398.483194 -61.270215 2006.554552
386 387 388 389 390
-5383.735718 1593.866401 85.075320 446.233013 364.247630
391 392 393 394 395
211.575315 -353.489569 -2062.862285 4234.482677 -6199.111883
396 397 398 399 400
3712.798981 -480.110115 -6949.088627 347.097654 5540.850501
401 402 403 404 405
-2058.527317 586.628153 815.992878 -3252.172767 -1639.129780
406 407 408 409 410
-1155.084059 -682.175458 161.294988 266.453118 -5000.528162
411 412 413 414 415
4958.071007 -937.482146 1675.626363 5109.128002 -10702.652116
416 417 418 419 420
-5180.946020 11155.125275 -443.565059 -3332.403928 -3668.241362
421 422 423 424 425
-2348.678631 5066.650137 3417.018967 -14401.817777 4095.498813
426 427 428 429 430
-2470.243562 24153.403695 -12733.619516 7888.505936 17046.328143
431
153933.352077
> postscript(file="/var/www/html/rcomp/tmp/6pell1291297725.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 18112.030820 NA
1 3189.149011 18112.030820
2 -5975.240229 3189.149011
3 -5896.762010 -5975.240229
4 -4963.274386 -5896.762010
5 -3586.218801 -4963.274386
6 -3639.496594 -3586.218801
7 -1771.217901 -3639.496594
8 -5120.670902 -1771.217901
9 -7276.418018 -5120.670902
10 3969.352621 -7276.418018
11 -4492.098869 3969.352621
12 -10688.605698 -4492.098869
13 -3692.210291 -10688.605698
14 -3658.988881 -3692.210291
15 -571.565104 -3658.988881
16 8580.303693 -571.565104
17 -2205.676265 8580.303693
18 -7792.585933 -2205.676265
19 749.071535 -7792.585933
20 -3340.109525 749.071535
21 -3344.907924 -3340.109525
22 1462.809853 -3344.907924
23 760.744286 1462.809853
24 2502.076015 760.744286
25 -3164.005247 2502.076015
26 -3829.109271 -3164.005247
27 -2410.225134 -3829.109271
28 -5503.404364 -2410.225134
29 -3004.098847 -5503.404364
30 -1916.099554 -3004.098847
31 -860.802328 -1916.099554
32 -2790.326038 -860.802328
33 -779.212044 -2790.326038
34 8916.934765 -779.212044
35 -407.005211 8916.934765
36 -2080.687204 -407.005211
37 -935.992092 -2080.687204
38 3988.236079 -935.992092
39 346.561017 3988.236079
40 1884.552701 346.561017
41 -220.002535 1884.552701
42 -705.251195 -220.002535
43 1014.684702 -705.251195
44 -1268.559636 1014.684702
45 446.715966 -1268.559636
46 5762.193305 446.715966
47 1607.948969 5762.193305
48 818.275166 1607.948969
49 329.161773 818.275166
50 -121.639573 329.161773
51 1192.597994 -121.639573
52 3825.734858 1192.597994
53 1898.113007 3825.734858
54 3883.404235 1898.113007
55 -1663.508514 3883.404235
56 -3703.968891 -1663.508514
57 -777.589508 -3703.968891
58 656.930989 -777.589508
59 -707.185995 656.930989
60 -745.130968 -707.185995
61 -164.659171 -745.130968
62 -120.864169 -164.659171
63 -715.556170 -120.864169
64 -1514.636577 -715.556170
65 -650.065442 -1514.636577
66 217.786117 -650.065442
67 318.539775 217.786117
68 1989.611071 318.539775
69 19.621278 1989.611071
70 -2214.551765 19.621278
71 -2783.762310 -2214.551765
72 1173.976909 -2783.762310
73 -213.043438 1173.976909
74 1246.830179 -213.043438
75 82.895705 1246.830179
76 1252.997086 82.895705
77 22.872089 1252.997086
78 -676.919991 22.872089
79 1592.900422 -676.919991
80 -1299.502738 1592.900422
81 -530.426761 -1299.502738
82 2538.618480 -530.426761
83 -2052.259268 2538.618480
84 635.543149 -2052.259268
85 -1470.195210 635.543149
86 -656.375725 -1470.195210
87 829.144472 -656.375725
88 1120.037178 829.144472
89 609.702605 1120.037178
90 91.538895 609.702605
91 -504.930538 91.538895
92 -1849.356020 -504.930538
93 1524.364018 -1849.356020
94 -2564.966168 1524.364018
95 -1618.168056 -2564.966168
96 849.475477 -1618.168056
97 -761.256688 849.475477
98 1699.459918 -761.256688
99 -533.611160 1699.459918
100 -328.546992 -533.611160
101 821.087835 -328.546992
102 1179.268373 821.087835
103 1533.402505 1179.268373
104 -1072.452550 1533.402505
105 -304.685718 -1072.452550
106 -199.410741 -304.685718
107 -481.158937 -199.410741
108 561.232612 -481.158937
109 -191.644522 561.232612
110 -1300.599370 -191.644522
111 -985.347261 -1300.599370
112 10.897599 -985.347261
113 -399.340499 10.897599
114 -1661.614585 -399.340499
115 -1532.042011 -1661.614585
116 11.808421 -1532.042011
117 -562.540244 11.808421
118 202.842125 -562.540244
119 -1088.325075 202.842125
120 375.001446 -1088.325075
121 -914.745069 375.001446
122 179.341066 -914.745069
123 185.809164 179.341066
124 -1335.232702 185.809164
125 47.161008 -1335.232702
126 -1061.735897 47.161008
127 -1061.158843 -1061.735897
128 -1245.133403 -1061.158843
129 -1061.735897 -1245.133403
130 -1063.834583 -1061.735897
131 -1725.494028 -1063.834583
132 -1061.735897 -1725.494028
133 -1045.872342 -1061.735897
134 -1082.659339 -1045.872342
135 -1090.258677 -1082.659339
136 -1080.264849 -1090.258677
137 -1059.405046 -1080.264849
138 -1114.409716 -1059.405046
139 -1061.735897 -1114.409716
140 -1231.006356 -1061.735897
141 -1061.735897 -1231.006356
142 -1223.175503 -1061.735897
143 -1058.869179 -1223.175503
144 -1047.115252 -1058.869179
145 -1274.122439 -1047.115252
146 -1061.735897 -1274.122439
147 -1055.892234 -1061.735897
148 -1135.815656 -1055.892234
149 -1050.887300 -1135.815656
150 -1059.353040 -1050.887300
151 -1138.519584 -1059.353040
152 -1062.843185 -1138.519584
153 -1092.243375 -1062.843185
154 -1051.874989 -1092.243375
155 -1061.735897 -1051.874989
156 -1061.735897 -1061.735897
157 -1029.374768 -1061.735897
158 -1203.209376 -1029.374768
159 -1061.735897 -1203.209376
160 -1061.735897 -1061.735897
161 -1061.735897 -1061.735897
162 -1061.735897 -1061.735897
163 -1034.423174 -1061.735897
164 -893.054439 -1034.423174
165 -1061.927612 -893.054439
166 -1061.735897 -1061.927612
167 -2242.074409 -1061.735897
168 -1508.435528 -2242.074409
169 -1061.735897 -1508.435528
170 -1061.735897 -1061.735897
171 -1005.825011 -1061.735897
172 -1061.735897 -1005.825011
173 -982.565463 -1061.735897
174 -1061.735897 -982.565463
175 -1060.752605 -1061.735897
176 -1098.819809 -1060.752605
177 -1069.429333 -1098.819809
178 -1061.735897 -1069.429333
179 -1061.735897 -1061.735897
180 -841.713596 -1061.735897
181 -643.649635 -841.713596
182 -1453.169244 -643.649635
183 -434.857855 -1453.169244
184 -970.539261 -434.857855
185 -1223.462113 -970.539261
186 -1214.789853 -1223.462113
187 -995.191188 -1214.789853
188 -692.406099 -995.191188
189 -1061.735897 -692.406099
190 -1159.573192 -1061.735897
191 -1061.735897 -1159.573192
192 -1026.952296 -1061.735897
193 -1061.735897 -1026.952296
194 -1057.512909 -1061.735897
195 -1061.735897 -1057.512909
196 670.444619 -1061.735897
197 -1061.735897 670.444619
198 -1061.735897 -1061.735897
199 122.450100 -1061.735897
200 -1061.735897 122.450100
201 -1061.735897 -1061.735897
202 -1061.735897 -1061.735897
203 -1481.573495 -1061.735897
204 160.496827 -1481.573495
205 -1021.841361 160.496827
206 -839.829927 -1021.841361
207 -1061.735897 -839.829927
208 -1900.832563 -1061.735897
209 -1029.650703 -1900.832563
210 -139.270215 -1029.650703
211 127.729785 -139.270215
212 1739.108604 127.729785
213 -1184.234238 1739.108604
214 -1096.934730 -1184.234238
215 -124.391717 -1096.934730
216 247.952847 -124.391717
217 -7.271879 247.952847
218 1726.608283 -7.271879
219 1484.825340 1726.608283
220 -1015.125087 1484.825340
221 -136.270215 -1015.125087
222 332.260570 -136.270215
223 -279.024713 332.260570
224 129.498368 -279.024713
225 -253.640739 129.498368
226 -180.646183 -253.640739
227 52.729785 -180.646183
228 -1298.125172 52.729785
229 177.729785 -1298.125172
230 -930.761739 177.729785
231 155.492980 -930.761739
232 228.729785 155.492980
233 -952.389572 228.729785
234 -137.270215 -952.389572
235 -139.093865 -137.270215
236 -86.270215 -139.093865
237 -37.133233 -86.270215
238 -139.270215 -37.133233
239 -139.270215 -139.270215
240 -45.270215 -139.270215
241 -403.968451 -45.270215
242 -115.270215 -403.968451
243 2404.716804 -115.270215
244 -2344.517538 2404.716804
245 -139.270215 -2344.517538
246 -8.270215 -139.270215
247 -695.891470 -8.270215
248 -139.270215 -695.891470
249 66.729785 -139.270215
250 625.561964 66.729785
251 27.729785 625.561964
252 -694.834735 27.729785
253 -1076.587638 -694.834735
254 -8060.784634 -1076.587638
255 225.729785 -8060.784634
256 1244.262432 225.729785
257 1908.386400 1244.262432
258 -139.270215 1908.386400
259 -139.270215 -139.270215
260 -139.270215 -139.270215
261 86.729785 -139.270215
262 -1109.947060 86.729785
263 1089.257151 -1109.947060
264 -139.270215 1089.257151
265 -139.270215 -139.270215
266 48.729785 -139.270215
267 -649.474595 48.729785
268 -1.270215 -649.474595
269 -1365.460633 -1.270215
270 -139.270215 -1365.460633
271 -139.270215 -139.270215
272 -14.270215 -139.270215
273 -412.769924 -14.270215
274 142.729785 -412.769924
275 39.764643 142.729785
276 -658.661062 39.764643
277 1184.729785 -658.661062
278 -2612.549884 1184.729785
279 -2014.058505 -2612.549884
280 -139.270215 -2014.058505
281 109.729785 -139.270215
282 650.616351 109.729785
283 193.729785 650.616351
284 -1589.334055 193.729785
285 461.729785 -1589.334055
286 79.057971 461.729785
287 -299.019094 79.057971
288 109.729785 -299.019094
289 162.247354 109.729785
290 25.729785 162.247354
291 325.753632 25.729785
292 -2168.191760 325.753632
293 -86.270215 -2168.191760
294 -127.865024 -86.270215
295 1420.496699 -127.865024
296 -139.270215 1420.496699
297 -139.270215 -139.270215
298 -139.270215 -139.270215
299 -109.816972 -139.270215
300 24.072522 -109.816972
301 -412.072342 24.072522
302 -139.330966 -412.072342
303 226.729785 -139.330966
304 1051.433365 226.729785
305 -142.631572 1051.433365
306 69.729785 -142.631572
307 -423.947404 69.729785
308 923.587931 -423.947404
309 -139.270215 923.587931
310 225.729785 -139.270215
311 -4465.476016 225.729785
312 -90.270215 -4465.476016
313 -559.628047 -90.270215
314 -2.242662 -559.628047
315 -89.365450 -2.242662
316 370.348242 -89.365450
317 -1136.768660 370.348242
318 -139.746521 -1136.768660
319 -139.270215 -139.746521
320 -139.270215 -139.270215
321 -139.270215 -139.270215
322 -139.270215 -139.270215
323 -139.270215 -139.270215
324 -117.270215 -139.270215
325 -281.518718 -117.270215
326 -139.270215 -281.518718
327 -139.270215 -139.270215
328 -139.270215 -139.270215
329 -139.270215 -139.270215
330 -139.270215 -139.270215
331 -139.270215 -139.270215
332 -43.270215 -139.270215
333 -144.309852 -43.270215
334 173.512840 -144.309852
335 704.043521 173.512840
336 -1133.525590 704.043521
337 -113.270215 -1133.525590
338 -127.184237 -113.270215
339 496.643627 -127.184237
340 -1613.230251 496.643627
341 -139.270215 -1613.230251
342 -139.270215 -139.270215
343 481.729785 -139.270215
344 -490.925235 481.729785
345 -20.270215 -490.925235
346 -361.063762 -20.270215
347 -139.270215 -361.063762
348 1455.729785 -139.270215
349 368.585152 1455.729785
350 553.501829 368.585152
351 292.282379 553.501829
352 752.484755 292.282379
353 -494.815251 752.484755
354 -139.270215 -494.815251
355 374.729785 -139.270215
356 -3795.045538 374.729785
357 -139.270215 -3795.045538
358 -139.270215 -139.270215
359 -138.270215 -139.270215
360 -139.069613 -138.270215
361 -139.270215 -139.069613
362 1623.729785 -139.270215
363 -502.779063 1623.729785
364 150.011155 -502.779063
365 -139.270215 150.011155
366 -139.270215 -139.270215
367 -139.270215 -139.270215
368 78.729785 -139.270215
369 -131.175020 78.729785
370 308.729785 -131.175020
371 688.326689 308.729785
372 -2445.155383 688.326689
373 -1439.270720 -2445.155383
374 -139.270215 -1439.270720
375 -18.270215 -139.270215
376 487.406100 -18.270215
377 2477.986548 487.406100
378 80.964476 2477.986548
379 -139.270215 80.964476
380 390.729785 -139.270215
381 -2290.634823 390.729785
382 2398.483194 -2290.634823
383 -61.270215 2398.483194
384 2006.554552 -61.270215
385 -5383.735718 2006.554552
386 1593.866401 -5383.735718
387 85.075320 1593.866401
388 446.233013 85.075320
389 364.247630 446.233013
390 211.575315 364.247630
391 -353.489569 211.575315
392 -2062.862285 -353.489569
393 4234.482677 -2062.862285
394 -6199.111883 4234.482677
395 3712.798981 -6199.111883
396 -480.110115 3712.798981
397 -6949.088627 -480.110115
398 347.097654 -6949.088627
399 5540.850501 347.097654
400 -2058.527317 5540.850501
401 586.628153 -2058.527317
402 815.992878 586.628153
403 -3252.172767 815.992878
404 -1639.129780 -3252.172767
405 -1155.084059 -1639.129780
406 -682.175458 -1155.084059
407 161.294988 -682.175458
408 266.453118 161.294988
409 -5000.528162 266.453118
410 4958.071007 -5000.528162
411 -937.482146 4958.071007
412 1675.626363 -937.482146
413 5109.128002 1675.626363
414 -10702.652116 5109.128002
415 -5180.946020 -10702.652116
416 11155.125275 -5180.946020
417 -443.565059 11155.125275
418 -3332.403928 -443.565059
419 -3668.241362 -3332.403928
420 -2348.678631 -3668.241362
421 5066.650137 -2348.678631
422 3417.018967 5066.650137
423 -14401.817777 3417.018967
424 4095.498813 -14401.817777
425 -2470.243562 4095.498813
426 24153.403695 -2470.243562
427 -12733.619516 24153.403695
428 7888.505936 -12733.619516
429 17046.328143 7888.505936
430 153933.352077 17046.328143
431 NA 153933.352077
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] 3189.149011 18112.030820
[2,] -5975.240229 3189.149011
[3,] -5896.762010 -5975.240229
[4,] -4963.274386 -5896.762010
[5,] -3586.218801 -4963.274386
[6,] -3639.496594 -3586.218801
[7,] -1771.217901 -3639.496594
[8,] -5120.670902 -1771.217901
[9,] -7276.418018 -5120.670902
[10,] 3969.352621 -7276.418018
[11,] -4492.098869 3969.352621
[12,] -10688.605698 -4492.098869
[13,] -3692.210291 -10688.605698
[14,] -3658.988881 -3692.210291
[15,] -571.565104 -3658.988881
[16,] 8580.303693 -571.565104
[17,] -2205.676265 8580.303693
[18,] -7792.585933 -2205.676265
[19,] 749.071535 -7792.585933
[20,] -3340.109525 749.071535
[21,] -3344.907924 -3340.109525
[22,] 1462.809853 -3344.907924
[23,] 760.744286 1462.809853
[24,] 2502.076015 760.744286
[25,] -3164.005247 2502.076015
[26,] -3829.109271 -3164.005247
[27,] -2410.225134 -3829.109271
[28,] -5503.404364 -2410.225134
[29,] -3004.098847 -5503.404364
[30,] -1916.099554 -3004.098847
[31,] -860.802328 -1916.099554
[32,] -2790.326038 -860.802328
[33,] -779.212044 -2790.326038
[34,] 8916.934765 -779.212044
[35,] -407.005211 8916.934765
[36,] -2080.687204 -407.005211
[37,] -935.992092 -2080.687204
[38,] 3988.236079 -935.992092
[39,] 346.561017 3988.236079
[40,] 1884.552701 346.561017
[41,] -220.002535 1884.552701
[42,] -705.251195 -220.002535
[43,] 1014.684702 -705.251195
[44,] -1268.559636 1014.684702
[45,] 446.715966 -1268.559636
[46,] 5762.193305 446.715966
[47,] 1607.948969 5762.193305
[48,] 818.275166 1607.948969
[49,] 329.161773 818.275166
[50,] -121.639573 329.161773
[51,] 1192.597994 -121.639573
[52,] 3825.734858 1192.597994
[53,] 1898.113007 3825.734858
[54,] 3883.404235 1898.113007
[55,] -1663.508514 3883.404235
[56,] -3703.968891 -1663.508514
[57,] -777.589508 -3703.968891
[58,] 656.930989 -777.589508
[59,] -707.185995 656.930989
[60,] -745.130968 -707.185995
[61,] -164.659171 -745.130968
[62,] -120.864169 -164.659171
[63,] -715.556170 -120.864169
[64,] -1514.636577 -715.556170
[65,] -650.065442 -1514.636577
[66,] 217.786117 -650.065442
[67,] 318.539775 217.786117
[68,] 1989.611071 318.539775
[69,] 19.621278 1989.611071
[70,] -2214.551765 19.621278
[71,] -2783.762310 -2214.551765
[72,] 1173.976909 -2783.762310
[73,] -213.043438 1173.976909
[74,] 1246.830179 -213.043438
[75,] 82.895705 1246.830179
[76,] 1252.997086 82.895705
[77,] 22.872089 1252.997086
[78,] -676.919991 22.872089
[79,] 1592.900422 -676.919991
[80,] -1299.502738 1592.900422
[81,] -530.426761 -1299.502738
[82,] 2538.618480 -530.426761
[83,] -2052.259268 2538.618480
[84,] 635.543149 -2052.259268
[85,] -1470.195210 635.543149
[86,] -656.375725 -1470.195210
[87,] 829.144472 -656.375725
[88,] 1120.037178 829.144472
[89,] 609.702605 1120.037178
[90,] 91.538895 609.702605
[91,] -504.930538 91.538895
[92,] -1849.356020 -504.930538
[93,] 1524.364018 -1849.356020
[94,] -2564.966168 1524.364018
[95,] -1618.168056 -2564.966168
[96,] 849.475477 -1618.168056
[97,] -761.256688 849.475477
[98,] 1699.459918 -761.256688
[99,] -533.611160 1699.459918
[100,] -328.546992 -533.611160
[101,] 821.087835 -328.546992
[102,] 1179.268373 821.087835
[103,] 1533.402505 1179.268373
[104,] -1072.452550 1533.402505
[105,] -304.685718 -1072.452550
[106,] -199.410741 -304.685718
[107,] -481.158937 -199.410741
[108,] 561.232612 -481.158937
[109,] -191.644522 561.232612
[110,] -1300.599370 -191.644522
[111,] -985.347261 -1300.599370
[112,] 10.897599 -985.347261
[113,] -399.340499 10.897599
[114,] -1661.614585 -399.340499
[115,] -1532.042011 -1661.614585
[116,] 11.808421 -1532.042011
[117,] -562.540244 11.808421
[118,] 202.842125 -562.540244
[119,] -1088.325075 202.842125
[120,] 375.001446 -1088.325075
[121,] -914.745069 375.001446
[122,] 179.341066 -914.745069
[123,] 185.809164 179.341066
[124,] -1335.232702 185.809164
[125,] 47.161008 -1335.232702
[126,] -1061.735897 47.161008
[127,] -1061.158843 -1061.735897
[128,] -1245.133403 -1061.158843
[129,] -1061.735897 -1245.133403
[130,] -1063.834583 -1061.735897
[131,] -1725.494028 -1063.834583
[132,] -1061.735897 -1725.494028
[133,] -1045.872342 -1061.735897
[134,] -1082.659339 -1045.872342
[135,] -1090.258677 -1082.659339
[136,] -1080.264849 -1090.258677
[137,] -1059.405046 -1080.264849
[138,] -1114.409716 -1059.405046
[139,] -1061.735897 -1114.409716
[140,] -1231.006356 -1061.735897
[141,] -1061.735897 -1231.006356
[142,] -1223.175503 -1061.735897
[143,] -1058.869179 -1223.175503
[144,] -1047.115252 -1058.869179
[145,] -1274.122439 -1047.115252
[146,] -1061.735897 -1274.122439
[147,] -1055.892234 -1061.735897
[148,] -1135.815656 -1055.892234
[149,] -1050.887300 -1135.815656
[150,] -1059.353040 -1050.887300
[151,] -1138.519584 -1059.353040
[152,] -1062.843185 -1138.519584
[153,] -1092.243375 -1062.843185
[154,] -1051.874989 -1092.243375
[155,] -1061.735897 -1051.874989
[156,] -1061.735897 -1061.735897
[157,] -1029.374768 -1061.735897
[158,] -1203.209376 -1029.374768
[159,] -1061.735897 -1203.209376
[160,] -1061.735897 -1061.735897
[161,] -1061.735897 -1061.735897
[162,] -1061.735897 -1061.735897
[163,] -1034.423174 -1061.735897
[164,] -893.054439 -1034.423174
[165,] -1061.927612 -893.054439
[166,] -1061.735897 -1061.927612
[167,] -2242.074409 -1061.735897
[168,] -1508.435528 -2242.074409
[169,] -1061.735897 -1508.435528
[170,] -1061.735897 -1061.735897
[171,] -1005.825011 -1061.735897
[172,] -1061.735897 -1005.825011
[173,] -982.565463 -1061.735897
[174,] -1061.735897 -982.565463
[175,] -1060.752605 -1061.735897
[176,] -1098.819809 -1060.752605
[177,] -1069.429333 -1098.819809
[178,] -1061.735897 -1069.429333
[179,] -1061.735897 -1061.735897
[180,] -841.713596 -1061.735897
[181,] -643.649635 -841.713596
[182,] -1453.169244 -643.649635
[183,] -434.857855 -1453.169244
[184,] -970.539261 -434.857855
[185,] -1223.462113 -970.539261
[186,] -1214.789853 -1223.462113
[187,] -995.191188 -1214.789853
[188,] -692.406099 -995.191188
[189,] -1061.735897 -692.406099
[190,] -1159.573192 -1061.735897
[191,] -1061.735897 -1159.573192
[192,] -1026.952296 -1061.735897
[193,] -1061.735897 -1026.952296
[194,] -1057.512909 -1061.735897
[195,] -1061.735897 -1057.512909
[196,] 670.444619 -1061.735897
[197,] -1061.735897 670.444619
[198,] -1061.735897 -1061.735897
[199,] 122.450100 -1061.735897
[200,] -1061.735897 122.450100
[201,] -1061.735897 -1061.735897
[202,] -1061.735897 -1061.735897
[203,] -1481.573495 -1061.735897
[204,] 160.496827 -1481.573495
[205,] -1021.841361 160.496827
[206,] -839.829927 -1021.841361
[207,] -1061.735897 -839.829927
[208,] -1900.832563 -1061.735897
[209,] -1029.650703 -1900.832563
[210,] -139.270215 -1029.650703
[211,] 127.729785 -139.270215
[212,] 1739.108604 127.729785
[213,] -1184.234238 1739.108604
[214,] -1096.934730 -1184.234238
[215,] -124.391717 -1096.934730
[216,] 247.952847 -124.391717
[217,] -7.271879 247.952847
[218,] 1726.608283 -7.271879
[219,] 1484.825340 1726.608283
[220,] -1015.125087 1484.825340
[221,] -136.270215 -1015.125087
[222,] 332.260570 -136.270215
[223,] -279.024713 332.260570
[224,] 129.498368 -279.024713
[225,] -253.640739 129.498368
[226,] -180.646183 -253.640739
[227,] 52.729785 -180.646183
[228,] -1298.125172 52.729785
[229,] 177.729785 -1298.125172
[230,] -930.761739 177.729785
[231,] 155.492980 -930.761739
[232,] 228.729785 155.492980
[233,] -952.389572 228.729785
[234,] -137.270215 -952.389572
[235,] -139.093865 -137.270215
[236,] -86.270215 -139.093865
[237,] -37.133233 -86.270215
[238,] -139.270215 -37.133233
[239,] -139.270215 -139.270215
[240,] -45.270215 -139.270215
[241,] -403.968451 -45.270215
[242,] -115.270215 -403.968451
[243,] 2404.716804 -115.270215
[244,] -2344.517538 2404.716804
[245,] -139.270215 -2344.517538
[246,] -8.270215 -139.270215
[247,] -695.891470 -8.270215
[248,] -139.270215 -695.891470
[249,] 66.729785 -139.270215
[250,] 625.561964 66.729785
[251,] 27.729785 625.561964
[252,] -694.834735 27.729785
[253,] -1076.587638 -694.834735
[254,] -8060.784634 -1076.587638
[255,] 225.729785 -8060.784634
[256,] 1244.262432 225.729785
[257,] 1908.386400 1244.262432
[258,] -139.270215 1908.386400
[259,] -139.270215 -139.270215
[260,] -139.270215 -139.270215
[261,] 86.729785 -139.270215
[262,] -1109.947060 86.729785
[263,] 1089.257151 -1109.947060
[264,] -139.270215 1089.257151
[265,] -139.270215 -139.270215
[266,] 48.729785 -139.270215
[267,] -649.474595 48.729785
[268,] -1.270215 -649.474595
[269,] -1365.460633 -1.270215
[270,] -139.270215 -1365.460633
[271,] -139.270215 -139.270215
[272,] -14.270215 -139.270215
[273,] -412.769924 -14.270215
[274,] 142.729785 -412.769924
[275,] 39.764643 142.729785
[276,] -658.661062 39.764643
[277,] 1184.729785 -658.661062
[278,] -2612.549884 1184.729785
[279,] -2014.058505 -2612.549884
[280,] -139.270215 -2014.058505
[281,] 109.729785 -139.270215
[282,] 650.616351 109.729785
[283,] 193.729785 650.616351
[284,] -1589.334055 193.729785
[285,] 461.729785 -1589.334055
[286,] 79.057971 461.729785
[287,] -299.019094 79.057971
[288,] 109.729785 -299.019094
[289,] 162.247354 109.729785
[290,] 25.729785 162.247354
[291,] 325.753632 25.729785
[292,] -2168.191760 325.753632
[293,] -86.270215 -2168.191760
[294,] -127.865024 -86.270215
[295,] 1420.496699 -127.865024
[296,] -139.270215 1420.496699
[297,] -139.270215 -139.270215
[298,] -139.270215 -139.270215
[299,] -109.816972 -139.270215
[300,] 24.072522 -109.816972
[301,] -412.072342 24.072522
[302,] -139.330966 -412.072342
[303,] 226.729785 -139.330966
[304,] 1051.433365 226.729785
[305,] -142.631572 1051.433365
[306,] 69.729785 -142.631572
[307,] -423.947404 69.729785
[308,] 923.587931 -423.947404
[309,] -139.270215 923.587931
[310,] 225.729785 -139.270215
[311,] -4465.476016 225.729785
[312,] -90.270215 -4465.476016
[313,] -559.628047 -90.270215
[314,] -2.242662 -559.628047
[315,] -89.365450 -2.242662
[316,] 370.348242 -89.365450
[317,] -1136.768660 370.348242
[318,] -139.746521 -1136.768660
[319,] -139.270215 -139.746521
[320,] -139.270215 -139.270215
[321,] -139.270215 -139.270215
[322,] -139.270215 -139.270215
[323,] -139.270215 -139.270215
[324,] -117.270215 -139.270215
[325,] -281.518718 -117.270215
[326,] -139.270215 -281.518718
[327,] -139.270215 -139.270215
[328,] -139.270215 -139.270215
[329,] -139.270215 -139.270215
[330,] -139.270215 -139.270215
[331,] -139.270215 -139.270215
[332,] -43.270215 -139.270215
[333,] -144.309852 -43.270215
[334,] 173.512840 -144.309852
[335,] 704.043521 173.512840
[336,] -1133.525590 704.043521
[337,] -113.270215 -1133.525590
[338,] -127.184237 -113.270215
[339,] 496.643627 -127.184237
[340,] -1613.230251 496.643627
[341,] -139.270215 -1613.230251
[342,] -139.270215 -139.270215
[343,] 481.729785 -139.270215
[344,] -490.925235 481.729785
[345,] -20.270215 -490.925235
[346,] -361.063762 -20.270215
[347,] -139.270215 -361.063762
[348,] 1455.729785 -139.270215
[349,] 368.585152 1455.729785
[350,] 553.501829 368.585152
[351,] 292.282379 553.501829
[352,] 752.484755 292.282379
[353,] -494.815251 752.484755
[354,] -139.270215 -494.815251
[355,] 374.729785 -139.270215
[356,] -3795.045538 374.729785
[357,] -139.270215 -3795.045538
[358,] -139.270215 -139.270215
[359,] -138.270215 -139.270215
[360,] -139.069613 -138.270215
[361,] -139.270215 -139.069613
[362,] 1623.729785 -139.270215
[363,] -502.779063 1623.729785
[364,] 150.011155 -502.779063
[365,] -139.270215 150.011155
[366,] -139.270215 -139.270215
[367,] -139.270215 -139.270215
[368,] 78.729785 -139.270215
[369,] -131.175020 78.729785
[370,] 308.729785 -131.175020
[371,] 688.326689 308.729785
[372,] -2445.155383 688.326689
[373,] -1439.270720 -2445.155383
[374,] -139.270215 -1439.270720
[375,] -18.270215 -139.270215
[376,] 487.406100 -18.270215
[377,] 2477.986548 487.406100
[378,] 80.964476 2477.986548
[379,] -139.270215 80.964476
[380,] 390.729785 -139.270215
[381,] -2290.634823 390.729785
[382,] 2398.483194 -2290.634823
[383,] -61.270215 2398.483194
[384,] 2006.554552 -61.270215
[385,] -5383.735718 2006.554552
[386,] 1593.866401 -5383.735718
[387,] 85.075320 1593.866401
[388,] 446.233013 85.075320
[389,] 364.247630 446.233013
[390,] 211.575315 364.247630
[391,] -353.489569 211.575315
[392,] -2062.862285 -353.489569
[393,] 4234.482677 -2062.862285
[394,] -6199.111883 4234.482677
[395,] 3712.798981 -6199.111883
[396,] -480.110115 3712.798981
[397,] -6949.088627 -480.110115
[398,] 347.097654 -6949.088627
[399,] 5540.850501 347.097654
[400,] -2058.527317 5540.850501
[401,] 586.628153 -2058.527317
[402,] 815.992878 586.628153
[403,] -3252.172767 815.992878
[404,] -1639.129780 -3252.172767
[405,] -1155.084059 -1639.129780
[406,] -682.175458 -1155.084059
[407,] 161.294988 -682.175458
[408,] 266.453118 161.294988
[409,] -5000.528162 266.453118
[410,] 4958.071007 -5000.528162
[411,] -937.482146 4958.071007
[412,] 1675.626363 -937.482146
[413,] 5109.128002 1675.626363
[414,] -10702.652116 5109.128002
[415,] -5180.946020 -10702.652116
[416,] 11155.125275 -5180.946020
[417,] -443.565059 11155.125275
[418,] -3332.403928 -443.565059
[419,] -3668.241362 -3332.403928
[420,] -2348.678631 -3668.241362
[421,] 5066.650137 -2348.678631
[422,] 3417.018967 5066.650137
[423,] -14401.817777 3417.018967
[424,] 4095.498813 -14401.817777
[425,] -2470.243562 4095.498813
[426,] 24153.403695 -2470.243562
[427,] -12733.619516 24153.403695
[428,] 7888.505936 -12733.619516
[429,] 17046.328143 7888.505936
[430,] 153933.352077 17046.328143
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 3189.149011 18112.030820
2 -5975.240229 3189.149011
3 -5896.762010 -5975.240229
4 -4963.274386 -5896.762010
5 -3586.218801 -4963.274386
6 -3639.496594 -3586.218801
7 -1771.217901 -3639.496594
8 -5120.670902 -1771.217901
9 -7276.418018 -5120.670902
10 3969.352621 -7276.418018
11 -4492.098869 3969.352621
12 -10688.605698 -4492.098869
13 -3692.210291 -10688.605698
14 -3658.988881 -3692.210291
15 -571.565104 -3658.988881
16 8580.303693 -571.565104
17 -2205.676265 8580.303693
18 -7792.585933 -2205.676265
19 749.071535 -7792.585933
20 -3340.109525 749.071535
21 -3344.907924 -3340.109525
22 1462.809853 -3344.907924
23 760.744286 1462.809853
24 2502.076015 760.744286
25 -3164.005247 2502.076015
26 -3829.109271 -3164.005247
27 -2410.225134 -3829.109271
28 -5503.404364 -2410.225134
29 -3004.098847 -5503.404364
30 -1916.099554 -3004.098847
31 -860.802328 -1916.099554
32 -2790.326038 -860.802328
33 -779.212044 -2790.326038
34 8916.934765 -779.212044
35 -407.005211 8916.934765
36 -2080.687204 -407.005211
37 -935.992092 -2080.687204
38 3988.236079 -935.992092
39 346.561017 3988.236079
40 1884.552701 346.561017
41 -220.002535 1884.552701
42 -705.251195 -220.002535
43 1014.684702 -705.251195
44 -1268.559636 1014.684702
45 446.715966 -1268.559636
46 5762.193305 446.715966
47 1607.948969 5762.193305
48 818.275166 1607.948969
49 329.161773 818.275166
50 -121.639573 329.161773
51 1192.597994 -121.639573
52 3825.734858 1192.597994
53 1898.113007 3825.734858
54 3883.404235 1898.113007
55 -1663.508514 3883.404235
56 -3703.968891 -1663.508514
57 -777.589508 -3703.968891
58 656.930989 -777.589508
59 -707.185995 656.930989
60 -745.130968 -707.185995
61 -164.659171 -745.130968
62 -120.864169 -164.659171
63 -715.556170 -120.864169
64 -1514.636577 -715.556170
65 -650.065442 -1514.636577
66 217.786117 -650.065442
67 318.539775 217.786117
68 1989.611071 318.539775
69 19.621278 1989.611071
70 -2214.551765 19.621278
71 -2783.762310 -2214.551765
72 1173.976909 -2783.762310
73 -213.043438 1173.976909
74 1246.830179 -213.043438
75 82.895705 1246.830179
76 1252.997086 82.895705
77 22.872089 1252.997086
78 -676.919991 22.872089
79 1592.900422 -676.919991
80 -1299.502738 1592.900422
81 -530.426761 -1299.502738
82 2538.618480 -530.426761
83 -2052.259268 2538.618480
84 635.543149 -2052.259268
85 -1470.195210 635.543149
86 -656.375725 -1470.195210
87 829.144472 -656.375725
88 1120.037178 829.144472
89 609.702605 1120.037178
90 91.538895 609.702605
91 -504.930538 91.538895
92 -1849.356020 -504.930538
93 1524.364018 -1849.356020
94 -2564.966168 1524.364018
95 -1618.168056 -2564.966168
96 849.475477 -1618.168056
97 -761.256688 849.475477
98 1699.459918 -761.256688
99 -533.611160 1699.459918
100 -328.546992 -533.611160
101 821.087835 -328.546992
102 1179.268373 821.087835
103 1533.402505 1179.268373
104 -1072.452550 1533.402505
105 -304.685718 -1072.452550
106 -199.410741 -304.685718
107 -481.158937 -199.410741
108 561.232612 -481.158937
109 -191.644522 561.232612
110 -1300.599370 -191.644522
111 -985.347261 -1300.599370
112 10.897599 -985.347261
113 -399.340499 10.897599
114 -1661.614585 -399.340499
115 -1532.042011 -1661.614585
116 11.808421 -1532.042011
117 -562.540244 11.808421
118 202.842125 -562.540244
119 -1088.325075 202.842125
120 375.001446 -1088.325075
121 -914.745069 375.001446
122 179.341066 -914.745069
123 185.809164 179.341066
124 -1335.232702 185.809164
125 47.161008 -1335.232702
126 -1061.735897 47.161008
127 -1061.158843 -1061.735897
128 -1245.133403 -1061.158843
129 -1061.735897 -1245.133403
130 -1063.834583 -1061.735897
131 -1725.494028 -1063.834583
132 -1061.735897 -1725.494028
133 -1045.872342 -1061.735897
134 -1082.659339 -1045.872342
135 -1090.258677 -1082.659339
136 -1080.264849 -1090.258677
137 -1059.405046 -1080.264849
138 -1114.409716 -1059.405046
139 -1061.735897 -1114.409716
140 -1231.006356 -1061.735897
141 -1061.735897 -1231.006356
142 -1223.175503 -1061.735897
143 -1058.869179 -1223.175503
144 -1047.115252 -1058.869179
145 -1274.122439 -1047.115252
146 -1061.735897 -1274.122439
147 -1055.892234 -1061.735897
148 -1135.815656 -1055.892234
149 -1050.887300 -1135.815656
150 -1059.353040 -1050.887300
151 -1138.519584 -1059.353040
152 -1062.843185 -1138.519584
153 -1092.243375 -1062.843185
154 -1051.874989 -1092.243375
155 -1061.735897 -1051.874989
156 -1061.735897 -1061.735897
157 -1029.374768 -1061.735897
158 -1203.209376 -1029.374768
159 -1061.735897 -1203.209376
160 -1061.735897 -1061.735897
161 -1061.735897 -1061.735897
162 -1061.735897 -1061.735897
163 -1034.423174 -1061.735897
164 -893.054439 -1034.423174
165 -1061.927612 -893.054439
166 -1061.735897 -1061.927612
167 -2242.074409 -1061.735897
168 -1508.435528 -2242.074409
169 -1061.735897 -1508.435528
170 -1061.735897 -1061.735897
171 -1005.825011 -1061.735897
172 -1061.735897 -1005.825011
173 -982.565463 -1061.735897
174 -1061.735897 -982.565463
175 -1060.752605 -1061.735897
176 -1098.819809 -1060.752605
177 -1069.429333 -1098.819809
178 -1061.735897 -1069.429333
179 -1061.735897 -1061.735897
180 -841.713596 -1061.735897
181 -643.649635 -841.713596
182 -1453.169244 -643.649635
183 -434.857855 -1453.169244
184 -970.539261 -434.857855
185 -1223.462113 -970.539261
186 -1214.789853 -1223.462113
187 -995.191188 -1214.789853
188 -692.406099 -995.191188
189 -1061.735897 -692.406099
190 -1159.573192 -1061.735897
191 -1061.735897 -1159.573192
192 -1026.952296 -1061.735897
193 -1061.735897 -1026.952296
194 -1057.512909 -1061.735897
195 -1061.735897 -1057.512909
196 670.444619 -1061.735897
197 -1061.735897 670.444619
198 -1061.735897 -1061.735897
199 122.450100 -1061.735897
200 -1061.735897 122.450100
201 -1061.735897 -1061.735897
202 -1061.735897 -1061.735897
203 -1481.573495 -1061.735897
204 160.496827 -1481.573495
205 -1021.841361 160.496827
206 -839.829927 -1021.841361
207 -1061.735897 -839.829927
208 -1900.832563 -1061.735897
209 -1029.650703 -1900.832563
210 -139.270215 -1029.650703
211 127.729785 -139.270215
212 1739.108604 127.729785
213 -1184.234238 1739.108604
214 -1096.934730 -1184.234238
215 -124.391717 -1096.934730
216 247.952847 -124.391717
217 -7.271879 247.952847
218 1726.608283 -7.271879
219 1484.825340 1726.608283
220 -1015.125087 1484.825340
221 -136.270215 -1015.125087
222 332.260570 -136.270215
223 -279.024713 332.260570
224 129.498368 -279.024713
225 -253.640739 129.498368
226 -180.646183 -253.640739
227 52.729785 -180.646183
228 -1298.125172 52.729785
229 177.729785 -1298.125172
230 -930.761739 177.729785
231 155.492980 -930.761739
232 228.729785 155.492980
233 -952.389572 228.729785
234 -137.270215 -952.389572
235 -139.093865 -137.270215
236 -86.270215 -139.093865
237 -37.133233 -86.270215
238 -139.270215 -37.133233
239 -139.270215 -139.270215
240 -45.270215 -139.270215
241 -403.968451 -45.270215
242 -115.270215 -403.968451
243 2404.716804 -115.270215
244 -2344.517538 2404.716804
245 -139.270215 -2344.517538
246 -8.270215 -139.270215
247 -695.891470 -8.270215
248 -139.270215 -695.891470
249 66.729785 -139.270215
250 625.561964 66.729785
251 27.729785 625.561964
252 -694.834735 27.729785
253 -1076.587638 -694.834735
254 -8060.784634 -1076.587638
255 225.729785 -8060.784634
256 1244.262432 225.729785
257 1908.386400 1244.262432
258 -139.270215 1908.386400
259 -139.270215 -139.270215
260 -139.270215 -139.270215
261 86.729785 -139.270215
262 -1109.947060 86.729785
263 1089.257151 -1109.947060
264 -139.270215 1089.257151
265 -139.270215 -139.270215
266 48.729785 -139.270215
267 -649.474595 48.729785
268 -1.270215 -649.474595
269 -1365.460633 -1.270215
270 -139.270215 -1365.460633
271 -139.270215 -139.270215
272 -14.270215 -139.270215
273 -412.769924 -14.270215
274 142.729785 -412.769924
275 39.764643 142.729785
276 -658.661062 39.764643
277 1184.729785 -658.661062
278 -2612.549884 1184.729785
279 -2014.058505 -2612.549884
280 -139.270215 -2014.058505
281 109.729785 -139.270215
282 650.616351 109.729785
283 193.729785 650.616351
284 -1589.334055 193.729785
285 461.729785 -1589.334055
286 79.057971 461.729785
287 -299.019094 79.057971
288 109.729785 -299.019094
289 162.247354 109.729785
290 25.729785 162.247354
291 325.753632 25.729785
292 -2168.191760 325.753632
293 -86.270215 -2168.191760
294 -127.865024 -86.270215
295 1420.496699 -127.865024
296 -139.270215 1420.496699
297 -139.270215 -139.270215
298 -139.270215 -139.270215
299 -109.816972 -139.270215
300 24.072522 -109.816972
301 -412.072342 24.072522
302 -139.330966 -412.072342
303 226.729785 -139.330966
304 1051.433365 226.729785
305 -142.631572 1051.433365
306 69.729785 -142.631572
307 -423.947404 69.729785
308 923.587931 -423.947404
309 -139.270215 923.587931
310 225.729785 -139.270215
311 -4465.476016 225.729785
312 -90.270215 -4465.476016
313 -559.628047 -90.270215
314 -2.242662 -559.628047
315 -89.365450 -2.242662
316 370.348242 -89.365450
317 -1136.768660 370.348242
318 -139.746521 -1136.768660
319 -139.270215 -139.746521
320 -139.270215 -139.270215
321 -139.270215 -139.270215
322 -139.270215 -139.270215
323 -139.270215 -139.270215
324 -117.270215 -139.270215
325 -281.518718 -117.270215
326 -139.270215 -281.518718
327 -139.270215 -139.270215
328 -139.270215 -139.270215
329 -139.270215 -139.270215
330 -139.270215 -139.270215
331 -139.270215 -139.270215
332 -43.270215 -139.270215
333 -144.309852 -43.270215
334 173.512840 -144.309852
335 704.043521 173.512840
336 -1133.525590 704.043521
337 -113.270215 -1133.525590
338 -127.184237 -113.270215
339 496.643627 -127.184237
340 -1613.230251 496.643627
341 -139.270215 -1613.230251
342 -139.270215 -139.270215
343 481.729785 -139.270215
344 -490.925235 481.729785
345 -20.270215 -490.925235
346 -361.063762 -20.270215
347 -139.270215 -361.063762
348 1455.729785 -139.270215
349 368.585152 1455.729785
350 553.501829 368.585152
351 292.282379 553.501829
352 752.484755 292.282379
353 -494.815251 752.484755
354 -139.270215 -494.815251
355 374.729785 -139.270215
356 -3795.045538 374.729785
357 -139.270215 -3795.045538
358 -139.270215 -139.270215
359 -138.270215 -139.270215
360 -139.069613 -138.270215
361 -139.270215 -139.069613
362 1623.729785 -139.270215
363 -502.779063 1623.729785
364 150.011155 -502.779063
365 -139.270215 150.011155
366 -139.270215 -139.270215
367 -139.270215 -139.270215
368 78.729785 -139.270215
369 -131.175020 78.729785
370 308.729785 -131.175020
371 688.326689 308.729785
372 -2445.155383 688.326689
373 -1439.270720 -2445.155383
374 -139.270215 -1439.270720
375 -18.270215 -139.270215
376 487.406100 -18.270215
377 2477.986548 487.406100
378 80.964476 2477.986548
379 -139.270215 80.964476
380 390.729785 -139.270215
381 -2290.634823 390.729785
382 2398.483194 -2290.634823
383 -61.270215 2398.483194
384 2006.554552 -61.270215
385 -5383.735718 2006.554552
386 1593.866401 -5383.735718
387 85.075320 1593.866401
388 446.233013 85.075320
389 364.247630 446.233013
390 211.575315 364.247630
391 -353.489569 211.575315
392 -2062.862285 -353.489569
393 4234.482677 -2062.862285
394 -6199.111883 4234.482677
395 3712.798981 -6199.111883
396 -480.110115 3712.798981
397 -6949.088627 -480.110115
398 347.097654 -6949.088627
399 5540.850501 347.097654
400 -2058.527317 5540.850501
401 586.628153 -2058.527317
402 815.992878 586.628153
403 -3252.172767 815.992878
404 -1639.129780 -3252.172767
405 -1155.084059 -1639.129780
406 -682.175458 -1155.084059
407 161.294988 -682.175458
408 266.453118 161.294988
409 -5000.528162 266.453118
410 4958.071007 -5000.528162
411 -937.482146 4958.071007
412 1675.626363 -937.482146
413 5109.128002 1675.626363
414 -10702.652116 5109.128002
415 -5180.946020 -10702.652116
416 11155.125275 -5180.946020
417 -443.565059 11155.125275
418 -3332.403928 -443.565059
419 -3668.241362 -3332.403928
420 -2348.678631 -3668.241362
421 5066.650137 -2348.678631
422 3417.018967 5066.650137
423 -14401.817777 3417.018967
424 4095.498813 -14401.817777
425 -2470.243562 4095.498813
426 24153.403695 -2470.243562
427 -12733.619516 24153.403695
428 7888.505936 -12733.619516
429 17046.328143 7888.505936
430 153933.352077 17046.328143
> plot(z,main=paste('Residual Lag plot, lowess, and regression line'), ylab='values of Residuals', xlab='lagged values of Residuals')
> lines(lowess(z))
> abline(lm(z))
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/70ok61291297725.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> acf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Autocorrelation Function')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/80ok61291297725.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> pacf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Partial Autocorrelation Function')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/90ok61291297725.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> opar <- par(mfrow = c(2,2), oma = c(0, 0, 1.1, 0))
> plot(mylm, las = 1, sub='Residual Diagnostics')
> par(opar)
> dev.off()
null device
1
> if (n > n25) {
+ postscript(file="/var/www/html/rcomp/tmp/10bf291291297725.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
+ plot(kp3:nmkm3,gqarr[,2], main='Goldfeld-Quandt test',ylab='2-sided p-value',xlab='breakpoint')
+ grid()
+ dev.off()
+ }
null device
1
>
> #Note: the /var/www/html/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/www/html/rcomp/createtable")
>
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Estimated Regression Equation', 1, TRUE)
> a<-table.row.end(a)
> myeq <- colnames(x)[1]
> myeq <- paste(myeq, '[t] = ', sep='')
> for (i in 1:k){
+ if (mysum$coefficients[i,1] > 0) myeq <- paste(myeq, '+', '')
+ myeq <- paste(myeq, mysum$coefficients[i,1], sep=' ')
+ if (rownames(mysum$coefficients)[i] != '(Intercept)') {
+ myeq <- paste(myeq, rownames(mysum$coefficients)[i], sep='')
+ if (rownames(mysum$coefficients)[i] != 't') myeq <- paste(myeq, '[t]', sep='')
+ }
+ }
> myeq <- paste(myeq, ' + e[t]')
> a<-table.row.start(a)
> a<-table.element(a, myeq)
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/11wx0x1291297725.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,hyperlink('http://www.xycoon.com/ols1.htm','Multiple Linear Regression - Ordinary Least Squares',''), 6, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Variable',header=TRUE)
> a<-table.element(a,'Parameter',header=TRUE)
> a<-table.element(a,'S.D.',header=TRUE)
> a<-table.element(a,'T-STAT
H0: parameter = 0',header=TRUE)
> a<-table.element(a,'2-tail p-value',header=TRUE)
> a<-table.element(a,'1-tail p-value',header=TRUE)
> a<-table.row.end(a)
> for (i in 1:k){
+ a<-table.row.start(a)
+ a<-table.element(a,rownames(mysum$coefficients)[i],header=TRUE)
+ a<-table.element(a,mysum$coefficients[i,1])
+ a<-table.element(a, round(mysum$coefficients[i,2],6))
+ a<-table.element(a, round(mysum$coefficients[i,3],4))
+ a<-table.element(a, round(mysum$coefficients[i,4],6))
+ a<-table.element(a, round(mysum$coefficients[i,4]/2,6))
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/12zgzl1291297725.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Regression Statistics', 2, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple R',1,TRUE)
> a<-table.element(a, sqrt(mysum$r.squared))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'R-squared',1,TRUE)
> a<-table.element(a, mysum$r.squared)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Adjusted R-squared',1,TRUE)
> a<-table.element(a, mysum$adj.r.squared)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (value)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[1])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (DF numerator)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[2])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (DF denominator)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[3])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'p-value',1,TRUE)
> a<-table.element(a, 1-pf(mysum$fstatistic[1],mysum$fstatistic[2],mysum$fstatistic[3]))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Residual Statistics', 2, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Residual Standard Deviation',1,TRUE)
> a<-table.element(a, mysum$sigma)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Sum Squared Residuals',1,TRUE)
> a<-table.element(a, sum(myerror*myerror))
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/13v7wc1291297725.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Actuals, Interpolation, and Residuals', 4, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Time or Index', 1, TRUE)
> a<-table.element(a, 'Actuals', 1, TRUE)
> a<-table.element(a, 'Interpolation
Forecast', 1, TRUE)
> a<-table.element(a, 'Residuals
Prediction Error', 1, TRUE)
> a<-table.row.end(a)
> for (i in 1:n) {
+ a<-table.row.start(a)
+ a<-table.element(a,i, 1, TRUE)
+ a<-table.element(a,x[i])
+ a<-table.element(a,x[i]-mysum$resid[i])
+ a<-table.element(a,mysum$resid[i])
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/14hqv01291297725.tab")
> if (n > n25) {
+ a<-table.start()
+ a<-table.row.start(a)
+ a<-table.element(a,'Goldfeld-Quandt test for Heteroskedasticity',4,TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'p-values',header=TRUE)
+ a<-table.element(a,'Alternative Hypothesis',3,header=TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'breakpoint index',header=TRUE)
+ a<-table.element(a,'greater',header=TRUE)
+ a<-table.element(a,'2-sided',header=TRUE)
+ a<-table.element(a,'less',header=TRUE)
+ a<-table.row.end(a)
+ for (mypoint in kp3:nmkm3) {
+ a<-table.row.start(a)
+ a<-table.element(a,mypoint,header=TRUE)
+ a<-table.element(a,gqarr[mypoint-kp3+1,1])
+ a<-table.element(a,gqarr[mypoint-kp3+1,2])
+ a<-table.element(a,gqarr[mypoint-kp3+1,3])
+ a<-table.row.end(a)
+ }
+ a<-table.end(a)
+ table.save(a,file="/var/www/html/rcomp/tmp/15kqb51291297725.tab")
+ a<-table.start()
+ a<-table.row.start(a)
+ a<-table.element(a,'Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity',4,TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'Description',header=TRUE)
+ a<-table.element(a,'# significant tests',header=TRUE)
+ a<-table.element(a,'% significant tests',header=TRUE)
+ a<-table.element(a,'OK/NOK',header=TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'1% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant1)
+ a<-table.element(a,numsignificant1/numgqtests)
+ if (numsignificant1/numgqtests < 0.01) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'5% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant5)
+ a<-table.element(a,numsignificant5/numgqtests)
+ if (numsignificant5/numgqtests < 0.05) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'10% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant10)
+ a<-table.element(a,numsignificant10/numgqtests)
+ if (numsignificant10/numgqtests < 0.1) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.end(a)
+ table.save(a,file="/var/www/html/rcomp/tmp/166rat1291297725.tab")
+ }
>
> try(system("convert tmp/1menx1291297725.ps tmp/1menx1291297725.png",intern=TRUE))
character(0)
> try(system("convert tmp/2menx1291297725.ps tmp/2menx1291297725.png",intern=TRUE))
character(0)
> try(system("convert tmp/3f5401291297725.ps tmp/3f5401291297725.png",intern=TRUE))
character(0)
> try(system("convert tmp/4f5401291297725.ps tmp/4f5401291297725.png",intern=TRUE))
character(0)
> try(system("convert tmp/5f5401291297725.ps tmp/5f5401291297725.png",intern=TRUE))
character(0)
> try(system("convert tmp/6pell1291297725.ps tmp/6pell1291297725.png",intern=TRUE))
character(0)
> try(system("convert tmp/70ok61291297725.ps tmp/70ok61291297725.png",intern=TRUE))
character(0)
> try(system("convert tmp/80ok61291297725.ps tmp/80ok61291297725.png",intern=TRUE))
character(0)
> try(system("convert tmp/90ok61291297725.ps tmp/90ok61291297725.png",intern=TRUE))
character(0)
> try(system("convert tmp/10bf291291297725.ps tmp/10bf291291297725.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
12.029 2.177 25.622