R version 2.12.0 (2010-10-15)
Copyright (C) 2010 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: x86_64-redhat-linux-gnu (64-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> x <- array(list(1
+ ,1595
+ ,17
+ ,60720
+ ,319369
+ ,0
+ ,0
+ ,5565
+ ,47
+ ,94355
+ ,493408
+ ,39
+ ,0
+ ,601
+ ,6
+ ,60720
+ ,319210
+ ,0
+ ,1
+ ,188
+ ,11
+ ,77655
+ ,381180
+ ,0
+ ,1
+ ,7146
+ ,105
+ ,134028
+ ,297978
+ ,305
+ ,0
+ ,1135
+ ,17
+ ,62285
+ ,290476
+ ,-9
+ ,1
+ ,450
+ ,8
+ ,59325
+ ,292136
+ ,-1
+ ,1
+ ,34
+ ,8
+ ,60630
+ ,314353
+ ,0
+ ,1
+ ,133
+ ,14
+ ,65990
+ ,339445
+ ,0
+ ,1
+ ,119
+ ,6
+ ,59118
+ ,303677
+ ,0
+ ,0
+ ,2053
+ ,53
+ ,100423
+ ,397144
+ ,33
+ ,0
+ ,4036
+ ,63
+ ,100269
+ ,424898
+ ,32
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,4655
+ ,393
+ ,60720
+ ,341570
+ ,0
+ ,1
+ ,131
+ ,11
+ ,61808
+ ,308989
+ ,0
+ ,1
+ ,1766
+ ,73
+ ,60438
+ ,305959
+ ,0
+ ,1
+ ,312
+ ,14
+ ,58598
+ ,318690
+ ,0
+ ,1
+ ,448
+ ,40
+ ,59781
+ ,323361
+ ,0
+ ,1
+ ,115
+ ,13
+ ,60945
+ ,318903
+ ,0
+ ,1
+ ,60
+ ,10
+ ,61124
+ ,314049
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,364
+ ,35
+ ,47705
+ ,298466
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,1442
+ ,118
+ ,121173
+ ,438493
+ ,-3
+ ,0
+ ,1389
+ ,9
+ ,57530
+ ,378049
+ ,6
+ ,1
+ ,149
+ ,5
+ ,62041
+ ,313332
+ ,0
+ ,1
+ ,2212
+ ,14
+ ,60720
+ ,319864
+ ,0
+ ,0
+ ,7489
+ ,322
+ ,136996
+ ,430866
+ ,70
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,402
+ ,26
+ ,84990
+ ,332743
+ ,5
+ ,0
+ ,7419
+ ,29
+ ,63255
+ ,286963
+ ,-5
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,307
+ ,15
+ ,58856
+ ,331955
+ ,0
+ ,0
+ ,1134
+ ,162
+ ,146216
+ ,527021
+ ,58
+ ,0
+ ,7561
+ ,87
+ ,82425
+ ,364304
+ ,37
+ ,1
+ ,5131
+ ,71
+ ,86111
+ ,292154
+ ,9
+ ,1
+ ,9
+ ,6
+ ,60835
+ ,314987
+ ,0
+ ,1
+ ,2264
+ ,24
+ ,55174
+ ,272713
+ ,-31
+ ,1
+ ,40949
+ ,481
+ ,129352
+ ,1398893
+ ,2865
+ ,0
+ ,4980
+ ,91
+ ,113521
+ ,409642
+ ,0
+ ,1
+ ,272
+ ,8
+ ,112995
+ ,429112
+ ,14
+ ,1
+ ,757
+ ,19
+ ,45689
+ ,382712
+ ,0
+ ,1
+ ,1424
+ ,94
+ ,131116
+ ,374943
+ ,32
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,378
+ ,52
+ ,71873
+ ,297413
+ ,0
+ ,1
+ ,1061
+ ,5
+ ,70415
+ ,304252
+ ,-4
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,203
+ ,9
+ ,61938
+ ,309596
+ ,0
+ ,0
+ ,3689
+ ,25
+ ,101792
+ ,377305
+ ,74
+ ,1
+ ,567
+ ,11
+ ,103345
+ ,386212
+ ,20
+ ,0
+ ,27664
+ ,282
+ ,74996
+ ,657954
+ ,247
+ ,1
+ ,1686
+ ,40
+ ,68696
+ ,368186
+ ,11
+ ,0
+ ,1164
+ ,17
+ ,57320
+ ,269753
+ ,-6
+ ,1
+ ,15824
+ ,64
+ ,63838
+ ,585715
+ ,284
+ ,1
+ ,6575
+ ,159
+ ,151352
+ ,1071292
+ ,819
+ ,0
+ ,37597
+ ,248
+ ,37884
+ ,992426
+ ,652
+ ,1
+ ,314
+ ,5
+ ,60720
+ ,315371
+ ,0
+ ,0
+ ,1932
+ ,36
+ ,72954
+ ,249898
+ ,10
+ ,0
+ ,5692
+ ,72
+ ,64239
+ ,357312
+ ,22
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,6533
+ ,75
+ ,126910
+ ,364839
+ ,16
+ ,0
+ ,2126
+ ,138
+ ,21001
+ ,230621
+ ,-3
+ ,1
+ ,14
+ ,1
+ ,61016
+ ,315877
+ ,0
+ ,0
+ ,3192
+ ,60
+ ,59831
+ ,125390
+ ,10
+ ,1
+ ,47
+ ,3
+ ,60720
+ ,314882
+ ,0
+ ,0
+ ,4891
+ ,217
+ ,210568
+ ,370837
+ ,73
+ ,0
+ ,531
+ ,10
+ ,71641
+ ,330068
+ ,0
+ ,1
+ ,2332
+ ,4
+ ,72680
+ ,324385
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,16275
+ ,109
+ ,47818
+ ,1443586
+ ,1974
+ ,1
+ ,9252
+ ,48
+ ,79510
+ ,253537
+ ,-129
+ ,0
+ ,29790
+ ,309
+ ,81767
+ ,4321023
+ ,28328
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,2300
+ ,29
+ ,64873
+ ,352108
+ ,0
+ ,1
+ ,267
+ ,32
+ ,57640
+ ,330059
+ ,0
+ ,0
+ ,382
+ ,21
+ ,59365
+ ,340968
+ ,0
+ ,1
+ ,397
+ ,28
+ ,60646
+ ,317736
+ ,0
+ ,1
+ ,11920
+ ,182
+ ,31080
+ ,209458
+ ,238
+ ,0
+ ,54660
+ ,302
+ ,126942
+ ,1491348
+ ,3040
+ ,1
+ ,5
+ ,2
+ ,60720
+ ,314887
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,2647
+ ,37
+ ,62920
+ ,353058
+ ,1
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,94
+ ,7
+ ,60793
+ ,314533
+ ,0
+ ,1
+ ,3
+ ,2
+ ,60698
+ ,315354
+ ,0
+ ,1
+ ,454
+ ,51
+ ,63261
+ ,302187
+ ,0
+ ,1
+ ,227
+ ,9
+ ,76644
+ ,336639
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,206
+ ,9
+ ,53110
+ ,296702
+ ,0
+ ,1
+ ,25830
+ ,115
+ ,245546
+ ,1073089
+ ,1910
+ ,1
+ ,1528
+ ,38
+ ,60326
+ ,146494
+ ,38
+ ,1
+ ,271
+ ,2
+ ,60720
+ ,325249
+ ,0
+ ,1
+ ,138
+ ,6
+ ,69817
+ ,331420
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,278
+ ,26
+ ,61404
+ ,314922
+ ,0
+ ,1
+ ,282
+ ,30
+ ,62452
+ ,320016
+ ,0
+ ,1
+ ,571
+ ,80
+ ,41477
+ ,280398
+ ,0
+ ,1
+ ,2253
+ ,99
+ ,63593
+ ,452469
+ ,17
+ ,1
+ ,290
+ ,2
+ ,58790
+ ,301164
+ ,0
+ ,1
+ ,78
+ ,16
+ ,62700
+ ,317330
+ ,0
+ ,1
+ ,20
+ ,2
+ ,60805
+ ,315576
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,18073
+ ,180
+ ,27284
+ ,-7170
+ ,-366
+ ,1
+ ,1866
+ ,29
+ ,56225
+ ,322331
+ ,0
+ ,0
+ ,42634
+ ,156
+ ,157214
+ ,1629616
+ ,2964
+ ,1
+ ,249
+ ,62
+ ,54323
+ ,292754
+ ,0
+ ,1
+ ,422
+ ,28
+ ,57935
+ ,318056
+ ,0
+ ,1
+ ,2675
+ ,79
+ ,59017
+ ,355178
+ ,36
+ ,1
+ ,965
+ ,3
+ ,73490
+ ,204325
+ ,-8
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,621
+ ,16
+ ,68005
+ ,317046
+ ,-5
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,365
+ ,8
+ ,54820
+ ,309560
+ ,0
+ ,0
+ ,3122
+ ,34
+ ,94670
+ ,414462
+ ,38
+ ,0
+ ,12988
+ ,119
+ ,82340
+ ,857217
+ ,601
+ ,0
+ ,5336
+ ,81
+ ,112477
+ ,697458
+ ,132
+ ,0
+ ,3160
+ ,108
+ ,108094
+ ,530670
+ ,66
+ ,1
+ ,2489
+ ,40
+ ,67804
+ ,238125
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,6984
+ ,45
+ ,80570
+ ,741409
+ ,195
+ ,1
+ ,2001
+ ,52
+ ,95551
+ ,393343
+ ,19
+ ,1
+ ,15236
+ ,90
+ ,77440
+ ,372631
+ ,-158
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,530
+ ,60
+ ,73433
+ ,317291
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,35624
+ ,223
+ ,157278
+ ,306275
+ ,-70
+ ,0
+ ,1383
+ ,11
+ ,73221
+ ,317892
+ ,13
+ ,1
+ ,875
+ ,14
+ ,67000
+ ,334280
+ ,3
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,72
+ ,3
+ ,60398
+ ,314210
+ ,0
+ ,0
+ ,265
+ ,31
+ ,60720
+ ,306948
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,335
+ ,33
+ ,64175
+ ,320398
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,4525
+ ,226
+ ,93811
+ ,501749
+ ,101
+ ,0
+ ,3045
+ ,58
+ ,27330
+ ,202055
+ ,31
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,638
+ ,14
+ ,60370
+ ,333210
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,607
+ ,11
+ ,60436
+ ,322340
+ ,0
+ ,1
+ ,1558
+ ,63
+ ,55637
+ ,369448
+ ,0
+ ,0
+ ,1324
+ ,117
+ ,67440
+ ,291841
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,611
+ ,8
+ ,59190
+ ,296919
+ ,-4
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,923
+ ,8
+ ,58620
+ ,309038
+ ,0
+ ,1
+ ,661
+ ,3
+ ,65920
+ ,246541
+ ,-12
+ ,0
+ ,2397
+ ,14
+ ,74020
+ ,289513
+ ,6
+ ,1
+ ,366
+ ,53
+ ,61808
+ ,344425
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,135
+ ,2
+ ,62065
+ ,314210
+ ,0
+ ,0
+ ,1659
+ ,29
+ ,107577
+ ,480382
+ ,2
+ ,1
+ ,316
+ ,9
+ ,60505
+ ,315009
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,309
+ ,11
+ ,56535
+ ,312878
+ ,0
+ ,0
+ ,49
+ ,8
+ ,64107
+ ,322031
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,4519
+ ,49
+ ,102129
+ ,597793
+ ,163
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,837
+ ,69
+ ,61262
+ ,315688
+ ,4
+ ,1
+ ,5119
+ ,49
+ ,39039
+ ,378525
+ ,31
+ ,1
+ ,1280
+ ,117
+ ,69465
+ ,312378
+ ,26
+ ,1
+ ,2564
+ ,22
+ ,130140
+ ,403560
+ ,13
+ ,0
+ ,2045
+ ,20
+ ,97890
+ ,510834
+ ,10
+ ,0
+ ,2234
+ ,16
+ ,77200
+ ,214215
+ ,5
+ ,1
+ ,975
+ ,23
+ ,90534
+ ,235133
+ ,13
+ ,0
+ ,1136
+ ,32
+ ,48522
+ ,343613
+ ,0
+ ,1
+ ,453
+ ,21
+ ,81125
+ ,365959
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,61
+ ,14
+ ,60720
+ ,314551
+ ,0
+ ,1
+ ,368
+ ,17
+ ,61686
+ ,303230
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,4901
+ ,84
+ ,121920
+ ,469107
+ ,67
+ ,1
+ ,540
+ ,14
+ ,103960
+ ,354228
+ ,22
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,2
+ ,9
+ ,60735
+ ,315394
+ ,0
+ ,1
+ ,36
+ ,4
+ ,61564
+ ,312412
+ ,0
+ ,1
+ ,776
+ ,9
+ ,64230
+ ,333505
+ ,-4
+ ,1
+ ,84738
+ ,588
+ ,-26007
+ ,223193
+ ,4145
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,3
+ ,4
+ ,60761
+ ,315656
+ ,0
+ ,1
+ ,529
+ ,9
+ ,63870
+ ,296261
+ ,-3
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,405
+ ,7
+ ,60845
+ ,336425
+ ,0
+ ,1
+ ,972
+ ,34
+ ,71642
+ ,359335
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,2099
+ ,64
+ ,106611
+ ,308636
+ ,20
+ ,0
+ ,3437
+ ,47
+ ,48022
+ ,158492
+ ,9
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,22330
+ ,84
+ ,79801
+ ,711969
+ ,446
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,483
+ ,27
+ ,60830
+ ,306268
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,2239
+ ,21
+ ,88590
+ ,442882
+ ,5
+ ,0
+ ,2949
+ ,41
+ ,82903
+ ,378509
+ ,18
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,365
+ ,83
+ ,87192
+ ,346611
+ ,0
+ ,1
+ ,2461
+ ,57
+ ,55792
+ ,314289
+ ,0
+ ,0
+ ,21950
+ ,519
+ ,114337
+ ,856956
+ ,745
+ ,0
+ ,3294
+ ,23
+ ,75832
+ ,217193
+ ,1
+ ,1
+ ,141
+ ,6
+ ,61630
+ ,315366
+ ,0
+ ,1
+ ,572
+ ,16
+ ,58580
+ ,307930
+ ,-1
+ ,1
+ ,13326
+ ,102
+ ,165548
+ ,702380
+ ,436
+ ,1
+ ,2284
+ ,33
+ ,76403
+ ,194493
+ ,66
+ ,0
+ ,10
+ ,2
+ ,61656
+ ,316155
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,1414
+ ,198
+ ,70184
+ ,330546
+ ,13
+ ,1
+ ,1975
+ ,35
+ ,118881
+ ,394510
+ ,30
+ ,1
+ ,43
+ ,3
+ ,60887
+ ,312846
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,844
+ ,78
+ ,60925
+ ,296139
+ ,0
+ ,1
+ ,304
+ ,15
+ ,62969
+ ,295580
+ ,0
+ ,1
+ ,458
+ ,11
+ ,58625
+ ,297765
+ ,0
+ ,1
+ ,18562
+ ,155
+ ,102313
+ ,377934
+ ,143
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,7123
+ ,109
+ ,288170
+ ,638830
+ ,769
+ ,1
+ ,622
+ ,73
+ ,73007
+ ,304376
+ ,0
+ ,1
+ ,174
+ ,1
+ ,64820
+ ,307424
+ ,0
+ ,1
+ ,2220
+ ,22
+ ,301670
+ ,644190
+ ,154
+ ,1
+ ,121
+ ,24
+ ,56178
+ ,295370
+ ,0
+ ,0
+ ,11819
+ ,85
+ ,106113
+ ,574339
+ ,199
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,125
+ ,12
+ ,60798
+ ,310201
+ ,0
+ ,1
+ ,1182
+ ,13
+ ,70694
+ ,327007
+ ,-5
+ ,1
+ ,1503
+ ,39
+ ,56364
+ ,343466
+ ,6
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,30
+ ,4
+ ,62045
+ ,318098
+ ,0
+ ,0
+ ,3310
+ ,33
+ ,75230
+ ,448243
+ ,-6
+ ,1
+ ,554
+ ,11
+ ,79285
+ ,325738
+ ,30
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,468
+ ,21
+ ,60720
+ ,312161
+ ,0
+ ,1
+ ,4917
+ ,42
+ ,52811
+ ,243650
+ ,-3
+ ,1
+ ,3256
+ ,66
+ ,35250
+ ,407159
+ ,-8
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,125
+ ,13
+ ,59734
+ ,317698
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,22
+ ,2
+ ,60722
+ ,312502
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,514
+ ,52
+ ,78780
+ ,322378
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,10327
+ ,91
+ ,88577
+ ,640273
+ ,135
+ ,1
+ ,13718
+ ,61
+ ,96448
+ ,345783
+ ,104
+ ,1
+ ,3748
+ ,9
+ ,50350
+ ,652925
+ ,53
+ ,1
+ ,14416
+ ,361
+ ,49857
+ ,439798
+ ,220
+ ,0
+ ,1526
+ ,25
+ ,69351
+ ,278990
+ ,-4
+ ,0
+ ,666
+ ,171
+ ,117869
+ ,339836
+ ,25
+ ,1
+ ,2844
+ ,67
+ ,72683
+ ,240897
+ ,-15
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,368
+ ,17
+ ,61167
+ ,297141
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,333
+ ,13
+ ,70811
+ ,331323
+ ,0
+ ,1
+ ,26
+ ,1
+ ,60896
+ ,313880
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,1303
+ ,61
+ ,69863
+ ,309422
+ ,-5
+ ,1
+ ,20
+ ,3
+ ,60938
+ ,315245
+ ,0
+ ,1
+ ,2384
+ ,97
+ ,61348
+ ,405972
+ ,54
+ ,0
+ ,203
+ ,27
+ ,50804
+ ,300962
+ ,0
+ ,0
+ ,71
+ ,1
+ ,60745
+ ,316176
+ ,0
+ ,1
+ ,53
+ ,14
+ ,59506
+ ,302409
+ ,0
+ ,1
+ ,562
+ ,13
+ ,58456
+ ,283587
+ ,-3
+ ,1
+ ,622
+ ,7
+ ,60950
+ ,263276
+ ,-2
+ ,1
+ ,645
+ ,5
+ ,60720
+ ,312075
+ ,0
+ ,0
+ ,1763
+ ,24
+ ,61600
+ ,308336
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,317
+ ,9
+ ,63915
+ ,298700
+ ,0
+ ,1
+ ,1
+ ,6
+ ,60719
+ ,315372
+ ,0
+ ,1
+ ,275
+ ,4
+ ,59500
+ ,318745
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,936
+ ,30
+ ,67939
+ ,408881
+ ,9
+ ,1
+ ,8568
+ ,37
+ ,32168
+ ,786690
+ ,223
+ ,0
+ ,11528
+ ,70
+ ,-14545
+ ,-83265
+ ,-25
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,738
+ ,10
+ ,60720
+ ,321376
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,592
+ ,7
+ ,64270
+ ,276898
+ ,-4
+ ,1
+ ,126
+ ,10
+ ,60951
+ ,315547
+ ,0
+ ,1
+ ,2
+ ,1
+ ,60743
+ ,315487
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,18014
+ ,239
+ ,-1710
+ ,1405225
+ ,943
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,37704
+ ,452
+ ,60448
+ ,983660
+ ,2167
+ ,0
+ ,63
+ ,4
+ ,65688
+ ,318574
+ ,0
+ ,1
+ ,1431
+ ,34
+ ,106885
+ ,310768
+ ,30
+ ,1
+ ,94
+ ,3
+ ,61360
+ ,312887
+ ,0
+ ,1
+ ,192
+ ,4
+ ,65276
+ ,312339
+ ,0
+ ,1
+ ,32
+ ,1
+ ,59988
+ ,314964
+ ,0
+ ,1
+ ,6869
+ ,105
+ ,117520
+ ,379983
+ ,90
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,1
+ ,4
+ ,60722
+ ,315398
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,2328
+ ,94
+ ,82732
+ ,253588
+ ,0
+ ,0
+ ,209
+ ,13
+ ,64016
+ ,316647
+ ,0
+ ,1
+ ,28
+ ,2
+ ,60890
+ ,315688
+ ,0
+ ,1
+ ,176
+ ,28
+ ,68136
+ ,310670
+ ,0
+ ,1
+ ,1920
+ ,15
+ ,79420
+ ,165404
+ ,-6
+ ,0
+ ,87550
+ ,458
+ ,153198
+ ,4111912
+ ,17936
+ ,1
+ ,520
+ ,10
+ ,58650
+ ,291650
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,1013
+ ,2
+ ,69770
+ ,253468
+ ,-14
+ ,1
+ ,15
+ ,8
+ ,60831
+ ,315688
+ ,0
+ ,1
+ ,587
+ ,12
+ ,59595
+ ,325699
+ ,0
+ ,1
+ ,5371
+ ,66
+ ,87720
+ ,446211
+ ,49
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,1012
+ ,47
+ ,114768
+ ,368078
+ ,26
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,876
+ ,59
+ ,138971
+ ,352850
+ ,24
+ ,0
+ ,162556
+ ,1081
+ ,213118
+ ,6282154
+ ,56289
+ ,1
+ ,43556
+ ,247
+ ,32648
+ ,227132
+ ,929
+ ,0
+ ,3425
+ ,12
+ ,83620
+ ,283910
+ ,-3
+ ,1
+ ,810
+ ,43
+ ,74015
+ ,236761
+ ,26
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,2365
+ ,79
+ ,191778
+ ,550608
+ ,120
+ ,0
+ ,1261
+ ,33
+ ,76114
+ ,307528
+ ,13
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,1585
+ ,107
+ ,93099
+ ,355864
+ ,2
+ ,1
+ ,16189
+ ,295
+ ,116384
+ ,358589
+ ,379
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,10579
+ ,92
+ ,110309
+ ,375195
+ ,4
+ ,1
+ ,474
+ ,20
+ ,61977
+ ,288985
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,3642
+ ,86
+ ,90262
+ ,458343
+ ,22
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,472
+ ,20
+ ,80045
+ ,269587
+ ,0
+ ,1
+ ,98
+ ,2
+ ,61490
+ ,315236
+ ,0
+ ,1
+ ,3999
+ ,30
+ ,51252
+ ,42754
+ ,-58
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,621
+ ,48
+ ,60720
+ ,308256
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,30
+ ,4
+ ,60798
+ ,313164
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,746
+ ,23
+ ,71561
+ ,269661
+ ,10
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,4150
+ ,150
+ ,134759
+ ,518365
+ ,66
+ ,0
+ ,4658
+ ,158
+ ,156608
+ ,233773
+ ,58
+ ,1
+ ,814
+ ,18
+ ,39625
+ ,301881
+ ,23
+ ,1
+ ,1002
+ ,15
+ ,56750
+ ,298568
+ ,0
+ ,0
+ ,496
+ ,17
+ ,87390
+ ,325479
+ ,0
+ ,1
+ ,389
+ ,14
+ ,58990
+ ,325506
+ ,0
+ ,1
+ ,12679
+ ,110
+ ,48020
+ ,984885
+ ,1871
+ ,1
+ ,400
+ ,152
+ ,60720
+ ,313267
+ ,0
+ ,1
+ ,53
+ ,6
+ ,60349
+ ,315793
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,5109
+ ,97
+ ,67038
+ ,215362
+ ,-31
+ ,1
+ ,576
+ ,18
+ ,113761
+ ,314073
+ ,12
+ ,1
+ ,438
+ ,11
+ ,58320
+ ,298096
+ ,-5
+ ,1
+ ,165
+ ,24
+ ,62841
+ ,325176
+ ,0
+ ,1
+ ,4069
+ ,145
+ ,79804
+ ,207393
+ ,119
+ ,0
+ ,23
+ ,5
+ ,62555
+ ,314806
+ ,0
+ ,1
+ ,1285
+ ,107
+ ,99489
+ ,341340
+ ,49
+ ,0
+ ,4677
+ ,75
+ ,90131
+ ,426280
+ ,24
+ ,0
+ ,24811
+ ,505
+ ,95350
+ ,929118
+ ,822
+ ,1
+ ,167
+ ,2
+ ,64245
+ ,307322
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,4628
+ ,62
+ ,69159
+ ,387475
+ ,85
+ ,1
+ ,226
+ ,5
+ ,65745
+ ,291787
+ ,0
+ ,0
+ ,1765
+ ,20
+ ,77623
+ ,247060
+ ,28
+ ,1
+ ,460
+ ,15
+ ,63346
+ ,329784
+ ,0
+ ,1
+ ,36
+ ,3
+ ,60894
+ ,315834
+ ,0
+ ,1
+ ,989
+ ,11
+ ,58930
+ ,304555
+ ,-5
+ ,0
+ ,3055
+ ,43
+ ,60247
+ ,376641
+ ,14
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,13253
+ ,947
+ ,90829
+ ,357760
+ ,35
+ ,1
+ ,24
+ ,8
+ ,59818
+ ,315637
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,131
+ ,6
+ ,59661
+ ,327071
+ ,0
+ ,1
+ ,514
+ ,15
+ ,102725
+ ,377516
+ ,23
+ ,1
+ ,3366
+ ,80
+ ,108479
+ ,299243
+ ,12
+ ,0
+ ,9327
+ ,101
+ ,87419
+ ,387699
+ ,44
+ ,1
+ ,384
+ ,22
+ ,54683
+ ,309836
+ ,0
+ ,1
+ ,17821
+ ,373
+ ,122844
+ ,444477
+ ,643
+ ,1
+ ,397
+ ,10
+ ,62710
+ ,322327
+ ,-1
+ ,1
+ ,897
+ ,5
+ ,60720
+ ,314913
+ ,0
+ ,1
+ ,218
+ ,21
+ ,60720
+ ,315553
+ ,0
+ ,0
+ ,3369
+ ,28
+ ,87161
+ ,688779
+ ,217
+ ,1
+ ,5702
+ ,61
+ ,101481
+ ,321896
+ ,10
+ ,1
+ ,8636
+ ,855
+ ,128294
+ ,301607
+ ,108
+ ,1
+ ,534
+ ,7
+ ,62620
+ ,304485
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,726
+ ,3
+ ,69980
+ ,231861
+ ,-5
+ ,1
+ ,1380
+ ,26
+ ,60982
+ ,347385
+ ,25
+ ,1
+ ,180
+ ,11
+ ,59635
+ ,316386
+ ,0
+ ,0
+ ,7285
+ ,115
+ ,188873
+ ,491303
+ ,59
+ ,1
+ ,880
+ ,74
+ ,80791
+ ,261216
+ ,34
+ ,1
+ ,1
+ ,2
+ ,60727
+ ,315388
+ ,0
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,96
+ ,5
+ ,60379
+ ,313729
+ ,0
+ ,1
+ ,1889
+ ,45
+ ,37527
+ ,358649
+ ,14
+ ,0
+ ,45187
+ ,353
+ ,234817
+ ,1926517
+ ,2854
+ ,1
+ ,288
+ ,4
+ ,60510
+ ,296656
+ ,0
+ ,1
+ ,1270
+ ,26
+ ,69206
+ ,275311
+ ,-10
+ ,1
+ ,6526
+ ,26
+ ,55830
+ ,-42143
+ ,-28
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,226
+ ,12
+ ,72835
+ ,343929
+ ,0
+ ,1
+ ,694
+ ,8
+ ,68060
+ ,367655
+ ,-2
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,249
+ ,72
+ ,56726
+ ,313491
+ ,0)
+ ,dim=c(6
+ ,431)
+ ,dimnames=list(c('Group'
+ ,'Costs'
+ ,'Trades'
+ ,'Dividends'
+ ,'Wealth'
+ ,'CCScore
')
+ ,1:431))
> y <- array(NA,dim=c(6,431),dimnames=list(c('Group','Costs','Trades','Dividends','Wealth','CCScore
'),1:431))
> for (i in 1:dim(x)[1])
+ {
+ for (j in 1:dim(x)[2])
+ {
+ y[i,j] <- as.numeric(x[i,j])
+ }
+ }
> par20 = ''
> par19 = ''
> par18 = ''
> par17 = ''
> par16 = ''
> par15 = ''
> par14 = ''
> par13 = ''
> par12 = ''
> par11 = ''
> par10 = ''
> par9 = ''
> par8 = ''
> par7 = ''
> par6 = ''
> par5 = ''
> par4 = ''
> par3 = 'No Linear Trend'
> par2 = 'Do not include Seasonal Dummies'
> par1 = '6'
> ylab = ''
> xlab = ''
> main = ''
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Prof. Dr. P. Wessa
> #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #Technical description: Write here your technical program description (don't use hard returns!)
> library(lattice)
> library(lmtest)
Loading required package: zoo
> n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test
> par1 <- as.numeric(par1)
> x <- t(y)
> k <- length(x[1,])
> n <- length(x[,1])
> x1 <- cbind(x[,par1], x[,1:k!=par1])
> mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1])
> colnames(x1) <- mycolnames #colnames(x)[par1]
> x <- x1
> if (par3 == 'First Differences'){
+ x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep='')))
+ for (i in 1:n-1) {
+ for (j in 1:k) {
+ x2[i,j] <- x[i+1,j] - x[i,j]
+ }
+ }
+ x <- x2
+ }
> if (par2 == 'Include Monthly Dummies'){
+ x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep ='')))
+ for (i in 1:11){
+ x2[seq(i,n,12),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> if (par2 == 'Include Quarterly Dummies'){
+ x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep ='')))
+ for (i in 1:3){
+ x2[seq(i,n,4),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> k <- length(x[1,])
> if (par3 == 'Linear Trend'){
+ x <- cbind(x, c(1:n))
+ colnames(x)[k+1] <- 't'
+ }
> x
CCScore\r Group Costs Trades Dividends Wealth
1 0 1 1595 17 60720 319369
2 39 0 5565 47 94355 493408
3 0 0 601 6 60720 319210
4 0 1 188 11 77655 381180
5 305 1 7146 105 134028 297978
6 -9 0 1135 17 62285 290476
7 -1 1 450 8 59325 292136
8 0 1 34 8 60630 314353
9 0 1 133 14 65990 339445
10 0 1 119 6 59118 303677
11 33 0 2053 53 100423 397144
12 32 0 4036 63 100269 424898
13 0 1 0 0 60720 315380
14 0 1 4655 393 60720 341570
15 0 1 131 11 61808 308989
16 0 1 1766 73 60438 305959
17 0 1 312 14 58598 318690
18 0 1 448 40 59781 323361
19 0 1 115 13 60945 318903
20 0 1 60 10 61124 314049
21 0 1 0 0 60720 315380
22 0 1 364 35 47705 298466
23 0 1 0 0 60720 315380
24 -3 0 1442 118 121173 438493
25 6 0 1389 9 57530 378049
26 0 1 149 5 62041 313332
27 0 1 2212 14 60720 319864
28 70 0 7489 322 136996 430866
29 0 1 0 0 60720 315380
30 5 1 402 26 84990 332743
31 -5 0 7419 29 63255 286963
32 0 1 0 0 60720 315380
33 0 1 307 15 58856 331955
34 58 0 1134 162 146216 527021
35 37 0 7561 87 82425 364304
36 9 1 5131 71 86111 292154
37 0 1 9 6 60835 314987
38 -31 1 2264 24 55174 272713
39 2865 1 40949 481 129352 1398893
40 0 0 4980 91 113521 409642
41 14 1 272 8 112995 429112
42 0 1 757 19 45689 382712
43 32 1 1424 94 131116 374943
44 0 1 0 0 60720 315380
45 0 0 378 52 71873 297413
46 -4 1 1061 5 70415 304252
47 0 1 0 0 60720 315380
48 0 1 203 9 61938 309596
49 74 0 3689 25 101792 377305
50 20 1 567 11 103345 386212
51 247 0 27664 282 74996 657954
52 11 1 1686 40 68696 368186
53 -6 0 1164 17 57320 269753
54 284 1 15824 64 63838 585715
55 819 1 6575 159 151352 1071292
56 652 0 37597 248 37884 992426
57 0 1 314 5 60720 315371
58 10 0 1932 36 72954 249898
59 22 0 5692 72 64239 357312
60 0 1 0 0 60720 315380
61 0 1 0 0 60720 315380
62 16 0 6533 75 126910 364839
63 -3 0 2126 138 21001 230621
64 0 1 14 1 61016 315877
65 10 0 3192 60 59831 125390
66 0 1 47 3 60720 314882
67 73 0 4891 217 210568 370837
68 0 0 531 10 71641 330068
69 0 1 2332 4 72680 324385
70 0 1 0 0 60720 315380
71 1974 0 16275 109 47818 1443586
72 -129 1 9252 48 79510 253537
73 28328 0 29790 309 81767 4321023
74 0 1 0 0 60720 315380
75 0 0 2300 29 64873 352108
76 0 1 267 32 57640 330059
77 0 0 382 21 59365 340968
78 0 1 397 28 60646 317736
79 238 1 11920 182 31080 209458
80 3040 0 54660 302 126942 1491348
81 0 1 5 2 60720 314887
82 0 1 0 0 60720 315380
83 0 1 0 0 60720 315380
84 1 0 2647 37 62920 353058
85 0 1 0 0 60720 315380
86 0 0 0 0 60720 315380
87 0 1 94 7 60793 314533
88 0 1 3 2 60698 315354
89 0 1 454 51 63261 302187
90 0 1 227 9 76644 336639
91 0 0 0 0 60720 315380
92 0 1 206 9 53110 296702
93 1910 1 25830 115 245546 1073089
94 38 1 1528 38 60326 146494
95 0 1 271 2 60720 325249
96 0 1 138 6 69817 331420
97 0 1 0 0 60720 315380
98 0 1 278 26 61404 314922
99 0 1 282 30 62452 320016
100 0 1 571 80 41477 280398
101 17 1 2253 99 63593 452469
102 0 1 290 2 58790 301164
103 0 1 78 16 62700 317330
104 0 1 20 2 60805 315576
105 0 1 0 0 60720 315380
106 -366 1 18073 180 27284 -7170
107 0 1 1866 29 56225 322331
108 2964 0 42634 156 157214 1629616
109 0 1 249 62 54323 292754
110 0 1 422 28 57935 318056
111 36 1 2675 79 59017 355178
112 -8 1 965 3 73490 204325
113 0 1 0 0 60720 315380
114 -5 1 621 16 68005 317046
115 0 0 0 0 60720 315380
116 0 1 365 8 54820 309560
117 38 0 3122 34 94670 414462
118 601 0 12988 119 82340 857217
119 132 0 5336 81 112477 697458
120 66 0 3160 108 108094 530670
121 0 1 2489 40 67804 238125
122 0 1 0 0 60720 315380
123 195 0 6984 45 80570 741409
124 19 1 2001 52 95551 393343
125 -158 1 15236 90 77440 372631
126 0 0 0 0 60720 315380
127 0 1 530 60 73433 317291
128 0 1 0 0 60720 315380
129 -70 1 35624 223 157278 306275
130 13 0 1383 11 73221 317892
131 3 1 875 14 67000 334280
132 0 0 0 0 60720 315380
133 0 0 0 0 60720 315380
134 0 0 0 0 60720 315380
135 0 1 72 3 60398 314210
136 1 0 265 31 60720 306948
137 1 0 0 60720 315380 0
138 0 335 33 64175 320398 0
139 1 0 0 60720 315380 0
140 1 0 0 60720 315380 0
141 0 4525 226 93811 501749 101
142 1 3045 58 27330 202055 31
143 1 0 0 60720 315380 0
144 0 0 0 60720 315380 0
145 0 638 14 60370 333210 0
146 1 0 0 60720 315380 0
147 1 607 11 60436 322340 0
148 0 1558 63 55637 369448 0
149 1 1324 117 67440 291841 0
150 0 0 0 60720 315380 0
151 1 0 0 60720 315380 0
152 1 611 8 59190 296919 -4
153 1 0 0 60720 315380 0
154 1 923 8 58620 309038 0
155 0 661 3 65920 246541 -12
156 1 2397 14 74020 289513 6
157 1 366 53 61808 344425 0
158 1 0 0 60720 315380 0
159 0 135 2 62065 314210 0
160 1 1659 29 107577 480382 2
161 1 316 9 60505 315009 0
162 1 0 0 60720 315380 0
163 0 309 11 56535 312878 0
164 1 49 8 64107 322031 0
165 0 0 0 60720 315380 0
166 1 4519 49 102129 597793 163
167 1 0 0 60720 315380 0
168 1 837 69 61262 315688 4
169 1 5119 49 39039 378525 31
170 1 1280 117 69465 312378 26
171 0 2564 22 130140 403560 13
172 0 2045 20 97890 510834 10
173 1 2234 16 77200 214215 5
174 0 975 23 90534 235133 13
175 1 1136 32 48522 343613 0
176 1 453 21 81125 365959 0
177 1 0 0 60720 315380 0
178 1 61 14 60720 314551 0
179 0 368 17 61686 303230 0
180 1 0 0 60720 315380 0
181 1 4901 84 121920 469107 67
182 1 540 14 103960 354228 22
183 1 0 0 60720 315380 0
184 0 0 0 60720 315380 0
185 1 0 0 60720 315380 0
186 1 2 9 60735 315394 0
187 1 36 4 61564 312412 0
188 1 776 9 64230 333505 -4
189 1 84738 588 -26007 223193 4145
190 1 0 0 60720 315380 0
191 1 3 4 60761 315656 0
192 1 529 9 63870 296261 -3
193 1 0 0 60720 315380 0
194 1 405 7 60845 336425 0
195 1 972 34 71642 359335 0
196 1 0 0 60720 315380 0
197 0 2099 64 106611 308636 20
198 0 3437 47 48022 158492 9
199 1 0 0 60720 315380 0
200 0 0 0 60720 315380 0
201 0 22330 84 79801 711969 446
202 0 0 0 60720 315380 0
203 1 0 0 60720 315380 0
204 1 483 27 60830 306268 0
205 1 0 0 60720 315380 0
206 0 2239 21 88590 442882 5
207 1 2949 41 82903 378509 18
208 1 0 0 60720 315380 0
209 1 365 83 87192 346611 0
210 0 2461 57 55792 314289 0
211 0 21950 519 114337 856956 745
212 1 3294 23 75832 217193 1
213 1 141 6 61630 315366 0
214 1 572 16 58580 307930 -1
215 1 13326 102 165548 702380 436
216 0 2284 33 76403 194493 66
217 1 10 2 61656 316155 0
218 0 0 0 60720 315380 0
219 1 1414 198 70184 330546 13
220 1 1975 35 118881 394510 30
221 1 43 3 60887 312846 0
222 1 0 0 60720 315380 0
223 1 844 78 60925 296139 0
224 1 304 15 62969 295580 0
225 1 458 11 58625 297765 0
226 1 18562 155 102313 377934 143
227 1 0 0 60720 315380 0
228 1 7123 109 288170 638830 769
229 1 622 73 73007 304376 0
230 1 174 1 64820 307424 0
231 1 2220 22 301670 644190 154
232 0 121 24 56178 295370 0
233 1 11819 85 106113 574339 199
234 1 0 0 60720 315380 0
235 1 125 12 60798 310201 0
236 1 1182 13 70694 327007 -5
237 1 1503 39 56364 343466 6
238 1 0 0 60720 315380 0
239 1 0 0 60720 315380 0
240 0 30 4 62045 318098 0
241 1 3310 33 75230 448243 -6
242 0 554 11 79285 325738 30
243 1 0 0 60720 315380 0
244 1 468 21 60720 312161 0
245 1 4917 42 52811 243650 -3
246 1 3256 66 35250 407159 -8
247 1 0 0 60720 315380 0
248 1 125 13 59734 317698 0
249 1 0 0 60720 315380 0
250 0 22 2 60722 312502 0
251 1 0 0 60720 315380 0
252 0 514 52 78780 322378 0
253 1 0 0 60720 315380 0
254 0 0 0 60720 315380 0
255 1 0 0 60720 315380 0
256 1 0 0 60720 315380 0
257 1 0 0 60720 315380 0
258 1 0 0 60720 315380 0
259 1 10327 91 88577 640273 135
260 1 13718 61 96448 345783 104
261 1 3748 9 50350 652925 53
262 0 14416 361 49857 439798 220
263 0 1526 25 69351 278990 -4
264 1 666 171 117869 339836 25
265 1 2844 67 72683 240897 -15
266 1 0 0 60720 315380 0
267 1 368 17 61167 297141 0
268 1 0 0 60720 315380 0
269 1 333 13 70811 331323 0
270 1 26 1 60896 313880 0
271 1 0 0 60720 315380 0
272 1 0 0 60720 315380 0
273 1 1303 61 69863 309422 -5
274 1 20 3 60938 315245 0
275 0 2384 97 61348 405972 54
276 0 203 27 50804 300962 0
277 1 71 1 60745 316176 0
278 1 53 14 59506 302409 0
279 1 562 13 58456 283587 -3
280 1 622 7 60950 263276 -2
281 0 645 5 60720 312075 0
282 1 1763 24 61600 308336 0
283 1 0 0 60720 315380 0
284 1 317 9 63915 298700 0
285 1 1 6 60719 315372 0
286 1 275 4 59500 318745 0
287 1 0 0 60720 315380 0
288 1 936 30 67939 408881 9
289 0 8568 37 32168 786690 223
290 0 11528 70 -14545 -83265 -25
291 1 0 0 60720 315380 0
292 1 738 10 60720 321376 0
293 1 0 0 60720 315380 0
294 1 592 7 64270 276898 -4
295 1 126 10 60951 315547 0
296 1 2 1 60743 315487 0
297 0 0 0 60720 315380 0
298 0 18014 239 -1710 1405225 943
299 0 0 0 60720 315380 0
300 0 0 0 60720 315380 0
301 0 37704 452 60448 983660 2167
302 1 63 4 65688 318574 0
303 1 1431 34 106885 310768 30
304 1 94 3 61360 312887 0
305 1 192 4 65276 312339 0
306 1 32 1 59988 314964 0
307 1 6869 105 117520 379983 90
308 1 0 0 60720 315380 0
309 1 0 0 60720 315380 0
310 0 1 4 60722 315398 0
311 1 0 0 60720 315380 0
312 0 2328 94 82732 253588 0
313 1 209 13 64016 316647 0
314 1 28 2 60890 315688 0
315 1 176 28 68136 310670 0
316 0 1920 15 79420 165404 -6
317 1 87550 458 153198 4111912 17936
318 0 520 10 58650 291650 0
319 1 0 0 60720 315380 0
320 1 1013 2 69770 253468 -14
321 1 15 8 60831 315688 0
322 1 587 12 59595 325699 0
323 1 5371 66 87720 446211 49
324 1 0 0 60720 315380 0
325 1 1012 47 114768 368078 26
326 1 0 0 60720 315380 0
327 0 876 59 138971 352850 24
328 1 162556 1081 213118 6282154 56289
329 0 43556 247 32648 227132 929
330 1 3425 12 83620 283910 -3
331 1 810 43 74015 236761 26
332 0 0 0 60720 315380 0
333 0 0 0 60720 315380 0
334 0 2365 79 191778 550608 120
335 1 1261 33 76114 307528 13
336 1 0 0 60720 315380 0
337 1 1585 107 93099 355864 2
338 1 16189 295 116384 358589 379
339 0 0 0 60720 315380 0
340 1 0 0 60720 315380 0
341 1 10579 92 110309 375195 4
342 1 474 20 61977 288985 0
343 1 0 0 60720 315380 0
344 1 3642 86 90262 458343 22
345 0 0 0 60720 315380 0
346 1 472 20 80045 269587 0
347 1 98 2 61490 315236 0
348 1 3999 30 51252 42754 -58
349 1 0 0 60720 315380 0
350 1 621 48 60720 308256 0
351 1 0 0 60720 315380 0
352 0 0 0 60720 315380 0
353 1 30 4 60798 313164 0
354 1 0 0 60720 315380 0
355 0 746 23 71561 269661 10
356 1 0 0 60720 315380 0
357 0 0 0 60720 315380 0
358 1 0 0 60720 315380 0
359 0 0 0 60720 315380 0
360 0 4150 150 134759 518365 66
361 1 4658 158 156608 233773 58
362 1 814 18 39625 301881 23
363 0 1002 15 56750 298568 0
364 1 496 17 87390 325479 0
365 1 389 14 58990 325506 0
366 1 12679 110 48020 984885 1871
367 1 400 152 60720 313267 0
368 1 53 6 60349 315793 0
369 0 0 0 60720 315380 0
370 1 5109 97 67038 215362 -31
371 1 576 18 113761 314073 12
372 1 438 11 58320 298096 -5
373 1 165 24 62841 325176 0
374 0 4069 145 79804 207393 119
375 1 23 5 62555 314806 0
376 0 1285 107 99489 341340 49
377 0 4677 75 90131 426280 24
378 1 24811 505 95350 929118 822
379 1 167 2 64245 307322 0
380 1 0 0 60720 315380 0
381 1 4628 62 69159 387475 85
382 0 226 5 65745 291787 0
383 1 1765 20 77623 247060 28
384 1 460 15 63346 329784 0
385 1 36 3 60894 315834 0
386 0 989 11 58930 304555 -5
387 0 3055 43 60247 376641 14
388 1 0 0 60720 315380 0
389 0 0 0 60720 315380 0
390 1 0 0 60720 315380 0
391 1 0 0 60720 315380 0
392 1 0 0 60720 315380 0
393 1 13253 947 90829 357760 35
394 1 24 8 59818 315637 0
395 0 0 0 60720 315380 0
396 1 0 0 60720 315380 0
397 1 0 0 60720 315380 0
398 1 0 0 60720 315380 0
399 1 131 6 59661 327071 0
400 1 514 15 102725 377516 23
401 0 3366 80 108479 299243 12
402 1 9327 101 87419 387699 44
403 1 384 22 54683 309836 0
404 1 17821 373 122844 444477 643
405 1 397 10 62710 322327 -1
406 1 897 5 60720 314913 0
407 0 218 21 60720 315553 0
408 1 3369 28 87161 688779 217
409 1 5702 61 101481 321896 10
410 1 8636 855 128294 301607 108
411 1 534 7 62620 304485 0
412 0 0 0 60720 315380 0
413 1 0 0 60720 315380 0
414 1 726 3 69980 231861 -5
415 1 1380 26 60982 347385 25
416 0 180 11 59635 316386 0
417 1 7285 115 188873 491303 59
418 1 880 74 80791 261216 34
419 1 1 2 60727 315388 0
420 1 0 0 60720 315380 0
421 1 96 5 60379 313729 0
422 0 1889 45 37527 358649 14
423 1 45187 353 234817 1926517 2854
424 1 288 4 60510 296656 0
425 1 1270 26 69206 275311 -10
426 0 6526 26 55830 -42143 -28
427 1 0 0 60720 315380 0
428 1 226 12 72835 343929 0
429 1 694 8 68060 367655 -2
430 1 0 0 60720 315380 0
431 1 249 72 56726 313491 0
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) Group Costs Trades Dividends Wealth
-1.185e+03 -1.050e-02 -6.235e-02 1.372e-02 1.202e-04 5.227e-03
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-3378.6 -428.6 247.2 315.4 8771.3
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -1.185e+03 7.036e+01 -16.846 < 2e-16 ***
Group -1.050e-02 6.360e-03 -1.651 0.0995 .
Costs -6.235e-02 8.697e-03 -7.170 3.35e-12 ***
Trades 1.372e-02 1.131e-03 12.128 < 2e-16 ***
Dividends 1.202e-04 1.966e-04 0.611 0.5413
Wealth 5.227e-03 1.755e-04 29.790 < 2e-16 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 729.1 on 425 degrees of freedom
Multiple R-squared: 0.729, Adjusted R-squared: 0.7258
F-statistic: 228.6 on 5 and 425 DF, p-value: < 2.2e-16
> if (n > n25) {
+ kp3 <- k + 3
+ nmkm3 <- n - k - 3
+ gqarr <- array(NA, dim=c(nmkm3-kp3+1,3))
+ numgqtests <- 0
+ numsignificant1 <- 0
+ numsignificant5 <- 0
+ numsignificant10 <- 0
+ for (mypoint in kp3:nmkm3) {
+ j <- 0
+ numgqtests <- numgqtests + 1
+ for (myalt in c('greater', 'two.sided', 'less')) {
+ j <- j + 1
+ gqarr[mypoint-kp3+1,j] <- gqtest(mylm, point=mypoint, alternative=myalt)$p.value
+ }
+ if (gqarr[mypoint-kp3+1,2] < 0.01) numsignificant1 <- numsignificant1 + 1
+ if (gqarr[mypoint-kp3+1,2] < 0.05) numsignificant5 <- numsignificant5 + 1
+ if (gqarr[mypoint-kp3+1,2] < 0.10) numsignificant10 <- numsignificant10 + 1
+ }
+ gqarr
+ }
[,1] [,2] [,3]
[1,] 1.677346e-05 3.354693e-05 9.999832e-01
[2,] 3.451887e-07 6.903774e-07 9.999997e-01
[3,] 1.078783e-08 2.157566e-08 1.000000e+00
[4,] 2.509298e-10 5.018595e-10 1.000000e+00
[5,] 4.963145e-12 9.926291e-12 1.000000e+00
[6,] 9.039728e-13 1.807946e-12 1.000000e+00
[7,] 1.938604e-14 3.877208e-14 1.000000e+00
[8,] 7.691098e-16 1.538220e-15 1.000000e+00
[9,] 2.015417e-17 4.030833e-17 1.000000e+00
[10,] 5.345730e-19 1.069146e-18 1.000000e+00
[11,] 1.163285e-20 2.326571e-20 1.000000e+00
[12,] 2.300441e-22 4.600882e-22 1.000000e+00
[13,] 4.488800e-24 8.977599e-24 1.000000e+00
[14,] 3.298070e-25 6.596140e-25 1.000000e+00
[15,] 6.454870e-27 1.290974e-26 1.000000e+00
[16,] 4.879198e-28 9.758395e-28 1.000000e+00
[17,] 1.189891e-27 2.379781e-27 1.000000e+00
[18,] 3.035039e-29 6.070078e-29 1.000000e+00
[19,] 8.414731e-30 1.682946e-29 1.000000e+00
[20,] 4.047575e-30 8.095150e-30 1.000000e+00
[21,] 1.212880e-31 2.425760e-31 1.000000e+00
[22,] 6.472157e-33 1.294431e-32 1.000000e+00
[23,] 6.714550e-31 1.342910e-30 1.000000e+00
[24,] 2.696599e-32 5.393198e-32 1.000000e+00
[25,] 1.118978e-33 2.237957e-33 1.000000e+00
[26,] 5.807141e-35 1.161428e-34 1.000000e+00
[27,] 2.688470e-36 5.376939e-36 1.000000e+00
[28,] 9.209374e-36 1.841875e-35 1.000000e+00
[29,] 4.032605e-37 8.065210e-37 1.000000e+00
[30,] 5.127136e-38 1.025427e-37 1.000000e+00
[31,] 9.240865e-28 1.848173e-27 1.000000e+00
[32,] 1.665227e-28 3.330454e-28 1.000000e+00
[33,] 5.400613e-29 1.080123e-28 1.000000e+00
[34,] 1.689924e-29 3.379849e-29 1.000000e+00
[35,] 1.781298e-30 3.562596e-30 1.000000e+00
[36,] 1.869708e-31 3.739416e-31 1.000000e+00
[37,] 7.357510e-31 1.471502e-30 1.000000e+00
[38,] 7.924913e-32 1.584983e-31 1.000000e+00
[39,] 8.775655e-33 1.755131e-32 1.000000e+00
[40,] 9.692653e-34 1.938531e-33 1.000000e+00
[41,] 1.234657e-34 2.469314e-34 1.000000e+00
[42,] 1.602581e-35 3.205162e-35 1.000000e+00
[43,] 2.748169e-25 5.496338e-25 1.000000e+00
[44,] 7.024019e-26 1.404804e-25 1.000000e+00
[45,] 6.052559e-26 1.210512e-25 1.000000e+00
[46,] 6.074398e-25 1.214880e-24 1.000000e+00
[47,] 5.207021e-23 1.041404e-22 1.000000e+00
[48,] 1.457187e-20 2.914374e-20 1.000000e+00
[49,] 3.944977e-21 7.889954e-21 1.000000e+00
[50,] 2.362866e-21 4.725732e-21 1.000000e+00
[51,] 6.645107e-22 1.329021e-21 1.000000e+00
[52,] 1.802758e-22 3.605516e-22 1.000000e+00
[53,] 4.863480e-23 9.726961e-23 1.000000e+00
[54,] 1.258580e-23 2.517160e-23 1.000000e+00
[55,] 6.884702e-24 1.376940e-23 1.000000e+00
[56,] 1.846280e-24 3.692560e-24 1.000000e+00
[57,] 2.616074e-24 5.232149e-24 1.000000e+00
[58,] 7.130378e-25 1.426076e-24 1.000000e+00
[59,] 1.928661e-25 3.857323e-25 1.000000e+00
[60,] 6.478062e-26 1.295612e-25 1.000000e+00
[61,] 1.721957e-26 3.443913e-26 1.000000e+00
[62,] 4.646390e-27 9.292781e-27 1.000000e+00
[63,] 1.300471e-25 2.600943e-25 1.000000e+00
[64,] 3.647360e-26 7.294719e-26 1.000000e+00
[65,] 1.000000e+00 2.762772e-146 1.381386e-146
[66,] 1.000000e+00 1.779885e-145 8.899427e-146
[67,] 1.000000e+00 5.727666e-145 2.863833e-145
[68,] 1.000000e+00 3.144854e-144 1.572427e-144
[69,] 1.000000e+00 1.257698e-143 6.288489e-144
[70,] 1.000000e+00 7.241033e-143 3.620516e-143
[71,] 1.000000e+00 4.231106e-146 2.115553e-146
[72,] 1.000000e+00 8.348916e-153 4.174458e-153
[73,] 1.000000e+00 5.101400e-152 2.550700e-152
[74,] 1.000000e+00 3.070059e-151 1.535030e-151
[75,] 1.000000e+00 1.821170e-150 9.105848e-151
[76,] 1.000000e+00 5.135912e-150 2.567956e-150
[77,] 1.000000e+00 2.943525e-149 1.471763e-149
[78,] 1.000000e+00 1.611638e-148 8.058191e-149
[79,] 1.000000e+00 8.874433e-148 4.437216e-148
[80,] 1.000000e+00 4.846123e-147 2.423062e-147
[81,] 1.000000e+00 2.551639e-146 1.275820e-146
[82,] 1.000000e+00 1.182848e-145 5.914241e-146
[83,] 1.000000e+00 5.960437e-145 2.980218e-145
[84,] 1.000000e+00 3.160320e-144 1.580160e-144
[85,] 1.000000e+00 1.810311e-153 9.051554e-154
[86,] 1.000000e+00 1.632360e-154 8.161801e-155
[87,] 1.000000e+00 8.211026e-154 4.105513e-154
[88,] 1.000000e+00 4.007385e-153 2.003693e-153
[89,] 1.000000e+00 2.120810e-152 1.060405e-152
[90,] 1.000000e+00 1.062634e-151 5.313168e-152
[91,] 1.000000e+00 5.031077e-151 2.515538e-151
[92,] 1.000000e+00 2.727164e-150 1.363582e-150
[93,] 1.000000e+00 5.755644e-151 2.877822e-151
[94,] 1.000000e+00 2.726361e-150 1.363180e-150
[95,] 1.000000e+00 1.203205e-149 6.016024e-150
[96,] 1.000000e+00 5.254052e-149 2.627026e-149
[97,] 1.000000e+00 2.231271e-148 1.115635e-148
[98,] 1.000000e+00 5.965881e-163 2.982940e-163
[99,] 1.000000e+00 1.995278e-162 9.976392e-163
[100,] 1.000000e+00 0.000000e+00 0.000000e+00
[101,] 1.000000e+00 0.000000e+00 0.000000e+00
[102,] 1.000000e+00 9.881313e-324 4.940656e-324
[103,] 1.000000e+00 1.284571e-322 6.422853e-323
[104,] 1.000000e+00 2.667954e-322 1.333977e-322
[105,] 1.000000e+00 1.333977e-321 6.669886e-322
[106,] 1.000000e+00 4.555285e-321 2.277643e-321
[107,] 1.000000e+00 1.747016e-320 8.735081e-321
[108,] 1.000000e+00 6.540441e-320 3.270221e-320
[109,] 1.000000e+00 5.744995e-319 2.872498e-319
[110,] 1.000000e+00 0.000000e+00 0.000000e+00
[111,] 1.000000e+00 0.000000e+00 0.000000e+00
[112,] 1.000000e+00 0.000000e+00 0.000000e+00
[113,] 1.000000e+00 0.000000e+00 0.000000e+00
[114,] 1.000000e+00 0.000000e+00 0.000000e+00
[115,] 1.000000e+00 0.000000e+00 0.000000e+00
[116,] 1.000000e+00 0.000000e+00 0.000000e+00
[117,] 1.000000e+00 0.000000e+00 0.000000e+00
[118,] 1.000000e+00 0.000000e+00 0.000000e+00
[119,] 1.000000e+00 0.000000e+00 0.000000e+00
[120,] 1.000000e+00 0.000000e+00 0.000000e+00
[121,] 1.000000e+00 0.000000e+00 0.000000e+00
[122,] 1.000000e+00 0.000000e+00 0.000000e+00
[123,] 1.000000e+00 0.000000e+00 0.000000e+00
[124,] 1.000000e+00 0.000000e+00 0.000000e+00
[125,] 1.000000e+00 0.000000e+00 0.000000e+00
[126,] 1.000000e+00 0.000000e+00 0.000000e+00
[127,] 1.000000e+00 0.000000e+00 0.000000e+00
[128,] 1.000000e+00 0.000000e+00 0.000000e+00
[129,] 1.000000e+00 0.000000e+00 0.000000e+00
[130,] 1.000000e+00 0.000000e+00 0.000000e+00
[131,] 1.000000e+00 0.000000e+00 0.000000e+00
[132,] 1.000000e+00 0.000000e+00 0.000000e+00
[133,] 1.000000e+00 0.000000e+00 0.000000e+00
[134,] 1.000000e+00 0.000000e+00 0.000000e+00
[135,] 1.000000e+00 0.000000e+00 0.000000e+00
[136,] 1.000000e+00 0.000000e+00 0.000000e+00
[137,] 1.000000e+00 0.000000e+00 0.000000e+00
[138,] 1.000000e+00 0.000000e+00 0.000000e+00
[139,] 1.000000e+00 0.000000e+00 0.000000e+00
[140,] 1.000000e+00 0.000000e+00 0.000000e+00
[141,] 1.000000e+00 0.000000e+00 0.000000e+00
[142,] 1.000000e+00 0.000000e+00 0.000000e+00
[143,] 1.000000e+00 0.000000e+00 0.000000e+00
[144,] 1.000000e+00 0.000000e+00 0.000000e+00
[145,] 1.000000e+00 0.000000e+00 0.000000e+00
[146,] 1.000000e+00 0.000000e+00 0.000000e+00
[147,] 1.000000e+00 0.000000e+00 0.000000e+00
[148,] 1.000000e+00 0.000000e+00 0.000000e+00
[149,] 1.000000e+00 0.000000e+00 0.000000e+00
[150,] 1.000000e+00 0.000000e+00 0.000000e+00
[151,] 1.000000e+00 0.000000e+00 0.000000e+00
[152,] 1.000000e+00 0.000000e+00 0.000000e+00
[153,] 1.000000e+00 0.000000e+00 0.000000e+00
[154,] 1.000000e+00 0.000000e+00 0.000000e+00
[155,] 1.000000e+00 0.000000e+00 0.000000e+00
[156,] 1.000000e+00 0.000000e+00 0.000000e+00
[157,] 1.000000e+00 0.000000e+00 0.000000e+00
[158,] 1.000000e+00 0.000000e+00 0.000000e+00
[159,] 1.000000e+00 0.000000e+00 0.000000e+00
[160,] 1.000000e+00 0.000000e+00 0.000000e+00
[161,] 1.000000e+00 0.000000e+00 0.000000e+00
[162,] 1.000000e+00 0.000000e+00 0.000000e+00
[163,] 1.000000e+00 0.000000e+00 0.000000e+00
[164,] 1.000000e+00 0.000000e+00 0.000000e+00
[165,] 1.000000e+00 0.000000e+00 0.000000e+00
[166,] 1.000000e+00 0.000000e+00 0.000000e+00
[167,] 1.000000e+00 0.000000e+00 0.000000e+00
[168,] 1.000000e+00 0.000000e+00 0.000000e+00
[169,] 1.000000e+00 0.000000e+00 0.000000e+00
[170,] 1.000000e+00 0.000000e+00 0.000000e+00
[171,] 1.000000e+00 0.000000e+00 0.000000e+00
[172,] 1.000000e+00 0.000000e+00 0.000000e+00
[173,] 1.000000e+00 0.000000e+00 0.000000e+00
[174,] 1.000000e+00 0.000000e+00 0.000000e+00
[175,] 1.000000e+00 0.000000e+00 0.000000e+00
[176,] 1.000000e+00 0.000000e+00 0.000000e+00
[177,] 1.000000e+00 0.000000e+00 0.000000e+00
[178,] 1.000000e+00 0.000000e+00 0.000000e+00
[179,] 1.000000e+00 0.000000e+00 0.000000e+00
[180,] 1.000000e+00 0.000000e+00 0.000000e+00
[181,] 1.000000e+00 0.000000e+00 0.000000e+00
[182,] 1.000000e+00 0.000000e+00 0.000000e+00
[183,] 1.000000e+00 0.000000e+00 0.000000e+00
[184,] 1.000000e+00 0.000000e+00 0.000000e+00
[185,] 1.000000e+00 0.000000e+00 0.000000e+00
[186,] 1.000000e+00 0.000000e+00 0.000000e+00
[187,] 1.000000e+00 0.000000e+00 0.000000e+00
[188,] 1.000000e+00 0.000000e+00 0.000000e+00
[189,] 1.000000e+00 0.000000e+00 0.000000e+00
[190,] 1.000000e+00 0.000000e+00 0.000000e+00
[191,] 1.000000e+00 0.000000e+00 0.000000e+00
[192,] 1.000000e+00 0.000000e+00 0.000000e+00
[193,] 1.000000e+00 0.000000e+00 0.000000e+00
[194,] 1.000000e+00 0.000000e+00 0.000000e+00
[195,] 1.000000e+00 0.000000e+00 0.000000e+00
[196,] 1.000000e+00 0.000000e+00 0.000000e+00
[197,] 1.000000e+00 0.000000e+00 0.000000e+00
[198,] 1.000000e+00 0.000000e+00 0.000000e+00
[199,] 1.000000e+00 0.000000e+00 0.000000e+00
[200,] 1.000000e+00 0.000000e+00 0.000000e+00
[201,] 1.000000e+00 0.000000e+00 0.000000e+00
[202,] 1.000000e+00 0.000000e+00 0.000000e+00
[203,] 1.000000e+00 0.000000e+00 0.000000e+00
[204,] 1.000000e+00 0.000000e+00 0.000000e+00
[205,] 1.000000e+00 0.000000e+00 0.000000e+00
[206,] 1.000000e+00 0.000000e+00 0.000000e+00
[207,] 1.000000e+00 0.000000e+00 0.000000e+00
[208,] 1.000000e+00 0.000000e+00 0.000000e+00
[209,] 1.000000e+00 0.000000e+00 0.000000e+00
[210,] 1.000000e+00 0.000000e+00 0.000000e+00
[211,] 1.000000e+00 0.000000e+00 0.000000e+00
[212,] 1.000000e+00 0.000000e+00 0.000000e+00
[213,] 1.000000e+00 0.000000e+00 0.000000e+00
[214,] 1.000000e+00 0.000000e+00 0.000000e+00
[215,] 1.000000e+00 0.000000e+00 0.000000e+00
[216,] 1.000000e+00 0.000000e+00 0.000000e+00
[217,] 1.000000e+00 0.000000e+00 0.000000e+00
[218,] 1.000000e+00 0.000000e+00 0.000000e+00
[219,] 1.000000e+00 0.000000e+00 0.000000e+00
[220,] 1.000000e+00 0.000000e+00 0.000000e+00
[221,] 1.000000e+00 0.000000e+00 0.000000e+00
[222,] 1.000000e+00 0.000000e+00 0.000000e+00
[223,] 1.000000e+00 0.000000e+00 0.000000e+00
[224,] 1.000000e+00 0.000000e+00 0.000000e+00
[225,] 1.000000e+00 0.000000e+00 0.000000e+00
[226,] 1.000000e+00 0.000000e+00 0.000000e+00
[227,] 1.000000e+00 0.000000e+00 0.000000e+00
[228,] 1.000000e+00 0.000000e+00 0.000000e+00
[229,] 1.000000e+00 0.000000e+00 0.000000e+00
[230,] 1.000000e+00 0.000000e+00 0.000000e+00
[231,] 1.000000e+00 0.000000e+00 0.000000e+00
[232,] 1.000000e+00 0.000000e+00 0.000000e+00
[233,] 1.000000e+00 0.000000e+00 0.000000e+00
[234,] 1.000000e+00 0.000000e+00 0.000000e+00
[235,] 1.000000e+00 0.000000e+00 0.000000e+00
[236,] 1.000000e+00 0.000000e+00 0.000000e+00
[237,] 1.000000e+00 0.000000e+00 0.000000e+00
[238,] 1.000000e+00 0.000000e+00 0.000000e+00
[239,] 1.000000e+00 0.000000e+00 0.000000e+00
[240,] 1.000000e+00 0.000000e+00 0.000000e+00
[241,] 1.000000e+00 0.000000e+00 0.000000e+00
[242,] 1.000000e+00 0.000000e+00 0.000000e+00
[243,] 1.000000e+00 0.000000e+00 0.000000e+00
[244,] 1.000000e+00 0.000000e+00 0.000000e+00
[245,] 1.000000e+00 0.000000e+00 0.000000e+00
[246,] 1.000000e+00 0.000000e+00 0.000000e+00
[247,] 1.000000e+00 0.000000e+00 0.000000e+00
[248,] 1.000000e+00 0.000000e+00 0.000000e+00
[249,] 1.000000e+00 0.000000e+00 0.000000e+00
[250,] 1.000000e+00 0.000000e+00 0.000000e+00
[251,] 1.000000e+00 0.000000e+00 0.000000e+00
[252,] 1.000000e+00 0.000000e+00 0.000000e+00
[253,] 1.000000e+00 0.000000e+00 0.000000e+00
[254,] 1.000000e+00 0.000000e+00 0.000000e+00
[255,] 1.000000e+00 0.000000e+00 0.000000e+00
[256,] 1.000000e+00 0.000000e+00 0.000000e+00
[257,] 1.000000e+00 0.000000e+00 0.000000e+00
[258,] 1.000000e+00 0.000000e+00 0.000000e+00
[259,] 1.000000e+00 0.000000e+00 0.000000e+00
[260,] 1.000000e+00 0.000000e+00 0.000000e+00
[261,] 1.000000e+00 0.000000e+00 0.000000e+00
[262,] 1.000000e+00 0.000000e+00 0.000000e+00
[263,] 1.000000e+00 0.000000e+00 0.000000e+00
[264,] 1.000000e+00 0.000000e+00 0.000000e+00
[265,] 1.000000e+00 0.000000e+00 0.000000e+00
[266,] 1.000000e+00 0.000000e+00 0.000000e+00
[267,] 1.000000e+00 0.000000e+00 0.000000e+00
[268,] 1.000000e+00 0.000000e+00 0.000000e+00
[269,] 1.000000e+00 0.000000e+00 0.000000e+00
[270,] 1.000000e+00 0.000000e+00 0.000000e+00
[271,] 1.000000e+00 0.000000e+00 0.000000e+00
[272,] 1.000000e+00 0.000000e+00 0.000000e+00
[273,] 1.000000e+00 0.000000e+00 0.000000e+00
[274,] 1.000000e+00 0.000000e+00 0.000000e+00
[275,] 1.000000e+00 0.000000e+00 0.000000e+00
[276,] 1.000000e+00 0.000000e+00 0.000000e+00
[277,] 1.000000e+00 0.000000e+00 0.000000e+00
[278,] 1.000000e+00 0.000000e+00 0.000000e+00
[279,] 1.000000e+00 0.000000e+00 0.000000e+00
[280,] 1.000000e+00 0.000000e+00 0.000000e+00
[281,] 1.000000e+00 0.000000e+00 0.000000e+00
[282,] 1.000000e+00 0.000000e+00 0.000000e+00
[283,] 1.000000e+00 0.000000e+00 0.000000e+00
[284,] 1.000000e+00 0.000000e+00 0.000000e+00
[285,] 1.000000e+00 0.000000e+00 0.000000e+00
[286,] 1.000000e+00 0.000000e+00 0.000000e+00
[287,] 1.000000e+00 0.000000e+00 0.000000e+00
[288,] 1.000000e+00 0.000000e+00 0.000000e+00
[289,] 1.000000e+00 0.000000e+00 0.000000e+00
[290,] 1.000000e+00 0.000000e+00 0.000000e+00
[291,] 1.000000e+00 0.000000e+00 0.000000e+00
[292,] 1.000000e+00 0.000000e+00 0.000000e+00
[293,] 1.000000e+00 0.000000e+00 0.000000e+00
[294,] 1.000000e+00 0.000000e+00 0.000000e+00
[295,] 1.000000e+00 0.000000e+00 0.000000e+00
[296,] 1.000000e+00 0.000000e+00 0.000000e+00
[297,] 1.000000e+00 0.000000e+00 0.000000e+00
[298,] 1.000000e+00 0.000000e+00 0.000000e+00
[299,] 1.000000e+00 0.000000e+00 0.000000e+00
[300,] 1.000000e+00 0.000000e+00 0.000000e+00
[301,] 1.000000e+00 0.000000e+00 0.000000e+00
[302,] 1.000000e+00 0.000000e+00 0.000000e+00
[303,] 1.000000e+00 0.000000e+00 0.000000e+00
[304,] 1.000000e+00 0.000000e+00 0.000000e+00
[305,] 1.000000e+00 0.000000e+00 0.000000e+00
[306,] 1.000000e+00 0.000000e+00 0.000000e+00
[307,] 1.000000e+00 0.000000e+00 0.000000e+00
[308,] 1.000000e+00 0.000000e+00 0.000000e+00
[309,] 1.000000e+00 0.000000e+00 0.000000e+00
[310,] 1.000000e+00 9.881313e-324 4.940656e-324
[311,] 1.000000e+00 1.077063e-320 5.385316e-321
[312,] 1.000000e+00 1.531626e-317 7.658131e-318
[313,] 1.000000e+00 2.291963e-314 1.145981e-314
[314,] 1.000000e+00 3.403938e-311 1.701969e-311
[315,] 1.000000e+00 5.208096e-308 2.604048e-308
[316,] 1.000000e+00 7.605405e-305 3.802703e-305
[317,] 1.000000e+00 1.036104e-301 5.180519e-302
[318,] 1.000000e+00 1.494995e-298 7.474976e-299
[319,] 1.000000e+00 6.476543e-296 3.238271e-296
[320,] 1.000000e+00 4.149540e-293 2.074770e-293
[321,] 1.000000e+00 2.511571e-290 1.255786e-290
[322,] 1.000000e+00 3.692146e-287 1.846073e-287
[323,] 1.000000e+00 5.378287e-284 2.689143e-284
[324,] 1.000000e+00 2.445436e-281 1.222718e-281
[325,] 1.000000e+00 1.054342e-278 5.271711e-279
[326,] 1.000000e+00 2.622121e-276 1.311060e-276
[327,] 1.000000e+00 4.068665e-273 2.034333e-273
[328,] 1.000000e+00 6.048531e-270 3.024265e-270
[329,] 1.000000e+00 9.656588e-267 4.828294e-267
[330,] 1.000000e+00 1.458358e-263 7.291788e-264
[331,] 1.000000e+00 6.126811e-261 3.063406e-261
[332,] 1.000000e+00 9.155059e-258 4.577530e-258
[333,] 1.000000e+00 1.261405e-254 6.307025e-255
[334,] 1.000000e+00 1.854940e-251 9.274701e-252
[335,] 1.000000e+00 2.721525e-248 1.360763e-248
[336,] 1.000000e+00 4.233826e-245 2.116913e-245
[337,] 1.000000e+00 1.816616e-242 9.083079e-243
[338,] 1.000000e+00 2.710570e-239 1.355285e-239
[339,] 1.000000e+00 3.972915e-236 1.986457e-236
[340,] 1.000000e+00 4.288895e-233 2.144448e-233
[341,] 1.000000e+00 6.207579e-230 3.103789e-230
[342,] 1.000000e+00 8.989598e-227 4.494799e-227
[343,] 1.000000e+00 1.278874e-223 6.394368e-224
[344,] 1.000000e+00 5.768274e-221 2.884137e-221
[345,] 1.000000e+00 8.321561e-218 4.160781e-218
[346,] 1.000000e+00 1.188932e-214 5.944661e-215
[347,] 1.000000e+00 5.773763e-212 2.886881e-212
[348,] 1.000000e+00 8.328352e-209 4.164176e-209
[349,] 1.000000e+00 3.568128e-206 1.784064e-206
[350,] 1.000000e+00 5.227315e-203 2.613657e-203
[351,] 1.000000e+00 2.129162e-200 1.064581e-200
[352,] 1.000000e+00 3.741876e-198 1.870938e-198
[353,] 1.000000e+00 4.974826e-195 2.487413e-195
[354,] 1.000000e+00 6.819293e-192 3.409647e-192
[355,] 1.000000e+00 2.788510e-189 1.394255e-189
[356,] 1.000000e+00 4.557584e-186 2.278792e-186
[357,] 1.000000e+00 7.057526e-183 3.528763e-183
[358,] 1.000000e+00 1.185826e-179 5.929129e-180
[359,] 1.000000e+00 1.911288e-176 9.556440e-177
[360,] 1.000000e+00 2.864863e-173 1.432432e-173
[361,] 1.000000e+00 1.079954e-170 5.399771e-171
[362,] 1.000000e+00 1.293567e-167 6.467834e-168
[363,] 1.000000e+00 2.026949e-164 1.013474e-164
[364,] 1.000000e+00 2.965199e-161 1.482599e-161
[365,] 1.000000e+00 4.538975e-158 2.269488e-158
[366,] 1.000000e+00 2.529751e-155 1.264875e-155
[367,] 1.000000e+00 3.833139e-152 1.916570e-152
[368,] 1.000000e+00 8.938501e-150 4.469250e-150
[369,] 1.000000e+00 2.452597e-147 1.226299e-147
[370,] 1.000000e+00 3.957988e-144 1.978994e-144
[371,] 1.000000e+00 6.325725e-141 3.162863e-141
[372,] 1.000000e+00 9.909567e-138 4.954784e-138
[373,] 1.000000e+00 1.410438e-134 7.052189e-135
[374,] 1.000000e+00 4.460902e-132 2.230451e-132
[375,] 1.000000e+00 7.006714e-129 3.503357e-129
[376,] 1.000000e+00 1.118503e-125 5.592516e-126
[377,] 1.000000e+00 1.747088e-122 8.735442e-123
[378,] 1.000000e+00 6.685272e-120 3.342636e-120
[379,] 1.000000e+00 2.138947e-117 1.069474e-117
[380,] 1.000000e+00 3.571440e-114 1.785720e-114
[381,] 1.000000e+00 8.163125e-112 4.081563e-112
[382,] 1.000000e+00 1.433081e-108 7.165406e-109
[383,] 1.000000e+00 2.474733e-105 1.237367e-105
[384,] 1.000000e+00 4.201036e-102 2.100518e-102
[385,] 1.000000e+00 7.982663e-99 3.991332e-99
[386,] 1.000000e+00 1.316488e-95 6.582439e-96
[387,] 1.000000e+00 3.080963e-93 1.540482e-93
[388,] 1.000000e+00 5.390582e-90 2.695291e-90
[389,] 1.000000e+00 9.248397e-87 4.624199e-87
[390,] 1.000000e+00 1.554471e-83 7.772355e-84
[391,] 1.000000e+00 2.573945e-80 1.286972e-80
[392,] 1.000000e+00 4.630680e-77 2.315340e-77
[393,] 1.000000e+00 7.223479e-75 3.611739e-75
[394,] 1.000000e+00 7.106914e-72 3.553457e-72
[395,] 1.000000e+00 1.148832e-68 5.744161e-69
[396,] 1.000000e+00 1.021952e-65 5.109762e-66
[397,] 1.000000e+00 1.780568e-62 8.902840e-63
[398,] 1.000000e+00 2.697331e-59 1.348666e-59
[399,] 1.000000e+00 4.678282e-57 2.339141e-57
[400,] 1.000000e+00 9.482406e-54 4.741203e-54
[401,] 1.000000e+00 6.771965e-51 3.385983e-51
[402,] 1.000000e+00 1.394498e-47 6.972492e-48
[403,] 1.000000e+00 2.547762e-44 1.273881e-44
[404,] 1.000000e+00 1.309126e-42 6.545628e-43
[405,] 1.000000e+00 3.321918e-39 1.660959e-39
[406,] 1.000000e+00 8.204901e-36 4.102451e-36
[407,] 1.000000e+00 1.161418e-32 5.807089e-33
[408,] 1.000000e+00 9.278476e-33 4.639238e-33
[409,] 1.000000e+00 4.304074e-29 2.152037e-29
[410,] 1.000000e+00 3.303227e-27 1.651613e-27
[411,] 1.000000e+00 3.025130e-23 1.512565e-23
[412,] 1.000000e+00 2.580880e-19 1.290440e-19
[413,] 1.000000e+00 1.989474e-15 9.947370e-16
[414,] 1.000000e+00 3.302574e-12 1.651287e-12
> postscript(file="/var/www/wessaorg/rcomp/tmp/1n3yc1293639134.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/wessaorg/rcomp/tmp/2n3yc1293639134.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/wessaorg/rcomp/tmp/3n3yc1293639134.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/wessaorg/rcomp/tmp/4gcff1293639134.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/wessaorg/rcomp/tmp/5gcff1293639134.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
-392.1615896 -1019.7563837 -453.1708724 -804.9260504 360.7594783
6 7 8 9 10
-279.0235781 -321.9235182 -463.1453074 -588.8512942 -401.8339358
11 12 13 14 15
-742.3992506 -764.9346651 -470.5344475 -322.5549404 -429.2427177
16 17 18 19 20
-312.1408758 -468.3176302 -484.7509013 -481.9832552 -460.0219261
21 22 23 24 25
-470.5344475 -358.3457368 -470.5344475 -1036.0095528 -705.2365682
26 27 28 29 30
-450.7663177 -356.2347827 -550.7845271 -470.5344475 -534.4957791
31 32 33 34 35
134.8952139 -470.5344475 -538.0085226 -1460.5526455 -221.6013291
36 37 38 39 40
-24.2183041 -468.0152199 -137.0108679 -730.3905299 -660.3075539
41 42 43 44 45
-1040.4291487 -773.7211062 -670.8199921 -470.5344475 -355.1175642
46 47 48 49 50
-351.4452806 -470.5344475 -427.9140673 -495.4713487 -790.6830912
51 52 53 54 55
-294.7326741 -631.9214648 -165.3021303 -614.0887094 -3205.6775652
56 57 58 59 60
-1013.6751239 -450.9765612 0.2260205 -314.2046138 -470.5344475
61 62 63 64 65
-470.5344475 -314.6810109 104.9268068 -472.3085284 730.8280805
66 67 68 69 70
-465.0419398 -403.4219727 -515.6566081 -373.6832878 -470.5344475
71 72 73 74 75
-3378.6463520 297.7009682 8771.2754479 -470.5344475 -519.9980963
76 77 78 79 80
-530.6796029 -580.5956424 -458.4692428 1065.4420485 -180.9583579
81 82 83 84 85
-467.6732617 -470.5344475 -470.5344475 -502.2014647 -470.5344475
86 87 88 89 90
-470.5449476 -460.3507269 -470.2362739 -374.2722757 -569.5353756
91 92 93 94 95
-470.5449476 -359.2706537 -934.1458012 545.0145513 -505.2476610
96 97 98 99 100
-546.9441833 -470.5344475 -451.2447467 -477.8018249 -250.8679047
101 102 103 104 105
-1031.2981959 -377.9418429 -476.3206260 -470.3494670 -470.5344475
106 107 108 109 110
1977.8810030 -390.3697034 -1731.1824846 -336.8265263 -458.2571058
111 112 113 114 115
-476.6329577 160.5315987 -470.5344475 -446.6151630 -470.5449476
116 117 118 119 120
-416.7551092 -760.3086408 -1896.0208159 -2010.2224842 -1339.9715274
121 122 123 124 125
87.0682504 -470.5344475 -2069.8586465 -739.1645505 19.0167286
126 127 128 129 130
-470.5449476 -449.8260321 -470.5344475 1713.7209097 -386.0915604
131 132 133 134 135
-512.7086149 -470.5449476 -470.5449476 -470.5449476 -459.9319090
136 137 138 139 140
-409.3731914 315.3994283 271.9793364 315.3994283 315.3994283
141 142 143 144 145
-100.8336478 822.4580612 315.3994283 314.3994283 324.6291894
146 147 148 149 150
315.3994283 325.5178793 397.9112373 247.2483103 314.3994283
151 152 153 154 155
315.3994283 345.5406820 315.3994283 355.1577545 258.5359260
156 157 158 159 160
162.0830619 304.1319653 315.3994283 297.6329451 -327.9527531
161 162 163 164 165
322.2723763 315.3994283 376.0361236 269.1539227 314.3994283
166 167 168 169 170
-236.8998305 315.3994283 320.9979119 661.8501481 216.4049847
171 172 173 174 175
-620.2062952 -196.2867470 125.9323451 -73.3090202 493.2488721
176 177 178 179 180
35.4909132 315.3994283 317.0125444 307.5333072 315.3994283
181 182 183 184 185
-486.2091934 -275.9641179 315.3994283 314.3994283 315.3994283
186 187 188 189 190
315.7741848 304.8064760 273.8043439 2420.8634789 315.3994283
191 192 193 194 195
315.0847766 280.6204554 315.3994283 315.8441910 172.6253245
196 197 198 199 200
315.3994283 -288.3499046 546.4077924 315.3994283 314.3994283
201 202 203 204 205
242.3681110 314.3994283 315.3994283 321.7409222 315.3994283
206 207 208 209 210
-58.4255405 36.9547142 315.3994283 -42.4623428 411.5226277
211 212 213 214 215
-227.2172816 155.9266370 304.7732975 352.6581697 -1025.0363920
216 217 218 219 220
139.5015577 302.6968945 314.3994283 210.8843809 -469.1418232
221 222 223 224 225
314.0518467 315.3994283 328.6259273 291.0572476 351.7487357
226 227 228 229 230
-58.8305535 315.3994283 -2765.8363590 159.2645634 262.0055034
231 232 233 234 235
-3005.3518062 381.8741556 -210.0222391 315.3994283 317.0127834
236 237 238 239 240
190.4365370 389.9566574 315.3994283 315.3994283 296.4621526
241 242 243 244 245
137.2401429 64.8448672 315.3994283 322.0098201 486.7721881
246 247 248 249 250
692.0842751 315.3994283 330.7688591 315.3994283 315.0736409
251 252 253 254 255
315.3994283 74.4691667 315.3994283 314.3994283 315.3994283
256 257 258 259 260
315.3994283 315.3994283 315.3994283 7.6359586 -31.0356061
261 262 263 264 265
456.7097840 621.1805299 217.9851292 -453.9265658 194.3743963
266 267 268 269 270
315.3994283 316.3843162 315.3994283 179.3722050 313.5008966
271 272 273 274 275
315.3994283 315.3994283 208.2126011 312.8224200 325.6942428
276 277 278 279 280
455.9650879 315.7686817 335.0404225 357.0034963 325.4853920
281 282 283 284 285
321.8809852 324.1832544 315.3994283 277.4683358 315.7987361
286 287 288 289 290
334.8665928 315.3994283 216.7893830 740.5007750 1520.2646473
291 292 293 294 295
315.3994283 323.0512807 315.3994283 278.0031205 314.1572864
296 297 298 299 300
315.1544314 314.3994283 1238.8720897 314.3994283 314.3994283
301 302 303 304 305
650.5540791 247.7804980 -300.3015591 308.0942909 255.5358441
306 307 308 309 310
325.8886133 -393.2884841 315.3994283 315.3994283 314.6297502
311 312 313 314 315
315.3994283 50.1942145 273.0411048 313.4492315 217.8344257
316 317 318 319 320
97.0458816 -555.3755167 351.7294826 315.3994283 209.5370554
321 322 323 324 325
314.4961618 336.5024572 -10.4298297 315.3994283 -418.8891389
326 327 328 329 330
315.3994283 -750.7199719 -1012.2136382 1177.9550887 41.7899478
331 332 333 334 335
153.5326488 314.3994283 314.3994283 -1482.4629570 120.4144274
336 337 338 339 340
315.3994283 -110.3046243 -266.9372740 314.3994283 315.3994283
341 342 343 344 345
-255.2058402 307.5539870 315.3994283 -63.5229823 314.3994283
346 347 348 349 350
62.0263233 306.0083792 522.2048761 315.3994283 325.7692883
351 352 353 354 355
315.3994283 314.3994283 315.1602883 315.3994283 180.4040504
356 357 358 359 360
315.3994283 314.3994283 315.3994283 314.3994283 -673.0071334
361 362 363 364 365
-931.6266442 615.9307969 382.3329805 -45.3781510 342.8701291
366 367 368 369 370
539.3398741 329.3313646 321.3694143 314.3994283 300.6128820
371 372 373 374 375
-404.8978638 355.7087606 288.3572654 116.7488226 290.8510489
376 377 378 379 380
-200.6059805 -48.7009041 54.3181962 269.8938791 315.3994283
381 382 383 384 385
242.9919123 250.9922927 111.3867762 283.4126263 313.5231695
386 387 388 389 390
351.3505507 348.2096572 315.3994283 314.3994283 315.3994283
391 392 393 394 395
315.3994283 315.3994283 95.3249937 328.4920881 314.3994283
396 397 398 399 400
315.3994283 315.3994283 315.3994283 330.2700782 -262.0387481
401 402 403 404 405
-298.5011399 44.4788159 404.2790165 -345.2528458 292.0649056
406 407 408 409 410
325.1858653 317.9770930 -56.1874205 -180.8788554 -466.4287385
411 412 413 414 415
296.6902579 314.3994283 315.3994283 206.2554596 323.9391740
416 417 418 419 420
331.7373343 -1380.2628987 60.2731184 315.4376576 315.3994283
421 422 423 424 425
321.5951448 649.9036455 -1784.7849797 323.8040500 218.8220148
426 427 428 429 430
494.7405549 315.3994283 148.9079456 216.2296820 315.3994283
431
377.5161438
> postscript(file="/var/www/wessaorg/rcomp/tmp/6gcff1293639134.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 -392.1615896 NA
1 -1019.7563837 -392.1615896
2 -453.1708724 -1019.7563837
3 -804.9260504 -453.1708724
4 360.7594783 -804.9260504
5 -279.0235781 360.7594783
6 -321.9235182 -279.0235781
7 -463.1453074 -321.9235182
8 -588.8512942 -463.1453074
9 -401.8339358 -588.8512942
10 -742.3992506 -401.8339358
11 -764.9346651 -742.3992506
12 -470.5344475 -764.9346651
13 -322.5549404 -470.5344475
14 -429.2427177 -322.5549404
15 -312.1408758 -429.2427177
16 -468.3176302 -312.1408758
17 -484.7509013 -468.3176302
18 -481.9832552 -484.7509013
19 -460.0219261 -481.9832552
20 -470.5344475 -460.0219261
21 -358.3457368 -470.5344475
22 -470.5344475 -358.3457368
23 -1036.0095528 -470.5344475
24 -705.2365682 -1036.0095528
25 -450.7663177 -705.2365682
26 -356.2347827 -450.7663177
27 -550.7845271 -356.2347827
28 -470.5344475 -550.7845271
29 -534.4957791 -470.5344475
30 134.8952139 -534.4957791
31 -470.5344475 134.8952139
32 -538.0085226 -470.5344475
33 -1460.5526455 -538.0085226
34 -221.6013291 -1460.5526455
35 -24.2183041 -221.6013291
36 -468.0152199 -24.2183041
37 -137.0108679 -468.0152199
38 -730.3905299 -137.0108679
39 -660.3075539 -730.3905299
40 -1040.4291487 -660.3075539
41 -773.7211062 -1040.4291487
42 -670.8199921 -773.7211062
43 -470.5344475 -670.8199921
44 -355.1175642 -470.5344475
45 -351.4452806 -355.1175642
46 -470.5344475 -351.4452806
47 -427.9140673 -470.5344475
48 -495.4713487 -427.9140673
49 -790.6830912 -495.4713487
50 -294.7326741 -790.6830912
51 -631.9214648 -294.7326741
52 -165.3021303 -631.9214648
53 -614.0887094 -165.3021303
54 -3205.6775652 -614.0887094
55 -1013.6751239 -3205.6775652
56 -450.9765612 -1013.6751239
57 0.2260205 -450.9765612
58 -314.2046138 0.2260205
59 -470.5344475 -314.2046138
60 -470.5344475 -470.5344475
61 -314.6810109 -470.5344475
62 104.9268068 -314.6810109
63 -472.3085284 104.9268068
64 730.8280805 -472.3085284
65 -465.0419398 730.8280805
66 -403.4219727 -465.0419398
67 -515.6566081 -403.4219727
68 -373.6832878 -515.6566081
69 -470.5344475 -373.6832878
70 -3378.6463520 -470.5344475
71 297.7009682 -3378.6463520
72 8771.2754479 297.7009682
73 -470.5344475 8771.2754479
74 -519.9980963 -470.5344475
75 -530.6796029 -519.9980963
76 -580.5956424 -530.6796029
77 -458.4692428 -580.5956424
78 1065.4420485 -458.4692428
79 -180.9583579 1065.4420485
80 -467.6732617 -180.9583579
81 -470.5344475 -467.6732617
82 -470.5344475 -470.5344475
83 -502.2014647 -470.5344475
84 -470.5344475 -502.2014647
85 -470.5449476 -470.5344475
86 -460.3507269 -470.5449476
87 -470.2362739 -460.3507269
88 -374.2722757 -470.2362739
89 -569.5353756 -374.2722757
90 -470.5449476 -569.5353756
91 -359.2706537 -470.5449476
92 -934.1458012 -359.2706537
93 545.0145513 -934.1458012
94 -505.2476610 545.0145513
95 -546.9441833 -505.2476610
96 -470.5344475 -546.9441833
97 -451.2447467 -470.5344475
98 -477.8018249 -451.2447467
99 -250.8679047 -477.8018249
100 -1031.2981959 -250.8679047
101 -377.9418429 -1031.2981959
102 -476.3206260 -377.9418429
103 -470.3494670 -476.3206260
104 -470.5344475 -470.3494670
105 1977.8810030 -470.5344475
106 -390.3697034 1977.8810030
107 -1731.1824846 -390.3697034
108 -336.8265263 -1731.1824846
109 -458.2571058 -336.8265263
110 -476.6329577 -458.2571058
111 160.5315987 -476.6329577
112 -470.5344475 160.5315987
113 -446.6151630 -470.5344475
114 -470.5449476 -446.6151630
115 -416.7551092 -470.5449476
116 -760.3086408 -416.7551092
117 -1896.0208159 -760.3086408
118 -2010.2224842 -1896.0208159
119 -1339.9715274 -2010.2224842
120 87.0682504 -1339.9715274
121 -470.5344475 87.0682504
122 -2069.8586465 -470.5344475
123 -739.1645505 -2069.8586465
124 19.0167286 -739.1645505
125 -470.5449476 19.0167286
126 -449.8260321 -470.5449476
127 -470.5344475 -449.8260321
128 1713.7209097 -470.5344475
129 -386.0915604 1713.7209097
130 -512.7086149 -386.0915604
131 -470.5449476 -512.7086149
132 -470.5449476 -470.5449476
133 -470.5449476 -470.5449476
134 -459.9319090 -470.5449476
135 -409.3731914 -459.9319090
136 315.3994283 -409.3731914
137 271.9793364 315.3994283
138 315.3994283 271.9793364
139 315.3994283 315.3994283
140 -100.8336478 315.3994283
141 822.4580612 -100.8336478
142 315.3994283 822.4580612
143 314.3994283 315.3994283
144 324.6291894 314.3994283
145 315.3994283 324.6291894
146 325.5178793 315.3994283
147 397.9112373 325.5178793
148 247.2483103 397.9112373
149 314.3994283 247.2483103
150 315.3994283 314.3994283
151 345.5406820 315.3994283
152 315.3994283 345.5406820
153 355.1577545 315.3994283
154 258.5359260 355.1577545
155 162.0830619 258.5359260
156 304.1319653 162.0830619
157 315.3994283 304.1319653
158 297.6329451 315.3994283
159 -327.9527531 297.6329451
160 322.2723763 -327.9527531
161 315.3994283 322.2723763
162 376.0361236 315.3994283
163 269.1539227 376.0361236
164 314.3994283 269.1539227
165 -236.8998305 314.3994283
166 315.3994283 -236.8998305
167 320.9979119 315.3994283
168 661.8501481 320.9979119
169 216.4049847 661.8501481
170 -620.2062952 216.4049847
171 -196.2867470 -620.2062952
172 125.9323451 -196.2867470
173 -73.3090202 125.9323451
174 493.2488721 -73.3090202
175 35.4909132 493.2488721
176 315.3994283 35.4909132
177 317.0125444 315.3994283
178 307.5333072 317.0125444
179 315.3994283 307.5333072
180 -486.2091934 315.3994283
181 -275.9641179 -486.2091934
182 315.3994283 -275.9641179
183 314.3994283 315.3994283
184 315.3994283 314.3994283
185 315.7741848 315.3994283
186 304.8064760 315.7741848
187 273.8043439 304.8064760
188 2420.8634789 273.8043439
189 315.3994283 2420.8634789
190 315.0847766 315.3994283
191 280.6204554 315.0847766
192 315.3994283 280.6204554
193 315.8441910 315.3994283
194 172.6253245 315.8441910
195 315.3994283 172.6253245
196 -288.3499046 315.3994283
197 546.4077924 -288.3499046
198 315.3994283 546.4077924
199 314.3994283 315.3994283
200 242.3681110 314.3994283
201 314.3994283 242.3681110
202 315.3994283 314.3994283
203 321.7409222 315.3994283
204 315.3994283 321.7409222
205 -58.4255405 315.3994283
206 36.9547142 -58.4255405
207 315.3994283 36.9547142
208 -42.4623428 315.3994283
209 411.5226277 -42.4623428
210 -227.2172816 411.5226277
211 155.9266370 -227.2172816
212 304.7732975 155.9266370
213 352.6581697 304.7732975
214 -1025.0363920 352.6581697
215 139.5015577 -1025.0363920
216 302.6968945 139.5015577
217 314.3994283 302.6968945
218 210.8843809 314.3994283
219 -469.1418232 210.8843809
220 314.0518467 -469.1418232
221 315.3994283 314.0518467
222 328.6259273 315.3994283
223 291.0572476 328.6259273
224 351.7487357 291.0572476
225 -58.8305535 351.7487357
226 315.3994283 -58.8305535
227 -2765.8363590 315.3994283
228 159.2645634 -2765.8363590
229 262.0055034 159.2645634
230 -3005.3518062 262.0055034
231 381.8741556 -3005.3518062
232 -210.0222391 381.8741556
233 315.3994283 -210.0222391
234 317.0127834 315.3994283
235 190.4365370 317.0127834
236 389.9566574 190.4365370
237 315.3994283 389.9566574
238 315.3994283 315.3994283
239 296.4621526 315.3994283
240 137.2401429 296.4621526
241 64.8448672 137.2401429
242 315.3994283 64.8448672
243 322.0098201 315.3994283
244 486.7721881 322.0098201
245 692.0842751 486.7721881
246 315.3994283 692.0842751
247 330.7688591 315.3994283
248 315.3994283 330.7688591
249 315.0736409 315.3994283
250 315.3994283 315.0736409
251 74.4691667 315.3994283
252 315.3994283 74.4691667
253 314.3994283 315.3994283
254 315.3994283 314.3994283
255 315.3994283 315.3994283
256 315.3994283 315.3994283
257 315.3994283 315.3994283
258 7.6359586 315.3994283
259 -31.0356061 7.6359586
260 456.7097840 -31.0356061
261 621.1805299 456.7097840
262 217.9851292 621.1805299
263 -453.9265658 217.9851292
264 194.3743963 -453.9265658
265 315.3994283 194.3743963
266 316.3843162 315.3994283
267 315.3994283 316.3843162
268 179.3722050 315.3994283
269 313.5008966 179.3722050
270 315.3994283 313.5008966
271 315.3994283 315.3994283
272 208.2126011 315.3994283
273 312.8224200 208.2126011
274 325.6942428 312.8224200
275 455.9650879 325.6942428
276 315.7686817 455.9650879
277 335.0404225 315.7686817
278 357.0034963 335.0404225
279 325.4853920 357.0034963
280 321.8809852 325.4853920
281 324.1832544 321.8809852
282 315.3994283 324.1832544
283 277.4683358 315.3994283
284 315.7987361 277.4683358
285 334.8665928 315.7987361
286 315.3994283 334.8665928
287 216.7893830 315.3994283
288 740.5007750 216.7893830
289 1520.2646473 740.5007750
290 315.3994283 1520.2646473
291 323.0512807 315.3994283
292 315.3994283 323.0512807
293 278.0031205 315.3994283
294 314.1572864 278.0031205
295 315.1544314 314.1572864
296 314.3994283 315.1544314
297 1238.8720897 314.3994283
298 314.3994283 1238.8720897
299 314.3994283 314.3994283
300 650.5540791 314.3994283
301 247.7804980 650.5540791
302 -300.3015591 247.7804980
303 308.0942909 -300.3015591
304 255.5358441 308.0942909
305 325.8886133 255.5358441
306 -393.2884841 325.8886133
307 315.3994283 -393.2884841
308 315.3994283 315.3994283
309 314.6297502 315.3994283
310 315.3994283 314.6297502
311 50.1942145 315.3994283
312 273.0411048 50.1942145
313 313.4492315 273.0411048
314 217.8344257 313.4492315
315 97.0458816 217.8344257
316 -555.3755167 97.0458816
317 351.7294826 -555.3755167
318 315.3994283 351.7294826
319 209.5370554 315.3994283
320 314.4961618 209.5370554
321 336.5024572 314.4961618
322 -10.4298297 336.5024572
323 315.3994283 -10.4298297
324 -418.8891389 315.3994283
325 315.3994283 -418.8891389
326 -750.7199719 315.3994283
327 -1012.2136382 -750.7199719
328 1177.9550887 -1012.2136382
329 41.7899478 1177.9550887
330 153.5326488 41.7899478
331 314.3994283 153.5326488
332 314.3994283 314.3994283
333 -1482.4629570 314.3994283
334 120.4144274 -1482.4629570
335 315.3994283 120.4144274
336 -110.3046243 315.3994283
337 -266.9372740 -110.3046243
338 314.3994283 -266.9372740
339 315.3994283 314.3994283
340 -255.2058402 315.3994283
341 307.5539870 -255.2058402
342 315.3994283 307.5539870
343 -63.5229823 315.3994283
344 314.3994283 -63.5229823
345 62.0263233 314.3994283
346 306.0083792 62.0263233
347 522.2048761 306.0083792
348 315.3994283 522.2048761
349 325.7692883 315.3994283
350 315.3994283 325.7692883
351 314.3994283 315.3994283
352 315.1602883 314.3994283
353 315.3994283 315.1602883
354 180.4040504 315.3994283
355 315.3994283 180.4040504
356 314.3994283 315.3994283
357 315.3994283 314.3994283
358 314.3994283 315.3994283
359 -673.0071334 314.3994283
360 -931.6266442 -673.0071334
361 615.9307969 -931.6266442
362 382.3329805 615.9307969
363 -45.3781510 382.3329805
364 342.8701291 -45.3781510
365 539.3398741 342.8701291
366 329.3313646 539.3398741
367 321.3694143 329.3313646
368 314.3994283 321.3694143
369 300.6128820 314.3994283
370 -404.8978638 300.6128820
371 355.7087606 -404.8978638
372 288.3572654 355.7087606
373 116.7488226 288.3572654
374 290.8510489 116.7488226
375 -200.6059805 290.8510489
376 -48.7009041 -200.6059805
377 54.3181962 -48.7009041
378 269.8938791 54.3181962
379 315.3994283 269.8938791
380 242.9919123 315.3994283
381 250.9922927 242.9919123
382 111.3867762 250.9922927
383 283.4126263 111.3867762
384 313.5231695 283.4126263
385 351.3505507 313.5231695
386 348.2096572 351.3505507
387 315.3994283 348.2096572
388 314.3994283 315.3994283
389 315.3994283 314.3994283
390 315.3994283 315.3994283
391 315.3994283 315.3994283
392 95.3249937 315.3994283
393 328.4920881 95.3249937
394 314.3994283 328.4920881
395 315.3994283 314.3994283
396 315.3994283 315.3994283
397 315.3994283 315.3994283
398 330.2700782 315.3994283
399 -262.0387481 330.2700782
400 -298.5011399 -262.0387481
401 44.4788159 -298.5011399
402 404.2790165 44.4788159
403 -345.2528458 404.2790165
404 292.0649056 -345.2528458
405 325.1858653 292.0649056
406 317.9770930 325.1858653
407 -56.1874205 317.9770930
408 -180.8788554 -56.1874205
409 -466.4287385 -180.8788554
410 296.6902579 -466.4287385
411 314.3994283 296.6902579
412 315.3994283 314.3994283
413 206.2554596 315.3994283
414 323.9391740 206.2554596
415 331.7373343 323.9391740
416 -1380.2628987 331.7373343
417 60.2731184 -1380.2628987
418 315.4376576 60.2731184
419 315.3994283 315.4376576
420 321.5951448 315.3994283
421 649.9036455 321.5951448
422 -1784.7849797 649.9036455
423 323.8040500 -1784.7849797
424 218.8220148 323.8040500
425 494.7405549 218.8220148
426 315.3994283 494.7405549
427 148.9079456 315.3994283
428 216.2296820 148.9079456
429 315.3994283 216.2296820
430 377.5161438 315.3994283
431 NA 377.5161438
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] -1019.7563837 -392.1615896
[2,] -453.1708724 -1019.7563837
[3,] -804.9260504 -453.1708724
[4,] 360.7594783 -804.9260504
[5,] -279.0235781 360.7594783
[6,] -321.9235182 -279.0235781
[7,] -463.1453074 -321.9235182
[8,] -588.8512942 -463.1453074
[9,] -401.8339358 -588.8512942
[10,] -742.3992506 -401.8339358
[11,] -764.9346651 -742.3992506
[12,] -470.5344475 -764.9346651
[13,] -322.5549404 -470.5344475
[14,] -429.2427177 -322.5549404
[15,] -312.1408758 -429.2427177
[16,] -468.3176302 -312.1408758
[17,] -484.7509013 -468.3176302
[18,] -481.9832552 -484.7509013
[19,] -460.0219261 -481.9832552
[20,] -470.5344475 -460.0219261
[21,] -358.3457368 -470.5344475
[22,] -470.5344475 -358.3457368
[23,] -1036.0095528 -470.5344475
[24,] -705.2365682 -1036.0095528
[25,] -450.7663177 -705.2365682
[26,] -356.2347827 -450.7663177
[27,] -550.7845271 -356.2347827
[28,] -470.5344475 -550.7845271
[29,] -534.4957791 -470.5344475
[30,] 134.8952139 -534.4957791
[31,] -470.5344475 134.8952139
[32,] -538.0085226 -470.5344475
[33,] -1460.5526455 -538.0085226
[34,] -221.6013291 -1460.5526455
[35,] -24.2183041 -221.6013291
[36,] -468.0152199 -24.2183041
[37,] -137.0108679 -468.0152199
[38,] -730.3905299 -137.0108679
[39,] -660.3075539 -730.3905299
[40,] -1040.4291487 -660.3075539
[41,] -773.7211062 -1040.4291487
[42,] -670.8199921 -773.7211062
[43,] -470.5344475 -670.8199921
[44,] -355.1175642 -470.5344475
[45,] -351.4452806 -355.1175642
[46,] -470.5344475 -351.4452806
[47,] -427.9140673 -470.5344475
[48,] -495.4713487 -427.9140673
[49,] -790.6830912 -495.4713487
[50,] -294.7326741 -790.6830912
[51,] -631.9214648 -294.7326741
[52,] -165.3021303 -631.9214648
[53,] -614.0887094 -165.3021303
[54,] -3205.6775652 -614.0887094
[55,] -1013.6751239 -3205.6775652
[56,] -450.9765612 -1013.6751239
[57,] 0.2260205 -450.9765612
[58,] -314.2046138 0.2260205
[59,] -470.5344475 -314.2046138
[60,] -470.5344475 -470.5344475
[61,] -314.6810109 -470.5344475
[62,] 104.9268068 -314.6810109
[63,] -472.3085284 104.9268068
[64,] 730.8280805 -472.3085284
[65,] -465.0419398 730.8280805
[66,] -403.4219727 -465.0419398
[67,] -515.6566081 -403.4219727
[68,] -373.6832878 -515.6566081
[69,] -470.5344475 -373.6832878
[70,] -3378.6463520 -470.5344475
[71,] 297.7009682 -3378.6463520
[72,] 8771.2754479 297.7009682
[73,] -470.5344475 8771.2754479
[74,] -519.9980963 -470.5344475
[75,] -530.6796029 -519.9980963
[76,] -580.5956424 -530.6796029
[77,] -458.4692428 -580.5956424
[78,] 1065.4420485 -458.4692428
[79,] -180.9583579 1065.4420485
[80,] -467.6732617 -180.9583579
[81,] -470.5344475 -467.6732617
[82,] -470.5344475 -470.5344475
[83,] -502.2014647 -470.5344475
[84,] -470.5344475 -502.2014647
[85,] -470.5449476 -470.5344475
[86,] -460.3507269 -470.5449476
[87,] -470.2362739 -460.3507269
[88,] -374.2722757 -470.2362739
[89,] -569.5353756 -374.2722757
[90,] -470.5449476 -569.5353756
[91,] -359.2706537 -470.5449476
[92,] -934.1458012 -359.2706537
[93,] 545.0145513 -934.1458012
[94,] -505.2476610 545.0145513
[95,] -546.9441833 -505.2476610
[96,] -470.5344475 -546.9441833
[97,] -451.2447467 -470.5344475
[98,] -477.8018249 -451.2447467
[99,] -250.8679047 -477.8018249
[100,] -1031.2981959 -250.8679047
[101,] -377.9418429 -1031.2981959
[102,] -476.3206260 -377.9418429
[103,] -470.3494670 -476.3206260
[104,] -470.5344475 -470.3494670
[105,] 1977.8810030 -470.5344475
[106,] -390.3697034 1977.8810030
[107,] -1731.1824846 -390.3697034
[108,] -336.8265263 -1731.1824846
[109,] -458.2571058 -336.8265263
[110,] -476.6329577 -458.2571058
[111,] 160.5315987 -476.6329577
[112,] -470.5344475 160.5315987
[113,] -446.6151630 -470.5344475
[114,] -470.5449476 -446.6151630
[115,] -416.7551092 -470.5449476
[116,] -760.3086408 -416.7551092
[117,] -1896.0208159 -760.3086408
[118,] -2010.2224842 -1896.0208159
[119,] -1339.9715274 -2010.2224842
[120,] 87.0682504 -1339.9715274
[121,] -470.5344475 87.0682504
[122,] -2069.8586465 -470.5344475
[123,] -739.1645505 -2069.8586465
[124,] 19.0167286 -739.1645505
[125,] -470.5449476 19.0167286
[126,] -449.8260321 -470.5449476
[127,] -470.5344475 -449.8260321
[128,] 1713.7209097 -470.5344475
[129,] -386.0915604 1713.7209097
[130,] -512.7086149 -386.0915604
[131,] -470.5449476 -512.7086149
[132,] -470.5449476 -470.5449476
[133,] -470.5449476 -470.5449476
[134,] -459.9319090 -470.5449476
[135,] -409.3731914 -459.9319090
[136,] 315.3994283 -409.3731914
[137,] 271.9793364 315.3994283
[138,] 315.3994283 271.9793364
[139,] 315.3994283 315.3994283
[140,] -100.8336478 315.3994283
[141,] 822.4580612 -100.8336478
[142,] 315.3994283 822.4580612
[143,] 314.3994283 315.3994283
[144,] 324.6291894 314.3994283
[145,] 315.3994283 324.6291894
[146,] 325.5178793 315.3994283
[147,] 397.9112373 325.5178793
[148,] 247.2483103 397.9112373
[149,] 314.3994283 247.2483103
[150,] 315.3994283 314.3994283
[151,] 345.5406820 315.3994283
[152,] 315.3994283 345.5406820
[153,] 355.1577545 315.3994283
[154,] 258.5359260 355.1577545
[155,] 162.0830619 258.5359260
[156,] 304.1319653 162.0830619
[157,] 315.3994283 304.1319653
[158,] 297.6329451 315.3994283
[159,] -327.9527531 297.6329451
[160,] 322.2723763 -327.9527531
[161,] 315.3994283 322.2723763
[162,] 376.0361236 315.3994283
[163,] 269.1539227 376.0361236
[164,] 314.3994283 269.1539227
[165,] -236.8998305 314.3994283
[166,] 315.3994283 -236.8998305
[167,] 320.9979119 315.3994283
[168,] 661.8501481 320.9979119
[169,] 216.4049847 661.8501481
[170,] -620.2062952 216.4049847
[171,] -196.2867470 -620.2062952
[172,] 125.9323451 -196.2867470
[173,] -73.3090202 125.9323451
[174,] 493.2488721 -73.3090202
[175,] 35.4909132 493.2488721
[176,] 315.3994283 35.4909132
[177,] 317.0125444 315.3994283
[178,] 307.5333072 317.0125444
[179,] 315.3994283 307.5333072
[180,] -486.2091934 315.3994283
[181,] -275.9641179 -486.2091934
[182,] 315.3994283 -275.9641179
[183,] 314.3994283 315.3994283
[184,] 315.3994283 314.3994283
[185,] 315.7741848 315.3994283
[186,] 304.8064760 315.7741848
[187,] 273.8043439 304.8064760
[188,] 2420.8634789 273.8043439
[189,] 315.3994283 2420.8634789
[190,] 315.0847766 315.3994283
[191,] 280.6204554 315.0847766
[192,] 315.3994283 280.6204554
[193,] 315.8441910 315.3994283
[194,] 172.6253245 315.8441910
[195,] 315.3994283 172.6253245
[196,] -288.3499046 315.3994283
[197,] 546.4077924 -288.3499046
[198,] 315.3994283 546.4077924
[199,] 314.3994283 315.3994283
[200,] 242.3681110 314.3994283
[201,] 314.3994283 242.3681110
[202,] 315.3994283 314.3994283
[203,] 321.7409222 315.3994283
[204,] 315.3994283 321.7409222
[205,] -58.4255405 315.3994283
[206,] 36.9547142 -58.4255405
[207,] 315.3994283 36.9547142
[208,] -42.4623428 315.3994283
[209,] 411.5226277 -42.4623428
[210,] -227.2172816 411.5226277
[211,] 155.9266370 -227.2172816
[212,] 304.7732975 155.9266370
[213,] 352.6581697 304.7732975
[214,] -1025.0363920 352.6581697
[215,] 139.5015577 -1025.0363920
[216,] 302.6968945 139.5015577
[217,] 314.3994283 302.6968945
[218,] 210.8843809 314.3994283
[219,] -469.1418232 210.8843809
[220,] 314.0518467 -469.1418232
[221,] 315.3994283 314.0518467
[222,] 328.6259273 315.3994283
[223,] 291.0572476 328.6259273
[224,] 351.7487357 291.0572476
[225,] -58.8305535 351.7487357
[226,] 315.3994283 -58.8305535
[227,] -2765.8363590 315.3994283
[228,] 159.2645634 -2765.8363590
[229,] 262.0055034 159.2645634
[230,] -3005.3518062 262.0055034
[231,] 381.8741556 -3005.3518062
[232,] -210.0222391 381.8741556
[233,] 315.3994283 -210.0222391
[234,] 317.0127834 315.3994283
[235,] 190.4365370 317.0127834
[236,] 389.9566574 190.4365370
[237,] 315.3994283 389.9566574
[238,] 315.3994283 315.3994283
[239,] 296.4621526 315.3994283
[240,] 137.2401429 296.4621526
[241,] 64.8448672 137.2401429
[242,] 315.3994283 64.8448672
[243,] 322.0098201 315.3994283
[244,] 486.7721881 322.0098201
[245,] 692.0842751 486.7721881
[246,] 315.3994283 692.0842751
[247,] 330.7688591 315.3994283
[248,] 315.3994283 330.7688591
[249,] 315.0736409 315.3994283
[250,] 315.3994283 315.0736409
[251,] 74.4691667 315.3994283
[252,] 315.3994283 74.4691667
[253,] 314.3994283 315.3994283
[254,] 315.3994283 314.3994283
[255,] 315.3994283 315.3994283
[256,] 315.3994283 315.3994283
[257,] 315.3994283 315.3994283
[258,] 7.6359586 315.3994283
[259,] -31.0356061 7.6359586
[260,] 456.7097840 -31.0356061
[261,] 621.1805299 456.7097840
[262,] 217.9851292 621.1805299
[263,] -453.9265658 217.9851292
[264,] 194.3743963 -453.9265658
[265,] 315.3994283 194.3743963
[266,] 316.3843162 315.3994283
[267,] 315.3994283 316.3843162
[268,] 179.3722050 315.3994283
[269,] 313.5008966 179.3722050
[270,] 315.3994283 313.5008966
[271,] 315.3994283 315.3994283
[272,] 208.2126011 315.3994283
[273,] 312.8224200 208.2126011
[274,] 325.6942428 312.8224200
[275,] 455.9650879 325.6942428
[276,] 315.7686817 455.9650879
[277,] 335.0404225 315.7686817
[278,] 357.0034963 335.0404225
[279,] 325.4853920 357.0034963
[280,] 321.8809852 325.4853920
[281,] 324.1832544 321.8809852
[282,] 315.3994283 324.1832544
[283,] 277.4683358 315.3994283
[284,] 315.7987361 277.4683358
[285,] 334.8665928 315.7987361
[286,] 315.3994283 334.8665928
[287,] 216.7893830 315.3994283
[288,] 740.5007750 216.7893830
[289,] 1520.2646473 740.5007750
[290,] 315.3994283 1520.2646473
[291,] 323.0512807 315.3994283
[292,] 315.3994283 323.0512807
[293,] 278.0031205 315.3994283
[294,] 314.1572864 278.0031205
[295,] 315.1544314 314.1572864
[296,] 314.3994283 315.1544314
[297,] 1238.8720897 314.3994283
[298,] 314.3994283 1238.8720897
[299,] 314.3994283 314.3994283
[300,] 650.5540791 314.3994283
[301,] 247.7804980 650.5540791
[302,] -300.3015591 247.7804980
[303,] 308.0942909 -300.3015591
[304,] 255.5358441 308.0942909
[305,] 325.8886133 255.5358441
[306,] -393.2884841 325.8886133
[307,] 315.3994283 -393.2884841
[308,] 315.3994283 315.3994283
[309,] 314.6297502 315.3994283
[310,] 315.3994283 314.6297502
[311,] 50.1942145 315.3994283
[312,] 273.0411048 50.1942145
[313,] 313.4492315 273.0411048
[314,] 217.8344257 313.4492315
[315,] 97.0458816 217.8344257
[316,] -555.3755167 97.0458816
[317,] 351.7294826 -555.3755167
[318,] 315.3994283 351.7294826
[319,] 209.5370554 315.3994283
[320,] 314.4961618 209.5370554
[321,] 336.5024572 314.4961618
[322,] -10.4298297 336.5024572
[323,] 315.3994283 -10.4298297
[324,] -418.8891389 315.3994283
[325,] 315.3994283 -418.8891389
[326,] -750.7199719 315.3994283
[327,] -1012.2136382 -750.7199719
[328,] 1177.9550887 -1012.2136382
[329,] 41.7899478 1177.9550887
[330,] 153.5326488 41.7899478
[331,] 314.3994283 153.5326488
[332,] 314.3994283 314.3994283
[333,] -1482.4629570 314.3994283
[334,] 120.4144274 -1482.4629570
[335,] 315.3994283 120.4144274
[336,] -110.3046243 315.3994283
[337,] -266.9372740 -110.3046243
[338,] 314.3994283 -266.9372740
[339,] 315.3994283 314.3994283
[340,] -255.2058402 315.3994283
[341,] 307.5539870 -255.2058402
[342,] 315.3994283 307.5539870
[343,] -63.5229823 315.3994283
[344,] 314.3994283 -63.5229823
[345,] 62.0263233 314.3994283
[346,] 306.0083792 62.0263233
[347,] 522.2048761 306.0083792
[348,] 315.3994283 522.2048761
[349,] 325.7692883 315.3994283
[350,] 315.3994283 325.7692883
[351,] 314.3994283 315.3994283
[352,] 315.1602883 314.3994283
[353,] 315.3994283 315.1602883
[354,] 180.4040504 315.3994283
[355,] 315.3994283 180.4040504
[356,] 314.3994283 315.3994283
[357,] 315.3994283 314.3994283
[358,] 314.3994283 315.3994283
[359,] -673.0071334 314.3994283
[360,] -931.6266442 -673.0071334
[361,] 615.9307969 -931.6266442
[362,] 382.3329805 615.9307969
[363,] -45.3781510 382.3329805
[364,] 342.8701291 -45.3781510
[365,] 539.3398741 342.8701291
[366,] 329.3313646 539.3398741
[367,] 321.3694143 329.3313646
[368,] 314.3994283 321.3694143
[369,] 300.6128820 314.3994283
[370,] -404.8978638 300.6128820
[371,] 355.7087606 -404.8978638
[372,] 288.3572654 355.7087606
[373,] 116.7488226 288.3572654
[374,] 290.8510489 116.7488226
[375,] -200.6059805 290.8510489
[376,] -48.7009041 -200.6059805
[377,] 54.3181962 -48.7009041
[378,] 269.8938791 54.3181962
[379,] 315.3994283 269.8938791
[380,] 242.9919123 315.3994283
[381,] 250.9922927 242.9919123
[382,] 111.3867762 250.9922927
[383,] 283.4126263 111.3867762
[384,] 313.5231695 283.4126263
[385,] 351.3505507 313.5231695
[386,] 348.2096572 351.3505507
[387,] 315.3994283 348.2096572
[388,] 314.3994283 315.3994283
[389,] 315.3994283 314.3994283
[390,] 315.3994283 315.3994283
[391,] 315.3994283 315.3994283
[392,] 95.3249937 315.3994283
[393,] 328.4920881 95.3249937
[394,] 314.3994283 328.4920881
[395,] 315.3994283 314.3994283
[396,] 315.3994283 315.3994283
[397,] 315.3994283 315.3994283
[398,] 330.2700782 315.3994283
[399,] -262.0387481 330.2700782
[400,] -298.5011399 -262.0387481
[401,] 44.4788159 -298.5011399
[402,] 404.2790165 44.4788159
[403,] -345.2528458 404.2790165
[404,] 292.0649056 -345.2528458
[405,] 325.1858653 292.0649056
[406,] 317.9770930 325.1858653
[407,] -56.1874205 317.9770930
[408,] -180.8788554 -56.1874205
[409,] -466.4287385 -180.8788554
[410,] 296.6902579 -466.4287385
[411,] 314.3994283 296.6902579
[412,] 315.3994283 314.3994283
[413,] 206.2554596 315.3994283
[414,] 323.9391740 206.2554596
[415,] 331.7373343 323.9391740
[416,] -1380.2628987 331.7373343
[417,] 60.2731184 -1380.2628987
[418,] 315.4376576 60.2731184
[419,] 315.3994283 315.4376576
[420,] 321.5951448 315.3994283
[421,] 649.9036455 321.5951448
[422,] -1784.7849797 649.9036455
[423,] 323.8040500 -1784.7849797
[424,] 218.8220148 323.8040500
[425,] 494.7405549 218.8220148
[426,] 315.3994283 494.7405549
[427,] 148.9079456 315.3994283
[428,] 216.2296820 148.9079456
[429,] 315.3994283 216.2296820
[430,] 377.5161438 315.3994283
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 -1019.7563837 -392.1615896
2 -453.1708724 -1019.7563837
3 -804.9260504 -453.1708724
4 360.7594783 -804.9260504
5 -279.0235781 360.7594783
6 -321.9235182 -279.0235781
7 -463.1453074 -321.9235182
8 -588.8512942 -463.1453074
9 -401.8339358 -588.8512942
10 -742.3992506 -401.8339358
11 -764.9346651 -742.3992506
12 -470.5344475 -764.9346651
13 -322.5549404 -470.5344475
14 -429.2427177 -322.5549404
15 -312.1408758 -429.2427177
16 -468.3176302 -312.1408758
17 -484.7509013 -468.3176302
18 -481.9832552 -484.7509013
19 -460.0219261 -481.9832552
20 -470.5344475 -460.0219261
21 -358.3457368 -470.5344475
22 -470.5344475 -358.3457368
23 -1036.0095528 -470.5344475
24 -705.2365682 -1036.0095528
25 -450.7663177 -705.2365682
26 -356.2347827 -450.7663177
27 -550.7845271 -356.2347827
28 -470.5344475 -550.7845271
29 -534.4957791 -470.5344475
30 134.8952139 -534.4957791
31 -470.5344475 134.8952139
32 -538.0085226 -470.5344475
33 -1460.5526455 -538.0085226
34 -221.6013291 -1460.5526455
35 -24.2183041 -221.6013291
36 -468.0152199 -24.2183041
37 -137.0108679 -468.0152199
38 -730.3905299 -137.0108679
39 -660.3075539 -730.3905299
40 -1040.4291487 -660.3075539
41 -773.7211062 -1040.4291487
42 -670.8199921 -773.7211062
43 -470.5344475 -670.8199921
44 -355.1175642 -470.5344475
45 -351.4452806 -355.1175642
46 -470.5344475 -351.4452806
47 -427.9140673 -470.5344475
48 -495.4713487 -427.9140673
49 -790.6830912 -495.4713487
50 -294.7326741 -790.6830912
51 -631.9214648 -294.7326741
52 -165.3021303 -631.9214648
53 -614.0887094 -165.3021303
54 -3205.6775652 -614.0887094
55 -1013.6751239 -3205.6775652
56 -450.9765612 -1013.6751239
57 0.2260205 -450.9765612
58 -314.2046138 0.2260205
59 -470.5344475 -314.2046138
60 -470.5344475 -470.5344475
61 -314.6810109 -470.5344475
62 104.9268068 -314.6810109
63 -472.3085284 104.9268068
64 730.8280805 -472.3085284
65 -465.0419398 730.8280805
66 -403.4219727 -465.0419398
67 -515.6566081 -403.4219727
68 -373.6832878 -515.6566081
69 -470.5344475 -373.6832878
70 -3378.6463520 -470.5344475
71 297.7009682 -3378.6463520
72 8771.2754479 297.7009682
73 -470.5344475 8771.2754479
74 -519.9980963 -470.5344475
75 -530.6796029 -519.9980963
76 -580.5956424 -530.6796029
77 -458.4692428 -580.5956424
78 1065.4420485 -458.4692428
79 -180.9583579 1065.4420485
80 -467.6732617 -180.9583579
81 -470.5344475 -467.6732617
82 -470.5344475 -470.5344475
83 -502.2014647 -470.5344475
84 -470.5344475 -502.2014647
85 -470.5449476 -470.5344475
86 -460.3507269 -470.5449476
87 -470.2362739 -460.3507269
88 -374.2722757 -470.2362739
89 -569.5353756 -374.2722757
90 -470.5449476 -569.5353756
91 -359.2706537 -470.5449476
92 -934.1458012 -359.2706537
93 545.0145513 -934.1458012
94 -505.2476610 545.0145513
95 -546.9441833 -505.2476610
96 -470.5344475 -546.9441833
97 -451.2447467 -470.5344475
98 -477.8018249 -451.2447467
99 -250.8679047 -477.8018249
100 -1031.2981959 -250.8679047
101 -377.9418429 -1031.2981959
102 -476.3206260 -377.9418429
103 -470.3494670 -476.3206260
104 -470.5344475 -470.3494670
105 1977.8810030 -470.5344475
106 -390.3697034 1977.8810030
107 -1731.1824846 -390.3697034
108 -336.8265263 -1731.1824846
109 -458.2571058 -336.8265263
110 -476.6329577 -458.2571058
111 160.5315987 -476.6329577
112 -470.5344475 160.5315987
113 -446.6151630 -470.5344475
114 -470.5449476 -446.6151630
115 -416.7551092 -470.5449476
116 -760.3086408 -416.7551092
117 -1896.0208159 -760.3086408
118 -2010.2224842 -1896.0208159
119 -1339.9715274 -2010.2224842
120 87.0682504 -1339.9715274
121 -470.5344475 87.0682504
122 -2069.8586465 -470.5344475
123 -739.1645505 -2069.8586465
124 19.0167286 -739.1645505
125 -470.5449476 19.0167286
126 -449.8260321 -470.5449476
127 -470.5344475 -449.8260321
128 1713.7209097 -470.5344475
129 -386.0915604 1713.7209097
130 -512.7086149 -386.0915604
131 -470.5449476 -512.7086149
132 -470.5449476 -470.5449476
133 -470.5449476 -470.5449476
134 -459.9319090 -470.5449476
135 -409.3731914 -459.9319090
136 315.3994283 -409.3731914
137 271.9793364 315.3994283
138 315.3994283 271.9793364
139 315.3994283 315.3994283
140 -100.8336478 315.3994283
141 822.4580612 -100.8336478
142 315.3994283 822.4580612
143 314.3994283 315.3994283
144 324.6291894 314.3994283
145 315.3994283 324.6291894
146 325.5178793 315.3994283
147 397.9112373 325.5178793
148 247.2483103 397.9112373
149 314.3994283 247.2483103
150 315.3994283 314.3994283
151 345.5406820 315.3994283
152 315.3994283 345.5406820
153 355.1577545 315.3994283
154 258.5359260 355.1577545
155 162.0830619 258.5359260
156 304.1319653 162.0830619
157 315.3994283 304.1319653
158 297.6329451 315.3994283
159 -327.9527531 297.6329451
160 322.2723763 -327.9527531
161 315.3994283 322.2723763
162 376.0361236 315.3994283
163 269.1539227 376.0361236
164 314.3994283 269.1539227
165 -236.8998305 314.3994283
166 315.3994283 -236.8998305
167 320.9979119 315.3994283
168 661.8501481 320.9979119
169 216.4049847 661.8501481
170 -620.2062952 216.4049847
171 -196.2867470 -620.2062952
172 125.9323451 -196.2867470
173 -73.3090202 125.9323451
174 493.2488721 -73.3090202
175 35.4909132 493.2488721
176 315.3994283 35.4909132
177 317.0125444 315.3994283
178 307.5333072 317.0125444
179 315.3994283 307.5333072
180 -486.2091934 315.3994283
181 -275.9641179 -486.2091934
182 315.3994283 -275.9641179
183 314.3994283 315.3994283
184 315.3994283 314.3994283
185 315.7741848 315.3994283
186 304.8064760 315.7741848
187 273.8043439 304.8064760
188 2420.8634789 273.8043439
189 315.3994283 2420.8634789
190 315.0847766 315.3994283
191 280.6204554 315.0847766
192 315.3994283 280.6204554
193 315.8441910 315.3994283
194 172.6253245 315.8441910
195 315.3994283 172.6253245
196 -288.3499046 315.3994283
197 546.4077924 -288.3499046
198 315.3994283 546.4077924
199 314.3994283 315.3994283
200 242.3681110 314.3994283
201 314.3994283 242.3681110
202 315.3994283 314.3994283
203 321.7409222 315.3994283
204 315.3994283 321.7409222
205 -58.4255405 315.3994283
206 36.9547142 -58.4255405
207 315.3994283 36.9547142
208 -42.4623428 315.3994283
209 411.5226277 -42.4623428
210 -227.2172816 411.5226277
211 155.9266370 -227.2172816
212 304.7732975 155.9266370
213 352.6581697 304.7732975
214 -1025.0363920 352.6581697
215 139.5015577 -1025.0363920
216 302.6968945 139.5015577
217 314.3994283 302.6968945
218 210.8843809 314.3994283
219 -469.1418232 210.8843809
220 314.0518467 -469.1418232
221 315.3994283 314.0518467
222 328.6259273 315.3994283
223 291.0572476 328.6259273
224 351.7487357 291.0572476
225 -58.8305535 351.7487357
226 315.3994283 -58.8305535
227 -2765.8363590 315.3994283
228 159.2645634 -2765.8363590
229 262.0055034 159.2645634
230 -3005.3518062 262.0055034
231 381.8741556 -3005.3518062
232 -210.0222391 381.8741556
233 315.3994283 -210.0222391
234 317.0127834 315.3994283
235 190.4365370 317.0127834
236 389.9566574 190.4365370
237 315.3994283 389.9566574
238 315.3994283 315.3994283
239 296.4621526 315.3994283
240 137.2401429 296.4621526
241 64.8448672 137.2401429
242 315.3994283 64.8448672
243 322.0098201 315.3994283
244 486.7721881 322.0098201
245 692.0842751 486.7721881
246 315.3994283 692.0842751
247 330.7688591 315.3994283
248 315.3994283 330.7688591
249 315.0736409 315.3994283
250 315.3994283 315.0736409
251 74.4691667 315.3994283
252 315.3994283 74.4691667
253 314.3994283 315.3994283
254 315.3994283 314.3994283
255 315.3994283 315.3994283
256 315.3994283 315.3994283
257 315.3994283 315.3994283
258 7.6359586 315.3994283
259 -31.0356061 7.6359586
260 456.7097840 -31.0356061
261 621.1805299 456.7097840
262 217.9851292 621.1805299
263 -453.9265658 217.9851292
264 194.3743963 -453.9265658
265 315.3994283 194.3743963
266 316.3843162 315.3994283
267 315.3994283 316.3843162
268 179.3722050 315.3994283
269 313.5008966 179.3722050
270 315.3994283 313.5008966
271 315.3994283 315.3994283
272 208.2126011 315.3994283
273 312.8224200 208.2126011
274 325.6942428 312.8224200
275 455.9650879 325.6942428
276 315.7686817 455.9650879
277 335.0404225 315.7686817
278 357.0034963 335.0404225
279 325.4853920 357.0034963
280 321.8809852 325.4853920
281 324.1832544 321.8809852
282 315.3994283 324.1832544
283 277.4683358 315.3994283
284 315.7987361 277.4683358
285 334.8665928 315.7987361
286 315.3994283 334.8665928
287 216.7893830 315.3994283
288 740.5007750 216.7893830
289 1520.2646473 740.5007750
290 315.3994283 1520.2646473
291 323.0512807 315.3994283
292 315.3994283 323.0512807
293 278.0031205 315.3994283
294 314.1572864 278.0031205
295 315.1544314 314.1572864
296 314.3994283 315.1544314
297 1238.8720897 314.3994283
298 314.3994283 1238.8720897
299 314.3994283 314.3994283
300 650.5540791 314.3994283
301 247.7804980 650.5540791
302 -300.3015591 247.7804980
303 308.0942909 -300.3015591
304 255.5358441 308.0942909
305 325.8886133 255.5358441
306 -393.2884841 325.8886133
307 315.3994283 -393.2884841
308 315.3994283 315.3994283
309 314.6297502 315.3994283
310 315.3994283 314.6297502
311 50.1942145 315.3994283
312 273.0411048 50.1942145
313 313.4492315 273.0411048
314 217.8344257 313.4492315
315 97.0458816 217.8344257
316 -555.3755167 97.0458816
317 351.7294826 -555.3755167
318 315.3994283 351.7294826
319 209.5370554 315.3994283
320 314.4961618 209.5370554
321 336.5024572 314.4961618
322 -10.4298297 336.5024572
323 315.3994283 -10.4298297
324 -418.8891389 315.3994283
325 315.3994283 -418.8891389
326 -750.7199719 315.3994283
327 -1012.2136382 -750.7199719
328 1177.9550887 -1012.2136382
329 41.7899478 1177.9550887
330 153.5326488 41.7899478
331 314.3994283 153.5326488
332 314.3994283 314.3994283
333 -1482.4629570 314.3994283
334 120.4144274 -1482.4629570
335 315.3994283 120.4144274
336 -110.3046243 315.3994283
337 -266.9372740 -110.3046243
338 314.3994283 -266.9372740
339 315.3994283 314.3994283
340 -255.2058402 315.3994283
341 307.5539870 -255.2058402
342 315.3994283 307.5539870
343 -63.5229823 315.3994283
344 314.3994283 -63.5229823
345 62.0263233 314.3994283
346 306.0083792 62.0263233
347 522.2048761 306.0083792
348 315.3994283 522.2048761
349 325.7692883 315.3994283
350 315.3994283 325.7692883
351 314.3994283 315.3994283
352 315.1602883 314.3994283
353 315.3994283 315.1602883
354 180.4040504 315.3994283
355 315.3994283 180.4040504
356 314.3994283 315.3994283
357 315.3994283 314.3994283
358 314.3994283 315.3994283
359 -673.0071334 314.3994283
360 -931.6266442 -673.0071334
361 615.9307969 -931.6266442
362 382.3329805 615.9307969
363 -45.3781510 382.3329805
364 342.8701291 -45.3781510
365 539.3398741 342.8701291
366 329.3313646 539.3398741
367 321.3694143 329.3313646
368 314.3994283 321.3694143
369 300.6128820 314.3994283
370 -404.8978638 300.6128820
371 355.7087606 -404.8978638
372 288.3572654 355.7087606
373 116.7488226 288.3572654
374 290.8510489 116.7488226
375 -200.6059805 290.8510489
376 -48.7009041 -200.6059805
377 54.3181962 -48.7009041
378 269.8938791 54.3181962
379 315.3994283 269.8938791
380 242.9919123 315.3994283
381 250.9922927 242.9919123
382 111.3867762 250.9922927
383 283.4126263 111.3867762
384 313.5231695 283.4126263
385 351.3505507 313.5231695
386 348.2096572 351.3505507
387 315.3994283 348.2096572
388 314.3994283 315.3994283
389 315.3994283 314.3994283
390 315.3994283 315.3994283
391 315.3994283 315.3994283
392 95.3249937 315.3994283
393 328.4920881 95.3249937
394 314.3994283 328.4920881
395 315.3994283 314.3994283
396 315.3994283 315.3994283
397 315.3994283 315.3994283
398 330.2700782 315.3994283
399 -262.0387481 330.2700782
400 -298.5011399 -262.0387481
401 44.4788159 -298.5011399
402 404.2790165 44.4788159
403 -345.2528458 404.2790165
404 292.0649056 -345.2528458
405 325.1858653 292.0649056
406 317.9770930 325.1858653
407 -56.1874205 317.9770930
408 -180.8788554 -56.1874205
409 -466.4287385 -180.8788554
410 296.6902579 -466.4287385
411 314.3994283 296.6902579
412 315.3994283 314.3994283
413 206.2554596 315.3994283
414 323.9391740 206.2554596
415 331.7373343 323.9391740
416 -1380.2628987 331.7373343
417 60.2731184 -1380.2628987
418 315.4376576 60.2731184
419 315.3994283 315.4376576
420 321.5951448 315.3994283
421 649.9036455 321.5951448
422 -1784.7849797 649.9036455
423 323.8040500 -1784.7849797
424 218.8220148 323.8040500
425 494.7405549 218.8220148
426 315.3994283 494.7405549
427 148.9079456 315.3994283
428 216.2296820 148.9079456
429 315.3994283 216.2296820
430 377.5161438 315.3994283
> 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/wessaorg/rcomp/tmp/7q3wi1293639134.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/wessaorg/rcomp/tmp/81ve31293639134.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/wessaorg/rcomp/tmp/91ve31293639134.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/wessaorg/rcomp/tmp/10c4dn1293639134.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/wessaorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/www/wessaorg/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/wessaorg/rcomp/tmp/11f4bb1293639134.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/wessaorg/rcomp/tmp/128ebw1293639134.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/wessaorg/rcomp/tmp/13fx8q1293639134.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/wessaorg/rcomp/tmp/14popt1293639134.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/wessaorg/rcomp/tmp/154yn21293639134.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/wessaorg/rcomp/tmp/167gl81293639134.tab")
+ }
>
> try(system("convert tmp/1n3yc1293639134.ps tmp/1n3yc1293639134.png",intern=TRUE))
character(0)
> try(system("convert tmp/2n3yc1293639134.ps tmp/2n3yc1293639134.png",intern=TRUE))
character(0)
> try(system("convert tmp/3n3yc1293639134.ps tmp/3n3yc1293639134.png",intern=TRUE))
character(0)
> try(system("convert tmp/4gcff1293639134.ps tmp/4gcff1293639134.png",intern=TRUE))
character(0)
> try(system("convert tmp/5gcff1293639134.ps tmp/5gcff1293639134.png",intern=TRUE))
character(0)
> try(system("convert tmp/6gcff1293639134.ps tmp/6gcff1293639134.png",intern=TRUE))
character(0)
> try(system("convert tmp/7q3wi1293639134.ps tmp/7q3wi1293639134.png",intern=TRUE))
character(0)
> try(system("convert tmp/81ve31293639134.ps tmp/81ve31293639134.png",intern=TRUE))
character(0)
> try(system("convert tmp/91ve31293639134.ps tmp/91ve31293639134.png",intern=TRUE))
character(0)
> try(system("convert tmp/10c4dn1293639134.ps tmp/10c4dn1293639134.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
16.48 0.50 17.27