R version 2.12.0 (2010-10-15)
Copyright (C) 2010 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i486-pc-linux-gnu (32-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> x <- array(list(1595
+ ,17
+ ,60720
+ ,319369
+ ,0
+ ,5565
+ ,47
+ ,94355
+ ,493408
+ ,39
+ ,601
+ ,6
+ ,60720
+ ,319210
+ ,0
+ ,188
+ ,11
+ ,77655
+ ,381180
+ ,0
+ ,7146
+ ,105
+ ,134028
+ ,297978
+ ,305
+ ,1135
+ ,17
+ ,62285
+ ,290476
+ ,-9
+ ,450
+ ,8
+ ,59325
+ ,292136
+ ,-1
+ ,34
+ ,8
+ ,60630
+ ,314353
+ ,0
+ ,133
+ ,14
+ ,65990
+ ,339445
+ ,0
+ ,119
+ ,6
+ ,59118
+ ,303677
+ ,0
+ ,2053
+ ,53
+ ,100423
+ ,397144
+ ,33
+ ,4036
+ ,63
+ ,100269
+ ,424898
+ ,32
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,4655
+ ,393
+ ,60720
+ ,341570
+ ,0
+ ,131
+ ,11
+ ,61808
+ ,308989
+ ,0
+ ,1766
+ ,73
+ ,60438
+ ,305959
+ ,0
+ ,312
+ ,14
+ ,58598
+ ,318690
+ ,0
+ ,448
+ ,40
+ ,59781
+ ,323361
+ ,0
+ ,115
+ ,13
+ ,60945
+ ,318903
+ ,0
+ ,60
+ ,10
+ ,61124
+ ,314049
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,364
+ ,35
+ ,47705
+ ,298466
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1442
+ ,118
+ ,121173
+ ,438493
+ ,-3
+ ,1389
+ ,9
+ ,57530
+ ,378049
+ ,6
+ ,149
+ ,5
+ ,62041
+ ,313332
+ ,0
+ ,2212
+ ,14
+ ,60720
+ ,319864
+ ,0
+ ,7489
+ ,322
+ ,136996
+ ,430866
+ ,70
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,402
+ ,26
+ ,84990
+ ,332743
+ ,5
+ ,7419
+ ,29
+ ,63255
+ ,286963
+ ,-5
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,307
+ ,15
+ ,58856
+ ,331955
+ ,0
+ ,1134
+ ,162
+ ,146216
+ ,527021
+ ,58
+ ,7561
+ ,87
+ ,82425
+ ,364304
+ ,37
+ ,5131
+ ,71
+ ,86111
+ ,292154
+ ,9
+ ,9
+ ,6
+ ,60835
+ ,314987
+ ,0
+ ,2264
+ ,24
+ ,55174
+ ,272713
+ ,-31
+ ,40949
+ ,481
+ ,129352
+ ,1398893
+ ,2865
+ ,4980
+ ,91
+ ,113521
+ ,409642
+ ,0
+ ,272
+ ,8
+ ,112995
+ ,429112
+ ,14
+ ,757
+ ,19
+ ,45689
+ ,382712
+ ,0
+ ,1424
+ ,94
+ ,131116
+ ,374943
+ ,32
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,378
+ ,52
+ ,71873
+ ,297413
+ ,0
+ ,1061
+ ,5
+ ,70415
+ ,304252
+ ,-4
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,203
+ ,9
+ ,61938
+ ,309596
+ ,0
+ ,3689
+ ,25
+ ,101792
+ ,377305
+ ,74
+ ,567
+ ,11
+ ,103345
+ ,386212
+ ,20
+ ,27664
+ ,282
+ ,74996
+ ,657954
+ ,247
+ ,1686
+ ,40
+ ,68696
+ ,368186
+ ,11
+ ,1164
+ ,17
+ ,57320
+ ,269753
+ ,-6
+ ,15824
+ ,64
+ ,63838
+ ,585715
+ ,284
+ ,6575
+ ,159
+ ,151352
+ ,1071292
+ ,819
+ ,37597
+ ,248
+ ,37884
+ ,992426
+ ,652
+ ,314
+ ,5
+ ,60720
+ ,315371
+ ,0
+ ,1932
+ ,36
+ ,72954
+ ,249898
+ ,10
+ ,5692
+ ,72
+ ,64239
+ ,357312
+ ,22
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,6533
+ ,75
+ ,126910
+ ,364839
+ ,16
+ ,2126
+ ,138
+ ,21001
+ ,230621
+ ,-3
+ ,14
+ ,1
+ ,61016
+ ,315877
+ ,0
+ ,3192
+ ,60
+ ,59831
+ ,125390
+ ,10
+ ,47
+ ,3
+ ,60720
+ ,314882
+ ,0
+ ,4891
+ ,217
+ ,210568
+ ,370837
+ ,73
+ ,531
+ ,10
+ ,71641
+ ,330068
+ ,0
+ ,2332
+ ,4
+ ,72680
+ ,324385
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,16275
+ ,109
+ ,47818
+ ,1443586
+ ,1974
+ ,9252
+ ,48
+ ,79510
+ ,253537
+ ,-129
+ ,29790
+ ,309
+ ,81767
+ ,4321023
+ ,28328
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,2300
+ ,29
+ ,64873
+ ,352108
+ ,0
+ ,267
+ ,32
+ ,57640
+ ,330059
+ ,0
+ ,382
+ ,21
+ ,59365
+ ,340968
+ ,0
+ ,397
+ ,28
+ ,60646
+ ,317736
+ ,0
+ ,11920
+ ,182
+ ,31080
+ ,209458
+ ,238
+ ,54660
+ ,302
+ ,126942
+ ,1491348
+ ,3040
+ ,5
+ ,2
+ ,60720
+ ,314887
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,2647
+ ,37
+ ,62920
+ ,353058
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,94
+ ,7
+ ,60793
+ ,314533
+ ,0
+ ,3
+ ,2
+ ,60698
+ ,315354
+ ,0
+ ,454
+ ,51
+ ,63261
+ ,302187
+ ,0
+ ,227
+ ,9
+ ,76644
+ ,336639
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,206
+ ,9
+ ,53110
+ ,296702
+ ,0
+ ,25830
+ ,115
+ ,245546
+ ,1073089
+ ,1910
+ ,1528
+ ,38
+ ,60326
+ ,146494
+ ,38
+ ,271
+ ,2
+ ,60720
+ ,325249
+ ,0
+ ,138
+ ,6
+ ,69817
+ ,331420
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,278
+ ,26
+ ,61404
+ ,314922
+ ,0
+ ,282
+ ,30
+ ,62452
+ ,320016
+ ,0
+ ,571
+ ,80
+ ,41477
+ ,280398
+ ,0
+ ,2253
+ ,99
+ ,63593
+ ,452469
+ ,17
+ ,290
+ ,2
+ ,58790
+ ,301164
+ ,0
+ ,78
+ ,16
+ ,62700
+ ,317330
+ ,0
+ ,20
+ ,2
+ ,60805
+ ,315576
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,18073
+ ,180
+ ,27284
+ ,-7170
+ ,-366
+ ,1866
+ ,29
+ ,56225
+ ,322331
+ ,0
+ ,42634
+ ,156
+ ,157214
+ ,1629616
+ ,2964
+ ,249
+ ,62
+ ,54323
+ ,292754
+ ,0
+ ,422
+ ,28
+ ,57935
+ ,318056
+ ,0
+ ,2675
+ ,79
+ ,59017
+ ,355178
+ ,36
+ ,965
+ ,3
+ ,73490
+ ,204325
+ ,-8
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,621
+ ,16
+ ,68005
+ ,317046
+ ,-5
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,365
+ ,8
+ ,54820
+ ,309560
+ ,0
+ ,3122
+ ,34
+ ,94670
+ ,414462
+ ,38
+ ,12988
+ ,119
+ ,82340
+ ,857217
+ ,601
+ ,5336
+ ,81
+ ,112477
+ ,697458
+ ,132
+ ,3160
+ ,108
+ ,108094
+ ,530670
+ ,66
+ ,2489
+ ,40
+ ,67804
+ ,238125
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,6984
+ ,45
+ ,80570
+ ,741409
+ ,195
+ ,2001
+ ,52
+ ,95551
+ ,393343
+ ,19
+ ,15236
+ ,90
+ ,77440
+ ,372631
+ ,-158
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,530
+ ,60
+ ,73433
+ ,317291
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,35624
+ ,223
+ ,157278
+ ,306275
+ ,-70
+ ,1383
+ ,11
+ ,73221
+ ,317892
+ ,13
+ ,875
+ ,14
+ ,67000
+ ,334280
+ ,3
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,72
+ ,3
+ ,60398
+ ,314210
+ ,0
+ ,265
+ ,31
+ ,60720
+ ,306948
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,335
+ ,33
+ ,64175
+ ,320398
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,4525
+ ,226
+ ,93811
+ ,501749
+ ,101
+ ,3045
+ ,58
+ ,27330
+ ,202055
+ ,31
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,638
+ ,14
+ ,60370
+ ,333210
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,607
+ ,11
+ ,60436
+ ,322340
+ ,0
+ ,1558
+ ,63
+ ,55637
+ ,369448
+ ,0
+ ,1324
+ ,117
+ ,67440
+ ,291841
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,611
+ ,8
+ ,59190
+ ,296919
+ ,-4
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,923
+ ,8
+ ,58620
+ ,309038
+ ,0
+ ,661
+ ,3
+ ,65920
+ ,246541
+ ,-12
+ ,2397
+ ,14
+ ,74020
+ ,289513
+ ,6
+ ,366
+ ,53
+ ,61808
+ ,344425
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,135
+ ,2
+ ,62065
+ ,314210
+ ,0
+ ,1659
+ ,29
+ ,107577
+ ,480382
+ ,2
+ ,316
+ ,9
+ ,60505
+ ,315009
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,309
+ ,11
+ ,56535
+ ,312878
+ ,0
+ ,49
+ ,8
+ ,64107
+ ,322031
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,4519
+ ,49
+ ,102129
+ ,597793
+ ,163
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,837
+ ,69
+ ,61262
+ ,315688
+ ,4
+ ,5119
+ ,49
+ ,39039
+ ,378525
+ ,31
+ ,1280
+ ,117
+ ,69465
+ ,312378
+ ,26
+ ,2564
+ ,22
+ ,130140
+ ,403560
+ ,13
+ ,2045
+ ,20
+ ,97890
+ ,510834
+ ,10
+ ,2234
+ ,16
+ ,77200
+ ,214215
+ ,5
+ ,975
+ ,23
+ ,90534
+ ,235133
+ ,13
+ ,1136
+ ,32
+ ,48522
+ ,343613
+ ,0
+ ,453
+ ,21
+ ,81125
+ ,365959
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,61
+ ,14
+ ,60720
+ ,314551
+ ,0
+ ,368
+ ,17
+ ,61686
+ ,303230
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,4901
+ ,84
+ ,121920
+ ,469107
+ ,67
+ ,540
+ ,14
+ ,103960
+ ,354228
+ ,22
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,2
+ ,9
+ ,60735
+ ,315394
+ ,0
+ ,36
+ ,4
+ ,61564
+ ,312412
+ ,0
+ ,776
+ ,9
+ ,64230
+ ,333505
+ ,-4
+ ,84738
+ ,588
+ ,-26007
+ ,223193
+ ,4145
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,3
+ ,4
+ ,60761
+ ,315656
+ ,0
+ ,529
+ ,9
+ ,63870
+ ,296261
+ ,-3
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,405
+ ,7
+ ,60845
+ ,336425
+ ,0
+ ,972
+ ,34
+ ,71642
+ ,359335
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,2099
+ ,64
+ ,106611
+ ,308636
+ ,20
+ ,3437
+ ,47
+ ,48022
+ ,158492
+ ,9
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,22330
+ ,84
+ ,79801
+ ,711969
+ ,446
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,483
+ ,27
+ ,60830
+ ,306268
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,2239
+ ,21
+ ,88590
+ ,442882
+ ,5
+ ,2949
+ ,41
+ ,82903
+ ,378509
+ ,18
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,365
+ ,83
+ ,87192
+ ,346611
+ ,0
+ ,2461
+ ,57
+ ,55792
+ ,314289
+ ,0
+ ,21950
+ ,519
+ ,114337
+ ,856956
+ ,745
+ ,3294
+ ,23
+ ,75832
+ ,217193
+ ,1
+ ,141
+ ,6
+ ,61630
+ ,315366
+ ,0
+ ,572
+ ,16
+ ,58580
+ ,307930
+ ,-1
+ ,13326
+ ,102
+ ,165548
+ ,702380
+ ,436
+ ,2284
+ ,33
+ ,76403
+ ,194493
+ ,66
+ ,10
+ ,2
+ ,61656
+ ,316155
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1414
+ ,198
+ ,70184
+ ,330546
+ ,13
+ ,1975
+ ,35
+ ,118881
+ ,394510
+ ,30
+ ,43
+ ,3
+ ,60887
+ ,312846
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,844
+ ,78
+ ,60925
+ ,296139
+ ,0
+ ,304
+ ,15
+ ,62969
+ ,295580
+ ,0
+ ,458
+ ,11
+ ,58625
+ ,297765
+ ,0
+ ,18562
+ ,155
+ ,102313
+ ,377934
+ ,143
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,7123
+ ,109
+ ,288170
+ ,638830
+ ,769
+ ,622
+ ,73
+ ,73007
+ ,304376
+ ,0
+ ,174
+ ,1
+ ,64820
+ ,307424
+ ,0
+ ,2220
+ ,22
+ ,301670
+ ,644190
+ ,154
+ ,121
+ ,24
+ ,56178
+ ,295370
+ ,0
+ ,11819
+ ,85
+ ,106113
+ ,574339
+ ,199
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,125
+ ,12
+ ,60798
+ ,310201
+ ,0
+ ,1182
+ ,13
+ ,70694
+ ,327007
+ ,-5
+ ,1503
+ ,39
+ ,56364
+ ,343466
+ ,6
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,30
+ ,4
+ ,62045
+ ,318098
+ ,0
+ ,3310
+ ,33
+ ,75230
+ ,448243
+ ,-6
+ ,554
+ ,11
+ ,79285
+ ,325738
+ ,30
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,468
+ ,21
+ ,60720
+ ,312161
+ ,0
+ ,4917
+ ,42
+ ,52811
+ ,243650
+ ,-3
+ ,3256
+ ,66
+ ,35250
+ ,407159
+ ,-8
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,125
+ ,13
+ ,59734
+ ,317698
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,22
+ ,2
+ ,60722
+ ,312502
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,514
+ ,52
+ ,78780
+ ,322378
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,10327
+ ,91
+ ,88577
+ ,640273
+ ,135
+ ,13718
+ ,61
+ ,96448
+ ,345783
+ ,104
+ ,3748
+ ,9
+ ,50350
+ ,652925
+ ,53
+ ,14416
+ ,361
+ ,49857
+ ,439798
+ ,220
+ ,1526
+ ,25
+ ,69351
+ ,278990
+ ,-4
+ ,666
+ ,171
+ ,117869
+ ,339836
+ ,25
+ ,2844
+ ,67
+ ,72683
+ ,240897
+ ,-15
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,368
+ ,17
+ ,61167
+ ,297141
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,333
+ ,13
+ ,70811
+ ,331323
+ ,0
+ ,26
+ ,1
+ ,60896
+ ,313880
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1303
+ ,61
+ ,69863
+ ,309422
+ ,-5
+ ,20
+ ,3
+ ,60938
+ ,315245
+ ,0
+ ,2384
+ ,97
+ ,61348
+ ,405972
+ ,54
+ ,203
+ ,27
+ ,50804
+ ,300962
+ ,0
+ ,71
+ ,1
+ ,60745
+ ,316176
+ ,0
+ ,53
+ ,14
+ ,59506
+ ,302409
+ ,0
+ ,562
+ ,13
+ ,58456
+ ,283587
+ ,-3
+ ,622
+ ,7
+ ,60950
+ ,263276
+ ,-2
+ ,645
+ ,5
+ ,60720
+ ,312075
+ ,0
+ ,1763
+ ,24
+ ,61600
+ ,308336
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,317
+ ,9
+ ,63915
+ ,298700
+ ,0
+ ,1
+ ,6
+ ,60719
+ ,315372
+ ,0
+ ,275
+ ,4
+ ,59500
+ ,318745
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,936
+ ,30
+ ,67939
+ ,408881
+ ,9
+ ,8568
+ ,37
+ ,32168
+ ,786690
+ ,223
+ ,11528
+ ,70
+ ,-14545
+ ,-83265
+ ,-25
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,738
+ ,10
+ ,60720
+ ,321376
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,592
+ ,7
+ ,64270
+ ,276898
+ ,-4
+ ,126
+ ,10
+ ,60951
+ ,315547
+ ,0
+ ,2
+ ,1
+ ,60743
+ ,315487
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,18014
+ ,239
+ ,-1710
+ ,1405225
+ ,943
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,37704
+ ,452
+ ,60448
+ ,983660
+ ,2167
+ ,63
+ ,4
+ ,65688
+ ,318574
+ ,0
+ ,1431
+ ,34
+ ,106885
+ ,310768
+ ,30
+ ,94
+ ,3
+ ,61360
+ ,312887
+ ,0
+ ,192
+ ,4
+ ,65276
+ ,312339
+ ,0
+ ,32
+ ,1
+ ,59988
+ ,314964
+ ,0
+ ,6869
+ ,105
+ ,117520
+ ,379983
+ ,90
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1
+ ,4
+ ,60722
+ ,315398
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,2328
+ ,94
+ ,82732
+ ,253588
+ ,0
+ ,209
+ ,13
+ ,64016
+ ,316647
+ ,0
+ ,28
+ ,2
+ ,60890
+ ,315688
+ ,0
+ ,176
+ ,28
+ ,68136
+ ,310670
+ ,0
+ ,1920
+ ,15
+ ,79420
+ ,165404
+ ,-6
+ ,87550
+ ,458
+ ,153198
+ ,4111912
+ ,17936
+ ,520
+ ,10
+ ,58650
+ ,291650
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1013
+ ,2
+ ,69770
+ ,253468
+ ,-14
+ ,15
+ ,8
+ ,60831
+ ,315688
+ ,0
+ ,587
+ ,12
+ ,59595
+ ,325699
+ ,0
+ ,5371
+ ,66
+ ,87720
+ ,446211
+ ,49
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1012
+ ,47
+ ,114768
+ ,368078
+ ,26
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,876
+ ,59
+ ,138971
+ ,352850
+ ,24
+ ,162556
+ ,1081
+ ,213118
+ ,6282154
+ ,56289
+ ,43556
+ ,247
+ ,32648
+ ,227132
+ ,929
+ ,3425
+ ,12
+ ,83620
+ ,283910
+ ,-3
+ ,810
+ ,43
+ ,74015
+ ,236761
+ ,26
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,2365
+ ,79
+ ,191778
+ ,550608
+ ,120
+ ,1261
+ ,33
+ ,76114
+ ,307528
+ ,13
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1585
+ ,107
+ ,93099
+ ,355864
+ ,2
+ ,16189
+ ,295
+ ,116384
+ ,358589
+ ,379
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,10579
+ ,92
+ ,110309
+ ,375195
+ ,4
+ ,474
+ ,20
+ ,61977
+ ,288985
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,3642
+ ,86
+ ,90262
+ ,458343
+ ,22
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,472
+ ,20
+ ,80045
+ ,269587
+ ,0
+ ,98
+ ,2
+ ,61490
+ ,315236
+ ,0
+ ,3999
+ ,30
+ ,51252
+ ,42754
+ ,-58
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,621
+ ,48
+ ,60720
+ ,308256
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,30
+ ,4
+ ,60798
+ ,313164
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,746
+ ,23
+ ,71561
+ ,269661
+ ,10
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,4150
+ ,150
+ ,134759
+ ,518365
+ ,66
+ ,4658
+ ,158
+ ,156608
+ ,233773
+ ,58
+ ,814
+ ,18
+ ,39625
+ ,301881
+ ,23
+ ,1002
+ ,15
+ ,56750
+ ,298568
+ ,0
+ ,496
+ ,17
+ ,87390
+ ,325479
+ ,0
+ ,389
+ ,14
+ ,58990
+ ,325506
+ ,0
+ ,12679
+ ,110
+ ,48020
+ ,984885
+ ,1871
+ ,400
+ ,152
+ ,60720
+ ,313267
+ ,0
+ ,53
+ ,6
+ ,60349
+ ,315793
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,5109
+ ,97
+ ,67038
+ ,215362
+ ,-31
+ ,576
+ ,18
+ ,113761
+ ,314073
+ ,12
+ ,438
+ ,11
+ ,58320
+ ,298096
+ ,-5
+ ,165
+ ,24
+ ,62841
+ ,325176
+ ,0
+ ,4069
+ ,145
+ ,79804
+ ,207393
+ ,119
+ ,23
+ ,5
+ ,62555
+ ,314806
+ ,0
+ ,1285
+ ,107
+ ,99489
+ ,341340
+ ,49
+ ,4677
+ ,75
+ ,90131
+ ,426280
+ ,24
+ ,24811
+ ,505
+ ,95350
+ ,929118
+ ,822
+ ,167
+ ,2
+ ,64245
+ ,307322
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,4628
+ ,62
+ ,69159
+ ,387475
+ ,85
+ ,226
+ ,5
+ ,65745
+ ,291787
+ ,0
+ ,1765
+ ,20
+ ,77623
+ ,247060
+ ,28
+ ,460
+ ,15
+ ,63346
+ ,329784
+ ,0
+ ,36
+ ,3
+ ,60894
+ ,315834
+ ,0
+ ,989
+ ,11
+ ,58930
+ ,304555
+ ,-5
+ ,3055
+ ,43
+ ,60247
+ ,376641
+ ,14
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,13253
+ ,947
+ ,90829
+ ,357760
+ ,35
+ ,24
+ ,8
+ ,59818
+ ,315637
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,131
+ ,6
+ ,59661
+ ,327071
+ ,0
+ ,514
+ ,15
+ ,102725
+ ,377516
+ ,23
+ ,3366
+ ,80
+ ,108479
+ ,299243
+ ,12
+ ,9327
+ ,101
+ ,87419
+ ,387699
+ ,44
+ ,384
+ ,22
+ ,54683
+ ,309836
+ ,0
+ ,17821
+ ,373
+ ,122844
+ ,444477
+ ,643
+ ,397
+ ,10
+ ,62710
+ ,322327
+ ,-1
+ ,897
+ ,5
+ ,60720
+ ,314913
+ ,0
+ ,218
+ ,21
+ ,60720
+ ,315553
+ ,0
+ ,3369
+ ,28
+ ,87161
+ ,688779
+ ,217
+ ,5702
+ ,61
+ ,101481
+ ,321896
+ ,10
+ ,8636
+ ,855
+ ,128294
+ ,301607
+ ,108
+ ,534
+ ,7
+ ,62620
+ ,304485
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,726
+ ,3
+ ,69980
+ ,231861
+ ,-5
+ ,1380
+ ,26
+ ,60982
+ ,347385
+ ,25
+ ,180
+ ,11
+ ,59635
+ ,316386
+ ,0
+ ,7285
+ ,115
+ ,188873
+ ,491303
+ ,59
+ ,880
+ ,74
+ ,80791
+ ,261216
+ ,34
+ ,1
+ ,2
+ ,60727
+ ,315388
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,96
+ ,5
+ ,60379
+ ,313729
+ ,0
+ ,1889
+ ,45
+ ,37527
+ ,358649
+ ,14
+ ,45187
+ ,353
+ ,234817
+ ,1926517
+ ,2854
+ ,288
+ ,4
+ ,60510
+ ,296656
+ ,0
+ ,1270
+ ,26
+ ,69206
+ ,275311
+ ,-10
+ ,6526
+ ,26
+ ,55830
+ ,-42143
+ ,-28
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,226
+ ,12
+ ,72835
+ ,343929
+ ,0
+ ,694
+ ,8
+ ,68060
+ ,367655
+ ,-2
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,249
+ ,72
+ ,56726
+ ,313491
+ ,0)
+ ,dim=c(5
+ ,431)
+ ,dimnames=list(c('Costs'
+ ,'Trades'
+ ,'Dividends'
+ ,'Wealth'
+ ,'CCScore
')
+ ,1:431))
> y <- array(NA,dim=c(5,431),dimnames=list(c('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])
+ }
+ }
> par3 = 'No Linear Trend'
> par2 = 'Do not include Seasonal Dummies'
> par1 = '5'
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Prof. Dr. P. Wessa
> #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #Technical description: Write here your technical program description (don't use hard returns!)
> library(lattice)
> library(lmtest)
Loading required package: zoo
> 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 Costs Trades Dividends Wealth
1 0 1595 17 60720 319369
2 39 5565 47 94355 493408
3 0 601 6 60720 319210
4 0 188 11 77655 381180
5 305 7146 105 134028 297978
6 -9 1135 17 62285 290476
7 -1 450 8 59325 292136
8 0 34 8 60630 314353
9 0 133 14 65990 339445
10 0 119 6 59118 303677
11 33 2053 53 100423 397144
12 32 4036 63 100269 424898
13 0 0 0 60720 315380
14 0 4655 393 60720 341570
15 0 131 11 61808 308989
16 0 1766 73 60438 305959
17 0 312 14 58598 318690
18 0 448 40 59781 323361
19 0 115 13 60945 318903
20 0 60 10 61124 314049
21 0 0 0 60720 315380
22 0 364 35 47705 298466
23 0 0 0 60720 315380
24 -3 1442 118 121173 438493
25 6 1389 9 57530 378049
26 0 149 5 62041 313332
27 0 2212 14 60720 319864
28 70 7489 322 136996 430866
29 0 0 0 60720 315380
30 5 402 26 84990 332743
31 -5 7419 29 63255 286963
32 0 0 0 60720 315380
33 0 307 15 58856 331955
34 58 1134 162 146216 527021
35 37 7561 87 82425 364304
36 9 5131 71 86111 292154
37 0 9 6 60835 314987
38 -31 2264 24 55174 272713
39 2865 40949 481 129352 1398893
40 0 4980 91 113521 409642
41 14 272 8 112995 429112
42 0 757 19 45689 382712
43 32 1424 94 131116 374943
44 0 0 0 60720 315380
45 0 378 52 71873 297413
46 -4 1061 5 70415 304252
47 0 0 0 60720 315380
48 0 203 9 61938 309596
49 74 3689 25 101792 377305
50 20 567 11 103345 386212
51 247 27664 282 74996 657954
52 11 1686 40 68696 368186
53 -6 1164 17 57320 269753
54 284 15824 64 63838 585715
55 819 6575 159 151352 1071292
56 652 37597 248 37884 992426
57 0 314 5 60720 315371
58 10 1932 36 72954 249898
59 22 5692 72 64239 357312
60 0 0 0 60720 315380
61 0 0 0 60720 315380
62 16 6533 75 126910 364839
63 -3 2126 138 21001 230621
64 0 14 1 61016 315877
65 10 3192 60 59831 125390
66 0 47 3 60720 314882
67 73 4891 217 210568 370837
68 0 531 10 71641 330068
69 0 2332 4 72680 324385
70 0 0 0 60720 315380
71 1974 16275 109 47818 1443586
72 -129 9252 48 79510 253537
73 28328 29790 309 81767 4321023
74 0 0 0 60720 315380
75 0 2300 29 64873 352108
76 0 267 32 57640 330059
77 0 382 21 59365 340968
78 0 397 28 60646 317736
79 238 11920 182 31080 209458
80 3040 54660 302 126942 1491348
81 0 5 2 60720 314887
82 0 0 0 60720 315380
83 0 0 0 60720 315380
84 1 2647 37 62920 353058
85 0 0 0 60720 315380
86 0 0 0 60720 315380
87 0 94 7 60793 314533
88 0 3 2 60698 315354
89 0 454 51 63261 302187
90 0 227 9 76644 336639
91 0 0 0 60720 315380
92 0 206 9 53110 296702
93 1910 25830 115 245546 1073089
94 38 1528 38 60326 146494
95 0 271 2 60720 325249
96 0 138 6 69817 331420
97 0 0 0 60720 315380
98 0 278 26 61404 314922
99 0 282 30 62452 320016
100 0 571 80 41477 280398
101 17 2253 99 63593 452469
102 0 290 2 58790 301164
103 0 78 16 62700 317330
104 0 20 2 60805 315576
105 0 0 0 60720 315380
106 -366 18073 180 27284 -7170
107 0 1866 29 56225 322331
108 2964 42634 156 157214 1629616
109 0 249 62 54323 292754
110 0 422 28 57935 318056
111 36 2675 79 59017 355178
112 -8 965 3 73490 204325
113 0 0 0 60720 315380
114 -5 621 16 68005 317046
115 0 0 0 60720 315380
116 0 365 8 54820 309560
117 38 3122 34 94670 414462
118 601 12988 119 82340 857217
119 132 5336 81 112477 697458
120 66 3160 108 108094 530670
121 0 2489 40 67804 238125
122 0 0 0 60720 315380
123 195 6984 45 80570 741409
124 19 2001 52 95551 393343
125 -158 15236 90 77440 372631
126 0 0 0 60720 315380
127 0 530 60 73433 317291
128 0 0 0 60720 315380
129 -70 35624 223 157278 306275
130 13 1383 11 73221 317892
131 3 875 14 67000 334280
132 0 0 0 60720 315380
133 0 0 0 60720 315380
134 0 0 0 60720 315380
135 0 72 3 60398 314210
136 0 265 31 60720 306948
137 335 0 60720 315380 0
138 0 33 64175 320398 0
139 0 0 60720 315380 0
140 4525 0 60720 315380 0
141 3045 226 93811 501749 101
142 0 58 27330 202055 31
143 0 0 60720 315380 0
144 638 0 60720 315380 0
145 0 14 60370 333210 0
146 607 0 60720 315380 0
147 1558 11 60436 322340 0
148 1324 63 55637 369448 0
149 0 117 67440 291841 0
150 0 0 60720 315380 0
151 611 0 60720 315380 0
152 0 8 59190 296919 -4
153 923 0 60720 315380 0
154 661 8 58620 309038 0
155 2397 3 65920 246541 -12
156 366 14 74020 289513 6
157 0 53 61808 344425 0
158 135 0 60720 315380 0
159 1659 2 62065 314210 0
160 316 29 107577 480382 2
161 0 9 60505 315009 0
162 309 0 60720 315380 0
163 49 11 56535 312878 0
164 0 8 64107 322031 0
165 4519 0 60720 315380 0
166 0 49 102129 597793 163
167 837 0 60720 315380 0
168 5119 69 61262 315688 4
169 1280 49 39039 378525 31
170 2564 117 69465 312378 26
171 2045 22 130140 403560 13
172 2234 20 97890 510834 10
173 975 16 77200 214215 5
174 1136 23 90534 235133 13
175 453 32 48522 343613 0
176 0 21 81125 365959 0
177 61 0 60720 315380 0
178 368 14 60720 314551 0
179 0 17 61686 303230 0
180 4901 0 60720 315380 0
181 540 84 121920 469107 67
182 0 14 103960 354228 22
183 0 0 60720 315380 0
184 0 0 60720 315380 0
185 2 0 60720 315380 0
186 36 9 60735 315394 0
187 776 4 61564 312412 0
188 84738 9 64230 333505 -4
189 0 588 -26007 223193 4145
190 3 0 60720 315380 0
191 529 4 60761 315656 0
192 0 9 63870 296261 -3
193 405 0 60720 315380 0
194 972 7 60845 336425 0
195 0 34 71642 359335 0
196 2099 0 60720 315380 0
197 3437 64 106611 308636 20
198 0 47 48022 158492 9
199 0 0 60720 315380 0
200 22330 0 60720 315380 0
201 0 84 79801 711969 446
202 0 0 60720 315380 0
203 483 0 60720 315380 0
204 0 27 60830 306268 0
205 2239 0 60720 315380 0
206 2949 21 88590 442882 5
207 0 41 82903 378509 18
208 365 0 60720 315380 0
209 2461 83 87192 346611 0
210 21950 57 55792 314289 0
211 3294 519 114337 856956 745
212 141 23 75832 217193 1
213 572 6 61630 315366 0
214 13326 16 58580 307930 -1
215 2284 102 165548 702380 436
216 10 33 76403 194493 66
217 0 2 61656 316155 0
218 1414 0 60720 315380 0
219 1975 198 70184 330546 13
220 43 35 118881 394510 30
221 0 3 60887 312846 0
222 844 0 60720 315380 0
223 304 78 60925 296139 0
224 458 15 62969 295580 0
225 18562 11 58625 297765 0
226 0 155 102313 377934 143
227 7123 0 60720 315380 0
228 622 109 288170 638830 769
229 174 73 73007 304376 0
230 2220 1 64820 307424 0
231 121 22 301670 644190 154
232 11819 24 56178 295370 0
233 0 85 106113 574339 199
234 125 0 60720 315380 0
235 1182 12 60798 310201 0
236 1503 13 70694 327007 -5
237 0 39 56364 343466 6
238 0 0 60720 315380 0
239 30 0 60720 315380 0
240 3310 4 62045 318098 0
241 554 33 75230 448243 -6
242 0 11 79285 325738 30
243 468 0 60720 315380 0
244 4917 21 60720 312161 0
245 3256 42 52811 243650 -3
246 0 66 35250 407159 -8
247 125 0 60720 315380 0
248 0 13 59734 317698 0
249 22 0 60720 315380 0
250 0 2 60722 312502 0
251 514 0 60720 315380 0
252 0 52 78780 322378 0
253 0 0 60720 315380 0
254 0 0 60720 315380 0
255 0 0 60720 315380 0
256 0 0 60720 315380 0
257 0 0 60720 315380 0
258 10327 0 60720 315380 0
259 13718 91 88577 640273 135
260 3748 61 96448 345783 104
261 14416 9 50350 652925 53
262 1526 361 49857 439798 220
263 666 25 69351 278990 -4
264 2844 171 117869 339836 25
265 0 67 72683 240897 -15
266 368 0 60720 315380 0
267 0 17 61167 297141 0
268 333 0 60720 315380 0
269 26 13 70811 331323 0
270 0 1 60896 313880 0
271 0 0 60720 315380 0
272 1303 0 60720 315380 0
273 20 61 69863 309422 -5
274 2384 3 60938 315245 0
275 203 97 61348 405972 54
276 71 27 50804 300962 0
277 53 1 60745 316176 0
278 562 14 59506 302409 0
279 622 13 58456 283587 -3
280 645 7 60950 263276 -2
281 1763 5 60720 312075 0
282 0 24 61600 308336 0
283 317 0 60720 315380 0
284 1 9 63915 298700 0
285 275 6 60719 315372 0
286 0 4 59500 318745 0
287 936 0 60720 315380 0
288 8568 30 67939 408881 9
289 11528 37 32168 786690 223
290 0 70 -14545 -83265 -25
291 738 0 60720 315380 0
292 0 10 60720 321376 0
293 592 0 60720 315380 0
294 126 7 64270 276898 -4
295 2 10 60951 315547 0
296 0 1 60743 315487 0
297 18014 0 60720 315380 0
298 0 239 -1710 1405225 943
299 0 0 60720 315380 0
300 37704 0 60720 315380 0
301 63 452 60448 983660 2167
302 1431 4 65688 318574 0
303 94 34 106885 310768 30
304 192 3 61360 312887 0
305 32 4 65276 312339 0
306 6869 1 59988 314964 0
307 0 105 117520 379983 90
308 0 0 60720 315380 0
309 1 0 60720 315380 0
310 0 4 60722 315398 0
311 2328 0 60720 315380 0
312 209 94 82732 253588 0
313 28 13 64016 316647 0
314 176 2 60890 315688 0
315 1920 28 68136 310670 0
316 87550 15 79420 165404 -6
317 520 458 153198 4111912 17936
318 0 10 58650 291650 0
319 1013 0 60720 315380 0
320 15 2 69770 253468 -14
321 587 8 60831 315688 0
322 5371 12 59595 325699 0
323 0 66 87720 446211 49
324 1012 0 60720 315380 0
325 0 47 114768 368078 26
326 876 0 60720 315380 0
327 162556 59 138971 352850 24
328 43556 1081 213118 6282154 56289
329 3425 247 32648 227132 929
330 810 12 83620 283910 -3
331 0 43 74015 236761 26
332 0 0 60720 315380 0
333 2365 0 60720 315380 0
334 1261 79 191778 550608 120
335 0 33 76114 307528 13
336 1585 0 60720 315380 0
337 16189 107 93099 355864 2
338 0 295 116384 358589 379
339 0 0 60720 315380 0
340 10579 0 60720 315380 0
341 474 92 110309 375195 4
342 0 20 61977 288985 0
343 3642 0 60720 315380 0
344 0 86 90262 458343 22
345 472 0 60720 315380 0
346 98 20 80045 269587 0
347 3999 2 61490 315236 0
348 0 30 51252 42754 -58
349 621 0 60720 315380 0
350 0 48 60720 308256 0
351 0 0 60720 315380 0
352 30 0 60720 315380 0
353 0 4 60798 313164 0
354 746 0 60720 315380 0
355 0 23 71561 269661 10
356 0 0 60720 315380 0
357 0 0 60720 315380 0
358 0 0 60720 315380 0
359 4150 0 60720 315380 0
360 4658 150 134759 518365 66
361 814 158 156608 233773 58
362 1002 18 39625 301881 23
363 496 15 56750 298568 0
364 389 17 87390 325479 0
365 12679 14 58990 325506 0
366 400 110 48020 984885 1871
367 53 152 60720 313267 0
368 0 6 60349 315793 0
369 5109 0 60720 315380 0
370 576 97 67038 215362 -31
371 438 18 113761 314073 12
372 165 11 58320 298096 -5
373 4069 24 62841 325176 0
374 23 145 79804 207393 119
375 1285 5 62555 314806 0
376 4677 107 99489 341340 49
377 24811 75 90131 426280 24
378 167 505 95350 929118 822
379 0 2 64245 307322 0
380 4628 0 60720 315380 0
381 226 62 69159 387475 85
382 1765 5 65745 291787 0
383 460 20 77623 247060 28
384 36 15 63346 329784 0
385 989 3 60894 315834 0
386 3055 11 58930 304555 -5
387 0 43 60247 376641 14
388 0 0 60720 315380 0
389 0 0 60720 315380 0
390 0 0 60720 315380 0
391 0 0 60720 315380 0
392 13253 0 60720 315380 0
393 24 947 90829 357760 35
394 0 8 59818 315637 0
395 0 0 60720 315380 0
396 0 0 60720 315380 0
397 0 0 60720 315380 0
398 131 0 60720 315380 0
399 514 6 59661 327071 0
400 3366 15 102725 377516 23
401 9327 80 108479 299243 12
402 384 101 87419 387699 44
403 17821 22 54683 309836 0
404 397 373 122844 444477 643
405 897 10 62710 322327 -1
406 218 5 60720 314913 0
407 3369 21 60720 315553 0
408 5702 28 87161 688779 217
409 8636 61 101481 321896 10
410 534 855 128294 301607 108
411 0 7 62620 304485 0
412 0 0 60720 315380 0
413 726 0 60720 315380 0
414 1380 3 69980 231861 -5
415 180 26 60982 347385 25
416 7285 11 59635 316386 0
417 880 115 188873 491303 59
418 1 74 80791 261216 34
419 0 2 60727 315388 0
420 96 0 60720 315380 0
421 1889 5 60379 313729 0
422 45187 45 37527 358649 14
423 288 353 234817 1926517 2854
424 1270 4 60510 296656 0
425 6526 26 69206 275311 -10
426 0 26 55830 -42143 -28
427 226 0 60720 315380 0
428 694 12 72835 343929 0
429 0 8 68060 367655 -2
430 249 0 60720 315380 0
431 1595 72 56726 313491 0
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) Costs Trades Dividends Wealth
-1.245e+03 -6.645e-02 4.547e-02 2.878e-03 5.080e-03
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-17095 -2424 -1405 -514 156470
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -1.245e+03 9.850e+02 -1.264 0.20697
Costs -6.645e-02 1.263e-01 -0.526 0.59912
Trades 4.547e-02 1.587e-02 2.866 0.00436 **
Dividends 2.878e-03 1.542e-03 1.866 0.06268 .
Wealth 5.080e-03 2.546e-03 1.995 0.04666 *
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 10590 on 426 degrees of freedom
Multiple R-squared: 0.04961, Adjusted R-squared: 0.04069
F-statistic: 5.56 on 4 and 426 DF, p-value: 0.000227
> 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,] 6.355009e-09 1.271002e-08 1.000000e+00
[2,] 9.135221e-12 1.827044e-11 1.000000e+00
[3,] 1.449072e-14 2.898144e-14 1.000000e+00
[4,] 1.080923e-16 2.161846e-16 1.000000e+00
[5,] 1.974032e-19 3.948064e-19 1.000000e+00
[6,] 3.066813e-22 6.133625e-22 1.000000e+00
[7,] 9.776785e-24 1.955357e-23 1.000000e+00
[8,] 1.855129e-26 3.710259e-26 1.000000e+00
[9,] 4.133112e-29 8.266224e-29 1.000000e+00
[10,] 1.039992e-31 2.079983e-31 1.000000e+00
[11,] 2.587061e-34 5.174123e-34 1.000000e+00
[12,] 5.256491e-37 1.051298e-36 1.000000e+00
[13,] 9.422029e-40 1.884406e-39 1.000000e+00
[14,] 1.695782e-42 3.391564e-42 1.000000e+00
[15,] 7.203974e-45 1.440795e-44 1.000000e+00
[16,] 1.245824e-47 2.491649e-47 1.000000e+00
[17,] 1.039243e-49 2.078485e-49 1.000000e+00
[18,] 9.860932e-52 1.972186e-51 1.000000e+00
[19,] 1.780336e-54 3.560673e-54 1.000000e+00
[20,] 1.703993e-56 3.407987e-56 1.000000e+00
[21,] 6.968396e-58 1.393679e-57 1.000000e+00
[22,] 1.471697e-60 2.943393e-60 1.000000e+00
[23,] 4.126217e-63 8.252433e-63 1.000000e+00
[24,] 1.463468e-61 2.926935e-61 1.000000e+00
[25,] 4.246066e-64 8.492133e-64 1.000000e+00
[26,] 1.377975e-66 2.755950e-66 1.000000e+00
[27,] 4.181158e-69 8.362316e-69 1.000000e+00
[28,] 1.429607e-71 2.859213e-71 1.000000e+00
[29,] 7.942876e-73 1.588575e-72 1.000000e+00
[30,] 2.354922e-75 4.709845e-75 1.000000e+00
[31,] 1.818977e-77 3.637954e-77 1.000000e+00
[32,] 1.026964e-63 2.053929e-63 1.000000e+00
[33,] 2.208274e-65 4.416549e-65 1.000000e+00
[34,] 1.950868e-67 3.901736e-67 1.000000e+00
[35,] 3.239172e-69 6.478344e-69 1.000000e+00
[36,] 5.639839e-71 1.127968e-70 1.000000e+00
[37,] 5.491757e-73 1.098351e-72 1.000000e+00
[38,] 7.647239e-75 1.529448e-74 1.000000e+00
[39,] 7.058672e-77 1.411734e-76 1.000000e+00
[40,] 6.598836e-79 1.319767e-78 1.000000e+00
[41,] 6.262627e-81 1.252525e-80 1.000000e+00
[42,] 5.401203e-83 1.080241e-82 1.000000e+00
[43,] 4.617989e-85 9.235977e-85 1.000000e+00
[44,] 2.738523e-77 5.477046e-77 1.000000e+00
[45,] 4.018980e-79 8.037961e-79 1.000000e+00
[46,] 7.058372e-81 1.411674e-80 1.000000e+00
[47,] 1.374041e-81 2.748083e-81 1.000000e+00
[48,] 1.255381e-81 2.510762e-81 1.000000e+00
[49,] 1.036583e-79 2.073166e-79 1.000000e+00
[50,] 1.838095e-81 3.676190e-81 1.000000e+00
[51,] 4.103905e-83 8.207810e-83 1.000000e+00
[52,] 7.080039e-85 1.416008e-84 1.000000e+00
[53,] 1.225312e-86 2.450624e-86 1.000000e+00
[54,] 2.097743e-88 4.195487e-88 1.000000e+00
[55,] 3.521624e-90 7.043248e-90 1.000000e+00
[56,] 6.550513e-92 1.310103e-91 1.000000e+00
[57,] 1.086812e-93 2.173624e-93 1.000000e+00
[58,] 5.969403e-95 1.193881e-94 1.000000e+00
[59,] 9.883463e-97 1.976693e-96 1.000000e+00
[60,] 1.585118e-98 3.170236e-98 1.000000e+00
[61,] 2.530176e-100 5.060352e-100 1.000000e+00
[62,] 4.003172e-102 8.006344e-102 1.000000e+00
[63,] 6.393274e-104 1.278655e-103 1.000000e+00
[64,] 1.699581e-105 3.399162e-105 1.000000e+00
[65,] 2.655375e-107 5.310749e-107 1.000000e+00
[66,] 1.621016e-61 3.242033e-61 1.000000e+00
[67,] 1.239384e-62 2.478768e-62 1.000000e+00
[68,] 9.404972e-64 1.880994e-63 1.000000e+00
[69,] 7.211589e-65 1.442318e-64 1.000000e+00
[70,] 5.541259e-66 1.108252e-65 1.000000e+00
[71,] 4.098929e-67 8.197857e-67 1.000000e+00
[72,] 2.242646e-67 4.485291e-67 1.000000e+00
[73,] 1.898967e-68 3.797935e-68 1.000000e+00
[74,] 1.411180e-69 2.822360e-69 1.000000e+00
[75,] 1.040439e-70 2.080878e-70 1.000000e+00
[76,] 7.610784e-72 1.522157e-71 1.000000e+00
[77,] 5.539715e-73 1.107943e-72 1.000000e+00
[78,] 3.989917e-74 7.979833e-74 1.000000e+00
[79,] 2.851788e-75 5.703576e-75 1.000000e+00
[80,] 2.022319e-76 4.044638e-76 1.000000e+00
[81,] 1.424042e-77 2.848084e-77 1.000000e+00
[82,] 9.952282e-79 1.990456e-78 1.000000e+00
[83,] 6.908787e-80 1.381757e-79 1.000000e+00
[84,] 4.757988e-81 9.515975e-81 1.000000e+00
[85,] 3.252465e-82 6.504930e-82 1.000000e+00
[86,] 2.458891e-83 4.917783e-83 1.000000e+00
[87,] 5.774142e-84 1.154828e-83 1.000000e+00
[88,] 3.948514e-85 7.897029e-85 1.000000e+00
[89,] 2.678534e-86 5.357068e-86 1.000000e+00
[90,] 1.796914e-87 3.593829e-87 1.000000e+00
[91,] 1.197050e-88 2.394100e-88 1.000000e+00
[92,] 7.946404e-90 1.589281e-89 1.000000e+00
[93,] 5.212504e-91 1.042501e-90 1.000000e+00
[94,] 6.265803e-92 1.253161e-91 1.000000e+00
[95,] 4.097152e-93 8.194303e-93 1.000000e+00
[96,] 2.657281e-94 5.314563e-94 1.000000e+00
[97,] 1.710868e-95 3.421737e-95 1.000000e+00
[98,] 1.094237e-96 2.188473e-96 1.000000e+00
[99,] 6.466360e-95 1.293272e-94 1.000000e+00
[100,] 4.365693e-96 8.731386e-96 1.000000e+00
[101,] 1.431757e-96 2.863514e-96 1.000000e+00
[102,] 9.743408e-98 1.948682e-97 1.000000e+00
[103,] 6.621878e-99 1.324376e-98 1.000000e+00
[104,] 4.603028e-100 9.206056e-100 1.000000e+00
[105,] 5.735959e-101 1.147192e-100 1.000000e+00
[106,] 3.817379e-102 7.634759e-102 1.000000e+00
[107,] 2.521017e-103 5.042035e-103 1.000000e+00
[108,] 1.657272e-104 3.314544e-104 1.000000e+00
[109,] 1.083129e-105 2.166257e-105 1.000000e+00
[110,] 7.098872e-107 1.419774e-106 1.000000e+00
[111,] 7.726397e-107 1.545279e-106 1.000000e+00
[112,] 3.583620e-107 7.167240e-107 1.000000e+00
[113,] 4.060135e-108 8.120271e-108 1.000000e+00
[114,] 3.907457e-109 7.814914e-109 1.000000e+00
[115,] 2.587052e-110 5.174105e-110 1.000000e+00
[116,] 2.422664e-110 4.845329e-110 1.000000e+00
[117,] 1.626297e-111 3.252594e-111 1.000000e+00
[118,] 1.830911e-112 3.661822e-112 1.000000e+00
[119,] 1.216073e-113 2.432145e-113 1.000000e+00
[120,] 8.029455e-115 1.605891e-114 1.000000e+00
[121,] 5.272320e-116 1.054464e-115 1.000000e+00
[122,] 1.435193e-113 2.870386e-113 1.000000e+00
[123,] 1.025980e-114 2.051960e-114 1.000000e+00
[124,] 7.101983e-116 1.420397e-115 1.000000e+00
[125,] 4.879584e-117 9.759168e-117 1.000000e+00
[126,] 3.335184e-118 6.670369e-118 1.000000e+00
[127,] 2.269188e-119 4.538377e-119 1.000000e+00
[128,] 1.539602e-120 3.079203e-120 1.000000e+00
[129,] 3.164801e-121 6.329602e-121 1.000000e+00
[130,] 5.549521e-122 1.109904e-121 1.000000e+00
[131,] 4.191066e-123 8.382132e-123 1.000000e+00
[132,] 2.858683e-124 5.717366e-124 1.000000e+00
[133,] 2.757753e-121 5.515506e-121 1.000000e+00
[134,] 2.272874e-122 4.545748e-122 1.000000e+00
[135,] 2.256224e-123 4.512448e-123 1.000000e+00
[136,] 4.026316e-124 8.052632e-124 1.000000e+00
[137,] 3.316278e-125 6.632557e-125 1.000000e+00
[138,] 3.722789e-126 7.445579e-126 1.000000e+00
[139,] 2.882557e-127 5.765115e-127 1.000000e+00
[140,] 2.603447e-128 5.206894e-128 1.000000e+00
[141,] 3.845593e-129 7.691187e-129 1.000000e+00
[142,] 1.552153e-129 3.104307e-129 1.000000e+00
[143,] 1.807332e-130 3.614665e-130 1.000000e+00
[144,] 1.335573e-131 2.671146e-131 1.000000e+00
[145,] 1.443281e-132 2.886563e-132 1.000000e+00
[146,] 1.022703e-133 2.045406e-133 1.000000e+00
[147,] 7.177108e-135 1.435422e-134 1.000000e+00
[148,] 6.724100e-136 1.344820e-135 1.000000e+00
[149,] 2.133752e-136 4.267505e-136 1.000000e+00
[150,] 1.917652e-137 3.835305e-137 1.000000e+00
[151,] 1.648861e-138 3.297722e-138 1.000000e+00
[152,] 1.718317e-139 3.436634e-139 1.000000e+00
[153,] 5.827789e-139 1.165558e-138 1.000000e+00
[154,] 4.927695e-140 9.855390e-140 1.000000e+00
[155,] 3.601478e-141 7.202955e-141 1.000000e+00
[156,] 2.614367e-142 5.228733e-142 1.000000e+00
[157,] 2.462159e-143 4.924318e-143 1.000000e+00
[158,] 9.304644e-141 1.860929e-140 1.000000e+00
[159,] 2.826516e-141 5.653031e-141 1.000000e+00
[160,] 2.159855e-142 4.319710e-142 1.000000e+00
[161,] 2.867877e-139 5.735753e-139 1.000000e+00
[162,] 3.459159e-139 6.918318e-139 1.000000e+00
[163,] 5.380062e-140 1.076012e-139 1.000000e+00
[164,] 6.236011e-140 1.247202e-139 1.000000e+00
[165,] 4.926119e-141 9.852237e-141 1.000000e+00
[166,] 5.098702e-142 1.019740e-141 1.000000e+00
[167,] 6.894762e-143 1.378952e-142 1.000000e+00
[168,] 5.921565e-144 1.184313e-143 1.000000e+00
[169,] 1.251246e-144 2.502493e-144 1.000000e+00
[170,] 1.104826e-145 2.209652e-145 1.000000e+00
[171,] 8.776002e-147 1.755200e-146 1.000000e+00
[172,] 8.116032e-148 1.623206e-147 1.000000e+00
[173,] 1.812266e-145 3.624531e-145 1.000000e+00
[174,] 2.603590e-145 5.207180e-145 1.000000e+00
[175,] 1.668377e-145 3.336754e-145 1.000000e+00
[176,] 1.523634e-146 3.047268e-146 1.000000e+00
[177,] 1.382216e-147 2.764433e-147 1.000000e+00
[178,] 1.244788e-148 2.489575e-148 1.000000e+00
[179,] 1.099044e-149 2.198088e-149 1.000000e+00
[180,] 9.141012e-151 1.828202e-150 1.000000e+00
[181,] 1.320131e-22 2.640262e-22 1.000000e+00
[182,] 6.539563e-23 1.307913e-22 1.000000e+00
[183,] 3.268426e-23 6.536852e-23 1.000000e+00
[184,] 1.595432e-23 3.190864e-23 1.000000e+00
[185,] 7.905886e-24 1.581177e-23 1.000000e+00
[186,] 3.838458e-24 7.676916e-24 1.000000e+00
[187,] 1.829378e-24 3.658756e-24 1.000000e+00
[188,] 9.145385e-25 1.829077e-24 1.000000e+00
[189,] 4.269519e-25 8.539039e-25 1.000000e+00
[190,] 1.989307e-25 3.978614e-25 1.000000e+00
[191,] 9.358103e-26 1.871621e-25 1.000000e+00
[192,] 4.490088e-26 8.980177e-26 1.000000e+00
[193,] 1.446296e-24 2.892591e-24 1.000000e+00
[194,] 8.391581e-25 1.678316e-24 1.000000e+00
[195,] 4.113225e-25 8.226450e-25 1.000000e+00
[196,] 1.971736e-25 3.943472e-25 1.000000e+00
[197,] 9.562663e-26 1.912533e-25 1.000000e+00
[198,] 4.434294e-26 8.868587e-26 1.000000e+00
[199,] 2.046302e-26 4.092604e-26 1.000000e+00
[200,] 1.034614e-26 2.069227e-26 1.000000e+00
[201,] 4.867161e-27 9.734322e-27 1.000000e+00
[202,] 2.222400e-27 4.444801e-27 1.000000e+00
[203,] 6.585074e-26 1.317015e-25 1.000000e+00
[204,] 3.335289e-26 6.670577e-26 1.000000e+00
[205,] 1.645543e-26 3.291085e-26 1.000000e+00
[206,] 7.774521e-27 1.554904e-26 1.000000e+00
[207,] 1.309976e-26 2.619953e-26 1.000000e+00
[208,] 7.863099e-27 1.572620e-26 1.000000e+00
[209,] 3.861431e-27 7.722861e-27 1.000000e+00
[210,] 1.849486e-27 3.698973e-27 1.000000e+00
[211,] 8.471627e-28 1.694325e-27 1.000000e+00
[212,] 3.840444e-28 7.680889e-28 1.000000e+00
[213,] 2.109667e-28 4.219334e-28 1.000000e+00
[214,] 9.911567e-29 1.982313e-28 1.000000e+00
[215,] 4.509720e-29 9.019439e-29 1.000000e+00
[216,] 2.065923e-29 4.131845e-29 1.000000e+00
[217,] 9.429323e-30 1.885865e-29 1.000000e+00
[218,] 7.720242e-29 1.544048e-28 1.000000e+00
[219,] 3.970073e-29 7.940146e-29 1.000000e+00
[220,] 2.313574e-29 4.627148e-29 1.000000e+00
[221,] 2.555871e-29 5.111742e-29 1.000000e+00
[222,] 1.188713e-29 2.377427e-29 1.000000e+00
[223,] 5.302085e-30 1.060417e-29 1.000000e+00
[224,] 7.120146e-30 1.424029e-29 1.000000e+00
[225,] 9.244596e-30 1.848919e-29 1.000000e+00
[226,] 5.250926e-30 1.050185e-29 1.000000e+00
[227,] 2.411001e-30 4.822002e-30 1.000000e+00
[228,] 1.070399e-30 2.140798e-30 1.000000e+00
[229,] 4.749662e-31 9.499323e-31 1.000000e+00
[230,] 2.164159e-31 4.328318e-31 1.000000e+00
[231,] 9.811165e-32 1.962233e-31 1.000000e+00
[232,] 4.423177e-32 8.846354e-32 1.000000e+00
[233,] 1.949868e-32 3.899737e-32 1.000000e+00
[234,] 9.007438e-33 1.801488e-32 1.000000e+00
[235,] 4.148043e-33 8.296087e-33 1.000000e+00
[236,] 1.809456e-33 3.618911e-33 1.000000e+00
[237,] 8.526697e-34 1.705339e-33 1.000000e+00
[238,] 3.793147e-34 7.586295e-34 1.000000e+00
[239,] 1.697779e-34 3.395557e-34 1.000000e+00
[240,] 7.371153e-35 1.474231e-34 1.000000e+00
[241,] 3.195337e-35 6.390675e-35 1.000000e+00
[242,] 1.381074e-35 2.762148e-35 1.000000e+00
[243,] 5.943969e-36 1.188794e-35 1.000000e+00
[244,] 2.500815e-36 5.001629e-36 1.000000e+00
[245,] 1.096252e-36 2.192504e-36 1.000000e+00
[246,] 4.662055e-37 9.324110e-37 1.000000e+00
[247,] 1.974129e-37 3.948258e-37 1.000000e+00
[248,] 8.323521e-38 1.664704e-37 1.000000e+00
[249,] 3.494415e-38 6.988830e-38 1.000000e+00
[250,] 1.460764e-38 2.921528e-38 1.000000e+00
[251,] 1.366975e-38 2.733951e-38 1.000000e+00
[252,] 1.365018e-38 2.730037e-38 1.000000e+00
[253,] 5.697156e-39 1.139431e-38 1.000000e+00
[254,] 6.882827e-39 1.376565e-38 1.000000e+00
[255,] 2.923155e-39 5.846311e-39 1.000000e+00
[256,] 1.187891e-39 2.375782e-39 1.000000e+00
[257,] 5.025793e-40 1.005159e-39 1.000000e+00
[258,] 2.066367e-40 4.132735e-40 1.000000e+00
[259,] 8.337316e-41 1.667463e-40 1.000000e+00
[260,] 3.379082e-41 6.758164e-41 1.000000e+00
[261,] 1.353180e-41 2.706361e-41 1.000000e+00
[262,] 5.582541e-42 1.116508e-41 1.000000e+00
[263,] 2.249545e-42 4.499091e-42 1.000000e+00
[264,] 9.028984e-43 1.805797e-42 1.000000e+00
[265,] 3.448727e-43 6.897454e-43 1.000000e+00
[266,] 1.374539e-43 2.749077e-43 1.000000e+00
[267,] 5.174104e-44 1.034821e-43 1.000000e+00
[268,] 2.095197e-44 4.190393e-44 1.000000e+00
[269,] 8.023418e-45 1.604684e-44 1.000000e+00
[270,] 3.125512e-45 6.251024e-45 1.000000e+00
[271,] 1.178087e-45 2.356174e-45 1.000000e+00
[272,] 4.397364e-46 8.794728e-46 1.000000e+00
[273,] 1.641481e-46 3.282962e-46 1.000000e+00
[274,] 5.972512e-47 1.194502e-46 1.000000e+00
[275,] 2.269851e-47 4.539702e-47 1.000000e+00
[276,] 8.497050e-48 1.699410e-47 1.000000e+00
[277,] 3.215841e-48 6.431682e-48 1.000000e+00
[278,] 1.194283e-48 2.388566e-48 1.000000e+00
[279,] 4.474488e-49 8.948977e-49 1.000000e+00
[280,] 1.607086e-49 3.214172e-49 1.000000e+00
[281,] 8.523482e-50 1.704696e-49 1.000000e+00
[282,] 5.994889e-50 1.198978e-49 1.000000e+00
[283,] 2.380717e-50 4.761434e-50 1.000000e+00
[284,] 8.464673e-51 1.692935e-50 1.000000e+00
[285,] 3.099942e-51 6.199884e-51 1.000000e+00
[286,] 1.097849e-51 2.195698e-51 1.000000e+00
[287,] 3.949921e-52 7.899843e-52 1.000000e+00
[288,] 1.424379e-52 2.848757e-52 1.000000e+00
[289,] 5.121351e-53 1.024270e-52 1.000000e+00
[290,] 5.812597e-52 1.162519e-51 1.000000e+00
[291,] 3.611454e-51 7.222909e-51 1.000000e+00
[292,] 1.303537e-51 2.607073e-51 1.000000e+00
[293,] 3.882168e-45 7.764337e-45 1.000000e+00
[294,] 2.444750e-45 4.889500e-45 1.000000e+00
[295,] 9.029282e-46 1.805856e-45 1.000000e+00
[296,] 4.445364e-46 8.890727e-46 1.000000e+00
[297,] 1.674072e-46 3.348143e-46 1.000000e+00
[298,] 6.403981e-47 1.280796e-46 1.000000e+00
[299,] 3.047112e-47 6.094225e-47 1.000000e+00
[300,] 1.549263e-47 3.098526e-47 1.000000e+00
[301,] 5.773357e-48 1.154671e-47 1.000000e+00
[302,] 2.140656e-48 4.281312e-48 1.000000e+00
[303,] 7.890016e-49 1.578003e-48 1.000000e+00
[304,] 2.759074e-49 5.518148e-49 1.000000e+00
[305,] 1.082943e-49 2.165887e-49 1.000000e+00
[306,] 3.963922e-50 7.927844e-50 1.000000e+00
[307,] 1.421729e-50 2.843457e-50 1.000000e+00
[308,] 4.907726e-51 9.815453e-51 1.000000e+00
[309,] 9.296435e-24 1.859287e-23 1.000000e+00
[310,] 5.790618e-24 1.158124e-23 1.000000e+00
[311,] 2.849208e-24 5.698415e-24 1.000000e+00
[312,] 1.366308e-24 2.732615e-24 1.000000e+00
[313,] 6.972594e-25 1.394519e-24 1.000000e+00
[314,] 3.330515e-25 6.661031e-25 1.000000e+00
[315,] 1.670830e-25 3.341661e-25 1.000000e+00
[316,] 8.429982e-26 1.685996e-25 1.000000e+00
[317,] 3.927732e-26 7.855463e-26 1.000000e+00
[318,] 2.456025e-26 4.912050e-26 1.000000e+00
[319,] 1.135525e-26 2.271050e-26 1.000000e+00
[320,] 1.000000e+00 1.686519e-12 8.432593e-13
[321,] 1.000000e+00 7.470124e-13 3.735062e-13
[322,] 1.000000e+00 9.618237e-13 4.809118e-13
[323,] 1.000000e+00 1.856430e-12 9.282149e-13
[324,] 1.000000e+00 3.547471e-12 1.773736e-12
[325,] 1.000000e+00 6.329640e-12 3.164820e-12
[326,] 1.000000e+00 1.220358e-11 6.101791e-12
[327,] 1.000000e+00 2.149302e-11 1.074651e-11
[328,] 1.000000e+00 3.863128e-11 1.931564e-11
[329,] 1.000000e+00 7.207940e-11 3.603970e-11
[330,] 1.000000e+00 2.378167e-11 1.189084e-11
[331,] 1.000000e+00 4.443872e-11 2.221936e-11
[332,] 1.000000e+00 7.841180e-11 3.920590e-11
[333,] 1.000000e+00 8.300183e-11 4.150091e-11
[334,] 1.000000e+00 1.505770e-10 7.528850e-11
[335,] 1.000000e+00 2.696785e-10 1.348393e-10
[336,] 1.000000e+00 5.057881e-10 2.528940e-10
[337,] 1.000000e+00 8.080858e-10 4.040429e-10
[338,] 1.000000e+00 1.427065e-09 7.135325e-10
[339,] 1.000000e+00 2.562581e-09 1.281291e-09
[340,] 1.000000e+00 4.665974e-09 2.332987e-09
[341,] 1.000000e+00 8.498698e-09 4.249349e-09
[342,] 1.000000e+00 1.471871e-08 7.359356e-09
[343,] 1.000000e+00 2.472924e-08 1.236462e-08
[344,] 1.000000e+00 4.104065e-08 2.052032e-08
[345,] 1.000000e+00 6.771729e-08 3.385864e-08
[346,] 9.999999e-01 1.109401e-07 5.547004e-08
[347,] 9.999999e-01 1.861573e-07 9.307863e-08
[348,] 9.999998e-01 3.126103e-07 1.563051e-07
[349,] 9.999997e-01 5.001518e-07 2.500759e-07
[350,] 9.999996e-01 7.941141e-07 3.970570e-07
[351,] 9.999994e-01 1.251163e-06 6.255816e-07
[352,] 9.999989e-01 2.113199e-06 1.056600e-06
[353,] 9.999982e-01 3.526443e-06 1.763222e-06
[354,] 9.999972e-01 5.561596e-06 2.780798e-06
[355,] 9.999956e-01 8.784262e-06 4.392131e-06
[356,] 9.999931e-01 1.372988e-05 6.864942e-06
[357,] 9.999892e-01 2.161115e-05 1.080557e-05
[358,] 9.999898e-01 2.047832e-05 1.023916e-05
[359,] 9.999836e-01 3.281802e-05 1.640901e-05
[360,] 9.999751e-01 4.970085e-05 2.485043e-05
[361,] 9.999629e-01 7.428604e-05 3.714302e-05
[362,] 9.999420e-01 1.160242e-04 5.801209e-05
[363,] 9.999094e-01 1.811884e-04 9.059422e-05
[364,] 9.998607e-01 2.785229e-04 1.392614e-04
[365,] 9.997954e-01 4.092862e-04 2.046431e-04
[366,] 9.996864e-01 6.272569e-04 3.136285e-04
[367,] 9.995247e-01 9.505048e-04 4.752524e-04
[368,] 9.993004e-01 1.399169e-03 6.995846e-04
[369,] 9.989851e-01 2.029748e-03 1.014874e-03
[370,] 9.998718e-01 2.564812e-04 1.282406e-04
[371,] 9.998283e-01 3.434906e-04 1.717453e-04
[372,] 9.997414e-01 5.171018e-04 2.585509e-04
[373,] 9.995979e-01 8.041495e-04 4.020747e-04
[374,] 9.994111e-01 1.177769e-03 5.888844e-04
[375,] 9.990972e-01 1.805523e-03 9.027615e-04
[376,] 9.986319e-01 2.736240e-03 1.368120e-03
[377,] 9.980594e-01 3.881248e-03 1.940624e-03
[378,] 9.971915e-01 5.617053e-03 2.808526e-03
[379,] 9.958603e-01 8.279328e-03 4.139664e-03
[380,] 9.944932e-01 1.101366e-02 5.506832e-03
[381,] 9.924994e-01 1.500127e-02 7.500634e-03
[382,] 9.898955e-01 2.020898e-02 1.010449e-02
[383,] 9.865383e-01 2.692345e-02 1.346172e-02
[384,] 9.822660e-01 3.546797e-02 1.773398e-02
[385,] 9.822820e-01 3.543609e-02 1.771805e-02
[386,] 9.787556e-01 4.248875e-02 2.124437e-02
[387,] 9.721992e-01 5.560150e-02 2.780075e-02
[388,] 9.638944e-01 7.221111e-02 3.610555e-02
[389,] 9.536783e-01 9.264345e-02 4.632173e-02
[390,] 9.412996e-01 1.174009e-01 5.870043e-02
[391,] 9.262661e-01 1.474677e-01 7.373386e-02
[392,] 9.085606e-01 1.828788e-01 9.143940e-02
[393,] 8.817873e-01 2.364254e-01 1.182127e-01
[394,] 8.809948e-01 2.380103e-01 1.190052e-01
[395,] 8.511728e-01 2.976543e-01 1.488272e-01
[396,] 8.754873e-01 2.490254e-01 1.245127e-01
[397,] 8.413464e-01 3.173071e-01 1.586536e-01
[398,] 8.029967e-01 3.940065e-01 1.970033e-01
[399,] 7.619043e-01 4.761915e-01 2.380957e-01
[400,] 7.064266e-01 5.871468e-01 2.935734e-01
[401,] 6.542459e-01 6.915081e-01 3.457541e-01
[402,] 6.338070e-01 7.323861e-01 3.661930e-01
[403,] 8.200569e-01 3.598863e-01 1.799431e-01
[404,] 7.666385e-01 4.667230e-01 2.333615e-01
[405,] 7.027285e-01 5.945430e-01 2.972715e-01
[406,] 6.285612e-01 7.428777e-01 3.714388e-01
[407,] 5.651762e-01 8.696475e-01 4.348238e-01
[408,] 5.252115e-01 9.495770e-01 4.747885e-01
[409,] 4.402295e-01 8.804589e-01 5.597705e-01
[410,] 7.995580e-01 4.008840e-01 2.004420e-01
[411,] 7.717690e-01 4.564620e-01 2.282310e-01
[412,] 6.967081e-01 6.065838e-01 3.032919e-01
[413,] 6.112946e-01 7.774108e-01 3.887054e-01
[414,] 5.061653e-01 9.876693e-01 4.938347e-01
[415,] 9.951747e-01 9.650621e-03 4.825311e-03
[416,] 9.811266e-01 3.774679e-02 1.887339e-02
> postscript(file="/var/www/rcomp/tmp/1uinc1292675864.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> plot(x[,1], type='l', main='Actuals and Interpolation', ylab='value of Actuals and Interpolation (dots)', xlab='time or index')
> points(x[,1]-mysum$resid)
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/rcomp/tmp/2nr4f1292675864.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> plot(mysum$resid, type='b', pch=19, main='Residuals', ylab='value of Residuals', xlab='time or index')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/rcomp/tmp/3nr4f1292675864.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> hist(mysum$resid, main='Residual Histogram', xlab='values of Residuals')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/rcomp/tmp/4nr4f1292675864.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> densityplot(~mysum$resid,col='black',main='Residual Density Plot', xlab='values of Residuals')
> dev.off()
null device
1
> postscript(file="/var/www/rcomp/tmp/5gi401292675864.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
-446.799014 -1126.182058 -511.543395 -902.728119 120.760443
6 7 8 9 10
-344.107520 -381.130939 -524.381751 -660.956188 -460.062309
11 12 13 14 15
-894.317669 -904.534116 -531.752996 -373.326204 -494.215890
16 17 18 19 20
-369.054268 -522.363459 -541.639199 -543.244901 -522.622415
21 22 23 24 25
-531.752996 -385.786972 -531.752996 -1243.621731 -743.010377
26 27 28 29 30
-515.477844 -408.176768 -784.862499 -531.752996 -659.261369
31 32 33 34 35
91.977173 -531.752996 -590.863508 -1726.837388 -307.246345
36 37 38 39 40
-140.114489 -529.762446 -180.711314 -668.761908 -835.720751
41 42 43 44 45
-1228.182183 -781.072067 -914.536718 -531.752996 -449.831058
46 47 48 49 50
-436.851090 -531.752996 -492.797868 -646.497911 -957.033521
51 52 53 54 55
-240.473167 -701.717963 -219.629339 -581.289588 -3383.558440
56 57 58 59 60
-766.028236 -511.069047 -97.596255 -357.909723 -531.752996
61 62 63 64 65
-531.752996 -526.747763 145.084534 -534.244498 655.246817
66 67 68 69 70
-526.236610 -856.532957 -602.958505 -457.130721 -531.752996
71 72 73 74 75
-3174.838918 211.927550 9354.448906 -531.752996 -578.746650
76 77 78 79 80
-581.164370 -633.399028 -518.399530 1113.396719 -37.206470
81 82 83 84 85
-529.007476 -531.752996 -531.752996 -554.257189 -531.752996
86 87 88 89 90
-531.752996 -521.732620 -531.449203 -444.201528 -670.889173
91 92 93 94 95
-531.752996 -401.698129 -1291.243244 465.050125 -563.965545
96 97 98 99 100
-630.510274 -531.752996 -514.103829 -542.910900 -264.378147
101 102 103 104 105
-1074.155039 -434.808579 -542.900405 -531.755114 -531.752996
106 107 108 109 110
2029.650731 -431.446382 -1695.157848 -384.687406 -510.562047
111 112 113 114 115
-518.842890 51.592784 -531.752996 -525.641602 -531.752996
116 117 118 119 120
-461.320462 -888.829048 -1887.590163 -2138.572831 -1490.581032
121 122 123 124 125
3.856788 -531.752996 -2095.853648 -878.399845 -20.323604
126 127 128 129 130
-531.752996 -545.554349 -531.752996 1523.730439 -476.086445
131 132 133 134 135
-585.320451 -531.752996 -531.752996 -531.752996 -520.235254
136 137 138 139 140
-472.722529 -2088.524181 -2592.868414 -2423.524181 2101.475819
141 142 143 144 145
-1404.972548 -575.483408 -2423.524181 -1785.524181 -2457.990409
146 147 148 149 150
-1816.524181 -871.909281 -1019.813082 -2653.563065 -2423.524181
151 152 153 154 155
-1812.524181 -2300.277558 -1500.524181 -1648.256178 -64.600804
156 157 158 159 160
-2586.926486 -2553.057962 -2288.524181 -822.180481 -4711.001339
161 162 163 164 165
-2412.082568 -2114.524181 -2176.304013 -2596.137187 2095.475819
166 167 168 169 170
-5116.656277 -1586.524181 2674.509898 -336.323074 -240.871400
171 172 173 174 175
-3787.369387 -2440.813076 -1905.689883 -2410.750338 -1495.011601
176 177 178 179 180
-3495.484878 -2362.524181 -2052.208190 -2431.352844 2477.475819
181 182 183 184 185
-5103.395872 -4500.593947 -2423.524181 -2423.524181 -2421.524181
186 187 188 189 190
-2387.648451 -1677.093216 82103.337225 1803.203792 -2420.524181
191 192 193 194 195
-1896.916885 -2511.118916 -2018.524181 -1517.305540 -3044.373271
196 197 198 199 200
-324.524181 -1049.596659 -1391.588625 -2423.524181 19906.475819
201 202 203 204 205
-4429.102815 -2423.524181 -1940.524181 -2400.509309 -184.524181
206 207 208 209 210
-1107.305949 -3611.207824 -2058.524181 -1250.547881 19757.475870
211 212 213 214 215
-3095.283585 -2685.572436 -1892.462246 11022.287859 -6015.116649
216 217 218 219 220
-2776.875487 -2468.180819 -1009.524181 -909.398752 -5250.609171
221 222 223 224 225
-2423.625909 -1579.524181 -2068.290876 -2009.807846 16285.157028
226 227 228 229 230
-4485.171685 4699.475819 -13070.996208 -2771.687434 -366.986238
231 232 233 234 235
-14203.913927 9661.176538 -5228.101162 -2298.524181 -1229.369359
236 237 238 239 240
-1406.605618 -2303.724111 -2423.524181 -2393.524181 818.673002
241 242 243 244 245
-2909.410634 -3296.891128 -1955.524181 2504.134858 1401.321731
246 247 248 249 250
-1525.114460 -2298.524181 -2384.498284 -2401.524181 -2415.199970
251 252 253 254 255
-1909.524181 -3261.382924 -2423.524181 -2423.524181 -2423.524181
256 257 258 259 260
-2423.524181 -2423.524181 7903.475819 9098.228752 -384.019079
261 262 263 264 265
11492.941439 -738.767412 -2043.565720 -2237.191848 -2748.599377
266 267 268 269 270
-2055.524181 -2390.231484 -2090.524181 -2901.371416 -2427.143658
271 272 273 274 275
-2423.524181 -1120.524181 -2798.025251 -48.848634 -2503.611189
276 277 278 279 280
-1858.364649 -2373.885174 -1768.066474 -1606.209373 -1638.563020
281 282 283 284 285
-650.680867 -2441.671267 -2106.524181 -2519.199175 -2148.056984
286 287 288 289 290
-2377.469561 -1487.524181 5549.105588 9047.715220 2150.713080
291 292 293 294 295
-1685.524181 -2440.114833 -1831.524181 -2347.712027 -2431.843668
296 297 298 299 300
-2424.811446 15590.475819 -2710.117007 -2423.524181 35280.475819
301 302 303 304 305
-4252.290035 -1227.341467 -4413.234993 -2253.250875 -2589.665167
306 307 308 309 310
4480.022951 -5185.573175 -2423.524181 -2422.524181 -2423.401116
311 312 313 314 315
-95.524181 -3031.324247 -2548.173266 -2256.007415 -825.309520
316 317 318 319 320
84708.824961 -17094.684688 -2260.448538 -1410.524181 -2641.648261
321 322 323 324 325
-1841.926021 2969.730390 -4023.560995 -1411.524181 -5029.710777
326 327 328 329 330
-1547.524181 156470.426224 16817.902759 2543.542624 -2563.394724
331 332 333 334 335
-2799.065069 -2423.524181 -58.524181 -7793.933302 -3098.755366
336 337 338 339 340
-838.524181 12183.821560 -5061.195168 -2423.524181 8155.475819
341 342 343 344 345
-4370.342296 -2403.391098 1218.475819 -4172.590994 -1951.524181
346 347 348 349 350
-3071.107195 1541.011763 -1206.175433 -1802.524181 -2399.833242
351 352 353 354 355
-2423.524181 -2393.524181 -2420.427824 -1677.524181 -2783.410082
356 357 358 359 360
-2423.524181 -2423.524181 -2423.524181 1726.475819 -1706.538957
361 362 363 364 365
-5724.432563 -422.422949 -1697.633094 -3275.123340 10305.927627
366 367 368 369 370
-3374.945564 -2354.342898 -2407.444889 2685.475819 -1840.366319
371 372 373 374 375
-4392.364717 -2098.901984 1522.439542 -2948.466625 -1219.976240
376 377 378 379 380
422.830876 20735.894586 -5567.945758 -2560.481408 2204.475819
381 382 383 384 385
-2785.024837 -818.779707 -2534.295064 -2547.381315 -1443.542997
386 387 388 389 390
744.774172 -2575.526532 -2423.524181 -2423.524181 -2423.524181
391 392 393 394 395
-2423.524181 10829.475819 -3827.768280 -2382.718857 -2423.524181
396 397 398 399 400
-2423.524181 -2423.524181 -2292.524181 -1894.617602 -1145.395706
401 402 403 404 405
4783.601521 -3455.138968 15689.390287 -5201.251309 -1636.330421
406 407 408 409 410
-2203.848003 946.373432 1002.420813 4344.352928 -4866.163850
411 412 413 414 415
-2478.097288 -2423.524181 -1697.524181 -1223.995912 -2345.939710
416 417 418 419 420
4908.645920 -7869.476244 -3174.521760 -2423.732586 -2327.524181
421 422 423 424 425
-513.935683 43696.445781 -14679.127902 -1089.826295 3833.711605
426 427 428 429 430
-1170.437279 -2197.524181 -2361.745002 -2907.162935 -2174.524181
431
-636.699222
> postscript(file="/var/www/rcomp/tmp/6gi401292675864.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 -446.799014 NA
1 -1126.182058 -446.799014
2 -511.543395 -1126.182058
3 -902.728119 -511.543395
4 120.760443 -902.728119
5 -344.107520 120.760443
6 -381.130939 -344.107520
7 -524.381751 -381.130939
8 -660.956188 -524.381751
9 -460.062309 -660.956188
10 -894.317669 -460.062309
11 -904.534116 -894.317669
12 -531.752996 -904.534116
13 -373.326204 -531.752996
14 -494.215890 -373.326204
15 -369.054268 -494.215890
16 -522.363459 -369.054268
17 -541.639199 -522.363459
18 -543.244901 -541.639199
19 -522.622415 -543.244901
20 -531.752996 -522.622415
21 -385.786972 -531.752996
22 -531.752996 -385.786972
23 -1243.621731 -531.752996
24 -743.010377 -1243.621731
25 -515.477844 -743.010377
26 -408.176768 -515.477844
27 -784.862499 -408.176768
28 -531.752996 -784.862499
29 -659.261369 -531.752996
30 91.977173 -659.261369
31 -531.752996 91.977173
32 -590.863508 -531.752996
33 -1726.837388 -590.863508
34 -307.246345 -1726.837388
35 -140.114489 -307.246345
36 -529.762446 -140.114489
37 -180.711314 -529.762446
38 -668.761908 -180.711314
39 -835.720751 -668.761908
40 -1228.182183 -835.720751
41 -781.072067 -1228.182183
42 -914.536718 -781.072067
43 -531.752996 -914.536718
44 -449.831058 -531.752996
45 -436.851090 -449.831058
46 -531.752996 -436.851090
47 -492.797868 -531.752996
48 -646.497911 -492.797868
49 -957.033521 -646.497911
50 -240.473167 -957.033521
51 -701.717963 -240.473167
52 -219.629339 -701.717963
53 -581.289588 -219.629339
54 -3383.558440 -581.289588
55 -766.028236 -3383.558440
56 -511.069047 -766.028236
57 -97.596255 -511.069047
58 -357.909723 -97.596255
59 -531.752996 -357.909723
60 -531.752996 -531.752996
61 -526.747763 -531.752996
62 145.084534 -526.747763
63 -534.244498 145.084534
64 655.246817 -534.244498
65 -526.236610 655.246817
66 -856.532957 -526.236610
67 -602.958505 -856.532957
68 -457.130721 -602.958505
69 -531.752996 -457.130721
70 -3174.838918 -531.752996
71 211.927550 -3174.838918
72 9354.448906 211.927550
73 -531.752996 9354.448906
74 -578.746650 -531.752996
75 -581.164370 -578.746650
76 -633.399028 -581.164370
77 -518.399530 -633.399028
78 1113.396719 -518.399530
79 -37.206470 1113.396719
80 -529.007476 -37.206470
81 -531.752996 -529.007476
82 -531.752996 -531.752996
83 -554.257189 -531.752996
84 -531.752996 -554.257189
85 -531.752996 -531.752996
86 -521.732620 -531.752996
87 -531.449203 -521.732620
88 -444.201528 -531.449203
89 -670.889173 -444.201528
90 -531.752996 -670.889173
91 -401.698129 -531.752996
92 -1291.243244 -401.698129
93 465.050125 -1291.243244
94 -563.965545 465.050125
95 -630.510274 -563.965545
96 -531.752996 -630.510274
97 -514.103829 -531.752996
98 -542.910900 -514.103829
99 -264.378147 -542.910900
100 -1074.155039 -264.378147
101 -434.808579 -1074.155039
102 -542.900405 -434.808579
103 -531.755114 -542.900405
104 -531.752996 -531.755114
105 2029.650731 -531.752996
106 -431.446382 2029.650731
107 -1695.157848 -431.446382
108 -384.687406 -1695.157848
109 -510.562047 -384.687406
110 -518.842890 -510.562047
111 51.592784 -518.842890
112 -531.752996 51.592784
113 -525.641602 -531.752996
114 -531.752996 -525.641602
115 -461.320462 -531.752996
116 -888.829048 -461.320462
117 -1887.590163 -888.829048
118 -2138.572831 -1887.590163
119 -1490.581032 -2138.572831
120 3.856788 -1490.581032
121 -531.752996 3.856788
122 -2095.853648 -531.752996
123 -878.399845 -2095.853648
124 -20.323604 -878.399845
125 -531.752996 -20.323604
126 -545.554349 -531.752996
127 -531.752996 -545.554349
128 1523.730439 -531.752996
129 -476.086445 1523.730439
130 -585.320451 -476.086445
131 -531.752996 -585.320451
132 -531.752996 -531.752996
133 -531.752996 -531.752996
134 -520.235254 -531.752996
135 -472.722529 -520.235254
136 -2088.524181 -472.722529
137 -2592.868414 -2088.524181
138 -2423.524181 -2592.868414
139 2101.475819 -2423.524181
140 -1404.972548 2101.475819
141 -575.483408 -1404.972548
142 -2423.524181 -575.483408
143 -1785.524181 -2423.524181
144 -2457.990409 -1785.524181
145 -1816.524181 -2457.990409
146 -871.909281 -1816.524181
147 -1019.813082 -871.909281
148 -2653.563065 -1019.813082
149 -2423.524181 -2653.563065
150 -1812.524181 -2423.524181
151 -2300.277558 -1812.524181
152 -1500.524181 -2300.277558
153 -1648.256178 -1500.524181
154 -64.600804 -1648.256178
155 -2586.926486 -64.600804
156 -2553.057962 -2586.926486
157 -2288.524181 -2553.057962
158 -822.180481 -2288.524181
159 -4711.001339 -822.180481
160 -2412.082568 -4711.001339
161 -2114.524181 -2412.082568
162 -2176.304013 -2114.524181
163 -2596.137187 -2176.304013
164 2095.475819 -2596.137187
165 -5116.656277 2095.475819
166 -1586.524181 -5116.656277
167 2674.509898 -1586.524181
168 -336.323074 2674.509898
169 -240.871400 -336.323074
170 -3787.369387 -240.871400
171 -2440.813076 -3787.369387
172 -1905.689883 -2440.813076
173 -2410.750338 -1905.689883
174 -1495.011601 -2410.750338
175 -3495.484878 -1495.011601
176 -2362.524181 -3495.484878
177 -2052.208190 -2362.524181
178 -2431.352844 -2052.208190
179 2477.475819 -2431.352844
180 -5103.395872 2477.475819
181 -4500.593947 -5103.395872
182 -2423.524181 -4500.593947
183 -2423.524181 -2423.524181
184 -2421.524181 -2423.524181
185 -2387.648451 -2421.524181
186 -1677.093216 -2387.648451
187 82103.337225 -1677.093216
188 1803.203792 82103.337225
189 -2420.524181 1803.203792
190 -1896.916885 -2420.524181
191 -2511.118916 -1896.916885
192 -2018.524181 -2511.118916
193 -1517.305540 -2018.524181
194 -3044.373271 -1517.305540
195 -324.524181 -3044.373271
196 -1049.596659 -324.524181
197 -1391.588625 -1049.596659
198 -2423.524181 -1391.588625
199 19906.475819 -2423.524181
200 -4429.102815 19906.475819
201 -2423.524181 -4429.102815
202 -1940.524181 -2423.524181
203 -2400.509309 -1940.524181
204 -184.524181 -2400.509309
205 -1107.305949 -184.524181
206 -3611.207824 -1107.305949
207 -2058.524181 -3611.207824
208 -1250.547881 -2058.524181
209 19757.475870 -1250.547881
210 -3095.283585 19757.475870
211 -2685.572436 -3095.283585
212 -1892.462246 -2685.572436
213 11022.287859 -1892.462246
214 -6015.116649 11022.287859
215 -2776.875487 -6015.116649
216 -2468.180819 -2776.875487
217 -1009.524181 -2468.180819
218 -909.398752 -1009.524181
219 -5250.609171 -909.398752
220 -2423.625909 -5250.609171
221 -1579.524181 -2423.625909
222 -2068.290876 -1579.524181
223 -2009.807846 -2068.290876
224 16285.157028 -2009.807846
225 -4485.171685 16285.157028
226 4699.475819 -4485.171685
227 -13070.996208 4699.475819
228 -2771.687434 -13070.996208
229 -366.986238 -2771.687434
230 -14203.913927 -366.986238
231 9661.176538 -14203.913927
232 -5228.101162 9661.176538
233 -2298.524181 -5228.101162
234 -1229.369359 -2298.524181
235 -1406.605618 -1229.369359
236 -2303.724111 -1406.605618
237 -2423.524181 -2303.724111
238 -2393.524181 -2423.524181
239 818.673002 -2393.524181
240 -2909.410634 818.673002
241 -3296.891128 -2909.410634
242 -1955.524181 -3296.891128
243 2504.134858 -1955.524181
244 1401.321731 2504.134858
245 -1525.114460 1401.321731
246 -2298.524181 -1525.114460
247 -2384.498284 -2298.524181
248 -2401.524181 -2384.498284
249 -2415.199970 -2401.524181
250 -1909.524181 -2415.199970
251 -3261.382924 -1909.524181
252 -2423.524181 -3261.382924
253 -2423.524181 -2423.524181
254 -2423.524181 -2423.524181
255 -2423.524181 -2423.524181
256 -2423.524181 -2423.524181
257 7903.475819 -2423.524181
258 9098.228752 7903.475819
259 -384.019079 9098.228752
260 11492.941439 -384.019079
261 -738.767412 11492.941439
262 -2043.565720 -738.767412
263 -2237.191848 -2043.565720
264 -2748.599377 -2237.191848
265 -2055.524181 -2748.599377
266 -2390.231484 -2055.524181
267 -2090.524181 -2390.231484
268 -2901.371416 -2090.524181
269 -2427.143658 -2901.371416
270 -2423.524181 -2427.143658
271 -1120.524181 -2423.524181
272 -2798.025251 -1120.524181
273 -48.848634 -2798.025251
274 -2503.611189 -48.848634
275 -1858.364649 -2503.611189
276 -2373.885174 -1858.364649
277 -1768.066474 -2373.885174
278 -1606.209373 -1768.066474
279 -1638.563020 -1606.209373
280 -650.680867 -1638.563020
281 -2441.671267 -650.680867
282 -2106.524181 -2441.671267
283 -2519.199175 -2106.524181
284 -2148.056984 -2519.199175
285 -2377.469561 -2148.056984
286 -1487.524181 -2377.469561
287 5549.105588 -1487.524181
288 9047.715220 5549.105588
289 2150.713080 9047.715220
290 -1685.524181 2150.713080
291 -2440.114833 -1685.524181
292 -1831.524181 -2440.114833
293 -2347.712027 -1831.524181
294 -2431.843668 -2347.712027
295 -2424.811446 -2431.843668
296 15590.475819 -2424.811446
297 -2710.117007 15590.475819
298 -2423.524181 -2710.117007
299 35280.475819 -2423.524181
300 -4252.290035 35280.475819
301 -1227.341467 -4252.290035
302 -4413.234993 -1227.341467
303 -2253.250875 -4413.234993
304 -2589.665167 -2253.250875
305 4480.022951 -2589.665167
306 -5185.573175 4480.022951
307 -2423.524181 -5185.573175
308 -2422.524181 -2423.524181
309 -2423.401116 -2422.524181
310 -95.524181 -2423.401116
311 -3031.324247 -95.524181
312 -2548.173266 -3031.324247
313 -2256.007415 -2548.173266
314 -825.309520 -2256.007415
315 84708.824961 -825.309520
316 -17094.684688 84708.824961
317 -2260.448538 -17094.684688
318 -1410.524181 -2260.448538
319 -2641.648261 -1410.524181
320 -1841.926021 -2641.648261
321 2969.730390 -1841.926021
322 -4023.560995 2969.730390
323 -1411.524181 -4023.560995
324 -5029.710777 -1411.524181
325 -1547.524181 -5029.710777
326 156470.426224 -1547.524181
327 16817.902759 156470.426224
328 2543.542624 16817.902759
329 -2563.394724 2543.542624
330 -2799.065069 -2563.394724
331 -2423.524181 -2799.065069
332 -58.524181 -2423.524181
333 -7793.933302 -58.524181
334 -3098.755366 -7793.933302
335 -838.524181 -3098.755366
336 12183.821560 -838.524181
337 -5061.195168 12183.821560
338 -2423.524181 -5061.195168
339 8155.475819 -2423.524181
340 -4370.342296 8155.475819
341 -2403.391098 -4370.342296
342 1218.475819 -2403.391098
343 -4172.590994 1218.475819
344 -1951.524181 -4172.590994
345 -3071.107195 -1951.524181
346 1541.011763 -3071.107195
347 -1206.175433 1541.011763
348 -1802.524181 -1206.175433
349 -2399.833242 -1802.524181
350 -2423.524181 -2399.833242
351 -2393.524181 -2423.524181
352 -2420.427824 -2393.524181
353 -1677.524181 -2420.427824
354 -2783.410082 -1677.524181
355 -2423.524181 -2783.410082
356 -2423.524181 -2423.524181
357 -2423.524181 -2423.524181
358 1726.475819 -2423.524181
359 -1706.538957 1726.475819
360 -5724.432563 -1706.538957
361 -422.422949 -5724.432563
362 -1697.633094 -422.422949
363 -3275.123340 -1697.633094
364 10305.927627 -3275.123340
365 -3374.945564 10305.927627
366 -2354.342898 -3374.945564
367 -2407.444889 -2354.342898
368 2685.475819 -2407.444889
369 -1840.366319 2685.475819
370 -4392.364717 -1840.366319
371 -2098.901984 -4392.364717
372 1522.439542 -2098.901984
373 -2948.466625 1522.439542
374 -1219.976240 -2948.466625
375 422.830876 -1219.976240
376 20735.894586 422.830876
377 -5567.945758 20735.894586
378 -2560.481408 -5567.945758
379 2204.475819 -2560.481408
380 -2785.024837 2204.475819
381 -818.779707 -2785.024837
382 -2534.295064 -818.779707
383 -2547.381315 -2534.295064
384 -1443.542997 -2547.381315
385 744.774172 -1443.542997
386 -2575.526532 744.774172
387 -2423.524181 -2575.526532
388 -2423.524181 -2423.524181
389 -2423.524181 -2423.524181
390 -2423.524181 -2423.524181
391 10829.475819 -2423.524181
392 -3827.768280 10829.475819
393 -2382.718857 -3827.768280
394 -2423.524181 -2382.718857
395 -2423.524181 -2423.524181
396 -2423.524181 -2423.524181
397 -2292.524181 -2423.524181
398 -1894.617602 -2292.524181
399 -1145.395706 -1894.617602
400 4783.601521 -1145.395706
401 -3455.138968 4783.601521
402 15689.390287 -3455.138968
403 -5201.251309 15689.390287
404 -1636.330421 -5201.251309
405 -2203.848003 -1636.330421
406 946.373432 -2203.848003
407 1002.420813 946.373432
408 4344.352928 1002.420813
409 -4866.163850 4344.352928
410 -2478.097288 -4866.163850
411 -2423.524181 -2478.097288
412 -1697.524181 -2423.524181
413 -1223.995912 -1697.524181
414 -2345.939710 -1223.995912
415 4908.645920 -2345.939710
416 -7869.476244 4908.645920
417 -3174.521760 -7869.476244
418 -2423.732586 -3174.521760
419 -2327.524181 -2423.732586
420 -513.935683 -2327.524181
421 43696.445781 -513.935683
422 -14679.127902 43696.445781
423 -1089.826295 -14679.127902
424 3833.711605 -1089.826295
425 -1170.437279 3833.711605
426 -2197.524181 -1170.437279
427 -2361.745002 -2197.524181
428 -2907.162935 -2361.745002
429 -2174.524181 -2907.162935
430 -636.699222 -2174.524181
431 NA -636.699222
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] -1126.182058 -446.799014
[2,] -511.543395 -1126.182058
[3,] -902.728119 -511.543395
[4,] 120.760443 -902.728119
[5,] -344.107520 120.760443
[6,] -381.130939 -344.107520
[7,] -524.381751 -381.130939
[8,] -660.956188 -524.381751
[9,] -460.062309 -660.956188
[10,] -894.317669 -460.062309
[11,] -904.534116 -894.317669
[12,] -531.752996 -904.534116
[13,] -373.326204 -531.752996
[14,] -494.215890 -373.326204
[15,] -369.054268 -494.215890
[16,] -522.363459 -369.054268
[17,] -541.639199 -522.363459
[18,] -543.244901 -541.639199
[19,] -522.622415 -543.244901
[20,] -531.752996 -522.622415
[21,] -385.786972 -531.752996
[22,] -531.752996 -385.786972
[23,] -1243.621731 -531.752996
[24,] -743.010377 -1243.621731
[25,] -515.477844 -743.010377
[26,] -408.176768 -515.477844
[27,] -784.862499 -408.176768
[28,] -531.752996 -784.862499
[29,] -659.261369 -531.752996
[30,] 91.977173 -659.261369
[31,] -531.752996 91.977173
[32,] -590.863508 -531.752996
[33,] -1726.837388 -590.863508
[34,] -307.246345 -1726.837388
[35,] -140.114489 -307.246345
[36,] -529.762446 -140.114489
[37,] -180.711314 -529.762446
[38,] -668.761908 -180.711314
[39,] -835.720751 -668.761908
[40,] -1228.182183 -835.720751
[41,] -781.072067 -1228.182183
[42,] -914.536718 -781.072067
[43,] -531.752996 -914.536718
[44,] -449.831058 -531.752996
[45,] -436.851090 -449.831058
[46,] -531.752996 -436.851090
[47,] -492.797868 -531.752996
[48,] -646.497911 -492.797868
[49,] -957.033521 -646.497911
[50,] -240.473167 -957.033521
[51,] -701.717963 -240.473167
[52,] -219.629339 -701.717963
[53,] -581.289588 -219.629339
[54,] -3383.558440 -581.289588
[55,] -766.028236 -3383.558440
[56,] -511.069047 -766.028236
[57,] -97.596255 -511.069047
[58,] -357.909723 -97.596255
[59,] -531.752996 -357.909723
[60,] -531.752996 -531.752996
[61,] -526.747763 -531.752996
[62,] 145.084534 -526.747763
[63,] -534.244498 145.084534
[64,] 655.246817 -534.244498
[65,] -526.236610 655.246817
[66,] -856.532957 -526.236610
[67,] -602.958505 -856.532957
[68,] -457.130721 -602.958505
[69,] -531.752996 -457.130721
[70,] -3174.838918 -531.752996
[71,] 211.927550 -3174.838918
[72,] 9354.448906 211.927550
[73,] -531.752996 9354.448906
[74,] -578.746650 -531.752996
[75,] -581.164370 -578.746650
[76,] -633.399028 -581.164370
[77,] -518.399530 -633.399028
[78,] 1113.396719 -518.399530
[79,] -37.206470 1113.396719
[80,] -529.007476 -37.206470
[81,] -531.752996 -529.007476
[82,] -531.752996 -531.752996
[83,] -554.257189 -531.752996
[84,] -531.752996 -554.257189
[85,] -531.752996 -531.752996
[86,] -521.732620 -531.752996
[87,] -531.449203 -521.732620
[88,] -444.201528 -531.449203
[89,] -670.889173 -444.201528
[90,] -531.752996 -670.889173
[91,] -401.698129 -531.752996
[92,] -1291.243244 -401.698129
[93,] 465.050125 -1291.243244
[94,] -563.965545 465.050125
[95,] -630.510274 -563.965545
[96,] -531.752996 -630.510274
[97,] -514.103829 -531.752996
[98,] -542.910900 -514.103829
[99,] -264.378147 -542.910900
[100,] -1074.155039 -264.378147
[101,] -434.808579 -1074.155039
[102,] -542.900405 -434.808579
[103,] -531.755114 -542.900405
[104,] -531.752996 -531.755114
[105,] 2029.650731 -531.752996
[106,] -431.446382 2029.650731
[107,] -1695.157848 -431.446382
[108,] -384.687406 -1695.157848
[109,] -510.562047 -384.687406
[110,] -518.842890 -510.562047
[111,] 51.592784 -518.842890
[112,] -531.752996 51.592784
[113,] -525.641602 -531.752996
[114,] -531.752996 -525.641602
[115,] -461.320462 -531.752996
[116,] -888.829048 -461.320462
[117,] -1887.590163 -888.829048
[118,] -2138.572831 -1887.590163
[119,] -1490.581032 -2138.572831
[120,] 3.856788 -1490.581032
[121,] -531.752996 3.856788
[122,] -2095.853648 -531.752996
[123,] -878.399845 -2095.853648
[124,] -20.323604 -878.399845
[125,] -531.752996 -20.323604
[126,] -545.554349 -531.752996
[127,] -531.752996 -545.554349
[128,] 1523.730439 -531.752996
[129,] -476.086445 1523.730439
[130,] -585.320451 -476.086445
[131,] -531.752996 -585.320451
[132,] -531.752996 -531.752996
[133,] -531.752996 -531.752996
[134,] -520.235254 -531.752996
[135,] -472.722529 -520.235254
[136,] -2088.524181 -472.722529
[137,] -2592.868414 -2088.524181
[138,] -2423.524181 -2592.868414
[139,] 2101.475819 -2423.524181
[140,] -1404.972548 2101.475819
[141,] -575.483408 -1404.972548
[142,] -2423.524181 -575.483408
[143,] -1785.524181 -2423.524181
[144,] -2457.990409 -1785.524181
[145,] -1816.524181 -2457.990409
[146,] -871.909281 -1816.524181
[147,] -1019.813082 -871.909281
[148,] -2653.563065 -1019.813082
[149,] -2423.524181 -2653.563065
[150,] -1812.524181 -2423.524181
[151,] -2300.277558 -1812.524181
[152,] -1500.524181 -2300.277558
[153,] -1648.256178 -1500.524181
[154,] -64.600804 -1648.256178
[155,] -2586.926486 -64.600804
[156,] -2553.057962 -2586.926486
[157,] -2288.524181 -2553.057962
[158,] -822.180481 -2288.524181
[159,] -4711.001339 -822.180481
[160,] -2412.082568 -4711.001339
[161,] -2114.524181 -2412.082568
[162,] -2176.304013 -2114.524181
[163,] -2596.137187 -2176.304013
[164,] 2095.475819 -2596.137187
[165,] -5116.656277 2095.475819
[166,] -1586.524181 -5116.656277
[167,] 2674.509898 -1586.524181
[168,] -336.323074 2674.509898
[169,] -240.871400 -336.323074
[170,] -3787.369387 -240.871400
[171,] -2440.813076 -3787.369387
[172,] -1905.689883 -2440.813076
[173,] -2410.750338 -1905.689883
[174,] -1495.011601 -2410.750338
[175,] -3495.484878 -1495.011601
[176,] -2362.524181 -3495.484878
[177,] -2052.208190 -2362.524181
[178,] -2431.352844 -2052.208190
[179,] 2477.475819 -2431.352844
[180,] -5103.395872 2477.475819
[181,] -4500.593947 -5103.395872
[182,] -2423.524181 -4500.593947
[183,] -2423.524181 -2423.524181
[184,] -2421.524181 -2423.524181
[185,] -2387.648451 -2421.524181
[186,] -1677.093216 -2387.648451
[187,] 82103.337225 -1677.093216
[188,] 1803.203792 82103.337225
[189,] -2420.524181 1803.203792
[190,] -1896.916885 -2420.524181
[191,] -2511.118916 -1896.916885
[192,] -2018.524181 -2511.118916
[193,] -1517.305540 -2018.524181
[194,] -3044.373271 -1517.305540
[195,] -324.524181 -3044.373271
[196,] -1049.596659 -324.524181
[197,] -1391.588625 -1049.596659
[198,] -2423.524181 -1391.588625
[199,] 19906.475819 -2423.524181
[200,] -4429.102815 19906.475819
[201,] -2423.524181 -4429.102815
[202,] -1940.524181 -2423.524181
[203,] -2400.509309 -1940.524181
[204,] -184.524181 -2400.509309
[205,] -1107.305949 -184.524181
[206,] -3611.207824 -1107.305949
[207,] -2058.524181 -3611.207824
[208,] -1250.547881 -2058.524181
[209,] 19757.475870 -1250.547881
[210,] -3095.283585 19757.475870
[211,] -2685.572436 -3095.283585
[212,] -1892.462246 -2685.572436
[213,] 11022.287859 -1892.462246
[214,] -6015.116649 11022.287859
[215,] -2776.875487 -6015.116649
[216,] -2468.180819 -2776.875487
[217,] -1009.524181 -2468.180819
[218,] -909.398752 -1009.524181
[219,] -5250.609171 -909.398752
[220,] -2423.625909 -5250.609171
[221,] -1579.524181 -2423.625909
[222,] -2068.290876 -1579.524181
[223,] -2009.807846 -2068.290876
[224,] 16285.157028 -2009.807846
[225,] -4485.171685 16285.157028
[226,] 4699.475819 -4485.171685
[227,] -13070.996208 4699.475819
[228,] -2771.687434 -13070.996208
[229,] -366.986238 -2771.687434
[230,] -14203.913927 -366.986238
[231,] 9661.176538 -14203.913927
[232,] -5228.101162 9661.176538
[233,] -2298.524181 -5228.101162
[234,] -1229.369359 -2298.524181
[235,] -1406.605618 -1229.369359
[236,] -2303.724111 -1406.605618
[237,] -2423.524181 -2303.724111
[238,] -2393.524181 -2423.524181
[239,] 818.673002 -2393.524181
[240,] -2909.410634 818.673002
[241,] -3296.891128 -2909.410634
[242,] -1955.524181 -3296.891128
[243,] 2504.134858 -1955.524181
[244,] 1401.321731 2504.134858
[245,] -1525.114460 1401.321731
[246,] -2298.524181 -1525.114460
[247,] -2384.498284 -2298.524181
[248,] -2401.524181 -2384.498284
[249,] -2415.199970 -2401.524181
[250,] -1909.524181 -2415.199970
[251,] -3261.382924 -1909.524181
[252,] -2423.524181 -3261.382924
[253,] -2423.524181 -2423.524181
[254,] -2423.524181 -2423.524181
[255,] -2423.524181 -2423.524181
[256,] -2423.524181 -2423.524181
[257,] 7903.475819 -2423.524181
[258,] 9098.228752 7903.475819
[259,] -384.019079 9098.228752
[260,] 11492.941439 -384.019079
[261,] -738.767412 11492.941439
[262,] -2043.565720 -738.767412
[263,] -2237.191848 -2043.565720
[264,] -2748.599377 -2237.191848
[265,] -2055.524181 -2748.599377
[266,] -2390.231484 -2055.524181
[267,] -2090.524181 -2390.231484
[268,] -2901.371416 -2090.524181
[269,] -2427.143658 -2901.371416
[270,] -2423.524181 -2427.143658
[271,] -1120.524181 -2423.524181
[272,] -2798.025251 -1120.524181
[273,] -48.848634 -2798.025251
[274,] -2503.611189 -48.848634
[275,] -1858.364649 -2503.611189
[276,] -2373.885174 -1858.364649
[277,] -1768.066474 -2373.885174
[278,] -1606.209373 -1768.066474
[279,] -1638.563020 -1606.209373
[280,] -650.680867 -1638.563020
[281,] -2441.671267 -650.680867
[282,] -2106.524181 -2441.671267
[283,] -2519.199175 -2106.524181
[284,] -2148.056984 -2519.199175
[285,] -2377.469561 -2148.056984
[286,] -1487.524181 -2377.469561
[287,] 5549.105588 -1487.524181
[288,] 9047.715220 5549.105588
[289,] 2150.713080 9047.715220
[290,] -1685.524181 2150.713080
[291,] -2440.114833 -1685.524181
[292,] -1831.524181 -2440.114833
[293,] -2347.712027 -1831.524181
[294,] -2431.843668 -2347.712027
[295,] -2424.811446 -2431.843668
[296,] 15590.475819 -2424.811446
[297,] -2710.117007 15590.475819
[298,] -2423.524181 -2710.117007
[299,] 35280.475819 -2423.524181
[300,] -4252.290035 35280.475819
[301,] -1227.341467 -4252.290035
[302,] -4413.234993 -1227.341467
[303,] -2253.250875 -4413.234993
[304,] -2589.665167 -2253.250875
[305,] 4480.022951 -2589.665167
[306,] -5185.573175 4480.022951
[307,] -2423.524181 -5185.573175
[308,] -2422.524181 -2423.524181
[309,] -2423.401116 -2422.524181
[310,] -95.524181 -2423.401116
[311,] -3031.324247 -95.524181
[312,] -2548.173266 -3031.324247
[313,] -2256.007415 -2548.173266
[314,] -825.309520 -2256.007415
[315,] 84708.824961 -825.309520
[316,] -17094.684688 84708.824961
[317,] -2260.448538 -17094.684688
[318,] -1410.524181 -2260.448538
[319,] -2641.648261 -1410.524181
[320,] -1841.926021 -2641.648261
[321,] 2969.730390 -1841.926021
[322,] -4023.560995 2969.730390
[323,] -1411.524181 -4023.560995
[324,] -5029.710777 -1411.524181
[325,] -1547.524181 -5029.710777
[326,] 156470.426224 -1547.524181
[327,] 16817.902759 156470.426224
[328,] 2543.542624 16817.902759
[329,] -2563.394724 2543.542624
[330,] -2799.065069 -2563.394724
[331,] -2423.524181 -2799.065069
[332,] -58.524181 -2423.524181
[333,] -7793.933302 -58.524181
[334,] -3098.755366 -7793.933302
[335,] -838.524181 -3098.755366
[336,] 12183.821560 -838.524181
[337,] -5061.195168 12183.821560
[338,] -2423.524181 -5061.195168
[339,] 8155.475819 -2423.524181
[340,] -4370.342296 8155.475819
[341,] -2403.391098 -4370.342296
[342,] 1218.475819 -2403.391098
[343,] -4172.590994 1218.475819
[344,] -1951.524181 -4172.590994
[345,] -3071.107195 -1951.524181
[346,] 1541.011763 -3071.107195
[347,] -1206.175433 1541.011763
[348,] -1802.524181 -1206.175433
[349,] -2399.833242 -1802.524181
[350,] -2423.524181 -2399.833242
[351,] -2393.524181 -2423.524181
[352,] -2420.427824 -2393.524181
[353,] -1677.524181 -2420.427824
[354,] -2783.410082 -1677.524181
[355,] -2423.524181 -2783.410082
[356,] -2423.524181 -2423.524181
[357,] -2423.524181 -2423.524181
[358,] 1726.475819 -2423.524181
[359,] -1706.538957 1726.475819
[360,] -5724.432563 -1706.538957
[361,] -422.422949 -5724.432563
[362,] -1697.633094 -422.422949
[363,] -3275.123340 -1697.633094
[364,] 10305.927627 -3275.123340
[365,] -3374.945564 10305.927627
[366,] -2354.342898 -3374.945564
[367,] -2407.444889 -2354.342898
[368,] 2685.475819 -2407.444889
[369,] -1840.366319 2685.475819
[370,] -4392.364717 -1840.366319
[371,] -2098.901984 -4392.364717
[372,] 1522.439542 -2098.901984
[373,] -2948.466625 1522.439542
[374,] -1219.976240 -2948.466625
[375,] 422.830876 -1219.976240
[376,] 20735.894586 422.830876
[377,] -5567.945758 20735.894586
[378,] -2560.481408 -5567.945758
[379,] 2204.475819 -2560.481408
[380,] -2785.024837 2204.475819
[381,] -818.779707 -2785.024837
[382,] -2534.295064 -818.779707
[383,] -2547.381315 -2534.295064
[384,] -1443.542997 -2547.381315
[385,] 744.774172 -1443.542997
[386,] -2575.526532 744.774172
[387,] -2423.524181 -2575.526532
[388,] -2423.524181 -2423.524181
[389,] -2423.524181 -2423.524181
[390,] -2423.524181 -2423.524181
[391,] 10829.475819 -2423.524181
[392,] -3827.768280 10829.475819
[393,] -2382.718857 -3827.768280
[394,] -2423.524181 -2382.718857
[395,] -2423.524181 -2423.524181
[396,] -2423.524181 -2423.524181
[397,] -2292.524181 -2423.524181
[398,] -1894.617602 -2292.524181
[399,] -1145.395706 -1894.617602
[400,] 4783.601521 -1145.395706
[401,] -3455.138968 4783.601521
[402,] 15689.390287 -3455.138968
[403,] -5201.251309 15689.390287
[404,] -1636.330421 -5201.251309
[405,] -2203.848003 -1636.330421
[406,] 946.373432 -2203.848003
[407,] 1002.420813 946.373432
[408,] 4344.352928 1002.420813
[409,] -4866.163850 4344.352928
[410,] -2478.097288 -4866.163850
[411,] -2423.524181 -2478.097288
[412,] -1697.524181 -2423.524181
[413,] -1223.995912 -1697.524181
[414,] -2345.939710 -1223.995912
[415,] 4908.645920 -2345.939710
[416,] -7869.476244 4908.645920
[417,] -3174.521760 -7869.476244
[418,] -2423.732586 -3174.521760
[419,] -2327.524181 -2423.732586
[420,] -513.935683 -2327.524181
[421,] 43696.445781 -513.935683
[422,] -14679.127902 43696.445781
[423,] -1089.826295 -14679.127902
[424,] 3833.711605 -1089.826295
[425,] -1170.437279 3833.711605
[426,] -2197.524181 -1170.437279
[427,] -2361.745002 -2197.524181
[428,] -2907.162935 -2361.745002
[429,] -2174.524181 -2907.162935
[430,] -636.699222 -2174.524181
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 -1126.182058 -446.799014
2 -511.543395 -1126.182058
3 -902.728119 -511.543395
4 120.760443 -902.728119
5 -344.107520 120.760443
6 -381.130939 -344.107520
7 -524.381751 -381.130939
8 -660.956188 -524.381751
9 -460.062309 -660.956188
10 -894.317669 -460.062309
11 -904.534116 -894.317669
12 -531.752996 -904.534116
13 -373.326204 -531.752996
14 -494.215890 -373.326204
15 -369.054268 -494.215890
16 -522.363459 -369.054268
17 -541.639199 -522.363459
18 -543.244901 -541.639199
19 -522.622415 -543.244901
20 -531.752996 -522.622415
21 -385.786972 -531.752996
22 -531.752996 -385.786972
23 -1243.621731 -531.752996
24 -743.010377 -1243.621731
25 -515.477844 -743.010377
26 -408.176768 -515.477844
27 -784.862499 -408.176768
28 -531.752996 -784.862499
29 -659.261369 -531.752996
30 91.977173 -659.261369
31 -531.752996 91.977173
32 -590.863508 -531.752996
33 -1726.837388 -590.863508
34 -307.246345 -1726.837388
35 -140.114489 -307.246345
36 -529.762446 -140.114489
37 -180.711314 -529.762446
38 -668.761908 -180.711314
39 -835.720751 -668.761908
40 -1228.182183 -835.720751
41 -781.072067 -1228.182183
42 -914.536718 -781.072067
43 -531.752996 -914.536718
44 -449.831058 -531.752996
45 -436.851090 -449.831058
46 -531.752996 -436.851090
47 -492.797868 -531.752996
48 -646.497911 -492.797868
49 -957.033521 -646.497911
50 -240.473167 -957.033521
51 -701.717963 -240.473167
52 -219.629339 -701.717963
53 -581.289588 -219.629339
54 -3383.558440 -581.289588
55 -766.028236 -3383.558440
56 -511.069047 -766.028236
57 -97.596255 -511.069047
58 -357.909723 -97.596255
59 -531.752996 -357.909723
60 -531.752996 -531.752996
61 -526.747763 -531.752996
62 145.084534 -526.747763
63 -534.244498 145.084534
64 655.246817 -534.244498
65 -526.236610 655.246817
66 -856.532957 -526.236610
67 -602.958505 -856.532957
68 -457.130721 -602.958505
69 -531.752996 -457.130721
70 -3174.838918 -531.752996
71 211.927550 -3174.838918
72 9354.448906 211.927550
73 -531.752996 9354.448906
74 -578.746650 -531.752996
75 -581.164370 -578.746650
76 -633.399028 -581.164370
77 -518.399530 -633.399028
78 1113.396719 -518.399530
79 -37.206470 1113.396719
80 -529.007476 -37.206470
81 -531.752996 -529.007476
82 -531.752996 -531.752996
83 -554.257189 -531.752996
84 -531.752996 -554.257189
85 -531.752996 -531.752996
86 -521.732620 -531.752996
87 -531.449203 -521.732620
88 -444.201528 -531.449203
89 -670.889173 -444.201528
90 -531.752996 -670.889173
91 -401.698129 -531.752996
92 -1291.243244 -401.698129
93 465.050125 -1291.243244
94 -563.965545 465.050125
95 -630.510274 -563.965545
96 -531.752996 -630.510274
97 -514.103829 -531.752996
98 -542.910900 -514.103829
99 -264.378147 -542.910900
100 -1074.155039 -264.378147
101 -434.808579 -1074.155039
102 -542.900405 -434.808579
103 -531.755114 -542.900405
104 -531.752996 -531.755114
105 2029.650731 -531.752996
106 -431.446382 2029.650731
107 -1695.157848 -431.446382
108 -384.687406 -1695.157848
109 -510.562047 -384.687406
110 -518.842890 -510.562047
111 51.592784 -518.842890
112 -531.752996 51.592784
113 -525.641602 -531.752996
114 -531.752996 -525.641602
115 -461.320462 -531.752996
116 -888.829048 -461.320462
117 -1887.590163 -888.829048
118 -2138.572831 -1887.590163
119 -1490.581032 -2138.572831
120 3.856788 -1490.581032
121 -531.752996 3.856788
122 -2095.853648 -531.752996
123 -878.399845 -2095.853648
124 -20.323604 -878.399845
125 -531.752996 -20.323604
126 -545.554349 -531.752996
127 -531.752996 -545.554349
128 1523.730439 -531.752996
129 -476.086445 1523.730439
130 -585.320451 -476.086445
131 -531.752996 -585.320451
132 -531.752996 -531.752996
133 -531.752996 -531.752996
134 -520.235254 -531.752996
135 -472.722529 -520.235254
136 -2088.524181 -472.722529
137 -2592.868414 -2088.524181
138 -2423.524181 -2592.868414
139 2101.475819 -2423.524181
140 -1404.972548 2101.475819
141 -575.483408 -1404.972548
142 -2423.524181 -575.483408
143 -1785.524181 -2423.524181
144 -2457.990409 -1785.524181
145 -1816.524181 -2457.990409
146 -871.909281 -1816.524181
147 -1019.813082 -871.909281
148 -2653.563065 -1019.813082
149 -2423.524181 -2653.563065
150 -1812.524181 -2423.524181
151 -2300.277558 -1812.524181
152 -1500.524181 -2300.277558
153 -1648.256178 -1500.524181
154 -64.600804 -1648.256178
155 -2586.926486 -64.600804
156 -2553.057962 -2586.926486
157 -2288.524181 -2553.057962
158 -822.180481 -2288.524181
159 -4711.001339 -822.180481
160 -2412.082568 -4711.001339
161 -2114.524181 -2412.082568
162 -2176.304013 -2114.524181
163 -2596.137187 -2176.304013
164 2095.475819 -2596.137187
165 -5116.656277 2095.475819
166 -1586.524181 -5116.656277
167 2674.509898 -1586.524181
168 -336.323074 2674.509898
169 -240.871400 -336.323074
170 -3787.369387 -240.871400
171 -2440.813076 -3787.369387
172 -1905.689883 -2440.813076
173 -2410.750338 -1905.689883
174 -1495.011601 -2410.750338
175 -3495.484878 -1495.011601
176 -2362.524181 -3495.484878
177 -2052.208190 -2362.524181
178 -2431.352844 -2052.208190
179 2477.475819 -2431.352844
180 -5103.395872 2477.475819
181 -4500.593947 -5103.395872
182 -2423.524181 -4500.593947
183 -2423.524181 -2423.524181
184 -2421.524181 -2423.524181
185 -2387.648451 -2421.524181
186 -1677.093216 -2387.648451
187 82103.337225 -1677.093216
188 1803.203792 82103.337225
189 -2420.524181 1803.203792
190 -1896.916885 -2420.524181
191 -2511.118916 -1896.916885
192 -2018.524181 -2511.118916
193 -1517.305540 -2018.524181
194 -3044.373271 -1517.305540
195 -324.524181 -3044.373271
196 -1049.596659 -324.524181
197 -1391.588625 -1049.596659
198 -2423.524181 -1391.588625
199 19906.475819 -2423.524181
200 -4429.102815 19906.475819
201 -2423.524181 -4429.102815
202 -1940.524181 -2423.524181
203 -2400.509309 -1940.524181
204 -184.524181 -2400.509309
205 -1107.305949 -184.524181
206 -3611.207824 -1107.305949
207 -2058.524181 -3611.207824
208 -1250.547881 -2058.524181
209 19757.475870 -1250.547881
210 -3095.283585 19757.475870
211 -2685.572436 -3095.283585
212 -1892.462246 -2685.572436
213 11022.287859 -1892.462246
214 -6015.116649 11022.287859
215 -2776.875487 -6015.116649
216 -2468.180819 -2776.875487
217 -1009.524181 -2468.180819
218 -909.398752 -1009.524181
219 -5250.609171 -909.398752
220 -2423.625909 -5250.609171
221 -1579.524181 -2423.625909
222 -2068.290876 -1579.524181
223 -2009.807846 -2068.290876
224 16285.157028 -2009.807846
225 -4485.171685 16285.157028
226 4699.475819 -4485.171685
227 -13070.996208 4699.475819
228 -2771.687434 -13070.996208
229 -366.986238 -2771.687434
230 -14203.913927 -366.986238
231 9661.176538 -14203.913927
232 -5228.101162 9661.176538
233 -2298.524181 -5228.101162
234 -1229.369359 -2298.524181
235 -1406.605618 -1229.369359
236 -2303.724111 -1406.605618
237 -2423.524181 -2303.724111
238 -2393.524181 -2423.524181
239 818.673002 -2393.524181
240 -2909.410634 818.673002
241 -3296.891128 -2909.410634
242 -1955.524181 -3296.891128
243 2504.134858 -1955.524181
244 1401.321731 2504.134858
245 -1525.114460 1401.321731
246 -2298.524181 -1525.114460
247 -2384.498284 -2298.524181
248 -2401.524181 -2384.498284
249 -2415.199970 -2401.524181
250 -1909.524181 -2415.199970
251 -3261.382924 -1909.524181
252 -2423.524181 -3261.382924
253 -2423.524181 -2423.524181
254 -2423.524181 -2423.524181
255 -2423.524181 -2423.524181
256 -2423.524181 -2423.524181
257 7903.475819 -2423.524181
258 9098.228752 7903.475819
259 -384.019079 9098.228752
260 11492.941439 -384.019079
261 -738.767412 11492.941439
262 -2043.565720 -738.767412
263 -2237.191848 -2043.565720
264 -2748.599377 -2237.191848
265 -2055.524181 -2748.599377
266 -2390.231484 -2055.524181
267 -2090.524181 -2390.231484
268 -2901.371416 -2090.524181
269 -2427.143658 -2901.371416
270 -2423.524181 -2427.143658
271 -1120.524181 -2423.524181
272 -2798.025251 -1120.524181
273 -48.848634 -2798.025251
274 -2503.611189 -48.848634
275 -1858.364649 -2503.611189
276 -2373.885174 -1858.364649
277 -1768.066474 -2373.885174
278 -1606.209373 -1768.066474
279 -1638.563020 -1606.209373
280 -650.680867 -1638.563020
281 -2441.671267 -650.680867
282 -2106.524181 -2441.671267
283 -2519.199175 -2106.524181
284 -2148.056984 -2519.199175
285 -2377.469561 -2148.056984
286 -1487.524181 -2377.469561
287 5549.105588 -1487.524181
288 9047.715220 5549.105588
289 2150.713080 9047.715220
290 -1685.524181 2150.713080
291 -2440.114833 -1685.524181
292 -1831.524181 -2440.114833
293 -2347.712027 -1831.524181
294 -2431.843668 -2347.712027
295 -2424.811446 -2431.843668
296 15590.475819 -2424.811446
297 -2710.117007 15590.475819
298 -2423.524181 -2710.117007
299 35280.475819 -2423.524181
300 -4252.290035 35280.475819
301 -1227.341467 -4252.290035
302 -4413.234993 -1227.341467
303 -2253.250875 -4413.234993
304 -2589.665167 -2253.250875
305 4480.022951 -2589.665167
306 -5185.573175 4480.022951
307 -2423.524181 -5185.573175
308 -2422.524181 -2423.524181
309 -2423.401116 -2422.524181
310 -95.524181 -2423.401116
311 -3031.324247 -95.524181
312 -2548.173266 -3031.324247
313 -2256.007415 -2548.173266
314 -825.309520 -2256.007415
315 84708.824961 -825.309520
316 -17094.684688 84708.824961
317 -2260.448538 -17094.684688
318 -1410.524181 -2260.448538
319 -2641.648261 -1410.524181
320 -1841.926021 -2641.648261
321 2969.730390 -1841.926021
322 -4023.560995 2969.730390
323 -1411.524181 -4023.560995
324 -5029.710777 -1411.524181
325 -1547.524181 -5029.710777
326 156470.426224 -1547.524181
327 16817.902759 156470.426224
328 2543.542624 16817.902759
329 -2563.394724 2543.542624
330 -2799.065069 -2563.394724
331 -2423.524181 -2799.065069
332 -58.524181 -2423.524181
333 -7793.933302 -58.524181
334 -3098.755366 -7793.933302
335 -838.524181 -3098.755366
336 12183.821560 -838.524181
337 -5061.195168 12183.821560
338 -2423.524181 -5061.195168
339 8155.475819 -2423.524181
340 -4370.342296 8155.475819
341 -2403.391098 -4370.342296
342 1218.475819 -2403.391098
343 -4172.590994 1218.475819
344 -1951.524181 -4172.590994
345 -3071.107195 -1951.524181
346 1541.011763 -3071.107195
347 -1206.175433 1541.011763
348 -1802.524181 -1206.175433
349 -2399.833242 -1802.524181
350 -2423.524181 -2399.833242
351 -2393.524181 -2423.524181
352 -2420.427824 -2393.524181
353 -1677.524181 -2420.427824
354 -2783.410082 -1677.524181
355 -2423.524181 -2783.410082
356 -2423.524181 -2423.524181
357 -2423.524181 -2423.524181
358 1726.475819 -2423.524181
359 -1706.538957 1726.475819
360 -5724.432563 -1706.538957
361 -422.422949 -5724.432563
362 -1697.633094 -422.422949
363 -3275.123340 -1697.633094
364 10305.927627 -3275.123340
365 -3374.945564 10305.927627
366 -2354.342898 -3374.945564
367 -2407.444889 -2354.342898
368 2685.475819 -2407.444889
369 -1840.366319 2685.475819
370 -4392.364717 -1840.366319
371 -2098.901984 -4392.364717
372 1522.439542 -2098.901984
373 -2948.466625 1522.439542
374 -1219.976240 -2948.466625
375 422.830876 -1219.976240
376 20735.894586 422.830876
377 -5567.945758 20735.894586
378 -2560.481408 -5567.945758
379 2204.475819 -2560.481408
380 -2785.024837 2204.475819
381 -818.779707 -2785.024837
382 -2534.295064 -818.779707
383 -2547.381315 -2534.295064
384 -1443.542997 -2547.381315
385 744.774172 -1443.542997
386 -2575.526532 744.774172
387 -2423.524181 -2575.526532
388 -2423.524181 -2423.524181
389 -2423.524181 -2423.524181
390 -2423.524181 -2423.524181
391 10829.475819 -2423.524181
392 -3827.768280 10829.475819
393 -2382.718857 -3827.768280
394 -2423.524181 -2382.718857
395 -2423.524181 -2423.524181
396 -2423.524181 -2423.524181
397 -2292.524181 -2423.524181
398 -1894.617602 -2292.524181
399 -1145.395706 -1894.617602
400 4783.601521 -1145.395706
401 -3455.138968 4783.601521
402 15689.390287 -3455.138968
403 -5201.251309 15689.390287
404 -1636.330421 -5201.251309
405 -2203.848003 -1636.330421
406 946.373432 -2203.848003
407 1002.420813 946.373432
408 4344.352928 1002.420813
409 -4866.163850 4344.352928
410 -2478.097288 -4866.163850
411 -2423.524181 -2478.097288
412 -1697.524181 -2423.524181
413 -1223.995912 -1697.524181
414 -2345.939710 -1223.995912
415 4908.645920 -2345.939710
416 -7869.476244 4908.645920
417 -3174.521760 -7869.476244
418 -2423.732586 -3174.521760
419 -2327.524181 -2423.732586
420 -513.935683 -2327.524181
421 43696.445781 -513.935683
422 -14679.127902 43696.445781
423 -1089.826295 -14679.127902
424 3833.711605 -1089.826295
425 -1170.437279 3833.711605
426 -2197.524181 -1170.437279
427 -2361.745002 -2197.524181
428 -2907.162935 -2361.745002
429 -2174.524181 -2907.162935
430 -636.699222 -2174.524181
> plot(z,main=paste('Residual Lag plot, lowess, and regression line'), ylab='values of Residuals', xlab='lagged values of Residuals')
> lines(lowess(z))
> abline(lm(z))
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/rcomp/tmp/7qr331292675864.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> acf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Autocorrelation Function')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/rcomp/tmp/8qr331292675864.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> pacf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Partial Autocorrelation Function')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/rcomp/tmp/9ctnj1292675865.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> opar <- par(mfrow = c(2,2), oma = c(0, 0, 1.1, 0))
> plot(mylm, las = 1, sub='Residual Diagnostics')
> par(opar)
> dev.off()
null device
1
> if (n > n25) {
+ postscript(file="/var/www/rcomp/tmp/10ctnj1292675865.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
+ plot(kp3:nmkm3,gqarr[,2], main='Goldfeld-Quandt test',ylab='2-sided p-value',xlab='breakpoint')
+ grid()
+ dev.off()
+ }
null device
1
>
> #Note: the /var/www/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/www/rcomp/createtable")
>
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Estimated Regression Equation', 1, TRUE)
> a<-table.row.end(a)
> myeq <- colnames(x)[1]
> myeq <- paste(myeq, '[t] = ', sep='')
> for (i in 1:k){
+ if (mysum$coefficients[i,1] > 0) myeq <- paste(myeq, '+', '')
+ myeq <- paste(myeq, mysum$coefficients[i,1], sep=' ')
+ if (rownames(mysum$coefficients)[i] != '(Intercept)') {
+ myeq <- paste(myeq, rownames(mysum$coefficients)[i], sep='')
+ if (rownames(mysum$coefficients)[i] != 't') myeq <- paste(myeq, '[t]', sep='')
+ }
+ }
> myeq <- paste(myeq, ' + e[t]')
> a<-table.row.start(a)
> a<-table.element(a, myeq)
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/www/rcomp/tmp/11ftl71292675865.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,hyperlink('http://www.xycoon.com/ols1.htm','Multiple Linear Regression - Ordinary Least Squares',''), 6, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Variable',header=TRUE)
> a<-table.element(a,'Parameter',header=TRUE)
> a<-table.element(a,'S.D.',header=TRUE)
> a<-table.element(a,'T-STAT
H0: parameter = 0',header=TRUE)
> a<-table.element(a,'2-tail p-value',header=TRUE)
> a<-table.element(a,'1-tail p-value',header=TRUE)
> a<-table.row.end(a)
> for (i in 1:k){
+ a<-table.row.start(a)
+ a<-table.element(a,rownames(mysum$coefficients)[i],header=TRUE)
+ a<-table.element(a,mysum$coefficients[i,1])
+ a<-table.element(a, round(mysum$coefficients[i,2],6))
+ a<-table.element(a, round(mysum$coefficients[i,3],4))
+ a<-table.element(a, round(mysum$coefficients[i,4],6))
+ a<-table.element(a, round(mysum$coefficients[i,4]/2,6))
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/www/rcomp/tmp/120c1v1292675865.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Regression Statistics', 2, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple R',1,TRUE)
> a<-table.element(a, sqrt(mysum$r.squared))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'R-squared',1,TRUE)
> a<-table.element(a, mysum$r.squared)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Adjusted R-squared',1,TRUE)
> a<-table.element(a, mysum$adj.r.squared)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (value)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[1])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (DF numerator)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[2])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (DF denominator)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[3])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'p-value',1,TRUE)
> a<-table.element(a, 1-pf(mysum$fstatistic[1],mysum$fstatistic[2],mysum$fstatistic[3]))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Residual Statistics', 2, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Residual Standard Deviation',1,TRUE)
> a<-table.element(a, mysum$sigma)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Sum Squared Residuals',1,TRUE)
> a<-table.element(a, sum(myerror*myerror))
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/www/rcomp/tmp/13pdhp1292675865.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Actuals, Interpolation, and Residuals', 4, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Time or Index', 1, TRUE)
> a<-table.element(a, 'Actuals', 1, TRUE)
> a<-table.element(a, 'Interpolation
Forecast', 1, TRUE)
> a<-table.element(a, 'Residuals
Prediction Error', 1, TRUE)
> a<-table.row.end(a)
> for (i in 1:n) {
+ a<-table.row.start(a)
+ a<-table.element(a,i, 1, TRUE)
+ a<-table.element(a,x[i])
+ a<-table.element(a,x[i]-mysum$resid[i])
+ a<-table.element(a,mysum$resid[i])
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/www/rcomp/tmp/14i4yr1292675865.tab")
> if (n > n25) {
+ a<-table.start()
+ a<-table.row.start(a)
+ a<-table.element(a,'Goldfeld-Quandt test for Heteroskedasticity',4,TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'p-values',header=TRUE)
+ a<-table.element(a,'Alternative Hypothesis',3,header=TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'breakpoint index',header=TRUE)
+ a<-table.element(a,'greater',header=TRUE)
+ a<-table.element(a,'2-sided',header=TRUE)
+ a<-table.element(a,'less',header=TRUE)
+ a<-table.row.end(a)
+ for (mypoint in kp3:nmkm3) {
+ a<-table.row.start(a)
+ a<-table.element(a,mypoint,header=TRUE)
+ a<-table.element(a,gqarr[mypoint-kp3+1,1])
+ a<-table.element(a,gqarr[mypoint-kp3+1,2])
+ a<-table.element(a,gqarr[mypoint-kp3+1,3])
+ a<-table.row.end(a)
+ }
+ a<-table.end(a)
+ table.save(a,file="/var/www/rcomp/tmp/15lnwx1292675865.tab")
+ a<-table.start()
+ a<-table.row.start(a)
+ a<-table.element(a,'Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity',4,TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'Description',header=TRUE)
+ a<-table.element(a,'# significant tests',header=TRUE)
+ a<-table.element(a,'% significant tests',header=TRUE)
+ a<-table.element(a,'OK/NOK',header=TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'1% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant1)
+ a<-table.element(a,numsignificant1/numgqtests)
+ if (numsignificant1/numgqtests < 0.01) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'5% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant5)
+ a<-table.element(a,numsignificant5/numgqtests)
+ if (numsignificant5/numgqtests < 0.05) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'10% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant10)
+ a<-table.element(a,numsignificant10/numgqtests)
+ if (numsignificant10/numgqtests < 0.1) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.end(a)
+ table.save(a,file="/var/www/rcomp/tmp/16zwuo1292675865.tab")
+ }
>
> try(system("convert tmp/1uinc1292675864.ps tmp/1uinc1292675864.png",intern=TRUE))
character(0)
> try(system("convert tmp/2nr4f1292675864.ps tmp/2nr4f1292675864.png",intern=TRUE))
character(0)
> try(system("convert tmp/3nr4f1292675864.ps tmp/3nr4f1292675864.png",intern=TRUE))
character(0)
> try(system("convert tmp/4nr4f1292675864.ps tmp/4nr4f1292675864.png",intern=TRUE))
character(0)
> try(system("convert tmp/5gi401292675864.ps tmp/5gi401292675864.png",intern=TRUE))
character(0)
> try(system("convert tmp/6gi401292675864.ps tmp/6gi401292675864.png",intern=TRUE))
character(0)
> try(system("convert tmp/7qr331292675864.ps tmp/7qr331292675864.png",intern=TRUE))
character(0)
> try(system("convert tmp/8qr331292675864.ps tmp/8qr331292675864.png",intern=TRUE))
character(0)
> try(system("convert tmp/9ctnj1292675865.ps tmp/9ctnj1292675865.png",intern=TRUE))
character(0)
> try(system("convert tmp/10ctnj1292675865.ps tmp/10ctnj1292675865.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
11.790 1.860 13.702