R version 2.8.0 (2008-10-20)
Copyright (C) 2008 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
Natural language support but running in an English locale
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> x <- array(list(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 = '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
Attaching package: 'zoo'
The following object(s) are masked from package:base :
as.Date.numeric
> n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test
> par1 <- as.numeric(par1)
> x <- t(y)
> k <- length(x[1,])
> n <- length(x[,1])
> x1 <- cbind(x[,par1], x[,1:k!=par1])
> mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1])
> colnames(x1) <- mycolnames #colnames(x)[par1]
> x <- x1
> if (par3 == 'First Differences'){
+ x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep='')))
+ for (i in 1:n-1) {
+ for (j in 1:k) {
+ x2[i,j] <- x[i+1,j] - x[i,j]
+ }
+ }
+ x <- x2
+ }
> if (par2 == 'Include Monthly Dummies'){
+ x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep ='')))
+ for (i in 1:11){
+ x2[seq(i,n,12),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> if (par2 == 'Include Quarterly Dummies'){
+ x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep ='')))
+ for (i in 1:3){
+ x2[seq(i,n,4),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> k <- length(x[1,])
> if (par3 == 'Linear Trend'){
+ x <- cbind(x, c(1:n))
+ colnames(x)[k+1] <- 't'
+ }
> x
CCScore\r Group Costs Trades Dividends Wealth t
1 0 1 1595 17 60720 319369 1
2 39 0 5565 47 94355 493408 2
3 0 0 601 6 60720 319210 3
4 0 1 188 11 77655 381180 4
5 305 1 7146 105 134028 297978 5
6 -9 0 1135 17 62285 290476 6
7 -1 1 450 8 59325 292136 7
8 0 1 34 8 60630 314353 8
9 0 1 133 14 65990 339445 9
10 0 1 119 6 59118 303677 10
11 33 0 2053 53 100423 397144 11
12 32 0 4036 63 100269 424898 12
13 0 1 0 0 60720 315380 13
14 0 1 4655 393 60720 341570 14
15 0 1 131 11 61808 308989 15
16 0 1 1766 73 60438 305959 16
17 0 1 312 14 58598 318690 17
18 0 1 448 40 59781 323361 18
19 0 1 115 13 60945 318903 19
20 0 1 60 10 61124 314049 20
21 0 1 0 0 60720 315380 21
22 0 1 364 35 47705 298466 22
23 0 1 0 0 60720 315380 23
24 -3 0 1442 118 121173 438493 24
25 6 0 1389 9 57530 378049 25
26 0 1 149 5 62041 313332 26
27 0 1 2212 14 60720 319864 27
28 70 0 7489 322 136996 430866 28
29 0 1 0 0 60720 315380 29
30 5 1 402 26 84990 332743 30
31 -5 0 7419 29 63255 286963 31
32 0 1 0 0 60720 315380 32
33 0 1 307 15 58856 331955 33
34 58 0 1134 162 146216 527021 34
35 37 0 7561 87 82425 364304 35
36 9 1 5131 71 86111 292154 36
37 0 1 9 6 60835 314987 37
38 -31 1 2264 24 55174 272713 38
39 2865 1 40949 481 129352 1398893 39
40 0 0 4980 91 113521 409642 40
41 14 1 272 8 112995 429112 41
42 0 1 757 19 45689 382712 42
43 32 1 1424 94 131116 374943 43
44 0 1 0 0 60720 315380 44
45 0 0 378 52 71873 297413 45
46 -4 1 1061 5 70415 304252 46
47 0 1 0 0 60720 315380 47
48 0 1 203 9 61938 309596 48
49 74 0 3689 25 101792 377305 49
50 20 1 567 11 103345 386212 50
51 247 0 27664 282 74996 657954 51
52 11 1 1686 40 68696 368186 52
53 -6 0 1164 17 57320 269753 53
54 284 1 15824 64 63838 585715 54
55 819 1 6575 159 151352 1071292 55
56 652 0 37597 248 37884 992426 56
57 0 1 314 5 60720 315371 57
58 10 0 1932 36 72954 249898 58
59 22 0 5692 72 64239 357312 59
60 0 1 0 0 60720 315380 60
61 0 1 0 0 60720 315380 61
62 16 0 6533 75 126910 364839 62
63 -3 0 2126 138 21001 230621 63
64 0 1 14 1 61016 315877 64
65 10 0 3192 60 59831 125390 65
66 0 1 47 3 60720 314882 66
67 73 0 4891 217 210568 370837 67
68 0 0 531 10 71641 330068 68
69 0 1 2332 4 72680 324385 69
70 0 1 0 0 60720 315380 70
71 1974 0 16275 109 47818 1443586 71
72 -129 1 9252 48 79510 253537 72
73 28328 0 29790 309 81767 4321023 73
74 0 1 0 0 60720 315380 74
75 0 0 2300 29 64873 352108 75
76 0 1 267 32 57640 330059 76
77 0 0 382 21 59365 340968 77
78 0 1 397 28 60646 317736 78
79 238 1 11920 182 31080 209458 79
80 3040 0 54660 302 126942 1491348 80
81 0 1 5 2 60720 314887 81
82 0 1 0 0 60720 315380 82
83 0 1 0 0 60720 315380 83
84 1 0 2647 37 62920 353058 84
85 0 1 0 0 60720 315380 85
86 0 0 0 0 60720 315380 86
87 0 1 94 7 60793 314533 87
88 0 1 3 2 60698 315354 88
89 0 1 454 51 63261 302187 89
90 0 1 227 9 76644 336639 90
91 0 0 0 0 60720 315380 91
92 0 1 206 9 53110 296702 92
93 1910 1 25830 115 245546 1073089 93
94 38 1 1528 38 60326 146494 94
95 0 1 271 2 60720 325249 95
96 0 1 138 6 69817 331420 96
97 0 1 0 0 60720 315380 97
98 0 1 278 26 61404 314922 98
99 0 1 282 30 62452 320016 99
100 0 1 571 80 41477 280398 100
101 17 1 2253 99 63593 452469 101
102 0 1 290 2 58790 301164 102
103 0 1 78 16 62700 317330 103
104 0 1 20 2 60805 315576 104
105 0 1 0 0 60720 315380 105
106 -366 1 18073 180 27284 -7170 106
107 0 1 1866 29 56225 322331 107
108 2964 0 42634 156 157214 1629616 108
109 0 1 249 62 54323 292754 109
110 0 1 422 28 57935 318056 110
111 36 1 2675 79 59017 355178 111
112 -8 1 965 3 73490 204325 112
113 0 1 0 0 60720 315380 113
114 -5 1 621 16 68005 317046 114
115 0 0 0 0 60720 315380 115
116 0 1 365 8 54820 309560 116
117 38 0 3122 34 94670 414462 117
118 601 0 12988 119 82340 857217 118
119 132 0 5336 81 112477 697458 119
120 66 0 3160 108 108094 530670 120
121 0 1 2489 40 67804 238125 121
122 0 1 0 0 60720 315380 122
123 195 0 6984 45 80570 741409 123
124 19 1 2001 52 95551 393343 124
125 -158 1 15236 90 77440 372631 125
126 0 0 0 0 60720 315380 126
127 0 1 530 60 73433 317291 127
128 0 1 0 0 60720 315380 128
129 -70 1 35624 223 157278 306275 129
130 13 0 1383 11 73221 317892 130
131 3 1 875 14 67000 334280 131
132 0 0 0 0 60720 315380 132
133 0 0 0 0 60720 315380 133
134 0 0 0 0 60720 315380 134
135 0 1 72 3 60398 314210 135
136 1 0 265 31 60720 306948 136
137 1 0 0 60720 315380 0 137
138 0 335 33 64175 320398 0 138
139 1 0 0 60720 315380 0 139
140 1 0 0 60720 315380 0 140
141 0 4525 226 93811 501749 101 141
142 1 3045 58 27330 202055 31 142
143 1 0 0 60720 315380 0 143
144 0 0 0 60720 315380 0 144
145 0 638 14 60370 333210 0 145
146 1 0 0 60720 315380 0 146
147 1 607 11 60436 322340 0 147
148 0 1558 63 55637 369448 0 148
149 1 1324 117 67440 291841 0 149
150 0 0 0 60720 315380 0 150
151 1 0 0 60720 315380 0 151
152 1 611 8 59190 296919 -4 152
153 1 0 0 60720 315380 0 153
154 1 923 8 58620 309038 0 154
155 0 661 3 65920 246541 -12 155
156 1 2397 14 74020 289513 6 156
157 1 366 53 61808 344425 0 157
158 1 0 0 60720 315380 0 158
159 0 135 2 62065 314210 0 159
160 1 1659 29 107577 480382 2 160
161 1 316 9 60505 315009 0 161
162 1 0 0 60720 315380 0 162
163 0 309 11 56535 312878 0 163
164 1 49 8 64107 322031 0 164
165 0 0 0 60720 315380 0 165
166 1 4519 49 102129 597793 163 166
167 1 0 0 60720 315380 0 167
168 1 837 69 61262 315688 4 168
169 1 5119 49 39039 378525 31 169
170 1 1280 117 69465 312378 26 170
171 0 2564 22 130140 403560 13 171
172 0 2045 20 97890 510834 10 172
173 1 2234 16 77200 214215 5 173
174 0 975 23 90534 235133 13 174
175 1 1136 32 48522 343613 0 175
176 1 453 21 81125 365959 0 176
177 1 0 0 60720 315380 0 177
178 1 61 14 60720 314551 0 178
179 0 368 17 61686 303230 0 179
180 1 0 0 60720 315380 0 180
181 1 4901 84 121920 469107 67 181
182 1 540 14 103960 354228 22 182
183 1 0 0 60720 315380 0 183
184 0 0 0 60720 315380 0 184
185 1 0 0 60720 315380 0 185
186 1 2 9 60735 315394 0 186
187 1 36 4 61564 312412 0 187
188 1 776 9 64230 333505 -4 188
189 1 84738 588 -26007 223193 4145 189
190 1 0 0 60720 315380 0 190
191 1 3 4 60761 315656 0 191
192 1 529 9 63870 296261 -3 192
193 1 0 0 60720 315380 0 193
194 1 405 7 60845 336425 0 194
195 1 972 34 71642 359335 0 195
196 1 0 0 60720 315380 0 196
197 0 2099 64 106611 308636 20 197
198 0 3437 47 48022 158492 9 198
199 1 0 0 60720 315380 0 199
200 0 0 0 60720 315380 0 200
201 0 22330 84 79801 711969 446 201
202 0 0 0 60720 315380 0 202
203 1 0 0 60720 315380 0 203
204 1 483 27 60830 306268 0 204
205 1 0 0 60720 315380 0 205
206 0 2239 21 88590 442882 5 206
207 1 2949 41 82903 378509 18 207
208 1 0 0 60720 315380 0 208
209 1 365 83 87192 346611 0 209
210 0 2461 57 55792 314289 0 210
211 0 21950 519 114337 856956 745 211
212 1 3294 23 75832 217193 1 212
213 1 141 6 61630 315366 0 213
214 1 572 16 58580 307930 -1 214
215 1 13326 102 165548 702380 436 215
216 0 2284 33 76403 194493 66 216
217 1 10 2 61656 316155 0 217
218 0 0 0 60720 315380 0 218
219 1 1414 198 70184 330546 13 219
220 1 1975 35 118881 394510 30 220
221 1 43 3 60887 312846 0 221
222 1 0 0 60720 315380 0 222
223 1 844 78 60925 296139 0 223
224 1 304 15 62969 295580 0 224
225 1 458 11 58625 297765 0 225
226 1 18562 155 102313 377934 143 226
227 1 0 0 60720 315380 0 227
228 1 7123 109 288170 638830 769 228
229 1 622 73 73007 304376 0 229
230 1 174 1 64820 307424 0 230
231 1 2220 22 301670 644190 154 231
232 0 121 24 56178 295370 0 232
233 1 11819 85 106113 574339 199 233
234 1 0 0 60720 315380 0 234
235 1 125 12 60798 310201 0 235
236 1 1182 13 70694 327007 -5 236
237 1 1503 39 56364 343466 6 237
238 1 0 0 60720 315380 0 238
239 1 0 0 60720 315380 0 239
240 0 30 4 62045 318098 0 240
241 1 3310 33 75230 448243 -6 241
242 0 554 11 79285 325738 30 242
243 1 0 0 60720 315380 0 243
244 1 468 21 60720 312161 0 244
245 1 4917 42 52811 243650 -3 245
246 1 3256 66 35250 407159 -8 246
247 1 0 0 60720 315380 0 247
248 1 125 13 59734 317698 0 248
249 1 0 0 60720 315380 0 249
250 0 22 2 60722 312502 0 250
251 1 0 0 60720 315380 0 251
252 0 514 52 78780 322378 0 252
253 1 0 0 60720 315380 0 253
254 0 0 0 60720 315380 0 254
255 1 0 0 60720 315380 0 255
256 1 0 0 60720 315380 0 256
257 1 0 0 60720 315380 0 257
258 1 0 0 60720 315380 0 258
259 1 10327 91 88577 640273 135 259
260 1 13718 61 96448 345783 104 260
261 1 3748 9 50350 652925 53 261
262 0 14416 361 49857 439798 220 262
263 0 1526 25 69351 278990 -4 263
264 1 666 171 117869 339836 25 264
265 1 2844 67 72683 240897 -15 265
266 1 0 0 60720 315380 0 266
267 1 368 17 61167 297141 0 267
268 1 0 0 60720 315380 0 268
269 1 333 13 70811 331323 0 269
270 1 26 1 60896 313880 0 270
271 1 0 0 60720 315380 0 271
272 1 0 0 60720 315380 0 272
273 1 1303 61 69863 309422 -5 273
274 1 20 3 60938 315245 0 274
275 0 2384 97 61348 405972 54 275
276 0 203 27 50804 300962 0 276
277 1 71 1 60745 316176 0 277
278 1 53 14 59506 302409 0 278
279 1 562 13 58456 283587 -3 279
280 1 622 7 60950 263276 -2 280
281 0 645 5 60720 312075 0 281
282 1 1763 24 61600 308336 0 282
283 1 0 0 60720 315380 0 283
284 1 317 9 63915 298700 0 284
285 1 1 6 60719 315372 0 285
286 1 275 4 59500 318745 0 286
287 1 0 0 60720 315380 0 287
288 1 936 30 67939 408881 9 288
289 0 8568 37 32168 786690 223 289
290 0 11528 70 -14545 -83265 -25 290
291 1 0 0 60720 315380 0 291
292 1 738 10 60720 321376 0 292
293 1 0 0 60720 315380 0 293
294 1 592 7 64270 276898 -4 294
295 1 126 10 60951 315547 0 295
296 1 2 1 60743 315487 0 296
297 0 0 0 60720 315380 0 297
298 0 18014 239 -1710 1405225 943 298
299 0 0 0 60720 315380 0 299
300 0 0 0 60720 315380 0 300
301 0 37704 452 60448 983660 2167 301
302 1 63 4 65688 318574 0 302
303 1 1431 34 106885 310768 30 303
304 1 94 3 61360 312887 0 304
305 1 192 4 65276 312339 0 305
306 1 32 1 59988 314964 0 306
307 1 6869 105 117520 379983 90 307
308 1 0 0 60720 315380 0 308
309 1 0 0 60720 315380 0 309
310 0 1 4 60722 315398 0 310
311 1 0 0 60720 315380 0 311
312 0 2328 94 82732 253588 0 312
313 1 209 13 64016 316647 0 313
314 1 28 2 60890 315688 0 314
315 1 176 28 68136 310670 0 315
316 0 1920 15 79420 165404 -6 316
317 1 87550 458 153198 4111912 17936 317
318 0 520 10 58650 291650 0 318
319 1 0 0 60720 315380 0 319
320 1 1013 2 69770 253468 -14 320
321 1 15 8 60831 315688 0 321
322 1 587 12 59595 325699 0 322
323 1 5371 66 87720 446211 49 323
324 1 0 0 60720 315380 0 324
325 1 1012 47 114768 368078 26 325
326 1 0 0 60720 315380 0 326
327 0 876 59 138971 352850 24 327
328 1 162556 1081 213118 6282154 56289 328
329 0 43556 247 32648 227132 929 329
330 1 3425 12 83620 283910 -3 330
331 1 810 43 74015 236761 26 331
332 0 0 0 60720 315380 0 332
333 0 0 0 60720 315380 0 333
334 0 2365 79 191778 550608 120 334
335 1 1261 33 76114 307528 13 335
336 1 0 0 60720 315380 0 336
337 1 1585 107 93099 355864 2 337
338 1 16189 295 116384 358589 379 338
339 0 0 0 60720 315380 0 339
340 1 0 0 60720 315380 0 340
341 1 10579 92 110309 375195 4 341
342 1 474 20 61977 288985 0 342
343 1 0 0 60720 315380 0 343
344 1 3642 86 90262 458343 22 344
345 0 0 0 60720 315380 0 345
346 1 472 20 80045 269587 0 346
347 1 98 2 61490 315236 0 347
348 1 3999 30 51252 42754 -58 348
349 1 0 0 60720 315380 0 349
350 1 621 48 60720 308256 0 350
351 1 0 0 60720 315380 0 351
352 0 0 0 60720 315380 0 352
353 1 30 4 60798 313164 0 353
354 1 0 0 60720 315380 0 354
355 0 746 23 71561 269661 10 355
356 1 0 0 60720 315380 0 356
357 0 0 0 60720 315380 0 357
358 1 0 0 60720 315380 0 358
359 0 0 0 60720 315380 0 359
360 0 4150 150 134759 518365 66 360
361 1 4658 158 156608 233773 58 361
362 1 814 18 39625 301881 23 362
363 0 1002 15 56750 298568 0 363
364 1 496 17 87390 325479 0 364
365 1 389 14 58990 325506 0 365
366 1 12679 110 48020 984885 1871 366
367 1 400 152 60720 313267 0 367
368 1 53 6 60349 315793 0 368
369 0 0 0 60720 315380 0 369
370 1 5109 97 67038 215362 -31 370
371 1 576 18 113761 314073 12 371
372 1 438 11 58320 298096 -5 372
373 1 165 24 62841 325176 0 373
374 0 4069 145 79804 207393 119 374
375 1 23 5 62555 314806 0 375
376 0 1285 107 99489 341340 49 376
377 0 4677 75 90131 426280 24 377
378 1 24811 505 95350 929118 822 378
379 1 167 2 64245 307322 0 379
380 1 0 0 60720 315380 0 380
381 1 4628 62 69159 387475 85 381
382 0 226 5 65745 291787 0 382
383 1 1765 20 77623 247060 28 383
384 1 460 15 63346 329784 0 384
385 1 36 3 60894 315834 0 385
386 0 989 11 58930 304555 -5 386
387 0 3055 43 60247 376641 14 387
388 1 0 0 60720 315380 0 388
389 0 0 0 60720 315380 0 389
390 1 0 0 60720 315380 0 390
391 1 0 0 60720 315380 0 391
392 1 0 0 60720 315380 0 392
393 1 13253 947 90829 357760 35 393
394 1 24 8 59818 315637 0 394
395 0 0 0 60720 315380 0 395
396 1 0 0 60720 315380 0 396
397 1 0 0 60720 315380 0 397
398 1 0 0 60720 315380 0 398
399 1 131 6 59661 327071 0 399
400 1 514 15 102725 377516 23 400
401 0 3366 80 108479 299243 12 401
402 1 9327 101 87419 387699 44 402
403 1 384 22 54683 309836 0 403
404 1 17821 373 122844 444477 643 404
405 1 397 10 62710 322327 -1 405
406 1 897 5 60720 314913 0 406
407 0 218 21 60720 315553 0 407
408 1 3369 28 87161 688779 217 408
409 1 5702 61 101481 321896 10 409
410 1 8636 855 128294 301607 108 410
411 1 534 7 62620 304485 0 411
412 0 0 0 60720 315380 0 412
413 1 0 0 60720 315380 0 413
414 1 726 3 69980 231861 -5 414
415 1 1380 26 60982 347385 25 415
416 0 180 11 59635 316386 0 416
417 1 7285 115 188873 491303 59 417
418 1 880 74 80791 261216 34 418
419 1 1 2 60727 315388 0 419
420 1 0 0 60720 315380 0 420
421 1 96 5 60379 313729 0 421
422 0 1889 45 37527 358649 14 422
423 1 45187 353 234817 1926517 2854 423
424 1 288 4 60510 296656 0 424
425 1 1270 26 69206 275311 -10 425
426 0 6526 26 55830 -42143 -28 426
427 1 0 0 60720 315380 0 427
428 1 226 12 72835 343929 0 428
429 1 694 8 68060 367655 -2 429
430 1 0 0 60720 315380 0 430
431 1 249 72 56726 313491 0 431
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) Group Costs Trades Dividends Wealth
-1.704e+03 -7.311e-03 -6.814e-02 8.654e-03 5.518e-05 5.647e-03
t
3.382e+00
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-3608.24 -282.40 -18.04 320.28 7407.69
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -1.704e+03 8.394e+01 -20.303 < 2e-16 ***
Group -7.311e-03 5.787e-03 -1.263 0.207
Costs -6.814e-02 7.924e-03 -8.600 < 2e-16 ***
Trades 8.654e-03 1.156e-03 7.483 4.23e-13 ***
Dividends 5.518e-05 1.788e-04 0.309 0.758
Wealth 5.647e-03 1.654e-04 34.149 < 2e-16 ***
t 3.382e+00 3.546e-01 9.539 < 2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 662.3 on 424 degrees of freedom
Multiple R-squared: 0.7769, Adjusted R-squared: 0.7737
F-statistic: 246 on 6 and 424 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,] 2.272059e-05 4.544119e-05 9.999773e-01
[2,] 7.578351e-07 1.515670e-06 9.999992e-01
[3,] 2.085640e-08 4.171279e-08 1.000000e+00
[4,] 4.619261e-10 9.238523e-10 1.000000e+00
[5,] 5.626539e-11 1.125308e-10 1.000000e+00
[6,] 1.397903e-12 2.795805e-12 1.000000e+00
[7,] 3.701185e-14 7.402369e-14 1.000000e+00
[8,] 1.655761e-15 3.311522e-15 1.000000e+00
[9,] 5.977628e-17 1.195526e-16 1.000000e+00
[10,] 1.560074e-18 3.120148e-18 1.000000e+00
[11,] 3.533265e-20 7.066530e-20 1.000000e+00
[12,] 7.835376e-22 1.567075e-21 1.000000e+00
[13,] 5.362653e-23 1.072531e-22 1.000000e+00
[14,] 1.146912e-24 2.293824e-24 1.000000e+00
[15,] 9.885209e-26 1.977042e-25 1.000000e+00
[16,] 7.726178e-26 1.545236e-25 1.000000e+00
[17,] 2.576333e-27 5.152665e-27 1.000000e+00
[18,] 1.556371e-27 3.112743e-27 1.000000e+00
[19,] 5.875372e-28 1.175074e-27 1.000000e+00
[20,] 1.967372e-29 3.934745e-29 1.000000e+00
[21,] 9.187748e-31 1.837550e-30 1.000000e+00
[22,] 9.944343e-30 1.988869e-29 1.000000e+00
[23,] 4.989103e-31 9.978206e-31 1.000000e+00
[24,] 3.170604e-32 6.341208e-32 1.000000e+00
[25,] 2.531664e-33 5.063328e-33 1.000000e+00
[26,] 9.910268e-35 1.982054e-34 1.000000e+00
[27,] 1.048047e-34 2.096094e-34 1.000000e+00
[28,] 6.064536e-36 1.212907e-35 1.000000e+00
[29,] 3.774455e-37 7.548910e-37 1.000000e+00
[30,] 6.168944e-26 1.233789e-25 1.000000e+00
[31,] 8.648461e-27 1.729692e-26 1.000000e+00
[32,] 1.787419e-27 3.574839e-27 1.000000e+00
[33,] 3.008830e-28 6.017659e-28 1.000000e+00
[34,] 3.431045e-29 6.862091e-29 1.000000e+00
[35,] 4.929986e-30 9.859972e-30 1.000000e+00
[36,] 7.298154e-29 1.459631e-28 1.000000e+00
[37,] 8.418735e-30 1.683747e-29 1.000000e+00
[38,] 1.037285e-30 2.074569e-30 1.000000e+00
[39,] 1.274431e-31 2.548861e-31 1.000000e+00
[40,] 1.844527e-32 3.689053e-32 1.000000e+00
[41,] 2.344464e-33 4.688928e-33 1.000000e+00
[42,] 9.692575e-24 1.938515e-23 1.000000e+00
[43,] 2.027655e-24 4.055309e-24 1.000000e+00
[44,] 3.674027e-24 7.348055e-24 1.000000e+00
[45,] 2.344964e-23 4.689928e-23 1.000000e+00
[46,] 1.218621e-21 2.437242e-21 1.000000e+00
[47,] 2.867984e-19 5.735968e-19 1.000000e+00
[48,] 8.993969e-20 1.798794e-19 1.000000e+00
[49,] 1.018241e-19 2.036481e-19 1.000000e+00
[50,] 3.376849e-20 6.753697e-20 1.000000e+00
[51,] 9.972092e-21 1.994418e-20 1.000000e+00
[52,] 2.902458e-21 5.804916e-21 1.000000e+00
[53,] 7.762575e-22 1.552515e-21 1.000000e+00
[54,] 6.043655e-22 1.208731e-21 1.000000e+00
[55,] 1.691754e-22 3.383508e-22 1.000000e+00
[56,] 3.105629e-22 6.211258e-22 1.000000e+00
[57,] 8.688141e-23 1.737628e-22 1.000000e+00
[58,] 2.397275e-23 4.794550e-23 1.000000e+00
[59,] 8.262567e-24 1.652513e-23 1.000000e+00
[60,] 2.242724e-24 4.485449e-24 1.000000e+00
[61,] 6.141088e-25 1.228218e-24 1.000000e+00
[62,] 2.408162e-23 4.816325e-23 1.000000e+00
[63,] 7.502575e-24 1.500515e-23 1.000000e+00
[64,] 1.000000e+00 7.910006e-151 3.955003e-151
[65,] 1.000000e+00 5.659588e-150 2.829794e-150
[66,] 1.000000e+00 2.189305e-149 1.094653e-149
[67,] 1.000000e+00 1.352871e-148 6.764356e-149
[68,] 1.000000e+00 6.197943e-148 3.098972e-148
[69,] 1.000000e+00 4.056313e-147 2.028156e-147
[70,] 1.000000e+00 4.305415e-150 2.152708e-150
[71,] 1.000000e+00 6.168338e-157 3.084169e-157
[72,] 1.000000e+00 4.165760e-156 2.082880e-156
[73,] 1.000000e+00 2.778689e-155 1.389344e-155
[74,] 1.000000e+00 1.832344e-154 9.161719e-155
[75,] 1.000000e+00 6.000646e-154 3.000323e-154
[76,] 1.000000e+00 3.851203e-153 1.925602e-153
[77,] 1.000000e+00 2.322593e-152 1.161297e-152
[78,] 1.000000e+00 1.445592e-151 7.227958e-152
[79,] 1.000000e+00 8.900082e-151 4.450041e-151
[80,] 1.000000e+00 5.462174e-150 2.731087e-150
[81,] 1.000000e+00 2.881241e-149 1.440621e-149
[82,] 1.000000e+00 1.622212e-148 8.111059e-149
[83,] 1.000000e+00 9.936588e-148 4.968294e-148
[84,] 1.000000e+00 5.649209e-157 2.824605e-157
[85,] 1.000000e+00 7.711590e-158 3.855795e-158
[86,] 1.000000e+00 4.172124e-157 2.086062e-157
[87,] 1.000000e+00 2.206284e-156 1.103142e-156
[88,] 1.000000e+00 1.289439e-155 6.447195e-156
[89,] 1.000000e+00 7.208794e-155 3.604397e-155
[90,] 1.000000e+00 3.777758e-154 1.888879e-154
[91,] 1.000000e+00 2.343686e-153 1.171843e-153
[92,] 1.000000e+00 4.053710e-154 2.026855e-154
[93,] 1.000000e+00 2.268261e-153 1.134131e-153
[94,] 1.000000e+00 1.151534e-152 5.757670e-153
[95,] 1.000000e+00 5.810277e-152 2.905138e-152
[96,] 1.000000e+00 2.861782e-151 1.430891e-151
[97,] 1.000000e+00 1.253665e-164 6.268326e-165
[98,] 1.000000e+00 4.733908e-164 2.366954e-164
[99,] 1.000000e+00 0.000000e+00 0.000000e+00
[100,] 1.000000e+00 0.000000e+00 0.000000e+00
[101,] 1.000000e+00 0.000000e+00 0.000000e+00
[102,] 1.000000e+00 6.916919e-323 3.458460e-323
[103,] 1.000000e+00 2.371515e-322 1.185758e-322
[104,] 1.000000e+00 1.324096e-321 6.620480e-322
[105,] 1.000000e+00 5.138283e-321 2.569141e-321
[106,] 1.000000e+00 2.176853e-320 1.088427e-320
[107,] 1.000000e+00 9.061164e-320 4.530582e-320
[108,] 1.000000e+00 7.862857e-319 3.931429e-319
[109,] 1.000000e+00 0.000000e+00 0.000000e+00
[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 1.976263e-323 9.881313e-324
[307,] 1.000000e+00 1.262832e-320 6.314159e-321
[308,] 1.000000e+00 3.978809e-318 1.989405e-318
[309,] 1.000000e+00 1.759635e-315 8.798177e-316
[310,] 1.000000e+00 2.267308e-312 1.133654e-312
[311,] 1.000000e+00 2.591487e-309 1.295744e-309
[312,] 1.000000e+00 3.229788e-306 1.614894e-306
[313,] 1.000000e+00 3.961815e-303 1.980908e-303
[314,] 1.000000e+00 4.419303e-300 2.209651e-300
[315,] 1.000000e+00 5.139431e-297 2.569715e-297
[316,] 1.000000e+00 4.682010e-294 2.341005e-294
[317,] 1.000000e+00 5.205358e-291 2.602679e-291
[318,] 1.000000e+00 2.296224e-288 1.148112e-288
[319,] 1.000000e+00 2.503991e-285 1.251995e-285
[320,] 1.000000e+00 1.984290e-282 9.921452e-283
[321,] 1.000000e+00 2.106762e-279 1.053381e-279
[322,] 1.000000e+00 2.261731e-276 1.130865e-276
[323,] 1.000000e+00 1.205015e-273 6.025074e-274
[324,] 1.000000e+00 5.741619e-271 2.870810e-271
[325,] 1.000000e+00 1.458253e-268 7.291264e-269
[326,] 1.000000e+00 1.819026e-265 9.095128e-266
[327,] 1.000000e+00 2.270729e-262 1.135364e-262
[328,] 1.000000e+00 2.745144e-259 1.372572e-259
[329,] 1.000000e+00 3.018687e-256 1.509343e-256
[330,] 1.000000e+00 1.376102e-253 6.880510e-254
[331,] 1.000000e+00 1.723811e-250 8.619056e-251
[332,] 1.000000e+00 1.689219e-247 8.446094e-248
[333,] 1.000000e+00 1.999757e-244 9.998783e-245
[334,] 1.000000e+00 2.339448e-241 1.169724e-241
[335,] 1.000000e+00 2.625487e-238 1.312744e-238
[336,] 1.000000e+00 1.351028e-235 6.755138e-236
[337,] 1.000000e+00 1.514450e-232 7.572250e-233
[338,] 1.000000e+00 1.740538e-229 8.702692e-230
[339,] 1.000000e+00 1.334936e-226 6.674680e-227
[340,] 1.000000e+00 1.438906e-223 7.194532e-224
[341,] 1.000000e+00 1.522195e-220 7.610974e-221
[342,] 1.000000e+00 1.532065e-217 7.660323e-218
[343,] 1.000000e+00 9.616569e-215 4.808285e-215
[344,] 1.000000e+00 1.002621e-211 5.013105e-212
[345,] 1.000000e+00 1.005556e-208 5.027779e-209
[346,] 1.000000e+00 6.611857e-206 3.305928e-206
[347,] 1.000000e+00 6.835684e-203 3.417842e-203
[348,] 1.000000e+00 4.053978e-200 2.026989e-200
[349,] 1.000000e+00 4.360112e-197 2.180056e-197
[350,] 1.000000e+00 2.399277e-194 1.199639e-194
[351,] 1.000000e+00 5.242699e-192 2.621349e-192
[352,] 1.000000e+00 5.089703e-189 2.544851e-189
[353,] 1.000000e+00 5.827155e-186 2.913578e-186
[354,] 1.000000e+00 2.918213e-183 1.459106e-183
[355,] 1.000000e+00 3.657216e-180 1.828608e-180
[356,] 1.000000e+00 4.531390e-177 2.265695e-177
[357,] 1.000000e+00 6.812222e-174 3.406111e-174
[358,] 1.000000e+00 8.737520e-171 4.368760e-171
[359,] 1.000000e+00 9.808639e-168 4.904320e-168
[360,] 1.000000e+00 4.934027e-165 2.467014e-165
[361,] 1.000000e+00 4.215718e-162 2.107859e-162
[362,] 1.000000e+00 4.571114e-159 2.285557e-159
[363,] 1.000000e+00 4.793523e-156 2.396761e-156
[364,] 1.000000e+00 5.176718e-153 2.588359e-153
[365,] 1.000000e+00 4.213046e-150 2.106523e-150
[366,] 1.000000e+00 4.572068e-147 2.286034e-147
[367,] 1.000000e+00 1.644322e-144 8.221609e-145
[368,] 1.000000e+00 5.221508e-142 2.610754e-142
[369,] 1.000000e+00 7.020048e-139 3.510024e-139
[370,] 1.000000e+00 9.032079e-136 4.516040e-136
[371,] 1.000000e+00 1.115760e-132 5.578801e-133
[372,] 1.000000e+00 1.169964e-129 5.849820e-130
[373,] 1.000000e+00 4.955620e-127 2.477810e-127
[374,] 1.000000e+00 5.873417e-124 2.936709e-124
[375,] 1.000000e+00 7.239121e-121 3.619561e-121
[376,] 1.000000e+00 8.381749e-118 4.190874e-118
[377,] 1.000000e+00 4.429685e-115 2.214842e-115
[378,] 1.000000e+00 1.677709e-112 8.388545e-113
[379,] 1.000000e+00 2.339444e-109 1.169722e-109
[380,] 1.000000e+00 5.796437e-107 2.898219e-107
[381,] 1.000000e+00 9.015827e-104 4.507913e-104
[382,] 1.000000e+00 1.351119e-100 6.755594e-101
[383,] 1.000000e+00 1.938316e-97 9.691582e-98
[384,] 1.000000e+00 3.204998e-94 1.602499e-94
[385,] 1.000000e+00 4.358726e-91 2.179363e-91
[386,] 1.000000e+00 1.146473e-88 5.732367e-89
[387,] 1.000000e+00 1.769303e-85 8.846517e-86
[388,] 1.000000e+00 2.611961e-82 1.305980e-82
[389,] 1.000000e+00 3.650428e-79 1.825214e-79
[390,] 1.000000e+00 4.825812e-76 2.412906e-76
[391,] 1.000000e+00 6.994267e-73 3.497134e-73
[392,] 1.000000e+00 1.382151e-70 6.910756e-71
[393,] 1.000000e+00 1.210544e-67 6.052718e-68
[394,] 1.000000e+00 1.656306e-64 8.281532e-65
[395,] 1.000000e+00 1.024960e-61 5.124800e-62
[396,] 1.000000e+00 1.293015e-58 6.465077e-59
[397,] 1.000000e+00 1.159636e-55 5.798181e-56
[398,] 1.000000e+00 4.194382e-53 2.097191e-53
[399,] 1.000000e+00 6.932129e-50 3.466065e-50
[400,] 1.000000e+00 2.027471e-47 1.013735e-47
[401,] 1.000000e+00 3.356941e-44 1.678471e-44
[402,] 1.000000e+00 3.241516e-41 1.620758e-41
[403,] 1.000000e+00 5.115093e-39 2.557547e-39
[404,] 1.000000e+00 9.622435e-36 4.811218e-36
[405,] 1.000000e+00 2.001148e-32 1.000574e-32
[406,] 1.000000e+00 5.298584e-30 2.649292e-30
[407,] 1.000000e+00 2.053586e-29 1.026793e-29
[408,] 1.000000e+00 1.154518e-25 5.772589e-26
[409,] 1.000000e+00 5.226577e-23 2.613288e-23
[410,] 1.000000e+00 4.068235e-19 2.034118e-19
[411,] 1.000000e+00 2.907728e-15 1.453864e-15
[412,] 1.000000e+00 1.783946e-11 8.919728e-12
> postscript(file="/var/www/html/freestat/rcomp/tmp/1gi6c1293543676.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> plot(x[,1], type='l', main='Actuals and Interpolation', ylab='value of Actuals and Interpolation (dots)', xlab='time or index')
> points(x[,1]-mysum$resid)
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/freestat/rcomp/tmp/299ox1293543676.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> plot(mysum$resid, type='b', pch=19, main='Residuals', ylab='value of Residuals', xlab='time or index')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/freestat/rcomp/tmp/399ox1293543676.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> hist(mysum$resid, main='Residual Histogram', xlab='values of Residuals')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/freestat/rcomp/tmp/4205z1293543676.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> densityplot(~mysum$resid,col='black',main='Residual Density Plot', xlab='values of Residuals')
> dev.off()
null device
1
> postscript(file="/var/www/html/freestat/rcomp/tmp/5205z1293543676.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> qqnorm(mysum$resid, main='Residual Normal Q-Q Plot')
> qqline(mysum$resid)
> grid()
> dev.off()
null device
1
> (myerror <- as.ts(mysum$resid))
Time Series:
Start = 1
End = 431
Frequency = 1
1 2 3 4 5 6
2.659729 -676.106104 -70.852260 -453.286768 788.370749 108.464831
7 8 9 10 11 12
57.279986 -98.978853 -237.655282 -39.564559 -408.654166 -434.714144
13 14 15 16 17 18
-123.943076 38.580355 -85.847247 38.831116 -134.908052 -155.690244
19 20 21 22 23 24
-156.420520 -136.124575 -151.002743 -33.654495 -157.767660 -765.462298
25 26 27 28 29 30
-417.679304 -146.313213 -46.010046 -253.508974 -178.062410 -248.663445
31 32 33 34 35 36
475.785799 -188.209786 -264.297271 -1260.946942 75.634570 286.034949
37 38 39 40 41 42
-202.347834 155.803306 -683.006057 -411.921502 -831.305677 -550.004027
43 44 45 46 47 48
-437.427953 -228.799287 -106.038849 -105.005296 -238.946662 -195.979717
49 50 51 52 53 54
-272.511980 -592.886976 -58.126186 -428.949200 71.760313 -427.634042
55 56 57 58 59 60
-3273.920292 -879.579239 -251.367211 234.273161 -107.284695 -282.918621
61 62 63 64 65 66
-286.301080 -112.113328 328.419908 -298.325974 1000.056573 -297.224513
67 68 69 70 71 72
-223.629406 -357.433368 -205.999392 -316.743205 -3608.238642 525.708299
73 74 75 76 77 78
7407.693803 -330.273039 -384.816886 -401.842120 -458.997753 -330.292979
79 80 81 82 83 84
1301.255990 -232.807091 -350.842919 -357.332706 -360.715165 -395.939886
85 86 87 88 89 90
-367.480081 -370.869851 -363.121352 -377.292306 -276.155443 -489.928471
91 92 93 94 95 96
-387.782143 -271.305059 -1014.350311 697.575507 -438.584956 -486.413928
97 98 99 100 101 102
-408.069582 -390.185113 -422.152817 -181.398286 -1026.221006 -324.855106
103 104 105 106 107 108
-434.308476 -431.512751 -435.129250 2248.711243 -353.996441 -2004.181683
109 110 111 112 113 114
-304.108333 -438.485528 -462.470532 225.337444 -462.188917 -438.203599
115 116 117 118 119 120
-468.961145 -414.343232 -786.662676 -2055.013187 -2148.004728 -1423.818104
121 122 123 124 125 126
115.871148 -492.631043 -2232.352396 -786.666624 52.435141 -506.168188
127 128 129 130 131 132
-485.440219 -512.925793 1885.327634 -427.428132 -567.643527 -526.462938
133 134 135 136 137 138
-529.845396 -533.227855 -525.098109 -473.588737 699.041048 669.180981
139 140 141 142 143 144
692.276132 688.893673 435.779774 1003.370068 678.746298 674.363840
145 146 147 148 149 150
678.644899 668.598923 672.477558 717.521359 619.249950 654.069089
151 152 153 154 155 156
651.686631 667.597937 644.921714 667.355453 601.060755 539.552096
157 158 159 160 161 162
626.661359 628.009422 613.175548 220.745869 622.666697 614.479588
163 164 165 166 167 168
649.459806 578.940853 603.332213 262.481864 597.567296 600.276210
169 170 171 172 173 174
815.529869 529.092780 -2.398903 263.468489 457.636493 327.937844
175 176 177 178 179 180
684.994135 392.497738 563.742712 561.805966 552.137655 553.595337
181 182 183 184 185 186
53.299615 175.277073 543.447962 539.065503 536.683045 533.797905
187 188 189 190 191 192
523.313917 501.470316 1914.953066 519.770753 516.312762 491.299394
193 194 195 196 197 198
509.623377 507.436001 415.340292 499.476002 117.931460 638.533724
199 200 201 202 203 204
489.328627 484.946169 461.023649 478.181252 475.798793 477.338368
205 206 207 208 209 210
469.033877 234.208937 291.072723 458.886501 233.023390 515.704428
211 212 213 214 215 216
145.508295 345.644094 435.539816 462.799772 -391.380967 320.355560
217 218 219 220 221 222
420.511110 424.061917 362.700378 -72.724049 414.128000 411.532084
223 224 225 226 227 228
418.923043 389.642200 424.584371 180.081343 394.619791 -1539.739817
229 230 231 232 233 234
291.655534 350.771398 -1705.310865 419.636964 58.296675 370.942583
235 236 237 238 239 240
368.902510 286.779688 410.553557 357.412749 354.030291 338.523558
241 242 243 244 245 246
240.851150 186.285327 340.500457 342.148254 444.964180 574.047385
247 248 249 250 251 252
326.970623 333.792590 320.205706 316.261876 313.440790 159.687047
253 254 255 256 257 258
306.675873 302.293414 299.910956 296.528498 293.146039 289.763581
259 260 261 262 263 264
108.328391 76.005266 378.446539 491.129629 212.051931 -210.052686
265 266 267 268 269 270
192.114977 262.703913 260.308589 255.938997 167.672446 247.992025
271 272 273 274 275 276
245.791621 242.409163 173.945766 234.115830 244.563256 317.809359
277 278 279 280 281 282
225.823844 234.677255 245.090019 221.270061 216.205860 215.883083
283 284 285 286 287 288
205.202120 178.022298 198.862460 207.709783 191.672287 129.495951
289 290 291 292 293 294
468.887067 943.040286 178.142453 180.506251 171.377536 144.225546
295 296 297 298 299 300
164.207024 161.107985 156.847703 776.246826 150.082786 146.700327
301 302 303 304 305 306
403.022668 98.500580 -249.082522 129.661345 93.205735 134.065168
307 308 309 310 311 312
-314.206314 120.640660 117.258202 113.137321 110.493285 -57.539722
313 314 315 316 317 318
77.549631 99.198777 36.241978 -45.874763 -349.546161 109.521970
319 320 321 322 323 324
83.433618 12.772713 76.345943 87.561699 -127.476932 66.521326
325 326 327 328 329 330
-397.031304 59.756409 -613.567675 -650.774152 626.439418 -124.331880
331 332 333 334 335 336
-59.163814 38.461658 35.079200 -1093.427693 -92.073337 25.931825
337 338 339 340 341 342
-241.015666 -328.595692 14.784449 12.401991 -339.819126 1.044138
343 344 345 346 347 348
2.254616 -232.301619 -5.510301 -167.785667 -17.077865 113.928559
349 350 351 352 353 354
-18.040135 -13.218420 -24.805052 -29.187510 -31.630779 -34.952427
355 356 357 358 359 360
-123.662340 -41.717344 -46.099802 -48.482260 -52.864719 -667.971518
361 362 363 364 365 366
-839.421527 128.331359 -22.763563 -295.344637 -53.949194 87.047247
367 368 369 370 371 372
-65.525766 -78.322752 -86.689303 -94.089263 -546.014883 -70.133881
373 374 375 376 377 378
-115.272492 -223.832900 -121.323121 -430.886276 -335.213622 -238.506941
379 380 381 382 383 384
-148.216416 -122.896345 -165.704622 -170.851401 -261.438943 -155.560397
385 386 387 388 389 390
-140.871809 -120.094988 -121.673798 -149.956013 -154.338471 -156.720930
391 392 393 394 395 396
-160.103388 -163.485846 -268.533814 -161.738671 -174.633222 -177.015680
397 398 399 400 401 402
-180.398138 -183.780597 -177.277225 -552.824104 -577.338159 -357.521288
403 404 405 406 407 408
-143.837737 -596.724102 -221.472446 -203.915608 -212.207444 -441.708296
409 410 411 412 413 414
-528.293298 -687.586616 -239.212288 -232.135015 -234.517473 -307.884414
415 416 417 418 419 420
-233.595658 -234.265473 -1305.991394 -410.845816 -254.729647 -258.194682
421 422 423 424 425 426
-257.492528 -50.842785 -1525.522093 -266.495855 -335.218189 -167.801944
427 428 429 430 431
-281.871891 -389.199636 -349.409295 -292.019266 -254.007770
> postscript(file="/var/www/html/freestat/rcomp/tmp/6205z1293543676.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 2.659729 NA
1 -676.106104 2.659729
2 -70.852260 -676.106104
3 -453.286768 -70.852260
4 788.370749 -453.286768
5 108.464831 788.370749
6 57.279986 108.464831
7 -98.978853 57.279986
8 -237.655282 -98.978853
9 -39.564559 -237.655282
10 -408.654166 -39.564559
11 -434.714144 -408.654166
12 -123.943076 -434.714144
13 38.580355 -123.943076
14 -85.847247 38.580355
15 38.831116 -85.847247
16 -134.908052 38.831116
17 -155.690244 -134.908052
18 -156.420520 -155.690244
19 -136.124575 -156.420520
20 -151.002743 -136.124575
21 -33.654495 -151.002743
22 -157.767660 -33.654495
23 -765.462298 -157.767660
24 -417.679304 -765.462298
25 -146.313213 -417.679304
26 -46.010046 -146.313213
27 -253.508974 -46.010046
28 -178.062410 -253.508974
29 -248.663445 -178.062410
30 475.785799 -248.663445
31 -188.209786 475.785799
32 -264.297271 -188.209786
33 -1260.946942 -264.297271
34 75.634570 -1260.946942
35 286.034949 75.634570
36 -202.347834 286.034949
37 155.803306 -202.347834
38 -683.006057 155.803306
39 -411.921502 -683.006057
40 -831.305677 -411.921502
41 -550.004027 -831.305677
42 -437.427953 -550.004027
43 -228.799287 -437.427953
44 -106.038849 -228.799287
45 -105.005296 -106.038849
46 -238.946662 -105.005296
47 -195.979717 -238.946662
48 -272.511980 -195.979717
49 -592.886976 -272.511980
50 -58.126186 -592.886976
51 -428.949200 -58.126186
52 71.760313 -428.949200
53 -427.634042 71.760313
54 -3273.920292 -427.634042
55 -879.579239 -3273.920292
56 -251.367211 -879.579239
57 234.273161 -251.367211
58 -107.284695 234.273161
59 -282.918621 -107.284695
60 -286.301080 -282.918621
61 -112.113328 -286.301080
62 328.419908 -112.113328
63 -298.325974 328.419908
64 1000.056573 -298.325974
65 -297.224513 1000.056573
66 -223.629406 -297.224513
67 -357.433368 -223.629406
68 -205.999392 -357.433368
69 -316.743205 -205.999392
70 -3608.238642 -316.743205
71 525.708299 -3608.238642
72 7407.693803 525.708299
73 -330.273039 7407.693803
74 -384.816886 -330.273039
75 -401.842120 -384.816886
76 -458.997753 -401.842120
77 -330.292979 -458.997753
78 1301.255990 -330.292979
79 -232.807091 1301.255990
80 -350.842919 -232.807091
81 -357.332706 -350.842919
82 -360.715165 -357.332706
83 -395.939886 -360.715165
84 -367.480081 -395.939886
85 -370.869851 -367.480081
86 -363.121352 -370.869851
87 -377.292306 -363.121352
88 -276.155443 -377.292306
89 -489.928471 -276.155443
90 -387.782143 -489.928471
91 -271.305059 -387.782143
92 -1014.350311 -271.305059
93 697.575507 -1014.350311
94 -438.584956 697.575507
95 -486.413928 -438.584956
96 -408.069582 -486.413928
97 -390.185113 -408.069582
98 -422.152817 -390.185113
99 -181.398286 -422.152817
100 -1026.221006 -181.398286
101 -324.855106 -1026.221006
102 -434.308476 -324.855106
103 -431.512751 -434.308476
104 -435.129250 -431.512751
105 2248.711243 -435.129250
106 -353.996441 2248.711243
107 -2004.181683 -353.996441
108 -304.108333 -2004.181683
109 -438.485528 -304.108333
110 -462.470532 -438.485528
111 225.337444 -462.470532
112 -462.188917 225.337444
113 -438.203599 -462.188917
114 -468.961145 -438.203599
115 -414.343232 -468.961145
116 -786.662676 -414.343232
117 -2055.013187 -786.662676
118 -2148.004728 -2055.013187
119 -1423.818104 -2148.004728
120 115.871148 -1423.818104
121 -492.631043 115.871148
122 -2232.352396 -492.631043
123 -786.666624 -2232.352396
124 52.435141 -786.666624
125 -506.168188 52.435141
126 -485.440219 -506.168188
127 -512.925793 -485.440219
128 1885.327634 -512.925793
129 -427.428132 1885.327634
130 -567.643527 -427.428132
131 -526.462938 -567.643527
132 -529.845396 -526.462938
133 -533.227855 -529.845396
134 -525.098109 -533.227855
135 -473.588737 -525.098109
136 699.041048 -473.588737
137 669.180981 699.041048
138 692.276132 669.180981
139 688.893673 692.276132
140 435.779774 688.893673
141 1003.370068 435.779774
142 678.746298 1003.370068
143 674.363840 678.746298
144 678.644899 674.363840
145 668.598923 678.644899
146 672.477558 668.598923
147 717.521359 672.477558
148 619.249950 717.521359
149 654.069089 619.249950
150 651.686631 654.069089
151 667.597937 651.686631
152 644.921714 667.597937
153 667.355453 644.921714
154 601.060755 667.355453
155 539.552096 601.060755
156 626.661359 539.552096
157 628.009422 626.661359
158 613.175548 628.009422
159 220.745869 613.175548
160 622.666697 220.745869
161 614.479588 622.666697
162 649.459806 614.479588
163 578.940853 649.459806
164 603.332213 578.940853
165 262.481864 603.332213
166 597.567296 262.481864
167 600.276210 597.567296
168 815.529869 600.276210
169 529.092780 815.529869
170 -2.398903 529.092780
171 263.468489 -2.398903
172 457.636493 263.468489
173 327.937844 457.636493
174 684.994135 327.937844
175 392.497738 684.994135
176 563.742712 392.497738
177 561.805966 563.742712
178 552.137655 561.805966
179 553.595337 552.137655
180 53.299615 553.595337
181 175.277073 53.299615
182 543.447962 175.277073
183 539.065503 543.447962
184 536.683045 539.065503
185 533.797905 536.683045
186 523.313917 533.797905
187 501.470316 523.313917
188 1914.953066 501.470316
189 519.770753 1914.953066
190 516.312762 519.770753
191 491.299394 516.312762
192 509.623377 491.299394
193 507.436001 509.623377
194 415.340292 507.436001
195 499.476002 415.340292
196 117.931460 499.476002
197 638.533724 117.931460
198 489.328627 638.533724
199 484.946169 489.328627
200 461.023649 484.946169
201 478.181252 461.023649
202 475.798793 478.181252
203 477.338368 475.798793
204 469.033877 477.338368
205 234.208937 469.033877
206 291.072723 234.208937
207 458.886501 291.072723
208 233.023390 458.886501
209 515.704428 233.023390
210 145.508295 515.704428
211 345.644094 145.508295
212 435.539816 345.644094
213 462.799772 435.539816
214 -391.380967 462.799772
215 320.355560 -391.380967
216 420.511110 320.355560
217 424.061917 420.511110
218 362.700378 424.061917
219 -72.724049 362.700378
220 414.128000 -72.724049
221 411.532084 414.128000
222 418.923043 411.532084
223 389.642200 418.923043
224 424.584371 389.642200
225 180.081343 424.584371
226 394.619791 180.081343
227 -1539.739817 394.619791
228 291.655534 -1539.739817
229 350.771398 291.655534
230 -1705.310865 350.771398
231 419.636964 -1705.310865
232 58.296675 419.636964
233 370.942583 58.296675
234 368.902510 370.942583
235 286.779688 368.902510
236 410.553557 286.779688
237 357.412749 410.553557
238 354.030291 357.412749
239 338.523558 354.030291
240 240.851150 338.523558
241 186.285327 240.851150
242 340.500457 186.285327
243 342.148254 340.500457
244 444.964180 342.148254
245 574.047385 444.964180
246 326.970623 574.047385
247 333.792590 326.970623
248 320.205706 333.792590
249 316.261876 320.205706
250 313.440790 316.261876
251 159.687047 313.440790
252 306.675873 159.687047
253 302.293414 306.675873
254 299.910956 302.293414
255 296.528498 299.910956
256 293.146039 296.528498
257 289.763581 293.146039
258 108.328391 289.763581
259 76.005266 108.328391
260 378.446539 76.005266
261 491.129629 378.446539
262 212.051931 491.129629
263 -210.052686 212.051931
264 192.114977 -210.052686
265 262.703913 192.114977
266 260.308589 262.703913
267 255.938997 260.308589
268 167.672446 255.938997
269 247.992025 167.672446
270 245.791621 247.992025
271 242.409163 245.791621
272 173.945766 242.409163
273 234.115830 173.945766
274 244.563256 234.115830
275 317.809359 244.563256
276 225.823844 317.809359
277 234.677255 225.823844
278 245.090019 234.677255
279 221.270061 245.090019
280 216.205860 221.270061
281 215.883083 216.205860
282 205.202120 215.883083
283 178.022298 205.202120
284 198.862460 178.022298
285 207.709783 198.862460
286 191.672287 207.709783
287 129.495951 191.672287
288 468.887067 129.495951
289 943.040286 468.887067
290 178.142453 943.040286
291 180.506251 178.142453
292 171.377536 180.506251
293 144.225546 171.377536
294 164.207024 144.225546
295 161.107985 164.207024
296 156.847703 161.107985
297 776.246826 156.847703
298 150.082786 776.246826
299 146.700327 150.082786
300 403.022668 146.700327
301 98.500580 403.022668
302 -249.082522 98.500580
303 129.661345 -249.082522
304 93.205735 129.661345
305 134.065168 93.205735
306 -314.206314 134.065168
307 120.640660 -314.206314
308 117.258202 120.640660
309 113.137321 117.258202
310 110.493285 113.137321
311 -57.539722 110.493285
312 77.549631 -57.539722
313 99.198777 77.549631
314 36.241978 99.198777
315 -45.874763 36.241978
316 -349.546161 -45.874763
317 109.521970 -349.546161
318 83.433618 109.521970
319 12.772713 83.433618
320 76.345943 12.772713
321 87.561699 76.345943
322 -127.476932 87.561699
323 66.521326 -127.476932
324 -397.031304 66.521326
325 59.756409 -397.031304
326 -613.567675 59.756409
327 -650.774152 -613.567675
328 626.439418 -650.774152
329 -124.331880 626.439418
330 -59.163814 -124.331880
331 38.461658 -59.163814
332 35.079200 38.461658
333 -1093.427693 35.079200
334 -92.073337 -1093.427693
335 25.931825 -92.073337
336 -241.015666 25.931825
337 -328.595692 -241.015666
338 14.784449 -328.595692
339 12.401991 14.784449
340 -339.819126 12.401991
341 1.044138 -339.819126
342 2.254616 1.044138
343 -232.301619 2.254616
344 -5.510301 -232.301619
345 -167.785667 -5.510301
346 -17.077865 -167.785667
347 113.928559 -17.077865
348 -18.040135 113.928559
349 -13.218420 -18.040135
350 -24.805052 -13.218420
351 -29.187510 -24.805052
352 -31.630779 -29.187510
353 -34.952427 -31.630779
354 -123.662340 -34.952427
355 -41.717344 -123.662340
356 -46.099802 -41.717344
357 -48.482260 -46.099802
358 -52.864719 -48.482260
359 -667.971518 -52.864719
360 -839.421527 -667.971518
361 128.331359 -839.421527
362 -22.763563 128.331359
363 -295.344637 -22.763563
364 -53.949194 -295.344637
365 87.047247 -53.949194
366 -65.525766 87.047247
367 -78.322752 -65.525766
368 -86.689303 -78.322752
369 -94.089263 -86.689303
370 -546.014883 -94.089263
371 -70.133881 -546.014883
372 -115.272492 -70.133881
373 -223.832900 -115.272492
374 -121.323121 -223.832900
375 -430.886276 -121.323121
376 -335.213622 -430.886276
377 -238.506941 -335.213622
378 -148.216416 -238.506941
379 -122.896345 -148.216416
380 -165.704622 -122.896345
381 -170.851401 -165.704622
382 -261.438943 -170.851401
383 -155.560397 -261.438943
384 -140.871809 -155.560397
385 -120.094988 -140.871809
386 -121.673798 -120.094988
387 -149.956013 -121.673798
388 -154.338471 -149.956013
389 -156.720930 -154.338471
390 -160.103388 -156.720930
391 -163.485846 -160.103388
392 -268.533814 -163.485846
393 -161.738671 -268.533814
394 -174.633222 -161.738671
395 -177.015680 -174.633222
396 -180.398138 -177.015680
397 -183.780597 -180.398138
398 -177.277225 -183.780597
399 -552.824104 -177.277225
400 -577.338159 -552.824104
401 -357.521288 -577.338159
402 -143.837737 -357.521288
403 -596.724102 -143.837737
404 -221.472446 -596.724102
405 -203.915608 -221.472446
406 -212.207444 -203.915608
407 -441.708296 -212.207444
408 -528.293298 -441.708296
409 -687.586616 -528.293298
410 -239.212288 -687.586616
411 -232.135015 -239.212288
412 -234.517473 -232.135015
413 -307.884414 -234.517473
414 -233.595658 -307.884414
415 -234.265473 -233.595658
416 -1305.991394 -234.265473
417 -410.845816 -1305.991394
418 -254.729647 -410.845816
419 -258.194682 -254.729647
420 -257.492528 -258.194682
421 -50.842785 -257.492528
422 -1525.522093 -50.842785
423 -266.495855 -1525.522093
424 -335.218189 -266.495855
425 -167.801944 -335.218189
426 -281.871891 -167.801944
427 -389.199636 -281.871891
428 -349.409295 -389.199636
429 -292.019266 -349.409295
430 -254.007770 -292.019266
431 NA -254.007770
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] -676.106104 2.659729
[2,] -70.852260 -676.106104
[3,] -453.286768 -70.852260
[4,] 788.370749 -453.286768
[5,] 108.464831 788.370749
[6,] 57.279986 108.464831
[7,] -98.978853 57.279986
[8,] -237.655282 -98.978853
[9,] -39.564559 -237.655282
[10,] -408.654166 -39.564559
[11,] -434.714144 -408.654166
[12,] -123.943076 -434.714144
[13,] 38.580355 -123.943076
[14,] -85.847247 38.580355
[15,] 38.831116 -85.847247
[16,] -134.908052 38.831116
[17,] -155.690244 -134.908052
[18,] -156.420520 -155.690244
[19,] -136.124575 -156.420520
[20,] -151.002743 -136.124575
[21,] -33.654495 -151.002743
[22,] -157.767660 -33.654495
[23,] -765.462298 -157.767660
[24,] -417.679304 -765.462298
[25,] -146.313213 -417.679304
[26,] -46.010046 -146.313213
[27,] -253.508974 -46.010046
[28,] -178.062410 -253.508974
[29,] -248.663445 -178.062410
[30,] 475.785799 -248.663445
[31,] -188.209786 475.785799
[32,] -264.297271 -188.209786
[33,] -1260.946942 -264.297271
[34,] 75.634570 -1260.946942
[35,] 286.034949 75.634570
[36,] -202.347834 286.034949
[37,] 155.803306 -202.347834
[38,] -683.006057 155.803306
[39,] -411.921502 -683.006057
[40,] -831.305677 -411.921502
[41,] -550.004027 -831.305677
[42,] -437.427953 -550.004027
[43,] -228.799287 -437.427953
[44,] -106.038849 -228.799287
[45,] -105.005296 -106.038849
[46,] -238.946662 -105.005296
[47,] -195.979717 -238.946662
[48,] -272.511980 -195.979717
[49,] -592.886976 -272.511980
[50,] -58.126186 -592.886976
[51,] -428.949200 -58.126186
[52,] 71.760313 -428.949200
[53,] -427.634042 71.760313
[54,] -3273.920292 -427.634042
[55,] -879.579239 -3273.920292
[56,] -251.367211 -879.579239
[57,] 234.273161 -251.367211
[58,] -107.284695 234.273161
[59,] -282.918621 -107.284695
[60,] -286.301080 -282.918621
[61,] -112.113328 -286.301080
[62,] 328.419908 -112.113328
[63,] -298.325974 328.419908
[64,] 1000.056573 -298.325974
[65,] -297.224513 1000.056573
[66,] -223.629406 -297.224513
[67,] -357.433368 -223.629406
[68,] -205.999392 -357.433368
[69,] -316.743205 -205.999392
[70,] -3608.238642 -316.743205
[71,] 525.708299 -3608.238642
[72,] 7407.693803 525.708299
[73,] -330.273039 7407.693803
[74,] -384.816886 -330.273039
[75,] -401.842120 -384.816886
[76,] -458.997753 -401.842120
[77,] -330.292979 -458.997753
[78,] 1301.255990 -330.292979
[79,] -232.807091 1301.255990
[80,] -350.842919 -232.807091
[81,] -357.332706 -350.842919
[82,] -360.715165 -357.332706
[83,] -395.939886 -360.715165
[84,] -367.480081 -395.939886
[85,] -370.869851 -367.480081
[86,] -363.121352 -370.869851
[87,] -377.292306 -363.121352
[88,] -276.155443 -377.292306
[89,] -489.928471 -276.155443
[90,] -387.782143 -489.928471
[91,] -271.305059 -387.782143
[92,] -1014.350311 -271.305059
[93,] 697.575507 -1014.350311
[94,] -438.584956 697.575507
[95,] -486.413928 -438.584956
[96,] -408.069582 -486.413928
[97,] -390.185113 -408.069582
[98,] -422.152817 -390.185113
[99,] -181.398286 -422.152817
[100,] -1026.221006 -181.398286
[101,] -324.855106 -1026.221006
[102,] -434.308476 -324.855106
[103,] -431.512751 -434.308476
[104,] -435.129250 -431.512751
[105,] 2248.711243 -435.129250
[106,] -353.996441 2248.711243
[107,] -2004.181683 -353.996441
[108,] -304.108333 -2004.181683
[109,] -438.485528 -304.108333
[110,] -462.470532 -438.485528
[111,] 225.337444 -462.470532
[112,] -462.188917 225.337444
[113,] -438.203599 -462.188917
[114,] -468.961145 -438.203599
[115,] -414.343232 -468.961145
[116,] -786.662676 -414.343232
[117,] -2055.013187 -786.662676
[118,] -2148.004728 -2055.013187
[119,] -1423.818104 -2148.004728
[120,] 115.871148 -1423.818104
[121,] -492.631043 115.871148
[122,] -2232.352396 -492.631043
[123,] -786.666624 -2232.352396
[124,] 52.435141 -786.666624
[125,] -506.168188 52.435141
[126,] -485.440219 -506.168188
[127,] -512.925793 -485.440219
[128,] 1885.327634 -512.925793
[129,] -427.428132 1885.327634
[130,] -567.643527 -427.428132
[131,] -526.462938 -567.643527
[132,] -529.845396 -526.462938
[133,] -533.227855 -529.845396
[134,] -525.098109 -533.227855
[135,] -473.588737 -525.098109
[136,] 699.041048 -473.588737
[137,] 669.180981 699.041048
[138,] 692.276132 669.180981
[139,] 688.893673 692.276132
[140,] 435.779774 688.893673
[141,] 1003.370068 435.779774
[142,] 678.746298 1003.370068
[143,] 674.363840 678.746298
[144,] 678.644899 674.363840
[145,] 668.598923 678.644899
[146,] 672.477558 668.598923
[147,] 717.521359 672.477558
[148,] 619.249950 717.521359
[149,] 654.069089 619.249950
[150,] 651.686631 654.069089
[151,] 667.597937 651.686631
[152,] 644.921714 667.597937
[153,] 667.355453 644.921714
[154,] 601.060755 667.355453
[155,] 539.552096 601.060755
[156,] 626.661359 539.552096
[157,] 628.009422 626.661359
[158,] 613.175548 628.009422
[159,] 220.745869 613.175548
[160,] 622.666697 220.745869
[161,] 614.479588 622.666697
[162,] 649.459806 614.479588
[163,] 578.940853 649.459806
[164,] 603.332213 578.940853
[165,] 262.481864 603.332213
[166,] 597.567296 262.481864
[167,] 600.276210 597.567296
[168,] 815.529869 600.276210
[169,] 529.092780 815.529869
[170,] -2.398903 529.092780
[171,] 263.468489 -2.398903
[172,] 457.636493 263.468489
[173,] 327.937844 457.636493
[174,] 684.994135 327.937844
[175,] 392.497738 684.994135
[176,] 563.742712 392.497738
[177,] 561.805966 563.742712
[178,] 552.137655 561.805966
[179,] 553.595337 552.137655
[180,] 53.299615 553.595337
[181,] 175.277073 53.299615
[182,] 543.447962 175.277073
[183,] 539.065503 543.447962
[184,] 536.683045 539.065503
[185,] 533.797905 536.683045
[186,] 523.313917 533.797905
[187,] 501.470316 523.313917
[188,] 1914.953066 501.470316
[189,] 519.770753 1914.953066
[190,] 516.312762 519.770753
[191,] 491.299394 516.312762
[192,] 509.623377 491.299394
[193,] 507.436001 509.623377
[194,] 415.340292 507.436001
[195,] 499.476002 415.340292
[196,] 117.931460 499.476002
[197,] 638.533724 117.931460
[198,] 489.328627 638.533724
[199,] 484.946169 489.328627
[200,] 461.023649 484.946169
[201,] 478.181252 461.023649
[202,] 475.798793 478.181252
[203,] 477.338368 475.798793
[204,] 469.033877 477.338368
[205,] 234.208937 469.033877
[206,] 291.072723 234.208937
[207,] 458.886501 291.072723
[208,] 233.023390 458.886501
[209,] 515.704428 233.023390
[210,] 145.508295 515.704428
[211,] 345.644094 145.508295
[212,] 435.539816 345.644094
[213,] 462.799772 435.539816
[214,] -391.380967 462.799772
[215,] 320.355560 -391.380967
[216,] 420.511110 320.355560
[217,] 424.061917 420.511110
[218,] 362.700378 424.061917
[219,] -72.724049 362.700378
[220,] 414.128000 -72.724049
[221,] 411.532084 414.128000
[222,] 418.923043 411.532084
[223,] 389.642200 418.923043
[224,] 424.584371 389.642200
[225,] 180.081343 424.584371
[226,] 394.619791 180.081343
[227,] -1539.739817 394.619791
[228,] 291.655534 -1539.739817
[229,] 350.771398 291.655534
[230,] -1705.310865 350.771398
[231,] 419.636964 -1705.310865
[232,] 58.296675 419.636964
[233,] 370.942583 58.296675
[234,] 368.902510 370.942583
[235,] 286.779688 368.902510
[236,] 410.553557 286.779688
[237,] 357.412749 410.553557
[238,] 354.030291 357.412749
[239,] 338.523558 354.030291
[240,] 240.851150 338.523558
[241,] 186.285327 240.851150
[242,] 340.500457 186.285327
[243,] 342.148254 340.500457
[244,] 444.964180 342.148254
[245,] 574.047385 444.964180
[246,] 326.970623 574.047385
[247,] 333.792590 326.970623
[248,] 320.205706 333.792590
[249,] 316.261876 320.205706
[250,] 313.440790 316.261876
[251,] 159.687047 313.440790
[252,] 306.675873 159.687047
[253,] 302.293414 306.675873
[254,] 299.910956 302.293414
[255,] 296.528498 299.910956
[256,] 293.146039 296.528498
[257,] 289.763581 293.146039
[258,] 108.328391 289.763581
[259,] 76.005266 108.328391
[260,] 378.446539 76.005266
[261,] 491.129629 378.446539
[262,] 212.051931 491.129629
[263,] -210.052686 212.051931
[264,] 192.114977 -210.052686
[265,] 262.703913 192.114977
[266,] 260.308589 262.703913
[267,] 255.938997 260.308589
[268,] 167.672446 255.938997
[269,] 247.992025 167.672446
[270,] 245.791621 247.992025
[271,] 242.409163 245.791621
[272,] 173.945766 242.409163
[273,] 234.115830 173.945766
[274,] 244.563256 234.115830
[275,] 317.809359 244.563256
[276,] 225.823844 317.809359
[277,] 234.677255 225.823844
[278,] 245.090019 234.677255
[279,] 221.270061 245.090019
[280,] 216.205860 221.270061
[281,] 215.883083 216.205860
[282,] 205.202120 215.883083
[283,] 178.022298 205.202120
[284,] 198.862460 178.022298
[285,] 207.709783 198.862460
[286,] 191.672287 207.709783
[287,] 129.495951 191.672287
[288,] 468.887067 129.495951
[289,] 943.040286 468.887067
[290,] 178.142453 943.040286
[291,] 180.506251 178.142453
[292,] 171.377536 180.506251
[293,] 144.225546 171.377536
[294,] 164.207024 144.225546
[295,] 161.107985 164.207024
[296,] 156.847703 161.107985
[297,] 776.246826 156.847703
[298,] 150.082786 776.246826
[299,] 146.700327 150.082786
[300,] 403.022668 146.700327
[301,] 98.500580 403.022668
[302,] -249.082522 98.500580
[303,] 129.661345 -249.082522
[304,] 93.205735 129.661345
[305,] 134.065168 93.205735
[306,] -314.206314 134.065168
[307,] 120.640660 -314.206314
[308,] 117.258202 120.640660
[309,] 113.137321 117.258202
[310,] 110.493285 113.137321
[311,] -57.539722 110.493285
[312,] 77.549631 -57.539722
[313,] 99.198777 77.549631
[314,] 36.241978 99.198777
[315,] -45.874763 36.241978
[316,] -349.546161 -45.874763
[317,] 109.521970 -349.546161
[318,] 83.433618 109.521970
[319,] 12.772713 83.433618
[320,] 76.345943 12.772713
[321,] 87.561699 76.345943
[322,] -127.476932 87.561699
[323,] 66.521326 -127.476932
[324,] -397.031304 66.521326
[325,] 59.756409 -397.031304
[326,] -613.567675 59.756409
[327,] -650.774152 -613.567675
[328,] 626.439418 -650.774152
[329,] -124.331880 626.439418
[330,] -59.163814 -124.331880
[331,] 38.461658 -59.163814
[332,] 35.079200 38.461658
[333,] -1093.427693 35.079200
[334,] -92.073337 -1093.427693
[335,] 25.931825 -92.073337
[336,] -241.015666 25.931825
[337,] -328.595692 -241.015666
[338,] 14.784449 -328.595692
[339,] 12.401991 14.784449
[340,] -339.819126 12.401991
[341,] 1.044138 -339.819126
[342,] 2.254616 1.044138
[343,] -232.301619 2.254616
[344,] -5.510301 -232.301619
[345,] -167.785667 -5.510301
[346,] -17.077865 -167.785667
[347,] 113.928559 -17.077865
[348,] -18.040135 113.928559
[349,] -13.218420 -18.040135
[350,] -24.805052 -13.218420
[351,] -29.187510 -24.805052
[352,] -31.630779 -29.187510
[353,] -34.952427 -31.630779
[354,] -123.662340 -34.952427
[355,] -41.717344 -123.662340
[356,] -46.099802 -41.717344
[357,] -48.482260 -46.099802
[358,] -52.864719 -48.482260
[359,] -667.971518 -52.864719
[360,] -839.421527 -667.971518
[361,] 128.331359 -839.421527
[362,] -22.763563 128.331359
[363,] -295.344637 -22.763563
[364,] -53.949194 -295.344637
[365,] 87.047247 -53.949194
[366,] -65.525766 87.047247
[367,] -78.322752 -65.525766
[368,] -86.689303 -78.322752
[369,] -94.089263 -86.689303
[370,] -546.014883 -94.089263
[371,] -70.133881 -546.014883
[372,] -115.272492 -70.133881
[373,] -223.832900 -115.272492
[374,] -121.323121 -223.832900
[375,] -430.886276 -121.323121
[376,] -335.213622 -430.886276
[377,] -238.506941 -335.213622
[378,] -148.216416 -238.506941
[379,] -122.896345 -148.216416
[380,] -165.704622 -122.896345
[381,] -170.851401 -165.704622
[382,] -261.438943 -170.851401
[383,] -155.560397 -261.438943
[384,] -140.871809 -155.560397
[385,] -120.094988 -140.871809
[386,] -121.673798 -120.094988
[387,] -149.956013 -121.673798
[388,] -154.338471 -149.956013
[389,] -156.720930 -154.338471
[390,] -160.103388 -156.720930
[391,] -163.485846 -160.103388
[392,] -268.533814 -163.485846
[393,] -161.738671 -268.533814
[394,] -174.633222 -161.738671
[395,] -177.015680 -174.633222
[396,] -180.398138 -177.015680
[397,] -183.780597 -180.398138
[398,] -177.277225 -183.780597
[399,] -552.824104 -177.277225
[400,] -577.338159 -552.824104
[401,] -357.521288 -577.338159
[402,] -143.837737 -357.521288
[403,] -596.724102 -143.837737
[404,] -221.472446 -596.724102
[405,] -203.915608 -221.472446
[406,] -212.207444 -203.915608
[407,] -441.708296 -212.207444
[408,] -528.293298 -441.708296
[409,] -687.586616 -528.293298
[410,] -239.212288 -687.586616
[411,] -232.135015 -239.212288
[412,] -234.517473 -232.135015
[413,] -307.884414 -234.517473
[414,] -233.595658 -307.884414
[415,] -234.265473 -233.595658
[416,] -1305.991394 -234.265473
[417,] -410.845816 -1305.991394
[418,] -254.729647 -410.845816
[419,] -258.194682 -254.729647
[420,] -257.492528 -258.194682
[421,] -50.842785 -257.492528
[422,] -1525.522093 -50.842785
[423,] -266.495855 -1525.522093
[424,] -335.218189 -266.495855
[425,] -167.801944 -335.218189
[426,] -281.871891 -167.801944
[427,] -389.199636 -281.871891
[428,] -349.409295 -389.199636
[429,] -292.019266 -349.409295
[430,] -254.007770 -292.019266
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 -676.106104 2.659729
2 -70.852260 -676.106104
3 -453.286768 -70.852260
4 788.370749 -453.286768
5 108.464831 788.370749
6 57.279986 108.464831
7 -98.978853 57.279986
8 -237.655282 -98.978853
9 -39.564559 -237.655282
10 -408.654166 -39.564559
11 -434.714144 -408.654166
12 -123.943076 -434.714144
13 38.580355 -123.943076
14 -85.847247 38.580355
15 38.831116 -85.847247
16 -134.908052 38.831116
17 -155.690244 -134.908052
18 -156.420520 -155.690244
19 -136.124575 -156.420520
20 -151.002743 -136.124575
21 -33.654495 -151.002743
22 -157.767660 -33.654495
23 -765.462298 -157.767660
24 -417.679304 -765.462298
25 -146.313213 -417.679304
26 -46.010046 -146.313213
27 -253.508974 -46.010046
28 -178.062410 -253.508974
29 -248.663445 -178.062410
30 475.785799 -248.663445
31 -188.209786 475.785799
32 -264.297271 -188.209786
33 -1260.946942 -264.297271
34 75.634570 -1260.946942
35 286.034949 75.634570
36 -202.347834 286.034949
37 155.803306 -202.347834
38 -683.006057 155.803306
39 -411.921502 -683.006057
40 -831.305677 -411.921502
41 -550.004027 -831.305677
42 -437.427953 -550.004027
43 -228.799287 -437.427953
44 -106.038849 -228.799287
45 -105.005296 -106.038849
46 -238.946662 -105.005296
47 -195.979717 -238.946662
48 -272.511980 -195.979717
49 -592.886976 -272.511980
50 -58.126186 -592.886976
51 -428.949200 -58.126186
52 71.760313 -428.949200
53 -427.634042 71.760313
54 -3273.920292 -427.634042
55 -879.579239 -3273.920292
56 -251.367211 -879.579239
57 234.273161 -251.367211
58 -107.284695 234.273161
59 -282.918621 -107.284695
60 -286.301080 -282.918621
61 -112.113328 -286.301080
62 328.419908 -112.113328
63 -298.325974 328.419908
64 1000.056573 -298.325974
65 -297.224513 1000.056573
66 -223.629406 -297.224513
67 -357.433368 -223.629406
68 -205.999392 -357.433368
69 -316.743205 -205.999392
70 -3608.238642 -316.743205
71 525.708299 -3608.238642
72 7407.693803 525.708299
73 -330.273039 7407.693803
74 -384.816886 -330.273039
75 -401.842120 -384.816886
76 -458.997753 -401.842120
77 -330.292979 -458.997753
78 1301.255990 -330.292979
79 -232.807091 1301.255990
80 -350.842919 -232.807091
81 -357.332706 -350.842919
82 -360.715165 -357.332706
83 -395.939886 -360.715165
84 -367.480081 -395.939886
85 -370.869851 -367.480081
86 -363.121352 -370.869851
87 -377.292306 -363.121352
88 -276.155443 -377.292306
89 -489.928471 -276.155443
90 -387.782143 -489.928471
91 -271.305059 -387.782143
92 -1014.350311 -271.305059
93 697.575507 -1014.350311
94 -438.584956 697.575507
95 -486.413928 -438.584956
96 -408.069582 -486.413928
97 -390.185113 -408.069582
98 -422.152817 -390.185113
99 -181.398286 -422.152817
100 -1026.221006 -181.398286
101 -324.855106 -1026.221006
102 -434.308476 -324.855106
103 -431.512751 -434.308476
104 -435.129250 -431.512751
105 2248.711243 -435.129250
106 -353.996441 2248.711243
107 -2004.181683 -353.996441
108 -304.108333 -2004.181683
109 -438.485528 -304.108333
110 -462.470532 -438.485528
111 225.337444 -462.470532
112 -462.188917 225.337444
113 -438.203599 -462.188917
114 -468.961145 -438.203599
115 -414.343232 -468.961145
116 -786.662676 -414.343232
117 -2055.013187 -786.662676
118 -2148.004728 -2055.013187
119 -1423.818104 -2148.004728
120 115.871148 -1423.818104
121 -492.631043 115.871148
122 -2232.352396 -492.631043
123 -786.666624 -2232.352396
124 52.435141 -786.666624
125 -506.168188 52.435141
126 -485.440219 -506.168188
127 -512.925793 -485.440219
128 1885.327634 -512.925793
129 -427.428132 1885.327634
130 -567.643527 -427.428132
131 -526.462938 -567.643527
132 -529.845396 -526.462938
133 -533.227855 -529.845396
134 -525.098109 -533.227855
135 -473.588737 -525.098109
136 699.041048 -473.588737
137 669.180981 699.041048
138 692.276132 669.180981
139 688.893673 692.276132
140 435.779774 688.893673
141 1003.370068 435.779774
142 678.746298 1003.370068
143 674.363840 678.746298
144 678.644899 674.363840
145 668.598923 678.644899
146 672.477558 668.598923
147 717.521359 672.477558
148 619.249950 717.521359
149 654.069089 619.249950
150 651.686631 654.069089
151 667.597937 651.686631
152 644.921714 667.597937
153 667.355453 644.921714
154 601.060755 667.355453
155 539.552096 601.060755
156 626.661359 539.552096
157 628.009422 626.661359
158 613.175548 628.009422
159 220.745869 613.175548
160 622.666697 220.745869
161 614.479588 622.666697
162 649.459806 614.479588
163 578.940853 649.459806
164 603.332213 578.940853
165 262.481864 603.332213
166 597.567296 262.481864
167 600.276210 597.567296
168 815.529869 600.276210
169 529.092780 815.529869
170 -2.398903 529.092780
171 263.468489 -2.398903
172 457.636493 263.468489
173 327.937844 457.636493
174 684.994135 327.937844
175 392.497738 684.994135
176 563.742712 392.497738
177 561.805966 563.742712
178 552.137655 561.805966
179 553.595337 552.137655
180 53.299615 553.595337
181 175.277073 53.299615
182 543.447962 175.277073
183 539.065503 543.447962
184 536.683045 539.065503
185 533.797905 536.683045
186 523.313917 533.797905
187 501.470316 523.313917
188 1914.953066 501.470316
189 519.770753 1914.953066
190 516.312762 519.770753
191 491.299394 516.312762
192 509.623377 491.299394
193 507.436001 509.623377
194 415.340292 507.436001
195 499.476002 415.340292
196 117.931460 499.476002
197 638.533724 117.931460
198 489.328627 638.533724
199 484.946169 489.328627
200 461.023649 484.946169
201 478.181252 461.023649
202 475.798793 478.181252
203 477.338368 475.798793
204 469.033877 477.338368
205 234.208937 469.033877
206 291.072723 234.208937
207 458.886501 291.072723
208 233.023390 458.886501
209 515.704428 233.023390
210 145.508295 515.704428
211 345.644094 145.508295
212 435.539816 345.644094
213 462.799772 435.539816
214 -391.380967 462.799772
215 320.355560 -391.380967
216 420.511110 320.355560
217 424.061917 420.511110
218 362.700378 424.061917
219 -72.724049 362.700378
220 414.128000 -72.724049
221 411.532084 414.128000
222 418.923043 411.532084
223 389.642200 418.923043
224 424.584371 389.642200
225 180.081343 424.584371
226 394.619791 180.081343
227 -1539.739817 394.619791
228 291.655534 -1539.739817
229 350.771398 291.655534
230 -1705.310865 350.771398
231 419.636964 -1705.310865
232 58.296675 419.636964
233 370.942583 58.296675
234 368.902510 370.942583
235 286.779688 368.902510
236 410.553557 286.779688
237 357.412749 410.553557
238 354.030291 357.412749
239 338.523558 354.030291
240 240.851150 338.523558
241 186.285327 240.851150
242 340.500457 186.285327
243 342.148254 340.500457
244 444.964180 342.148254
245 574.047385 444.964180
246 326.970623 574.047385
247 333.792590 326.970623
248 320.205706 333.792590
249 316.261876 320.205706
250 313.440790 316.261876
251 159.687047 313.440790
252 306.675873 159.687047
253 302.293414 306.675873
254 299.910956 302.293414
255 296.528498 299.910956
256 293.146039 296.528498
257 289.763581 293.146039
258 108.328391 289.763581
259 76.005266 108.328391
260 378.446539 76.005266
261 491.129629 378.446539
262 212.051931 491.129629
263 -210.052686 212.051931
264 192.114977 -210.052686
265 262.703913 192.114977
266 260.308589 262.703913
267 255.938997 260.308589
268 167.672446 255.938997
269 247.992025 167.672446
270 245.791621 247.992025
271 242.409163 245.791621
272 173.945766 242.409163
273 234.115830 173.945766
274 244.563256 234.115830
275 317.809359 244.563256
276 225.823844 317.809359
277 234.677255 225.823844
278 245.090019 234.677255
279 221.270061 245.090019
280 216.205860 221.270061
281 215.883083 216.205860
282 205.202120 215.883083
283 178.022298 205.202120
284 198.862460 178.022298
285 207.709783 198.862460
286 191.672287 207.709783
287 129.495951 191.672287
288 468.887067 129.495951
289 943.040286 468.887067
290 178.142453 943.040286
291 180.506251 178.142453
292 171.377536 180.506251
293 144.225546 171.377536
294 164.207024 144.225546
295 161.107985 164.207024
296 156.847703 161.107985
297 776.246826 156.847703
298 150.082786 776.246826
299 146.700327 150.082786
300 403.022668 146.700327
301 98.500580 403.022668
302 -249.082522 98.500580
303 129.661345 -249.082522
304 93.205735 129.661345
305 134.065168 93.205735
306 -314.206314 134.065168
307 120.640660 -314.206314
308 117.258202 120.640660
309 113.137321 117.258202
310 110.493285 113.137321
311 -57.539722 110.493285
312 77.549631 -57.539722
313 99.198777 77.549631
314 36.241978 99.198777
315 -45.874763 36.241978
316 -349.546161 -45.874763
317 109.521970 -349.546161
318 83.433618 109.521970
319 12.772713 83.433618
320 76.345943 12.772713
321 87.561699 76.345943
322 -127.476932 87.561699
323 66.521326 -127.476932
324 -397.031304 66.521326
325 59.756409 -397.031304
326 -613.567675 59.756409
327 -650.774152 -613.567675
328 626.439418 -650.774152
329 -124.331880 626.439418
330 -59.163814 -124.331880
331 38.461658 -59.163814
332 35.079200 38.461658
333 -1093.427693 35.079200
334 -92.073337 -1093.427693
335 25.931825 -92.073337
336 -241.015666 25.931825
337 -328.595692 -241.015666
338 14.784449 -328.595692
339 12.401991 14.784449
340 -339.819126 12.401991
341 1.044138 -339.819126
342 2.254616 1.044138
343 -232.301619 2.254616
344 -5.510301 -232.301619
345 -167.785667 -5.510301
346 -17.077865 -167.785667
347 113.928559 -17.077865
348 -18.040135 113.928559
349 -13.218420 -18.040135
350 -24.805052 -13.218420
351 -29.187510 -24.805052
352 -31.630779 -29.187510
353 -34.952427 -31.630779
354 -123.662340 -34.952427
355 -41.717344 -123.662340
356 -46.099802 -41.717344
357 -48.482260 -46.099802
358 -52.864719 -48.482260
359 -667.971518 -52.864719
360 -839.421527 -667.971518
361 128.331359 -839.421527
362 -22.763563 128.331359
363 -295.344637 -22.763563
364 -53.949194 -295.344637
365 87.047247 -53.949194
366 -65.525766 87.047247
367 -78.322752 -65.525766
368 -86.689303 -78.322752
369 -94.089263 -86.689303
370 -546.014883 -94.089263
371 -70.133881 -546.014883
372 -115.272492 -70.133881
373 -223.832900 -115.272492
374 -121.323121 -223.832900
375 -430.886276 -121.323121
376 -335.213622 -430.886276
377 -238.506941 -335.213622
378 -148.216416 -238.506941
379 -122.896345 -148.216416
380 -165.704622 -122.896345
381 -170.851401 -165.704622
382 -261.438943 -170.851401
383 -155.560397 -261.438943
384 -140.871809 -155.560397
385 -120.094988 -140.871809
386 -121.673798 -120.094988
387 -149.956013 -121.673798
388 -154.338471 -149.956013
389 -156.720930 -154.338471
390 -160.103388 -156.720930
391 -163.485846 -160.103388
392 -268.533814 -163.485846
393 -161.738671 -268.533814
394 -174.633222 -161.738671
395 -177.015680 -174.633222
396 -180.398138 -177.015680
397 -183.780597 -180.398138
398 -177.277225 -183.780597
399 -552.824104 -177.277225
400 -577.338159 -552.824104
401 -357.521288 -577.338159
402 -143.837737 -357.521288
403 -596.724102 -143.837737
404 -221.472446 -596.724102
405 -203.915608 -221.472446
406 -212.207444 -203.915608
407 -441.708296 -212.207444
408 -528.293298 -441.708296
409 -687.586616 -528.293298
410 -239.212288 -687.586616
411 -232.135015 -239.212288
412 -234.517473 -232.135015
413 -307.884414 -234.517473
414 -233.595658 -307.884414
415 -234.265473 -233.595658
416 -1305.991394 -234.265473
417 -410.845816 -1305.991394
418 -254.729647 -410.845816
419 -258.194682 -254.729647
420 -257.492528 -258.194682
421 -50.842785 -257.492528
422 -1525.522093 -50.842785
423 -266.495855 -1525.522093
424 -335.218189 -266.495855
425 -167.801944 -335.218189
426 -281.871891 -167.801944
427 -389.199636 -281.871891
428 -349.409295 -389.199636
429 -292.019266 -349.409295
430 -254.007770 -292.019266
> plot(z,main=paste('Residual Lag plot, lowess, and regression line'), ylab='values of Residuals', xlab='lagged values of Residuals')
> lines(lowess(z))
> abline(lm(z))
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/freestat/rcomp/tmp/7ds431293543676.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> acf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Autocorrelation Function')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/freestat/rcomp/tmp/8nj351293543676.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> pacf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Partial Autocorrelation Function')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/freestat/rcomp/tmp/9nj351293543676.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> opar <- par(mfrow = c(2,2), oma = c(0, 0, 1.1, 0))
> plot(mylm, las = 1, sub='Residual Diagnostics')
> par(opar)
> dev.off()
null device
1
> if (n > n25) {
+ postscript(file="/var/www/html/freestat/rcomp/tmp/10nj351293543676.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
+ plot(kp3:nmkm3,gqarr[,2], main='Goldfeld-Quandt test',ylab='2-sided p-value',xlab='breakpoint')
+ grid()
+ dev.off()
+ }
null device
1
>
> #Note: the /var/www/html/freestat/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/www/html/freestat/rcomp/createtable")
>
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Estimated Regression Equation', 1, TRUE)
> a<-table.row.end(a)
> myeq <- colnames(x)[1]
> myeq <- paste(myeq, '[t] = ', sep='')
> for (i in 1:k){
+ if (mysum$coefficients[i,1] > 0) myeq <- paste(myeq, '+', '')
+ myeq <- paste(myeq, mysum$coefficients[i,1], sep=' ')
+ if (rownames(mysum$coefficients)[i] != '(Intercept)') {
+ myeq <- paste(myeq, rownames(mysum$coefficients)[i], sep='')
+ if (rownames(mysum$coefficients)[i] != 't') myeq <- paste(myeq, '[t]', sep='')
+ }
+ }
> myeq <- paste(myeq, ' + e[t]')
> a<-table.row.start(a)
> a<-table.element(a, myeq)
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/www/html/freestat/rcomp/tmp/11r2kb1293543676.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,hyperlink('http://www.xycoon.com/ols1.htm','Multiple Linear Regression - Ordinary Least Squares',''), 6, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Variable',header=TRUE)
> a<-table.element(a,'Parameter',header=TRUE)
> a<-table.element(a,'S.D.',header=TRUE)
> a<-table.element(a,'T-STAT
H0: parameter = 0',header=TRUE)
> a<-table.element(a,'2-tail p-value',header=TRUE)
> a<-table.element(a,'1-tail p-value',header=TRUE)
> a<-table.row.end(a)
> for (i in 1:k){
+ a<-table.row.start(a)
+ a<-table.element(a,rownames(mysum$coefficients)[i],header=TRUE)
+ a<-table.element(a,mysum$coefficients[i,1])
+ a<-table.element(a, round(mysum$coefficients[i,2],6))
+ a<-table.element(a, round(mysum$coefficients[i,3],4))
+ a<-table.element(a, round(mysum$coefficients[i,4],6))
+ a<-table.element(a, round(mysum$coefficients[i,4]/2,6))
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/www/html/freestat/rcomp/tmp/12u20h1293543676.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Regression Statistics', 2, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple R',1,TRUE)
> a<-table.element(a, sqrt(mysum$r.squared))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'R-squared',1,TRUE)
> a<-table.element(a, mysum$r.squared)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Adjusted R-squared',1,TRUE)
> a<-table.element(a, mysum$adj.r.squared)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (value)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[1])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (DF numerator)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[2])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (DF denominator)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[3])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'p-value',1,TRUE)
> a<-table.element(a, 1-pf(mysum$fstatistic[1],mysum$fstatistic[2],mysum$fstatistic[3]))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Residual Statistics', 2, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Residual Standard Deviation',1,TRUE)
> a<-table.element(a, mysum$sigma)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Sum Squared Residuals',1,TRUE)
> a<-table.element(a, sum(myerror*myerror))
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/www/html/freestat/rcomp/tmp/1313fb1293543676.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Actuals, Interpolation, and Residuals', 4, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Time or Index', 1, TRUE)
> a<-table.element(a, 'Actuals', 1, TRUE)
> a<-table.element(a, 'Interpolation
Forecast', 1, TRUE)
> a<-table.element(a, 'Residuals
Prediction Error', 1, TRUE)
> a<-table.row.end(a)
> for (i in 1:n) {
+ a<-table.row.start(a)
+ a<-table.element(a,i, 1, TRUE)
+ a<-table.element(a,x[i])
+ a<-table.element(a,x[i]-mysum$resid[i])
+ a<-table.element(a,mysum$resid[i])
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/www/html/freestat/rcomp/tmp/14tcfe1293543676.tab")
> if (n > n25) {
+ a<-table.start()
+ a<-table.row.start(a)
+ a<-table.element(a,'Goldfeld-Quandt test for Heteroskedasticity',4,TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'p-values',header=TRUE)
+ a<-table.element(a,'Alternative Hypothesis',3,header=TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'breakpoint index',header=TRUE)
+ a<-table.element(a,'greater',header=TRUE)
+ a<-table.element(a,'2-sided',header=TRUE)
+ a<-table.element(a,'less',header=TRUE)
+ a<-table.row.end(a)
+ for (mypoint in kp3:nmkm3) {
+ a<-table.row.start(a)
+ a<-table.element(a,mypoint,header=TRUE)
+ a<-table.element(a,gqarr[mypoint-kp3+1,1])
+ a<-table.element(a,gqarr[mypoint-kp3+1,2])
+ a<-table.element(a,gqarr[mypoint-kp3+1,3])
+ a<-table.row.end(a)
+ }
+ a<-table.end(a)
+ table.save(a,file="/var/www/html/freestat/rcomp/tmp/15p4u41293543676.tab")
+ a<-table.start()
+ a<-table.row.start(a)
+ a<-table.element(a,'Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity',4,TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'Description',header=TRUE)
+ a<-table.element(a,'# significant tests',header=TRUE)
+ a<-table.element(a,'% significant tests',header=TRUE)
+ a<-table.element(a,'OK/NOK',header=TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'1% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant1)
+ a<-table.element(a,numsignificant1/numgqtests)
+ if (numsignificant1/numgqtests < 0.01) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'5% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant5)
+ a<-table.element(a,numsignificant5/numgqtests)
+ if (numsignificant5/numgqtests < 0.05) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'10% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant10)
+ a<-table.element(a,numsignificant10/numgqtests)
+ if (numsignificant10/numgqtests < 0.1) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.end(a)
+ table.save(a,file="/var/www/html/freestat/rcomp/tmp/16bnbs1293543676.tab")
+ }
>
> try(system("convert tmp/1gi6c1293543676.ps tmp/1gi6c1293543676.png",intern=TRUE))
character(0)
> try(system("convert tmp/299ox1293543676.ps tmp/299ox1293543676.png",intern=TRUE))
character(0)
> try(system("convert tmp/399ox1293543676.ps tmp/399ox1293543676.png",intern=TRUE))
character(0)
> try(system("convert tmp/4205z1293543676.ps tmp/4205z1293543676.png",intern=TRUE))
character(0)
> try(system("convert tmp/5205z1293543676.ps tmp/5205z1293543676.png",intern=TRUE))
character(0)
> try(system("convert tmp/6205z1293543676.ps tmp/6205z1293543676.png",intern=TRUE))
character(0)
> try(system("convert tmp/7ds431293543676.ps tmp/7ds431293543676.png",intern=TRUE))
character(0)
> try(system("convert tmp/8nj351293543676.ps tmp/8nj351293543676.png",intern=TRUE))
character(0)
> try(system("convert tmp/9nj351293543676.ps tmp/9nj351293543676.png",intern=TRUE))
character(0)
> try(system("convert tmp/10nj351293543676.ps tmp/10nj351293543676.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
14.836 3.143 15.517