R version 2.9.0 (2009-04-17)
Copyright (C) 2009 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> x <- array(list(162556
+ ,807
+ ,213118
+ ,6282154
+ ,1
+ ,5626
+ ,37
+ ,29790
+ ,444
+ ,81767
+ ,4321023
+ ,2
+ ,13337
+ ,138
+ ,87550
+ ,412
+ ,153198
+ ,4111912
+ ,3
+ ,8541
+ ,45
+ ,84738
+ ,428
+ ,-26007
+ ,223193
+ ,415
+ ,39
+ ,0
+ ,54660
+ ,315
+ ,126942
+ ,1491348
+ ,6
+ ,4276
+ ,24
+ ,42634
+ ,168
+ ,157214
+ ,1629616
+ ,5
+ ,9164
+ ,34
+ ,40949
+ ,263
+ ,129352
+ ,1398893
+ ,9
+ ,2493
+ ,29
+ ,45187
+ ,267
+ ,234817
+ ,1926517
+ ,4
+ ,4891
+ ,38
+ ,37704
+ ,228
+ ,60448
+ ,983660
+ ,14
+ ,1734
+ ,21
+ ,16275
+ ,129
+ ,47818
+ ,1443586
+ ,7
+ ,11409
+ ,76
+ ,25830
+ ,104
+ ,245546
+ ,1073089
+ ,10
+ ,7592
+ ,34
+ ,12679
+ ,122
+ ,48020
+ ,984885
+ ,13
+ ,7135
+ ,62
+ ,18014
+ ,393
+ ,-1710
+ ,1405225
+ ,8
+ ,5043
+ ,67
+ ,43556
+ ,190
+ ,32648
+ ,227132
+ ,414
+ ,110
+ ,1
+ ,24811
+ ,280
+ ,95350
+ ,929118
+ ,15
+ ,1444
+ ,29
+ ,6575
+ ,63
+ ,151352
+ ,1071292
+ ,11
+ ,5480
+ ,133
+ ,7123
+ ,102
+ ,288170
+ ,638830
+ ,28
+ ,4026
+ ,62
+ ,21950
+ ,265
+ ,114337
+ ,856956
+ ,17
+ ,1266
+ ,30
+ ,37597
+ ,234
+ ,37884
+ ,992426
+ ,12
+ ,3195
+ ,21
+ ,17821
+ ,277
+ ,122844
+ ,444477
+ ,46
+ ,655
+ ,14
+ ,12988
+ ,73
+ ,82340
+ ,857217
+ ,16
+ ,5523
+ ,51
+ ,22330
+ ,67
+ ,79801
+ ,711969
+ ,20
+ ,6095
+ ,23
+ ,13326
+ ,103
+ ,165548
+ ,702380
+ ,21
+ ,4925
+ ,38
+ ,16189
+ ,290
+ ,116384
+ ,358589
+ ,89
+ ,538
+ ,10
+ ,7146
+ ,83
+ ,134028
+ ,297978
+ ,367
+ ,933
+ ,14
+ ,15824
+ ,56
+ ,63838
+ ,585715
+ ,30
+ ,6027
+ ,24
+ ,27664
+ ,236
+ ,74996
+ ,657954
+ ,24
+ ,1624
+ ,17
+ ,11920
+ ,73
+ ,31080
+ ,209458
+ ,419
+ ,52
+ ,1
+ ,8568
+ ,34
+ ,32168
+ ,786690
+ ,18
+ ,15856
+ ,68
+ ,14416
+ ,139
+ ,49857
+ ,439798
+ ,48
+ ,664
+ ,17
+ ,3369
+ ,26
+ ,87161
+ ,688779
+ ,23
+ ,17456
+ ,145
+ ,11819
+ ,70
+ ,106113
+ ,574339
+ ,31
+ ,4404
+ ,32
+ ,6984
+ ,40
+ ,80570
+ ,741409
+ ,19
+ ,12031
+ ,78
+ ,4519
+ ,42
+ ,102129
+ ,597793
+ ,29
+ ,8118
+ ,88
+ ,2220
+ ,12
+ ,301670
+ ,644190
+ ,26
+ ,20190
+ ,200
+ ,18562
+ ,211
+ ,102313
+ ,377934
+ ,72
+ ,1148
+ ,10
+ ,10327
+ ,74
+ ,88577
+ ,640273
+ ,27
+ ,4838
+ ,43
+ ,5336
+ ,80
+ ,112477
+ ,697458
+ ,22
+ ,6141
+ ,93
+ ,2365
+ ,83
+ ,191778
+ ,550608
+ ,32
+ ,4438
+ ,148
+ ,4069
+ ,131
+ ,79804
+ ,207393
+ ,420
+ ,51
+ ,2
+ ,8636
+ ,203
+ ,128294
+ ,301607
+ ,359
+ ,119
+ ,12
+ ,13718
+ ,56
+ ,96448
+ ,345783
+ ,100
+ ,2390
+ ,11
+ ,4525
+ ,89
+ ,93811
+ ,501749
+ ,37
+ ,1335
+ ,67
+ ,6869
+ ,88
+ ,117520
+ ,379983
+ ,68
+ ,1714
+ ,26
+ ,4628
+ ,39
+ ,69159
+ ,387475
+ ,64
+ ,3024
+ ,41
+ ,3689
+ ,25
+ ,101792
+ ,377305
+ ,74
+ ,7092
+ ,48
+ ,4891
+ ,49
+ ,210568
+ ,370837
+ ,79
+ ,787
+ ,35
+ ,7489
+ ,149
+ ,136996
+ ,430866
+ ,50
+ ,717
+ ,31
+ ,4901
+ ,58
+ ,121920
+ ,469107
+ ,41
+ ,3204
+ ,55
+ ,2284
+ ,41
+ ,76403
+ ,194493
+ ,423
+ ,-167
+ ,-2
+ ,3160
+ ,90
+ ,108094
+ ,530670
+ ,33
+ ,3062
+ ,105
+ ,4150
+ ,136
+ ,134759
+ ,518365
+ ,35
+ ,2122
+ ,77
+ ,7285
+ ,97
+ ,188873
+ ,491303
+ ,39
+ ,2533
+ ,40
+ ,1134
+ ,63
+ ,146216
+ ,527021
+ ,34
+ ,2019
+ ,288
+ ,4658
+ ,114
+ ,156608
+ ,233773
+ ,411
+ ,214
+ ,7
+ ,2384
+ ,77
+ ,61348
+ ,405972
+ ,58
+ ,2123
+ ,86
+ ,3748
+ ,6
+ ,50350
+ ,652925
+ ,25
+ ,50325
+ ,121
+ ,5371
+ ,47
+ ,87720
+ ,446211
+ ,45
+ ,3730
+ ,46
+ ,1285
+ ,51
+ ,99489
+ ,341340
+ ,106
+ ,1321
+ ,110
+ ,9327
+ ,85
+ ,87419
+ ,387699
+ ,63
+ ,1858
+ ,20
+ ,5565
+ ,43
+ ,94355
+ ,493408
+ ,38
+ ,6243
+ ,53
+ ,1528
+ ,32
+ ,60326
+ ,146494
+ ,426
+ ,-1408
+ ,-35
+ ,3122
+ ,25
+ ,94670
+ ,414462
+ ,54
+ ,6308
+ ,69
+ ,7561
+ ,77
+ ,82425
+ ,364304
+ ,86
+ ,1889
+ ,22
+ ,2675
+ ,54
+ ,59017
+ ,355178
+ ,93
+ ,1964
+ ,58
+ ,13253
+ ,251
+ ,90829
+ ,357760
+ ,90
+ ,167
+ ,12
+ ,880
+ ,15
+ ,80791
+ ,261216
+ ,399
+ ,827
+ ,70
+ ,2053
+ ,44
+ ,100423
+ ,397144
+ ,60
+ ,3720
+ ,96
+ ,1424
+ ,73
+ ,131116
+ ,374943
+ ,77
+ ,1861
+ ,123
+ ,4036
+ ,85
+ ,100269
+ ,424898
+ ,53
+ ,3570
+ ,56
+ ,3045
+ ,49
+ ,27330
+ ,202055
+ ,422
+ ,35
+ ,1
+ ,5119
+ ,38
+ ,39039
+ ,378525
+ ,69
+ ,3643
+ ,35
+ ,1431
+ ,35
+ ,106885
+ ,310768
+ ,330
+ ,3258
+ ,77
+ ,554
+ ,9
+ ,79285
+ ,325738
+ ,125
+ ,11431
+ ,227
+ ,1975
+ ,34
+ ,118881
+ ,394510
+ ,61
+ ,5557
+ ,98
+ ,1765
+ ,20
+ ,77623
+ ,247060
+ ,404
+ ,2353
+ ,27
+ ,1012
+ ,29
+ ,114768
+ ,368078
+ ,82
+ ,3576
+ ,166
+ ,810
+ ,11
+ ,74015
+ ,236761
+ ,409
+ ,855
+ ,45
+ ,1280
+ ,52
+ ,69465
+ ,312378
+ ,326
+ ,960
+ ,88
+ ,666
+ ,13
+ ,117869
+ ,339836
+ ,108
+ ,818
+ ,210
+ ,1380
+ ,29
+ ,60982
+ ,347385
+ ,98
+ ,5669
+ ,107
+ ,4677
+ ,66
+ ,90131
+ ,426280
+ ,52
+ ,3017
+ ,48
+ ,876
+ ,33
+ ,138971
+ ,352850
+ ,96
+ ,2591
+ ,174
+ ,814
+ ,15
+ ,39625
+ ,301881
+ ,358
+ ,5660
+ ,125
+ ,514
+ ,15
+ ,102725
+ ,377516
+ ,73
+ ,11834
+ ,346
+ ,5692
+ ,68
+ ,64239
+ ,357312
+ ,91
+ ,2185
+ ,28
+ ,3642
+ ,100
+ ,90262
+ ,458343
+ ,42
+ ,3004
+ ,71
+ ,540
+ ,13
+ ,103960
+ ,354228
+ ,94
+ ,11016
+ ,286
+ ,2099
+ ,45
+ ,106611
+ ,308636
+ ,339
+ ,1697
+ ,52
+ ,567
+ ,14
+ ,103345
+ ,386212
+ ,65
+ ,16928
+ ,329
+ ,2001
+ ,36
+ ,95551
+ ,393343
+ ,62
+ ,3718
+ ,97
+ ,2949
+ ,40
+ ,82903
+ ,378509
+ ,70
+ ,4354
+ ,61
+ ,2253
+ ,68
+ ,63593
+ ,452469
+ ,43
+ ,2550
+ ,112
+ ,6533
+ ,29
+ ,126910
+ ,364839
+ ,85
+ ,2198
+ ,25
+ ,1889
+ ,43
+ ,37527
+ ,358649
+ ,88
+ ,3526
+ ,84
+ ,3055
+ ,30
+ ,60247
+ ,376641
+ ,75
+ ,4108
+ ,58
+ ,272
+ ,9
+ ,112995
+ ,429112
+ ,51
+ ,28639
+ ,844
+ ,1414
+ ,22
+ ,70184
+ ,330546
+ ,119
+ ,659
+ ,92
+ ,2564
+ ,19
+ ,130140
+ ,403560
+ ,59
+ ,9253
+ ,79
+ ,1383
+ ,9
+ ,73221
+ ,317892
+ ,151
+ ,10717
+ ,85
+ ,1261
+ ,31
+ ,76114
+ ,307528
+ ,343
+ ,3258
+ ,85
+ ,975
+ ,19
+ ,90534
+ ,235133
+ ,410
+ ,1528
+ ,36
+ ,3366
+ ,55
+ ,108479
+ ,299243
+ ,362
+ ,1241
+ ,29
+ ,576
+ ,8
+ ,113761
+ ,314073
+ ,313
+ ,6337
+ ,198
+ ,1686
+ ,28
+ ,68696
+ ,368186
+ ,81
+ ,4205
+ ,100
+ ,746
+ ,29
+ ,71561
+ ,269661
+ ,396
+ ,3029
+ ,93
+ ,3192
+ ,48
+ ,59831
+ ,125390
+ ,427
+ ,-1243
+ ,-23
+ ,2045
+ ,16
+ ,97890
+ ,510834
+ ,36
+ ,15542
+ ,152
+ ,5702
+ ,47
+ ,101481
+ ,321896
+ ,138
+ ,1998
+ ,21
+ ,1932
+ ,20
+ ,72954
+ ,249898
+ ,403
+ ,1386
+ ,26
+ ,936
+ ,22
+ ,67939
+ ,408881
+ ,56
+ ,6963
+ ,223
+ ,3437
+ ,33
+ ,48022
+ ,158492
+ ,425
+ ,-883
+ ,-12
+ ,5131
+ ,44
+ ,86111
+ ,292154
+ ,379
+ ,1298
+ ,18
+ ,2397
+ ,13
+ ,74020
+ ,289513
+ ,385
+ ,6394
+ ,37
+ ,1389
+ ,6
+ ,57530
+ ,378049
+ ,71
+ ,19783
+ ,128
+ ,1503
+ ,35
+ ,56364
+ ,343466
+ ,104
+ ,3679
+ ,95
+ ,402
+ ,8
+ ,84990
+ ,332743
+ ,115
+ ,5106
+ ,330
+ ,2239
+ ,17
+ ,88590
+ ,442882
+ ,47
+ ,11566
+ ,108
+ ,2234
+ ,11
+ ,77200
+ ,214215
+ ,418
+ ,888
+ ,6
+ ,837
+ ,21
+ ,61262
+ ,315688
+ ,165
+ ,1677
+ ,138
+ ,10579
+ ,92
+ ,110309
+ ,375195
+ ,76
+ ,1904
+ ,17
+ ,875
+ ,12
+ ,67000
+ ,334280
+ ,112
+ ,9591
+ ,153
+ ,1585
+ ,112
+ ,93099
+ ,355864
+ ,92
+ ,1457
+ ,98
+ ,1659
+ ,25
+ ,107577
+ ,480382
+ ,40
+ ,9668
+ ,169
+ ,2647
+ ,17
+ ,62920
+ ,353058
+ ,95
+ ,4137
+ ,58
+ ,3294
+ ,23
+ ,75832
+ ,217193
+ ,416
+ ,748
+ ,5
+ ,94
+ ,10
+ ,60793
+ ,314533
+ ,308
+ ,16362
+ ,1212
+ ,422
+ ,23
+ ,57935
+ ,318056
+ ,150
+ ,4216
+ ,280
+ ,34
+ ,7
+ ,60630
+ ,314353
+ ,309
+ ,14294
+ ,3340
+ ,1558
+ ,25
+ ,55637
+ ,369448
+ ,80
+ ,2690
+ ,109
+ ,43
+ ,20
+ ,60887
+ ,312846
+ ,323
+ ,37615
+ ,2613
+ ,645
+ ,4
+ ,60720
+ ,312075
+ ,329
+ ,22415
+ ,174
+ ,316
+ ,4
+ ,60505
+ ,315009
+ ,299
+ ,12779
+ ,364
+ ,115
+ ,10
+ ,60945
+ ,318903
+ ,145
+ ,9146
+ ,1037
+ ,5
+ ,1
+ ,60720
+ ,314887
+ ,304
+ ,57443
+ ,24814
+ ,897
+ ,4
+ ,60720
+ ,314913
+ ,303
+ ,22983
+ ,128
+ ,389
+ ,8
+ ,58990
+ ,325506
+ ,127
+ ,8965
+ ,322
+ ,1002
+ ,11
+ ,56750
+ ,298568
+ ,364
+ ,6571
+ ,98
+ ,36
+ ,4
+ ,60894
+ ,315834
+ ,162
+ ,38611
+ ,3215
+ ,460
+ ,15
+ ,63346
+ ,329784
+ ,122
+ ,8652
+ ,282
+ ,309
+ ,9
+ ,56535
+ ,312878
+ ,322
+ ,10262
+ ,365
+ ,9
+ ,7
+ ,60835
+ ,314987
+ ,300
+ ,19165
+ ,12458
+ ,271
+ ,2
+ ,60720
+ ,325249
+ ,129
+ ,62625
+ ,463
+ ,14
+ ,0
+ ,61016
+ ,315877
+ ,161
+ ,115877
+ ,8041
+ ,520
+ ,7
+ ,58650
+ ,291650
+ ,383
+ ,9165
+ ,176
+ ,1766
+ ,46
+ ,60438
+ ,305959
+ ,349
+ ,1451
+ ,60
+ ,458
+ ,7
+ ,58625
+ ,297765
+ ,368
+ ,8888
+ ,213
+ ,20
+ ,2
+ ,60938
+ ,315245
+ ,297
+ ,38415
+ ,5674
+ ,98
+ ,2
+ ,61490
+ ,315236
+ ,298
+ ,57618
+ ,1179
+ ,405
+ ,5
+ ,60845
+ ,336425
+ ,111
+ ,19489
+ ,337
+ ,483
+ ,7
+ ,60830
+ ,306268
+ ,348
+ ,3936
+ ,220
+ ,454
+ ,24
+ ,63261
+ ,302187
+ ,357
+ ,2004
+ ,225
+ ,47
+ ,1
+ ,60720
+ ,314882
+ ,305
+ ,38294
+ ,2432
+ ,757
+ ,18
+ ,45689
+ ,382712
+ ,66
+ ,9616
+ ,242
+ ,4655
+ ,55
+ ,60720
+ ,341570
+ ,105
+ ,360
+ ,30
+ ,36
+ ,3
+ ,61564
+ ,312412
+ ,325
+ ,28103
+ ,3146
+ ,203
+ ,9
+ ,61938
+ ,309596
+ ,334
+ ,12177
+ ,541
+ ,126
+ ,8
+ ,60951
+ ,315547
+ ,171
+ ,11555
+ ,921
+ ,400
+ ,113
+ ,60720
+ ,313267
+ ,319
+ ,745
+ ,283
+ ,71
+ ,0
+ ,60745
+ ,316176
+ ,159
+ ,116176
+ ,1646
+ ,972
+ ,19
+ ,71642
+ ,359335
+ ,87
+ ,4686
+ ,164
+ ,531
+ ,11
+ ,71641
+ ,330068
+ ,120
+ ,13007
+ ,245
+ ,2461
+ ,25
+ ,55792
+ ,314289
+ ,310
+ ,2005
+ ,46
+ ,378
+ ,16
+ ,71873
+ ,297413
+ ,369
+ ,1873
+ ,258
+ ,23
+ ,5
+ ,62555
+ ,314806
+ ,306
+ ,22961
+ ,4974
+ ,638
+ ,11
+ ,60370
+ ,333210
+ ,114
+ ,9515
+ ,209
+ ,2300
+ ,23
+ ,64873
+ ,352108
+ ,97
+ ,5245
+ ,66
+ ,149
+ ,6
+ ,62041
+ ,313332
+ ,318
+ ,22666
+ ,761
+ ,226
+ ,5
+ ,65745
+ ,291787
+ ,382
+ ,18357
+ ,406
+ ,275
+ ,7
+ ,59500
+ ,318745
+ ,146
+ ,29686
+ ,431
+ ,141
+ ,7
+ ,61630
+ ,315366
+ ,295
+ ,19228
+ ,819
+ ,28
+ ,3
+ ,60890
+ ,315688
+ ,166
+ ,57844
+ ,4148
+ ,4980
+ ,89
+ ,113521
+ ,409642
+ ,55
+ ,2304
+ ,42
+ ,472
+ ,19
+ ,80045
+ ,269587
+ ,397
+ ,3479
+ ,148
+ ,203
+ ,12
+ ,50804
+ ,300962
+ ,361
+ ,3739
+ ,498
+ ,496
+ ,12
+ ,87390
+ ,325479
+ ,128
+ ,7381
+ ,253
+ ,10
+ ,5
+ ,61656
+ ,316155
+ ,160
+ ,58077
+ ,11639
+ ,63
+ ,2
+ ,65688
+ ,318574
+ ,148
+ ,29643
+ ,1891
+ ,1136
+ ,26
+ ,48522
+ ,343613
+ ,103
+ ,4488
+ ,126
+ ,265
+ ,3
+ ,60720
+ ,306948
+ ,346
+ ,3450
+ ,404
+ ,267
+ ,11
+ ,57640
+ ,330059
+ ,121
+ ,4064
+ ,487
+ ,474
+ ,10
+ ,61977
+ ,288985
+ ,386
+ ,4449
+ ,188
+ ,534
+ ,5
+ ,62620
+ ,304485
+ ,351
+ ,14926
+ ,196
+ ,15
+ ,6
+ ,60831
+ ,315688
+ ,164
+ ,14461
+ ,7687
+ ,397
+ ,7
+ ,60646
+ ,317736
+ ,152
+ ,4205
+ ,297
+ ,1866
+ ,28
+ ,56225
+ ,322331
+ ,135
+ ,4218
+ ,66
+ ,288
+ ,3
+ ,60510
+ ,296656
+ ,373
+ ,24164
+ ,335
+ ,3
+ ,1
+ ,60698
+ ,315354
+ ,296
+ ,57677
+ ,36505
+ ,468
+ ,20
+ ,60720
+ ,312161
+ ,328
+ ,5341
+ ,239
+ ,20
+ ,1
+ ,60805
+ ,315576
+ ,169
+ ,57788
+ ,5924
+ ,278
+ ,22
+ ,61404
+ ,314922
+ ,302
+ ,4420
+ ,413
+ ,61
+ ,9
+ ,60720
+ ,314551
+ ,307
+ ,8182
+ ,1874
+ ,192
+ ,2
+ ,65276
+ ,312339
+ ,327
+ ,28085
+ ,586
+ ,317
+ ,7
+ ,63915
+ ,298700
+ ,363
+ ,10967
+ ,311
+ ,738
+ ,9
+ ,60720
+ ,321376
+ ,139
+ ,12138
+ ,164
+ ,368
+ ,13
+ ,61686
+ ,303230
+ ,355
+ ,6072
+ ,281
+ ,2
+ ,0
+ ,60743
+ ,315487
+ ,172
+ ,115487
+ ,52257
+ ,53
+ ,6
+ ,60349
+ ,315793
+ ,163
+ ,19299
+ ,2197
+ ,94
+ ,3
+ ,61360
+ ,312887
+ ,321
+ ,37629
+ ,1200
+ ,24
+ ,7
+ ,59818
+ ,315637
+ ,168
+ ,14455
+ ,4772
+ ,2332
+ ,2
+ ,72680
+ ,324385
+ ,131
+ ,31096
+ ,53
+ ,131
+ ,15
+ ,61808
+ ,308989
+ ,338
+ ,9908
+ ,834
+ ,206
+ ,9
+ ,53110
+ ,296702
+ ,372
+ ,10745
+ ,469
+ ,167
+ ,1
+ ,64245
+ ,307322
+ ,345
+ ,53661
+ ,642
+ ,622
+ ,38
+ ,73007
+ ,304376
+ ,352
+ ,1430
+ ,168
+ ,2328
+ ,57
+ ,82732
+ ,253588
+ ,400
+ ,570
+ ,23
+ ,365
+ ,7
+ ,54820
+ ,309560
+ ,335
+ ,13695
+ ,300
+ ,364
+ ,26
+ ,47705
+ ,298466
+ ,365
+ ,2813
+ ,271
+ ,226
+ ,13
+ ,72835
+ ,343929
+ ,102
+ ,11994
+ ,636
+ ,307
+ ,10
+ ,58856
+ ,331955
+ ,116
+ ,8797
+ ,430
+ ,188
+ ,9
+ ,77655
+ ,381180
+ ,67
+ ,16471
+ ,964
+ ,138
+ ,26
+ ,69817
+ ,331420
+ ,117
+ ,21903
+ ,954
+ ,125
+ ,19
+ ,60798
+ ,310201
+ ,332
+ ,9183
+ ,883
+ ,282
+ ,12
+ ,62452
+ ,320016
+ ,141
+ ,4001
+ ,425
+ ,335
+ ,23
+ ,64175
+ ,320398
+ ,140
+ ,3648
+ ,359
+ ,1324
+ ,29
+ ,67440
+ ,291841
+ ,381
+ ,785
+ ,69
+ ,176
+ ,8
+ ,68136
+ ,310670
+ ,331
+ ,3952
+ ,629
+ ,249
+ ,26
+ ,56726
+ ,313491
+ ,317
+ ,1576
+ ,455
+ ,333
+ ,9
+ ,70811
+ ,331323
+ ,118
+ ,10102
+ ,394
+ ,601
+ ,5
+ ,60720
+ ,319210
+ ,144
+ ,19868
+ ,198
+ ,30
+ ,3
+ ,62045
+ ,318098
+ ,149
+ ,29525
+ ,3912
+ ,249
+ ,13
+ ,54323
+ ,292754
+ ,378
+ ,1496
+ ,373
+ ,165
+ ,12
+ ,62841
+ ,325176
+ ,130
+ ,5216
+ ,759
+ ,453
+ ,19
+ ,81125
+ ,365959
+ ,84
+ ,7903
+ ,366
+ ,53
+ ,10
+ ,59506
+ ,302409
+ ,356
+ ,7315
+ ,1932
+ ,382
+ ,9
+ ,59365
+ ,340968
+ ,107
+ ,6713
+ ,369
+ ,30
+ ,4
+ ,60798
+ ,313164
+ ,320
+ ,28291
+ ,3808
+ ,290
+ ,1
+ ,58790
+ ,301164
+ ,360
+ ,50582
+ ,348
+ ,366
+ ,14
+ ,61808
+ ,344425
+ ,101
+ ,2725
+ ,394
+ ,2
+ ,12
+ ,60735
+ ,315394
+ ,174
+ ,12822
+ ,63057
+ ,209
+ ,19
+ ,64016
+ ,316647
+ ,157
+ ,8973
+ ,558
+ ,384
+ ,17
+ ,54683
+ ,309836
+ ,333
+ ,4993
+ ,286
+ ,365
+ ,32
+ ,87192
+ ,346611
+ ,99
+ ,1766
+ ,402
+ ,49
+ ,14
+ ,64107
+ ,322031
+ ,137
+ ,15254
+ ,2500
+ ,3
+ ,8
+ ,60761
+ ,315656
+ ,167
+ ,28914
+ ,37920
+ ,133
+ ,4
+ ,65990
+ ,339445
+ ,109
+ ,9960
+ ,1052
+ ,32
+ ,0
+ ,59988
+ ,314964
+ ,301
+ ,114964
+ ,3603
+ ,368
+ ,20
+ ,61167
+ ,297141
+ ,370
+ ,5714
+ ,264
+ ,1
+ ,5
+ ,60719
+ ,315372
+ ,293
+ ,19229
+ ,84213
+ ,22
+ ,1
+ ,60722
+ ,312502
+ ,324
+ ,56251
+ ,5154
+ ,96
+ ,4
+ ,60379
+ ,313729
+ ,316
+ ,22746
+ ,1188
+ ,1
+ ,1
+ ,60727
+ ,315388
+ ,175
+ ,57694
+ ,117743
+ ,314
+ ,4
+ ,60720
+ ,315371
+ ,294
+ ,23074
+ ,367
+ ,844
+ ,20
+ ,60925
+ ,296139
+ ,375
+ ,1233
+ ,114
+ ,26
+ ,1
+ ,60896
+ ,313880
+ ,315
+ ,113880
+ ,4328
+ ,125
+ ,10
+ ,59734
+ ,317698
+ ,153
+ ,9054
+ ,944
+ ,304
+ ,12
+ ,62969
+ ,295580
+ ,376
+ ,6372
+ ,314
+ ,621
+ ,13
+ ,60720
+ ,308256
+ ,341
+ ,2255
+ ,174
+ ,119
+ ,3
+ ,59118
+ ,303677
+ ,354
+ ,17280
+ ,872
+ ,1595
+ ,10
+ ,60720
+ ,319369
+ ,143
+ ,7022
+ ,75
+ ,312
+ ,3
+ ,58598
+ ,318690
+ ,147
+ ,8478
+ ,381
+ ,60
+ ,7
+ ,61124
+ ,314049
+ ,314
+ ,11405
+ ,1916
+ ,587
+ ,10
+ ,59595
+ ,325699
+ ,126
+ ,10475
+ ,214
+ ,135
+ ,1
+ ,62065
+ ,314210
+ ,311
+ ,57105
+ ,847
+ ,514
+ ,15
+ ,78780
+ ,322378
+ ,133
+ ,2353
+ ,238
+ ,1
+ ,4
+ ,60722
+ ,315398
+ ,173
+ ,28849
+ ,172235
+ ,1763
+ ,28
+ ,61600
+ ,308336
+ ,340
+ ,4514
+ ,61
+ ,180
+ ,9
+ ,59635
+ ,316386
+ ,158
+ ,10581
+ ,647
+ ,218
+ ,7
+ ,60720
+ ,315553
+ ,170
+ ,5503
+ ,531
+ ,448
+ ,7
+ ,59781
+ ,323361
+ ,132
+ ,3084
+ ,275
+ ,227
+ ,7
+ ,76644
+ ,336639
+ ,110
+ ,15182
+ ,603
+ ,174
+ ,3
+ ,64820
+ ,307424
+ ,344
+ ,107424
+ ,618
+ ,121
+ ,11
+ ,56178
+ ,295370
+ ,377
+ ,3974
+ ,788
+ ,607
+ ,7
+ ,60436
+ ,322340
+ ,134
+ ,11122
+ ,201
+ ,2212
+ ,10
+ ,60720
+ ,319864
+ ,142
+ ,8562
+ ,54
+ ,530
+ ,18
+ ,73433
+ ,317291
+ ,155
+ ,1955
+ ,221
+ ,571
+ ,14
+ ,41477
+ ,280398
+ ,390
+ ,1005
+ ,141
+ ,78
+ ,12
+ ,62700
+ ,317330
+ ,154
+ ,7333
+ ,1511
+ ,2489
+ ,29
+ ,67804
+ ,238125
+ ,408
+ ,953
+ ,15
+ ,131
+ ,3
+ ,59661
+ ,327071
+ ,123
+ ,21179
+ ,969
+ ,923
+ ,6
+ ,58620
+ ,309038
+ ,337
+ ,13630
+ ,118
+ ,72
+ ,3
+ ,60398
+ ,314210
+ ,312
+ ,38070
+ ,1591
+ ,572
+ ,8
+ ,58580
+ ,307930
+ ,342
+ ,6746
+ ,189
+ ,397
+ ,10
+ ,62710
+ ,322327
+ ,136
+ ,12233
+ ,308
+ ,450
+ ,6
+ ,59325
+ ,292136
+ ,380
+ ,11517
+ ,205
+ ,622
+ ,8
+ ,60950
+ ,263276
+ ,398
+ ,9039
+ ,102
+ ,694
+ ,6
+ ,68060
+ ,367655
+ ,83
+ ,20957
+ ,241
+ ,3425
+ ,9
+ ,83620
+ ,283910
+ ,388
+ ,6993
+ ,25
+ ,562
+ ,8
+ ,58456
+ ,283587
+ ,389
+ ,6430
+ ,149
+ ,4917
+ ,26
+ ,52811
+ ,243650
+ ,406
+ ,1039
+ ,9
+ ,1442
+ ,239
+ ,121173
+ ,438493
+ ,49
+ ,2021
+ ,165
+ ,529
+ ,7
+ ,63870
+ ,296261
+ ,374
+ ,10696
+ ,182
+ ,2126
+ ,41
+ ,21001
+ ,230621
+ ,413
+ ,222
+ ,14
+ ,1061
+ ,3
+ ,70415
+ ,304252
+ ,353
+ ,20850
+ ,98
+ ,776
+ ,8
+ ,64230
+ ,333505
+ ,113
+ ,14834
+ ,172
+ ,611
+ ,6
+ ,59190
+ ,296919
+ ,371
+ ,12115
+ ,159
+ ,1526
+ ,21
+ ,69351
+ ,278990
+ ,391
+ ,3160
+ ,52
+ ,592
+ ,7
+ ,64270
+ ,276898
+ ,392
+ ,10985
+ ,130
+ ,1182
+ ,11
+ ,70694
+ ,327007
+ ,124
+ ,9770
+ ,107
+ ,621
+ ,11
+ ,68005
+ ,317046
+ ,156
+ ,7315
+ ,188
+ ,989
+ ,12
+ ,58930
+ ,304555
+ ,350
+ ,9505
+ ,106
+ ,438
+ ,9
+ ,58320
+ ,298096
+ ,366
+ ,8918
+ ,224
+ ,726
+ ,3
+ ,69980
+ ,231861
+ ,412
+ ,10620
+ ,44
+ ,1303
+ ,57
+ ,69863
+ ,309422
+ ,336
+ ,1794
+ ,84
+ ,7419
+ ,21
+ ,63255
+ ,286963
+ ,387
+ ,2999
+ ,12
+ ,1164
+ ,15
+ ,57320
+ ,269753
+ ,395
+ ,4103
+ ,60
+ ,3310
+ ,32
+ ,75230
+ ,448243
+ ,44
+ ,7523
+ ,75
+ ,1920
+ ,11
+ ,79420
+ ,165404
+ ,424
+ ,-2306
+ ,-18
+ ,965
+ ,2
+ ,73490
+ ,204325
+ ,421
+ ,1442
+ ,4
+ ,3256
+ ,23
+ ,35250
+ ,407159
+ ,57
+ ,3139
+ ,64
+ ,1135
+ ,20
+ ,62285
+ ,290476
+ ,384
+ ,5322
+ ,80
+ ,1270
+ ,24
+ ,69206
+ ,275311
+ ,393
+ ,2897
+ ,59
+ ,661
+ ,1
+ ,65920
+ ,246541
+ ,405
+ ,15514
+ ,70
+ ,1013
+ ,1
+ ,69770
+ ,253468
+ ,402
+ ,26734
+ ,53
+ ,2844
+ ,74
+ ,72683
+ ,240897
+ ,407
+ ,610
+ ,14
+ ,11528
+ ,68
+ ,-14545
+ ,-83265
+ ,431
+ ,-4047
+ ,-25
+ ,6526
+ ,20
+ ,55830
+ ,-42143
+ ,430
+ ,-9313
+ ,-37
+ ,2264
+ ,20
+ ,55174
+ ,272713
+ ,394
+ ,3030
+ ,32
+ ,5109
+ ,82
+ ,67038
+ ,215362
+ ,417
+ ,158
+ ,3
+ ,3999
+ ,21
+ ,51252
+ ,42754
+ ,428
+ ,-5242
+ ,-39
+ ,35624
+ ,244
+ ,157278
+ ,306275
+ ,347
+ ,477
+ ,3
+ ,9252
+ ,32
+ ,79510
+ ,253537
+ ,401
+ ,1115
+ ,6
+ ,15236
+ ,86
+ ,77440
+ ,372631
+ ,78
+ ,1918
+ ,11
+ ,18073
+ ,69
+ ,27284
+ ,-7170
+ ,429
+ ,-1151
+ ,-11)
+ ,dim=c(7
+ ,314)
+ ,dimnames=list(c('Costs'
+ ,'Orders'
+ ,'Dividends'
+ ,'Wealth'
+ ,'Wrank'
+ ,'Profit/Trades'
+ ,'Profit/Cost')
+ ,1:314))
> y <- array(NA,dim=c(7,314),dimnames=list(c('Costs','Orders','Dividends','Wealth','Wrank','Profit/Trades','Profit/Cost'),1:314))
> for (i in 1:dim(x)[1])
+ {
+ for (j in 1:dim(x)[2])
+ {
+ y[i,j] <- as.numeric(x[i,j])
+ }
+ }
> par3 = 'No Linear Trend'
> par2 = 'Do not include Seasonal Dummies'
> par1 = '4'
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Prof. Dr. P. Wessa
> #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #Technical description: Write here your technical program description (don't use hard returns!)
> library(lattice)
> library(lmtest)
Loading required package: zoo
Attaching package: 'zoo'
The following object(s) are masked from package:base :
as.Date.numeric
> n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test
> par1 <- as.numeric(par1)
> x <- t(y)
> k <- length(x[1,])
> n <- length(x[,1])
> x1 <- cbind(x[,par1], x[,1:k!=par1])
> mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1])
> colnames(x1) <- mycolnames #colnames(x)[par1]
> x <- x1
> if (par3 == 'First Differences'){
+ x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep='')))
+ for (i in 1:n-1) {
+ for (j in 1:k) {
+ x2[i,j] <- x[i+1,j] - x[i,j]
+ }
+ }
+ x <- x2
+ }
> if (par2 == 'Include Monthly Dummies'){
+ x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep ='')))
+ for (i in 1:11){
+ x2[seq(i,n,12),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> if (par2 == 'Include Quarterly Dummies'){
+ x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep ='')))
+ for (i in 1:3){
+ x2[seq(i,n,4),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> k <- length(x[1,])
> if (par3 == 'Linear Trend'){
+ x <- cbind(x, c(1:n))
+ colnames(x)[k+1] <- 't'
+ }
> x
Wealth Costs Orders Dividends Wrank Profit/Trades Profit/Cost
1 6282154 162556 807 213118 1 5626 37
2 4321023 29790 444 81767 2 13337 138
3 4111912 87550 412 153198 3 8541 45
4 223193 84738 428 -26007 415 39 0
5 1491348 54660 315 126942 6 4276 24
6 1629616 42634 168 157214 5 9164 34
7 1398893 40949 263 129352 9 2493 29
8 1926517 45187 267 234817 4 4891 38
9 983660 37704 228 60448 14 1734 21
10 1443586 16275 129 47818 7 11409 76
11 1073089 25830 104 245546 10 7592 34
12 984885 12679 122 48020 13 7135 62
13 1405225 18014 393 -1710 8 5043 67
14 227132 43556 190 32648 414 110 1
15 929118 24811 280 95350 15 1444 29
16 1071292 6575 63 151352 11 5480 133
17 638830 7123 102 288170 28 4026 62
18 856956 21950 265 114337 17 1266 30
19 992426 37597 234 37884 12 3195 21
20 444477 17821 277 122844 46 655 14
21 857217 12988 73 82340 16 5523 51
22 711969 22330 67 79801 20 6095 23
23 702380 13326 103 165548 21 4925 38
24 358589 16189 290 116384 89 538 10
25 297978 7146 83 134028 367 933 14
26 585715 15824 56 63838 30 6027 24
27 657954 27664 236 74996 24 1624 17
28 209458 11920 73 31080 419 52 1
29 786690 8568 34 32168 18 15856 68
30 439798 14416 139 49857 48 664 17
31 688779 3369 26 87161 23 17456 145
32 574339 11819 70 106113 31 4404 32
33 741409 6984 40 80570 19 12031 78
34 597793 4519 42 102129 29 8118 88
35 644190 2220 12 301670 26 20190 200
36 377934 18562 211 102313 72 1148 10
37 640273 10327 74 88577 27 4838 43
38 697458 5336 80 112477 22 6141 93
39 550608 2365 83 191778 32 4438 148
40 207393 4069 131 79804 420 51 2
41 301607 8636 203 128294 359 119 12
42 345783 13718 56 96448 100 2390 11
43 501749 4525 89 93811 37 1335 67
44 379983 6869 88 117520 68 1714 26
45 387475 4628 39 69159 64 3024 41
46 377305 3689 25 101792 74 7092 48
47 370837 4891 49 210568 79 787 35
48 430866 7489 149 136996 50 717 31
49 469107 4901 58 121920 41 3204 55
50 194493 2284 41 76403 423 -167 -2
51 530670 3160 90 108094 33 3062 105
52 518365 4150 136 134759 35 2122 77
53 491303 7285 97 188873 39 2533 40
54 527021 1134 63 146216 34 2019 288
55 233773 4658 114 156608 411 214 7
56 405972 2384 77 61348 58 2123 86
57 652925 3748 6 50350 25 50325 121
58 446211 5371 47 87720 45 3730 46
59 341340 1285 51 99489 106 1321 110
60 387699 9327 85 87419 63 1858 20
61 493408 5565 43 94355 38 6243 53
62 146494 1528 32 60326 426 -1408 -35
63 414462 3122 25 94670 54 6308 69
64 364304 7561 77 82425 86 1889 22
65 355178 2675 54 59017 93 1964 58
66 357760 13253 251 90829 90 167 12
67 261216 880 15 80791 399 827 70
68 397144 2053 44 100423 60 3720 96
69 374943 1424 73 131116 77 1861 123
70 424898 4036 85 100269 53 3570 56
71 202055 3045 49 27330 422 35 1
72 378525 5119 38 39039 69 3643 35
73 310768 1431 35 106885 330 3258 77
74 325738 554 9 79285 125 11431 227
75 394510 1975 34 118881 61 5557 98
76 247060 1765 20 77623 404 2353 27
77 368078 1012 29 114768 82 3576 166
78 236761 810 11 74015 409 855 45
79 312378 1280 52 69465 326 960 88
80 339836 666 13 117869 108 818 210
81 347385 1380 29 60982 98 5669 107
82 426280 4677 66 90131 52 3017 48
83 352850 876 33 138971 96 2591 174
84 301881 814 15 39625 358 5660 125
85 377516 514 15 102725 73 11834 346
86 357312 5692 68 64239 91 2185 28
87 458343 3642 100 90262 42 3004 71
88 354228 540 13 103960 94 11016 286
89 308636 2099 45 106611 339 1697 52
90 386212 567 14 103345 65 16928 329
91 393343 2001 36 95551 62 3718 97
92 378509 2949 40 82903 70 4354 61
93 452469 2253 68 63593 43 2550 112
94 364839 6533 29 126910 85 2198 25
95 358649 1889 43 37527 88 3526 84
96 376641 3055 30 60247 75 4108 58
97 429112 272 9 112995 51 28639 844
98 330546 1414 22 70184 119 659 92
99 403560 2564 19 130140 59 9253 79
100 317892 1383 9 73221 151 10717 85
101 307528 1261 31 76114 343 3258 85
102 235133 975 19 90534 410 1528 36
103 299243 3366 55 108479 362 1241 29
104 314073 576 8 113761 313 6337 198
105 368186 1686 28 68696 81 4205 100
106 269661 746 29 71561 396 3029 93
107 125390 3192 48 59831 427 -1243 -23
108 510834 2045 16 97890 36 15542 152
109 321896 5702 47 101481 138 1998 21
110 249898 1932 20 72954 403 1386 26
111 408881 936 22 67939 56 6963 223
112 158492 3437 33 48022 425 -883 -12
113 292154 5131 44 86111 379 1298 18
114 289513 2397 13 74020 385 6394 37
115 378049 1389 6 57530 71 19783 128
116 343466 1503 35 56364 104 3679 95
117 332743 402 8 84990 115 5106 330
118 442882 2239 17 88590 47 11566 108
119 214215 2234 11 77200 418 888 6
120 315688 837 21 61262 165 1677 138
121 375195 10579 92 110309 76 1904 17
122 334280 875 12 67000 112 9591 153
123 355864 1585 112 93099 92 1457 98
124 480382 1659 25 107577 40 9668 169
125 353058 2647 17 62920 95 4137 58
126 217193 3294 23 75832 416 748 5
127 314533 94 10 60793 308 16362 1212
128 318056 422 23 57935 150 4216 280
129 314353 34 7 60630 309 14294 3340
130 369448 1558 25 55637 80 2690 109
131 312846 43 20 60887 323 37615 2613
132 312075 645 4 60720 329 22415 174
133 315009 316 4 60505 299 12779 364
134 318903 115 10 60945 145 9146 1037
135 314887 5 1 60720 304 57443 24814
136 314913 897 4 60720 303 22983 128
137 325506 389 8 58990 127 8965 322
138 298568 1002 11 56750 364 6571 98
139 315834 36 4 60894 162 38611 3215
140 329784 460 15 63346 122 8652 282
141 312878 309 9 56535 322 10262 365
142 314987 9 7 60835 300 19165 12458
143 325249 271 2 60720 129 62625 463
144 315877 14 0 61016 161 115877 8041
145 291650 520 7 58650 383 9165 176
146 305959 1766 46 60438 349 1451 60
147 297765 458 7 58625 368 8888 213
148 315245 20 2 60938 297 38415 5674
149 315236 98 2 61490 298 57618 1179
150 336425 405 5 60845 111 19489 337
151 306268 483 7 60830 348 3936 220
152 302187 454 24 63261 357 2004 225
153 314882 47 1 60720 305 38294 2432
154 382712 757 18 45689 66 9616 242
155 341570 4655 55 60720 105 360 30
156 312412 36 3 61564 325 28103 3146
157 309596 203 9 61938 334 12177 541
158 315547 126 8 60951 171 11555 921
159 313267 400 113 60720 319 745 283
160 316176 71 0 60745 159 116176 1646
161 359335 972 19 71642 87 4686 164
162 330068 531 11 71641 120 13007 245
163 314289 2461 25 55792 310 2005 46
164 297413 378 16 71873 369 1873 258
165 314806 23 5 62555 306 22961 4974
166 333210 638 11 60370 114 9515 209
167 352108 2300 23 64873 97 5245 66
168 313332 149 6 62041 318 22666 761
169 291787 226 5 65745 382 18357 406
170 318745 275 7 59500 146 29686 431
171 315366 141 7 61630 295 19228 819
172 315688 28 3 60890 166 57844 4148
173 409642 4980 89 113521 55 2304 42
174 269587 472 19 80045 397 3479 148
175 300962 203 12 50804 361 3739 498
176 325479 496 12 87390 128 7381 253
177 316155 10 5 61656 160 58077 11639
178 318574 63 2 65688 148 29643 1891
179 343613 1136 26 48522 103 4488 126
180 306948 265 3 60720 346 3450 404
181 330059 267 11 57640 121 4064 487
182 288985 474 10 61977 386 4449 188
183 304485 534 5 62620 351 14926 196
184 315688 15 6 60831 164 14461 7687
185 317736 397 7 60646 152 4205 297
186 322331 1866 28 56225 135 4218 66
187 296656 288 3 60510 373 24164 335
188 315354 3 1 60698 296 57677 36505
189 312161 468 20 60720 328 5341 239
190 315576 20 1 60805 169 57788 5924
191 314922 278 22 61404 302 4420 413
192 314551 61 9 60720 307 8182 1874
193 312339 192 2 65276 327 28085 586
194 298700 317 7 63915 363 10967 311
195 321376 738 9 60720 139 12138 164
196 303230 368 13 61686 355 6072 281
197 315487 2 0 60743 172 115487 52257
198 315793 53 6 60349 163 19299 2197
199 312887 94 3 61360 321 37629 1200
200 315637 24 7 59818 168 14455 4772
201 324385 2332 2 72680 131 31096 53
202 308989 131 15 61808 338 9908 834
203 296702 206 9 53110 372 10745 469
204 307322 167 1 64245 345 53661 642
205 304376 622 38 73007 352 1430 168
206 253588 2328 57 82732 400 570 23
207 309560 365 7 54820 335 13695 300
208 298466 364 26 47705 365 2813 271
209 343929 226 13 72835 102 11994 636
210 331955 307 10 58856 116 8797 430
211 381180 188 9 77655 67 16471 964
212 331420 138 26 69817 117 21903 954
213 310201 125 19 60798 332 9183 883
214 320016 282 12 62452 141 4001 425
215 320398 335 23 64175 140 3648 359
216 291841 1324 29 67440 381 785 69
217 310670 176 8 68136 331 3952 629
218 313491 249 26 56726 317 1576 455
219 331323 333 9 70811 118 10102 394
220 319210 601 5 60720 144 19868 198
221 318098 30 3 62045 149 29525 3912
222 292754 249 13 54323 378 1496 373
223 325176 165 12 62841 130 5216 759
224 365959 453 19 81125 84 7903 366
225 302409 53 10 59506 356 7315 1932
226 340968 382 9 59365 107 6713 369
227 313164 30 4 60798 320 28291 3808
228 301164 290 1 58790 360 50582 348
229 344425 366 14 61808 101 2725 394
230 315394 2 12 60735 174 12822 63057
231 316647 209 19 64016 157 8973 558
232 309836 384 17 54683 333 4993 286
233 346611 365 32 87192 99 1766 402
234 322031 49 14 64107 137 15254 2500
235 315656 3 8 60761 167 28914 37920
236 339445 133 4 65990 109 9960 1052
237 314964 32 0 59988 301 114964 3603
238 297141 368 20 61167 370 5714 264
239 315372 1 5 60719 293 19229 84213
240 312502 22 1 60722 324 56251 5154
241 313729 96 4 60379 316 22746 1188
242 315388 1 1 60727 175 57694 117743
243 315371 314 4 60720 294 23074 367
244 296139 844 20 60925 375 1233 114
245 313880 26 1 60896 315 113880 4328
246 317698 125 10 59734 153 9054 944
247 295580 304 12 62969 376 6372 314
248 308256 621 13 60720 341 2255 174
249 303677 119 3 59118 354 17280 872
250 319369 1595 10 60720 143 7022 75
251 318690 312 3 58598 147 8478 381
252 314049 60 7 61124 314 11405 1916
253 325699 587 10 59595 126 10475 214
254 314210 135 1 62065 311 57105 847
255 322378 514 15 78780 133 2353 238
256 315398 1 4 60722 173 28849 172235
257 308336 1763 28 61600 340 4514 61
258 316386 180 9 59635 158 10581 647
259 315553 218 7 60720 170 5503 531
260 323361 448 7 59781 132 3084 275
261 336639 227 7 76644 110 15182 603
262 307424 174 3 64820 344 107424 618
263 295370 121 11 56178 377 3974 788
264 322340 607 7 60436 134 11122 201
265 319864 2212 10 60720 142 8562 54
266 317291 530 18 73433 155 1955 221
267 280398 571 14 41477 390 1005 141
268 317330 78 12 62700 154 7333 1511
269 238125 2489 29 67804 408 953 15
270 327071 131 3 59661 123 21179 969
271 309038 923 6 58620 337 13630 118
272 314210 72 3 60398 312 38070 1591
273 307930 572 8 58580 342 6746 189
274 322327 397 10 62710 136 12233 308
275 292136 450 6 59325 380 11517 205
276 263276 622 8 60950 398 9039 102
277 367655 694 6 68060 83 20957 241
278 283910 3425 9 83620 388 6993 25
279 283587 562 8 58456 389 6430 149
280 243650 4917 26 52811 406 1039 9
281 438493 1442 239 121173 49 2021 165
282 296261 529 7 63870 374 10696 182
283 230621 2126 41 21001 413 222 14
284 304252 1061 3 70415 353 20850 98
285 333505 776 8 64230 113 14834 172
286 296919 611 6 59190 371 12115 159
287 278990 1526 21 69351 391 3160 52
288 276898 592 7 64270 392 10985 130
289 327007 1182 11 70694 124 9770 107
290 317046 621 11 68005 156 7315 188
291 304555 989 12 58930 350 9505 106
292 298096 438 9 58320 366 8918 224
293 231861 726 3 69980 412 10620 44
294 309422 1303 57 69863 336 1794 84
295 286963 7419 21 63255 387 2999 12
296 269753 1164 15 57320 395 4103 60
297 448243 3310 32 75230 44 7523 75
298 165404 1920 11 79420 424 -2306 -18
299 204325 965 2 73490 421 1442 4
300 407159 3256 23 35250 57 3139 64
301 290476 1135 20 62285 384 5322 80
302 275311 1270 24 69206 393 2897 59
303 246541 661 1 65920 405 15514 70
304 253468 1013 1 69770 402 26734 53
305 240897 2844 74 72683 407 610 14
306 -83265 11528 68 -14545 431 -4047 -25
307 -42143 6526 20 55830 430 -9313 -37
308 272713 2264 20 55174 394 3030 32
309 215362 5109 82 67038 417 158 3
310 42754 3999 21 51252 428 -5242 -39
311 306275 35624 244 157278 347 477 3
312 253537 9252 32 79510 401 1115 6
313 372631 15236 86 77440 78 1918 11
314 -7170 18073 69 27284 429 -1151 -11
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) Costs Orders Dividends
2.439e+05 2.118e+01 1.199e+03 1.173e+00
Wrank `Profit/Trades` `Profit/Cost`
-5.148e+02 2.872e+00 -1.753e-01
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-2084202 -39915 20558 84946 2780931
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 2.439e+05 5.786e+04 4.216 3.28e-05 ***
Costs 2.118e+01 2.450e+00 8.645 3.04e-16 ***
Orders 1.199e+03 4.155e+02 2.885 0.00419 **
Dividends 1.173e+00 4.901e-01 2.393 0.01733 *
Wrank -5.148e+02 1.240e+02 -4.153 4.26e-05 ***
`Profit/Trades` 2.872e+00 8.971e-01 3.201 0.00151 **
`Profit/Cost` -1.753e-01 1.201e+00 -0.146 0.88406
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 282400 on 307 degrees of freedom
Multiple R-squared: 0.6892, Adjusted R-squared: 0.6831
F-statistic: 113.4 on 6 and 307 DF, p-value: < 2.2e-16
> if (n > n25) {
+ kp3 <- k + 3
+ nmkm3 <- n - k - 3
+ gqarr <- array(NA, dim=c(nmkm3-kp3+1,3))
+ numgqtests <- 0
+ numsignificant1 <- 0
+ numsignificant5 <- 0
+ numsignificant10 <- 0
+ for (mypoint in kp3:nmkm3) {
+ j <- 0
+ numgqtests <- numgqtests + 1
+ for (myalt in c('greater', 'two.sided', 'less')) {
+ j <- j + 1
+ gqarr[mypoint-kp3+1,j] <- gqtest(mylm, point=mypoint, alternative=myalt)$p.value
+ }
+ if (gqarr[mypoint-kp3+1,2] < 0.01) numsignificant1 <- numsignificant1 + 1
+ if (gqarr[mypoint-kp3+1,2] < 0.05) numsignificant5 <- numsignificant5 + 1
+ if (gqarr[mypoint-kp3+1,2] < 0.10) numsignificant10 <- numsignificant10 + 1
+ }
+ gqarr
+ }
[,1] [,2] [,3]
[1,] 0.9999987 2.525134e-06 1.262567e-06
[2,] 0.9999980 4.047515e-06 2.023757e-06
[3,] 0.9999994 1.234505e-06 6.172525e-07
[4,] 1.0000000 4.468203e-08 2.234102e-08
[5,] 1.0000000 1.624734e-26 8.123671e-27
[6,] 1.0000000 4.974526e-28 2.487263e-28
[7,] 1.0000000 2.211599e-43 1.105800e-43
[8,] 1.0000000 2.054460e-43 1.027230e-43
[9,] 1.0000000 6.347633e-45 3.173817e-45
[10,] 1.0000000 4.439635e-48 2.219818e-48
[11,] 1.0000000 2.534908e-47 1.267454e-47
[12,] 1.0000000 1.212679e-52 6.063396e-53
[13,] 1.0000000 4.139878e-54 2.069939e-54
[14,] 1.0000000 1.093743e-54 5.468717e-55
[15,] 1.0000000 3.229075e-54 1.614537e-54
[16,] 1.0000000 9.551320e-71 4.775660e-71
[17,] 1.0000000 1.068486e-71 5.342432e-72
[18,] 1.0000000 6.487460e-73 3.243730e-73
[19,] 1.0000000 3.894814e-85 1.947407e-85
[20,] 1.0000000 2.358233e-98 1.179116e-98
[21,] 1.0000000 2.833226e-100 1.416613e-100
[22,] 1.0000000 1.875534e-116 9.377668e-117
[23,] 1.0000000 1.394238e-117 6.971188e-118
[24,] 1.0000000 1.467401e-124 7.337006e-125
[25,] 1.0000000 5.326689e-126 2.663345e-126
[26,] 1.0000000 8.010114e-140 4.005057e-140
[27,] 1.0000000 3.721455e-139 1.860727e-139
[28,] 1.0000000 7.983380e-144 3.991690e-144
[29,] 1.0000000 1.862002e-149 9.310012e-150
[30,] 1.0000000 9.199134e-149 4.599567e-149
[31,] 1.0000000 8.508442e-151 4.254221e-151
[32,] 1.0000000 6.722279e-150 3.361140e-150
[33,] 1.0000000 2.405855e-149 1.202927e-149
[34,] 1.0000000 1.062030e-149 5.310149e-150
[35,] 1.0000000 6.707223e-149 3.353611e-149
[36,] 1.0000000 1.275244e-148 6.376221e-149
[37,] 1.0000000 1.261501e-147 6.307504e-148
[38,] 1.0000000 6.520955e-148 3.260477e-148
[39,] 1.0000000 7.271175e-147 3.635588e-147
[40,] 1.0000000 2.985270e-146 1.492635e-146
[41,] 1.0000000 5.040632e-149 2.520316e-149
[42,] 1.0000000 4.652739e-149 2.326369e-149
[43,] 1.0000000 1.764392e-148 8.821959e-149
[44,] 1.0000000 1.696967e-147 8.484836e-148
[45,] 1.0000000 4.615535e-147 2.307768e-147
[46,] 1.0000000 2.609989e-147 1.304994e-147
[47,] 1.0000000 1.328123e-146 6.640614e-147
[48,] 1.0000000 4.026839e-162 2.013420e-162
[49,] 1.0000000 7.677348e-162 3.838674e-162
[50,] 1.0000000 4.532324e-161 2.266162e-161
[51,] 1.0000000 3.644451e-160 1.822226e-160
[52,] 1.0000000 7.495641e-161 3.747821e-161
[53,] 1.0000000 1.057744e-163 5.288720e-164
[54,] 1.0000000 8.540673e-163 4.270337e-163
[55,] 1.0000000 8.701117e-162 4.350558e-162
[56,] 1.0000000 7.614139e-161 3.807070e-161
[57,] 1.0000000 1.531478e-160 7.657388e-161
[58,] 1.0000000 2.443651e-160 1.221826e-160
[59,] 1.0000000 3.007943e-159 1.503971e-159
[60,] 1.0000000 9.726816e-159 4.863408e-159
[61,] 1.0000000 1.004477e-157 5.022384e-158
[62,] 1.0000000 8.747778e-158 4.373889e-158
[63,] 1.0000000 1.808770e-157 9.043852e-158
[64,] 1.0000000 1.362661e-156 6.813305e-157
[65,] 1.0000000 3.746790e-156 1.873395e-156
[66,] 1.0000000 3.874700e-155 1.937350e-155
[67,] 1.0000000 8.795495e-155 4.397748e-155
[68,] 1.0000000 5.219534e-154 2.609767e-154
[69,] 1.0000000 1.072739e-153 5.363697e-154
[70,] 1.0000000 8.587405e-153 4.293702e-153
[71,] 1.0000000 2.218810e-152 1.109405e-152
[72,] 1.0000000 2.530985e-151 1.265493e-151
[73,] 1.0000000 1.485332e-150 7.426659e-151
[74,] 1.0000000 1.118108e-150 5.590542e-151
[75,] 1.0000000 4.070565e-150 2.035282e-150
[76,] 1.0000000 2.878832e-150 1.439416e-150
[77,] 1.0000000 2.438560e-149 1.219280e-149
[78,] 1.0000000 1.033172e-148 5.165862e-149
[79,] 1.0000000 2.114796e-148 1.057398e-148
[80,] 1.0000000 1.293381e-147 6.466905e-148
[81,] 1.0000000 2.279747e-147 1.139873e-147
[82,] 1.0000000 2.393802e-146 1.196901e-146
[83,] 1.0000000 2.472528e-145 1.236264e-145
[84,] 1.0000000 2.520119e-145 1.260059e-145
[85,] 1.0000000 1.659042e-144 8.295212e-145
[86,] 1.0000000 7.326679e-144 3.663339e-144
[87,] 1.0000000 4.011709e-143 2.005854e-143
[88,] 1.0000000 4.826902e-144 2.413451e-144
[89,] 1.0000000 3.485461e-143 1.742731e-143
[90,] 1.0000000 1.822082e-142 9.110411e-143
[91,] 1.0000000 1.540531e-141 7.702656e-142
[92,] 1.0000000 1.009244e-140 5.046220e-141
[93,] 1.0000000 1.322643e-140 6.613214e-141
[94,] 1.0000000 7.566089e-140 3.783044e-140
[95,] 1.0000000 1.790938e-139 8.954690e-140
[96,] 1.0000000 1.676282e-138 8.381411e-139
[97,] 1.0000000 1.105802e-137 5.529011e-138
[98,] 1.0000000 1.597249e-138 7.986244e-139
[99,] 1.0000000 9.764353e-139 4.882176e-139
[100,] 1.0000000 5.009629e-138 2.504814e-138
[101,] 1.0000000 2.667949e-137 1.333974e-137
[102,] 1.0000000 1.472540e-136 7.362702e-137
[103,] 1.0000000 3.348855e-136 1.674427e-136
[104,] 1.0000000 3.137742e-135 1.568871e-135
[105,] 1.0000000 2.795249e-134 1.397624e-134
[106,] 1.0000000 1.622690e-133 8.113450e-134
[107,] 1.0000000 1.460384e-132 7.301922e-133
[108,] 1.0000000 6.984676e-132 3.492338e-132
[109,] 1.0000000 3.277830e-131 1.638915e-131
[110,] 1.0000000 5.770063e-131 2.885032e-131
[111,] 1.0000000 4.842780e-130 2.421390e-130
[112,] 1.0000000 4.248938e-129 2.124469e-129
[113,] 1.0000000 4.171305e-128 2.085653e-128
[114,] 1.0000000 1.814105e-127 9.070525e-128
[115,] 1.0000000 8.950501e-127 4.475251e-127
[116,] 1.0000000 6.997107e-126 3.498553e-126
[117,] 1.0000000 2.107255e-125 1.053628e-125
[118,] 1.0000000 1.615920e-124 8.079598e-125
[119,] 1.0000000 1.474160e-123 7.370802e-124
[120,] 1.0000000 1.230256e-122 6.151282e-123
[121,] 1.0000000 6.219312e-122 3.109656e-122
[122,] 1.0000000 2.037865e-121 1.018932e-121
[123,] 1.0000000 1.613755e-120 8.068773e-121
[124,] 1.0000000 1.352853e-119 6.764266e-120
[125,] 1.0000000 1.207318e-118 6.036591e-119
[126,] 1.0000000 1.117061e-117 5.585305e-118
[127,] 1.0000000 8.549778e-117 4.274889e-117
[128,] 1.0000000 7.996446e-116 3.998223e-116
[129,] 1.0000000 4.670720e-115 2.335360e-115
[130,] 1.0000000 2.103810e-114 1.051905e-114
[131,] 1.0000000 1.909116e-113 9.545579e-114
[132,] 1.0000000 1.223256e-112 6.116280e-113
[133,] 1.0000000 6.849419e-112 3.424710e-112
[134,] 1.0000000 6.198288e-112 3.099144e-112
[135,] 1.0000000 2.293453e-112 1.146727e-112
[136,] 1.0000000 1.525454e-111 7.627268e-112
[137,] 1.0000000 9.807707e-111 4.903854e-111
[138,] 1.0000000 6.364236e-110 3.182118e-110
[139,] 1.0000000 5.771455e-109 2.885727e-109
[140,] 1.0000000 5.013340e-108 2.506670e-108
[141,] 1.0000000 4.714155e-107 2.357078e-107
[142,] 1.0000000 2.825978e-106 1.412989e-106
[143,] 1.0000000 1.978219e-105 9.891097e-106
[144,] 1.0000000 1.713485e-104 8.567423e-105
[145,] 1.0000000 4.855005e-104 2.427503e-104
[146,] 1.0000000 3.442397e-103 1.721199e-103
[147,] 1.0000000 2.778161e-102 1.389081e-102
[148,] 1.0000000 1.996460e-101 9.982298e-102
[149,] 1.0000000 1.721276e-100 8.606379e-101
[150,] 1.0000000 1.469188e-99 7.345941e-100
[151,] 1.0000000 3.525508e-99 1.762754e-99
[152,] 1.0000000 3.136015e-98 1.568007e-98
[153,] 1.0000000 2.335471e-97 1.167735e-97
[154,] 1.0000000 9.610108e-97 4.805054e-97
[155,] 1.0000000 7.355859e-96 3.677929e-96
[156,] 1.0000000 5.794852e-95 2.897426e-95
[157,] 1.0000000 5.033115e-94 2.516557e-94
[158,] 1.0000000 4.027470e-93 2.013735e-93
[159,] 1.0000000 3.033541e-92 1.516771e-92
[160,] 1.0000000 2.394133e-91 1.197066e-91
[161,] 1.0000000 1.960397e-90 9.801986e-91
[162,] 1.0000000 1.479679e-89 7.398393e-90
[163,] 1.0000000 1.020288e-88 5.101439e-89
[164,] 1.0000000 7.581519e-88 3.790760e-88
[165,] 1.0000000 4.304522e-87 2.152261e-87
[166,] 1.0000000 1.872037e-86 9.360185e-87
[167,] 1.0000000 6.152988e-86 3.076494e-86
[168,] 1.0000000 4.040559e-85 2.020280e-85
[169,] 1.0000000 2.721877e-84 1.360939e-84
[170,] 1.0000000 1.652293e-83 8.261466e-84
[171,] 1.0000000 9.578775e-83 4.789388e-83
[172,] 1.0000000 7.638108e-82 3.819054e-82
[173,] 1.0000000 5.101305e-81 2.550653e-81
[174,] 1.0000000 3.446527e-80 1.723263e-80
[175,] 1.0000000 2.632402e-79 1.316201e-79
[176,] 1.0000000 2.023325e-78 1.011662e-78
[177,] 1.0000000 1.569457e-77 7.847283e-78
[178,] 1.0000000 1.091540e-76 5.457700e-77
[179,] 1.0000000 8.372230e-76 4.186115e-76
[180,] 1.0000000 4.932789e-75 2.466395e-75
[181,] 1.0000000 3.175961e-74 1.587981e-74
[182,] 1.0000000 2.024464e-73 1.012232e-73
[183,] 1.0000000 1.256402e-72 6.282008e-73
[184,] 1.0000000 9.113828e-72 4.556914e-72
[185,] 1.0000000 6.068084e-71 3.034042e-71
[186,] 1.0000000 4.512743e-70 2.256371e-70
[187,] 1.0000000 2.686034e-69 1.343017e-69
[188,] 1.0000000 1.161197e-68 5.805984e-69
[189,] 1.0000000 8.239817e-68 4.119909e-68
[190,] 1.0000000 5.824635e-67 2.912318e-67
[191,] 1.0000000 4.130838e-66 2.065419e-66
[192,] 1.0000000 2.486618e-65 1.243309e-65
[193,] 1.0000000 1.494189e-64 7.470943e-65
[194,] 1.0000000 7.053364e-64 3.526682e-64
[195,] 1.0000000 5.071761e-63 2.535880e-63
[196,] 1.0000000 3.476921e-62 1.738461e-62
[197,] 1.0000000 1.728765e-61 8.643827e-62
[198,] 1.0000000 8.237878e-61 4.118939e-61
[199,] 1.0000000 2.773566e-60 1.386783e-60
[200,] 1.0000000 1.713319e-59 8.566597e-60
[201,] 1.0000000 1.193645e-58 5.968223e-59
[202,] 1.0000000 8.157256e-58 4.078628e-58
[203,] 1.0000000 4.408161e-57 2.204080e-57
[204,] 1.0000000 2.466463e-56 1.233231e-56
[205,] 1.0000000 1.604788e-55 8.023942e-56
[206,] 1.0000000 1.014029e-54 5.070144e-55
[207,] 1.0000000 6.070931e-54 3.035466e-54
[208,] 1.0000000 3.647412e-53 1.823706e-53
[209,] 1.0000000 1.661000e-52 8.305000e-53
[210,] 1.0000000 9.256780e-52 4.628390e-52
[211,] 1.0000000 5.809445e-51 2.904723e-51
[212,] 1.0000000 3.302658e-50 1.651329e-50
[213,] 1.0000000 1.378480e-49 6.892400e-50
[214,] 1.0000000 8.456451e-49 4.228225e-49
[215,] 1.0000000 4.688401e-48 2.344201e-48
[216,] 1.0000000 2.313660e-47 1.156830e-47
[217,] 1.0000000 1.464439e-46 7.322193e-47
[218,] 1.0000000 8.421116e-46 4.210558e-46
[219,] 1.0000000 5.016848e-45 2.508424e-45
[220,] 1.0000000 3.096248e-44 1.548124e-44
[221,] 1.0000000 1.742650e-43 8.713248e-44
[222,] 1.0000000 9.791182e-43 4.895591e-43
[223,] 1.0000000 3.824263e-42 1.912131e-42
[224,] 1.0000000 1.308045e-41 6.540227e-42
[225,] 1.0000000 6.854726e-41 3.427363e-41
[226,] 1.0000000 3.831456e-40 1.915728e-40
[227,] 1.0000000 2.081104e-39 1.040552e-39
[228,] 1.0000000 1.038433e-38 5.192163e-39
[229,] 1.0000000 4.964494e-38 2.482247e-38
[230,] 1.0000000 2.339444e-37 1.169722e-37
[231,] 1.0000000 1.352187e-36 6.760937e-37
[232,] 1.0000000 6.990720e-36 3.495360e-36
[233,] 1.0000000 3.771872e-35 1.885936e-35
[234,] 1.0000000 2.002004e-34 1.001002e-34
[235,] 1.0000000 8.357375e-34 4.178687e-34
[236,] 1.0000000 3.724562e-33 1.862281e-33
[237,] 1.0000000 1.982140e-32 9.910701e-33
[238,] 1.0000000 9.103563e-32 4.551781e-32
[239,] 1.0000000 3.639339e-31 1.819670e-31
[240,] 1.0000000 1.638972e-30 8.194860e-31
[241,] 1.0000000 8.745947e-30 4.372974e-30
[242,] 1.0000000 4.546064e-29 2.273032e-29
[243,] 1.0000000 2.065137e-28 1.032569e-28
[244,] 1.0000000 1.046477e-27 5.232383e-28
[245,] 1.0000000 5.419533e-27 2.709766e-27
[246,] 1.0000000 1.877749e-26 9.388743e-27
[247,] 1.0000000 9.353558e-26 4.676779e-26
[248,] 1.0000000 3.512926e-25 1.756463e-25
[249,] 1.0000000 1.697589e-24 8.487944e-25
[250,] 1.0000000 8.225187e-24 4.112593e-24
[251,] 1.0000000 3.898850e-23 1.949425e-23
[252,] 1.0000000 1.209013e-22 6.045066e-23
[253,] 1.0000000 1.641860e-22 8.209300e-23
[254,] 1.0000000 4.007077e-22 2.003538e-22
[255,] 1.0000000 1.734643e-21 8.673215e-22
[256,] 1.0000000 7.966901e-21 3.983451e-21
[257,] 1.0000000 3.424878e-20 1.712439e-20
[258,] 1.0000000 5.871199e-20 2.935600e-20
[259,] 1.0000000 2.551862e-19 1.275931e-19
[260,] 1.0000000 1.214332e-18 6.071661e-19
[261,] 1.0000000 4.907813e-18 2.453906e-18
[262,] 1.0000000 2.134379e-17 1.067190e-17
[263,] 1.0000000 7.015405e-17 3.507703e-17
[264,] 1.0000000 2.716006e-16 1.358003e-16
[265,] 1.0000000 7.075495e-16 3.537747e-16
[266,] 1.0000000 3.170924e-15 1.585462e-15
[267,] 1.0000000 1.401752e-14 7.008761e-15
[268,] 1.0000000 2.666502e-14 1.333251e-14
[269,] 1.0000000 1.170024e-13 5.850121e-14
[270,] 1.0000000 4.729221e-13 2.364610e-13
[271,] 1.0000000 1.251287e-12 6.256436e-13
[272,] 1.0000000 5.418683e-13 2.709341e-13
[273,] 1.0000000 2.485644e-12 1.242822e-12
[274,] 1.0000000 3.650719e-12 1.825359e-12
[275,] 1.0000000 1.666906e-11 8.334528e-12
[276,] 1.0000000 2.225956e-11 1.112978e-11
[277,] 1.0000000 1.042511e-10 5.212554e-11
[278,] 1.0000000 4.001823e-10 2.000912e-10
[279,] 1.0000000 1.822071e-09 9.110355e-10
[280,] 1.0000000 2.804989e-09 1.402495e-09
[281,] 1.0000000 1.323714e-09 6.618569e-10
[282,] 1.0000000 6.717223e-09 3.358611e-09
[283,] 1.0000000 7.521235e-09 3.760618e-09
[284,] 1.0000000 3.724654e-08 1.862327e-08
[285,] 0.9999999 1.360604e-07 6.803022e-08
[286,] 1.0000000 6.579974e-08 3.289987e-08
[287,] 0.9999998 3.836734e-07 1.918367e-07
[288,] 0.9999995 9.775037e-07 4.887518e-07
[289,] 0.9999971 5.846424e-06 2.923212e-06
[290,] 0.9999851 2.985313e-05 1.492657e-05
[291,] 0.9999203 1.594143e-04 7.970714e-05
[292,] 0.9996357 7.286242e-04 3.643121e-04
[293,] 0.9985047 2.990554e-03 1.495277e-03
[294,] 0.9977156 4.568877e-03 2.284438e-03
[295,] 0.9993467 1.306557e-03 6.532784e-04
> postscript(file="/var/www/html/rcomp/tmp/1bbai1291232200.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> plot(x[,1], type='l', main='Actuals and Interpolation', ylab='value of Actuals and Interpolation (dots)', xlab='time or index')
> points(x[,1]-mysum$resid)
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/2bbai1291232200.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> plot(mysum$resid, type='b', pch=19, main='Residuals', ylab='value of Residuals', xlab='time or index')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/3mkal1291232200.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> hist(mysum$resid, main='Residual Histogram', xlab='values of Residuals')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/4mkal1291232200.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> densityplot(~mysum$resid,col='black',main='Residual Density Plot', xlab='values of Residuals')
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/5mkal1291232200.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 = 314
Frequency = 1
1 2 3 4 5
1362967.1736 2780931.1541 1317495.4894 -2084201.5018 -445713.3330
6 7 8 9 10
73375.2662 -181641.9261 118317.2713 -400659.6024 615160.3371
11 12 13 14 15
-147065.6848 256126.7917 300419.3313 -992377.4418 -284064.4725
16 17 18 19 20
425089.7262 -213230.7692 -298381.9459 -375587.3613 -631095.5993
21 22 23 24 25
146575.4045 -185929.8003 -144652.0777 -667955.7109 -167653.6538
26 27 28 29 30
-137159.5009 -534929.5535 -195277.5256 246581.8931 -311679.0206
31 32 33 34 35
201870.8099 -124892.6240 182402.4323 79700.1927 -59423.9376
36 37 38 39 40
-598182.5702 -14898.6672 106459.0671 -64003.9842 -157210.3445
41 42 43 44 45
-334464.9916 -324247.1903 -39455.2658 -222597.4140 -58031.0557
46 47 48 49 50
-76332.5003 -243886.2669 -287196.3969 -79173.5418 -18280.0960
51 52 53 54 55
-6585.8409 -122531.9144 -231818.6310 23874.2287 -218084.7311
56 57 58 59 60
-28892.8951 131753.3868 -58190.4800 -56790.4135 -231033.4100
61 62 63 64 65
-28907.1432 -15528.2888 -26858.7533 -189832.5784 -37078.0150
66 67 68 69 70
-528319.2771 88995.1516 -40530.5083 -106062.4759 -106911.7828
71 72 73 74 75
20023.1890 -40068.0969 29806.4332 -2114.3821 -55930.2579
76 77 78 79 80
52003.2369 -34640.2824 83825.5884 62657.2124 -18699.6284
81 82 83 84 85
2156.0473 -83372.4103 -70123.4819 144347.9482 -12073.3311
86 87 88 89 90
-123408.2490 -75403.0055 -21812.5294 10975.8713 -22779.3065
91 92 93 94 95
-26893.5165 -49481.5892 19588.7263 -163558.3325 14364.8869
96 97 98 99 100
-11763.4129 -19704.5514 7394.5410 -66224.2995 -4994.4797
101 102 103 104 105
77735.4451 48326.0561 -26277.7460 57943.2431 4081.3474
106 107 108 109 110
86456.0662 -90424.7736 63560.9019 -152800.5120 59041.3533
111 112 113 114 115
47968.1096 -32747.7456 -22746.1688 72301.3139 9813.1066
116 117 118 119 120
2657.6602 15653.1987 18280.8119 31920.6663 37182.9258
121 122 123 124 125
-298720.4283 9017.0950 -121840.2138 38069.6663 -4045.3824
126 127 128 129 130
-955.7540 97133.4249 34853.8901 108837.8583 30799.6250
131 132 133 134 135
31362.2250 83529.8754 105944.9960 37654.3602 -5657.6658
136 137 138 139 140
66006.0499 34271.6014 122236.1478 -31976.9224 21866.1329
141 142 143 144 145
121694.4761 92736.6837 -111375.4345 -248386.0315 130437.9464
146 147 148 149 150
74148.3704 130974.6342 40615.9476 -17114.3772 7813.6856
151 152 153 154 155
140288.8159 123778.7897 45033.2574 54005.7392 -85032.5698
156 157 158 159 160
79098.9801 115031.2701 42902.6339 16371.0585 -251985.7870
161 162 163 164 165
19402.5665 2171.2951 76705.5087 126664.9797 83514.9431
166 167 168 169 170
23196.3749 -9279.2649 85063.5039 104005.9545 -19182.2198
171 172 173 174 175
84592.6901 -83762.6200 -157828.8037 93454.5965 153982.9185
176 177 178 179 180
-1063.9769 -88654.7749 -14717.3841 27726.1118 150907.0199
181 182 183 184 185
50410.5136 136340.8610 107702.1801 37168.5810 52125.9209
186 187 188 189 190
-3204.6641 94774.2108 -7863.8557 116718.5415 -79191.7184
191 192 193 194 195
109593.7106 122227.8025 73201.0440 120165.1651 16561.5414
196 197 198 199 200
128971.2808 -233667.5167 21662.5213 48826.6673 38481.8023
201 202 203 204 205
-78403.2373 117537.1159 136088.8357 6946.1219 93266.7355
206 207 208 209 210
-668.5685 118420.8990 139610.0784 12404.8292 35055.2748
211 212 213 214 215
18783.7467 -30962.5639 114267.8417 43675.8982 28218.3265
216 217 218 219 220
99945.4123 132700.5090 125365.3889 18329.0421 2471.4575
221 222 223 224 225
-10213.5335 154648.0856 41763.4498 15150.3021 138207.4323
226 227 228 229 230
44422.5117 76679.5657 21093.2038 47729.4636 49632.7602
231 232 233 234 235
25610.6683 130428.6334 324.4753 12274.3581 412.3646
236 237 238 239 240
38221.7388 -174553.0343 123848.6213 104615.0985 1868.7387
241 242 243 244 245
89744.2163 -55925.2029 73957.3433 128466.9310 -167325.2007
246 247 248 249 250
42024.1763 132326.6958 133507.6425 117088.6655 11943.5225
251 252 253 254 255
47244.5077 118025.5208 22294.6438 -10291.5560 18965.2072
256 257 258 259 260
31857.6418 83370.6252 38998.8280 59227.3013 50603.8573
261 262 263 264 265
2776.1763 -151092.5860 152638.1761 23380.0287 -5569.0632
266 267 268 269 270
28679.7829 156884.4521 42333.3868 34539.9273 9483.0331
271 272 273 274 275
104009.3825 45912.0144 130341.1689 19409.9300 124519.3013
276 277 278 279 280
104079.9597 4616.6885 38288.1777 131453.3728 8544.5835
281 282 283 284 285
-245065.2680 119708.2855 79892.2718 73565.3265 3841.6443
286 287 288 289 290
119692.0345 88489.1070 106970.1481 -2234.1475 26377.8586
291 292 293 294 295
109111.2094 128575.7542 68524.4566 55507.1617 -22791.3372
296 297 298 299 300
107569.0796 8708.4356 -586.0401 63994.9277 45714.7608
301 302 303 304 305
107931.1100 88592.0011 74084.2139 35272.3096 -29392.3792
306 307 308 309 310
-402263.5485 -265592.9731 86319.3857 -99418.7196 -135733.4010
311 312 313 314
-991661.4399 -114661.6884 -353186.5569 -524362.4068
> postscript(file="/var/www/html/rcomp/tmp/6xcr61291232200.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 = 314
Frequency = 1
lag(myerror, k = 1) myerror
0 1362967.1736 NA
1 2780931.1541 1362967.1736
2 1317495.4894 2780931.1541
3 -2084201.5018 1317495.4894
4 -445713.3330 -2084201.5018
5 73375.2662 -445713.3330
6 -181641.9261 73375.2662
7 118317.2713 -181641.9261
8 -400659.6024 118317.2713
9 615160.3371 -400659.6024
10 -147065.6848 615160.3371
11 256126.7917 -147065.6848
12 300419.3313 256126.7917
13 -992377.4418 300419.3313
14 -284064.4725 -992377.4418
15 425089.7262 -284064.4725
16 -213230.7692 425089.7262
17 -298381.9459 -213230.7692
18 -375587.3613 -298381.9459
19 -631095.5993 -375587.3613
20 146575.4045 -631095.5993
21 -185929.8003 146575.4045
22 -144652.0777 -185929.8003
23 -667955.7109 -144652.0777
24 -167653.6538 -667955.7109
25 -137159.5009 -167653.6538
26 -534929.5535 -137159.5009
27 -195277.5256 -534929.5535
28 246581.8931 -195277.5256
29 -311679.0206 246581.8931
30 201870.8099 -311679.0206
31 -124892.6240 201870.8099
32 182402.4323 -124892.6240
33 79700.1927 182402.4323
34 -59423.9376 79700.1927
35 -598182.5702 -59423.9376
36 -14898.6672 -598182.5702
37 106459.0671 -14898.6672
38 -64003.9842 106459.0671
39 -157210.3445 -64003.9842
40 -334464.9916 -157210.3445
41 -324247.1903 -334464.9916
42 -39455.2658 -324247.1903
43 -222597.4140 -39455.2658
44 -58031.0557 -222597.4140
45 -76332.5003 -58031.0557
46 -243886.2669 -76332.5003
47 -287196.3969 -243886.2669
48 -79173.5418 -287196.3969
49 -18280.0960 -79173.5418
50 -6585.8409 -18280.0960
51 -122531.9144 -6585.8409
52 -231818.6310 -122531.9144
53 23874.2287 -231818.6310
54 -218084.7311 23874.2287
55 -28892.8951 -218084.7311
56 131753.3868 -28892.8951
57 -58190.4800 131753.3868
58 -56790.4135 -58190.4800
59 -231033.4100 -56790.4135
60 -28907.1432 -231033.4100
61 -15528.2888 -28907.1432
62 -26858.7533 -15528.2888
63 -189832.5784 -26858.7533
64 -37078.0150 -189832.5784
65 -528319.2771 -37078.0150
66 88995.1516 -528319.2771
67 -40530.5083 88995.1516
68 -106062.4759 -40530.5083
69 -106911.7828 -106062.4759
70 20023.1890 -106911.7828
71 -40068.0969 20023.1890
72 29806.4332 -40068.0969
73 -2114.3821 29806.4332
74 -55930.2579 -2114.3821
75 52003.2369 -55930.2579
76 -34640.2824 52003.2369
77 83825.5884 -34640.2824
78 62657.2124 83825.5884
79 -18699.6284 62657.2124
80 2156.0473 -18699.6284
81 -83372.4103 2156.0473
82 -70123.4819 -83372.4103
83 144347.9482 -70123.4819
84 -12073.3311 144347.9482
85 -123408.2490 -12073.3311
86 -75403.0055 -123408.2490
87 -21812.5294 -75403.0055
88 10975.8713 -21812.5294
89 -22779.3065 10975.8713
90 -26893.5165 -22779.3065
91 -49481.5892 -26893.5165
92 19588.7263 -49481.5892
93 -163558.3325 19588.7263
94 14364.8869 -163558.3325
95 -11763.4129 14364.8869
96 -19704.5514 -11763.4129
97 7394.5410 -19704.5514
98 -66224.2995 7394.5410
99 -4994.4797 -66224.2995
100 77735.4451 -4994.4797
101 48326.0561 77735.4451
102 -26277.7460 48326.0561
103 57943.2431 -26277.7460
104 4081.3474 57943.2431
105 86456.0662 4081.3474
106 -90424.7736 86456.0662
107 63560.9019 -90424.7736
108 -152800.5120 63560.9019
109 59041.3533 -152800.5120
110 47968.1096 59041.3533
111 -32747.7456 47968.1096
112 -22746.1688 -32747.7456
113 72301.3139 -22746.1688
114 9813.1066 72301.3139
115 2657.6602 9813.1066
116 15653.1987 2657.6602
117 18280.8119 15653.1987
118 31920.6663 18280.8119
119 37182.9258 31920.6663
120 -298720.4283 37182.9258
121 9017.0950 -298720.4283
122 -121840.2138 9017.0950
123 38069.6663 -121840.2138
124 -4045.3824 38069.6663
125 -955.7540 -4045.3824
126 97133.4249 -955.7540
127 34853.8901 97133.4249
128 108837.8583 34853.8901
129 30799.6250 108837.8583
130 31362.2250 30799.6250
131 83529.8754 31362.2250
132 105944.9960 83529.8754
133 37654.3602 105944.9960
134 -5657.6658 37654.3602
135 66006.0499 -5657.6658
136 34271.6014 66006.0499
137 122236.1478 34271.6014
138 -31976.9224 122236.1478
139 21866.1329 -31976.9224
140 121694.4761 21866.1329
141 92736.6837 121694.4761
142 -111375.4345 92736.6837
143 -248386.0315 -111375.4345
144 130437.9464 -248386.0315
145 74148.3704 130437.9464
146 130974.6342 74148.3704
147 40615.9476 130974.6342
148 -17114.3772 40615.9476
149 7813.6856 -17114.3772
150 140288.8159 7813.6856
151 123778.7897 140288.8159
152 45033.2574 123778.7897
153 54005.7392 45033.2574
154 -85032.5698 54005.7392
155 79098.9801 -85032.5698
156 115031.2701 79098.9801
157 42902.6339 115031.2701
158 16371.0585 42902.6339
159 -251985.7870 16371.0585
160 19402.5665 -251985.7870
161 2171.2951 19402.5665
162 76705.5087 2171.2951
163 126664.9797 76705.5087
164 83514.9431 126664.9797
165 23196.3749 83514.9431
166 -9279.2649 23196.3749
167 85063.5039 -9279.2649
168 104005.9545 85063.5039
169 -19182.2198 104005.9545
170 84592.6901 -19182.2198
171 -83762.6200 84592.6901
172 -157828.8037 -83762.6200
173 93454.5965 -157828.8037
174 153982.9185 93454.5965
175 -1063.9769 153982.9185
176 -88654.7749 -1063.9769
177 -14717.3841 -88654.7749
178 27726.1118 -14717.3841
179 150907.0199 27726.1118
180 50410.5136 150907.0199
181 136340.8610 50410.5136
182 107702.1801 136340.8610
183 37168.5810 107702.1801
184 52125.9209 37168.5810
185 -3204.6641 52125.9209
186 94774.2108 -3204.6641
187 -7863.8557 94774.2108
188 116718.5415 -7863.8557
189 -79191.7184 116718.5415
190 109593.7106 -79191.7184
191 122227.8025 109593.7106
192 73201.0440 122227.8025
193 120165.1651 73201.0440
194 16561.5414 120165.1651
195 128971.2808 16561.5414
196 -233667.5167 128971.2808
197 21662.5213 -233667.5167
198 48826.6673 21662.5213
199 38481.8023 48826.6673
200 -78403.2373 38481.8023
201 117537.1159 -78403.2373
202 136088.8357 117537.1159
203 6946.1219 136088.8357
204 93266.7355 6946.1219
205 -668.5685 93266.7355
206 118420.8990 -668.5685
207 139610.0784 118420.8990
208 12404.8292 139610.0784
209 35055.2748 12404.8292
210 18783.7467 35055.2748
211 -30962.5639 18783.7467
212 114267.8417 -30962.5639
213 43675.8982 114267.8417
214 28218.3265 43675.8982
215 99945.4123 28218.3265
216 132700.5090 99945.4123
217 125365.3889 132700.5090
218 18329.0421 125365.3889
219 2471.4575 18329.0421
220 -10213.5335 2471.4575
221 154648.0856 -10213.5335
222 41763.4498 154648.0856
223 15150.3021 41763.4498
224 138207.4323 15150.3021
225 44422.5117 138207.4323
226 76679.5657 44422.5117
227 21093.2038 76679.5657
228 47729.4636 21093.2038
229 49632.7602 47729.4636
230 25610.6683 49632.7602
231 130428.6334 25610.6683
232 324.4753 130428.6334
233 12274.3581 324.4753
234 412.3646 12274.3581
235 38221.7388 412.3646
236 -174553.0343 38221.7388
237 123848.6213 -174553.0343
238 104615.0985 123848.6213
239 1868.7387 104615.0985
240 89744.2163 1868.7387
241 -55925.2029 89744.2163
242 73957.3433 -55925.2029
243 128466.9310 73957.3433
244 -167325.2007 128466.9310
245 42024.1763 -167325.2007
246 132326.6958 42024.1763
247 133507.6425 132326.6958
248 117088.6655 133507.6425
249 11943.5225 117088.6655
250 47244.5077 11943.5225
251 118025.5208 47244.5077
252 22294.6438 118025.5208
253 -10291.5560 22294.6438
254 18965.2072 -10291.5560
255 31857.6418 18965.2072
256 83370.6252 31857.6418
257 38998.8280 83370.6252
258 59227.3013 38998.8280
259 50603.8573 59227.3013
260 2776.1763 50603.8573
261 -151092.5860 2776.1763
262 152638.1761 -151092.5860
263 23380.0287 152638.1761
264 -5569.0632 23380.0287
265 28679.7829 -5569.0632
266 156884.4521 28679.7829
267 42333.3868 156884.4521
268 34539.9273 42333.3868
269 9483.0331 34539.9273
270 104009.3825 9483.0331
271 45912.0144 104009.3825
272 130341.1689 45912.0144
273 19409.9300 130341.1689
274 124519.3013 19409.9300
275 104079.9597 124519.3013
276 4616.6885 104079.9597
277 38288.1777 4616.6885
278 131453.3728 38288.1777
279 8544.5835 131453.3728
280 -245065.2680 8544.5835
281 119708.2855 -245065.2680
282 79892.2718 119708.2855
283 73565.3265 79892.2718
284 3841.6443 73565.3265
285 119692.0345 3841.6443
286 88489.1070 119692.0345
287 106970.1481 88489.1070
288 -2234.1475 106970.1481
289 26377.8586 -2234.1475
290 109111.2094 26377.8586
291 128575.7542 109111.2094
292 68524.4566 128575.7542
293 55507.1617 68524.4566
294 -22791.3372 55507.1617
295 107569.0796 -22791.3372
296 8708.4356 107569.0796
297 -586.0401 8708.4356
298 63994.9277 -586.0401
299 45714.7608 63994.9277
300 107931.1100 45714.7608
301 88592.0011 107931.1100
302 74084.2139 88592.0011
303 35272.3096 74084.2139
304 -29392.3792 35272.3096
305 -402263.5485 -29392.3792
306 -265592.9731 -402263.5485
307 86319.3857 -265592.9731
308 -99418.7196 86319.3857
309 -135733.4010 -99418.7196
310 -991661.4399 -135733.4010
311 -114661.6884 -991661.4399
312 -353186.5569 -114661.6884
313 -524362.4068 -353186.5569
314 NA -524362.4068
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] 2780931.1541 1362967.1736
[2,] 1317495.4894 2780931.1541
[3,] -2084201.5018 1317495.4894
[4,] -445713.3330 -2084201.5018
[5,] 73375.2662 -445713.3330
[6,] -181641.9261 73375.2662
[7,] 118317.2713 -181641.9261
[8,] -400659.6024 118317.2713
[9,] 615160.3371 -400659.6024
[10,] -147065.6848 615160.3371
[11,] 256126.7917 -147065.6848
[12,] 300419.3313 256126.7917
[13,] -992377.4418 300419.3313
[14,] -284064.4725 -992377.4418
[15,] 425089.7262 -284064.4725
[16,] -213230.7692 425089.7262
[17,] -298381.9459 -213230.7692
[18,] -375587.3613 -298381.9459
[19,] -631095.5993 -375587.3613
[20,] 146575.4045 -631095.5993
[21,] -185929.8003 146575.4045
[22,] -144652.0777 -185929.8003
[23,] -667955.7109 -144652.0777
[24,] -167653.6538 -667955.7109
[25,] -137159.5009 -167653.6538
[26,] -534929.5535 -137159.5009
[27,] -195277.5256 -534929.5535
[28,] 246581.8931 -195277.5256
[29,] -311679.0206 246581.8931
[30,] 201870.8099 -311679.0206
[31,] -124892.6240 201870.8099
[32,] 182402.4323 -124892.6240
[33,] 79700.1927 182402.4323
[34,] -59423.9376 79700.1927
[35,] -598182.5702 -59423.9376
[36,] -14898.6672 -598182.5702
[37,] 106459.0671 -14898.6672
[38,] -64003.9842 106459.0671
[39,] -157210.3445 -64003.9842
[40,] -334464.9916 -157210.3445
[41,] -324247.1903 -334464.9916
[42,] -39455.2658 -324247.1903
[43,] -222597.4140 -39455.2658
[44,] -58031.0557 -222597.4140
[45,] -76332.5003 -58031.0557
[46,] -243886.2669 -76332.5003
[47,] -287196.3969 -243886.2669
[48,] -79173.5418 -287196.3969
[49,] -18280.0960 -79173.5418
[50,] -6585.8409 -18280.0960
[51,] -122531.9144 -6585.8409
[52,] -231818.6310 -122531.9144
[53,] 23874.2287 -231818.6310
[54,] -218084.7311 23874.2287
[55,] -28892.8951 -218084.7311
[56,] 131753.3868 -28892.8951
[57,] -58190.4800 131753.3868
[58,] -56790.4135 -58190.4800
[59,] -231033.4100 -56790.4135
[60,] -28907.1432 -231033.4100
[61,] -15528.2888 -28907.1432
[62,] -26858.7533 -15528.2888
[63,] -189832.5784 -26858.7533
[64,] -37078.0150 -189832.5784
[65,] -528319.2771 -37078.0150
[66,] 88995.1516 -528319.2771
[67,] -40530.5083 88995.1516
[68,] -106062.4759 -40530.5083
[69,] -106911.7828 -106062.4759
[70,] 20023.1890 -106911.7828
[71,] -40068.0969 20023.1890
[72,] 29806.4332 -40068.0969
[73,] -2114.3821 29806.4332
[74,] -55930.2579 -2114.3821
[75,] 52003.2369 -55930.2579
[76,] -34640.2824 52003.2369
[77,] 83825.5884 -34640.2824
[78,] 62657.2124 83825.5884
[79,] -18699.6284 62657.2124
[80,] 2156.0473 -18699.6284
[81,] -83372.4103 2156.0473
[82,] -70123.4819 -83372.4103
[83,] 144347.9482 -70123.4819
[84,] -12073.3311 144347.9482
[85,] -123408.2490 -12073.3311
[86,] -75403.0055 -123408.2490
[87,] -21812.5294 -75403.0055
[88,] 10975.8713 -21812.5294
[89,] -22779.3065 10975.8713
[90,] -26893.5165 -22779.3065
[91,] -49481.5892 -26893.5165
[92,] 19588.7263 -49481.5892
[93,] -163558.3325 19588.7263
[94,] 14364.8869 -163558.3325
[95,] -11763.4129 14364.8869
[96,] -19704.5514 -11763.4129
[97,] 7394.5410 -19704.5514
[98,] -66224.2995 7394.5410
[99,] -4994.4797 -66224.2995
[100,] 77735.4451 -4994.4797
[101,] 48326.0561 77735.4451
[102,] -26277.7460 48326.0561
[103,] 57943.2431 -26277.7460
[104,] 4081.3474 57943.2431
[105,] 86456.0662 4081.3474
[106,] -90424.7736 86456.0662
[107,] 63560.9019 -90424.7736
[108,] -152800.5120 63560.9019
[109,] 59041.3533 -152800.5120
[110,] 47968.1096 59041.3533
[111,] -32747.7456 47968.1096
[112,] -22746.1688 -32747.7456
[113,] 72301.3139 -22746.1688
[114,] 9813.1066 72301.3139
[115,] 2657.6602 9813.1066
[116,] 15653.1987 2657.6602
[117,] 18280.8119 15653.1987
[118,] 31920.6663 18280.8119
[119,] 37182.9258 31920.6663
[120,] -298720.4283 37182.9258
[121,] 9017.0950 -298720.4283
[122,] -121840.2138 9017.0950
[123,] 38069.6663 -121840.2138
[124,] -4045.3824 38069.6663
[125,] -955.7540 -4045.3824
[126,] 97133.4249 -955.7540
[127,] 34853.8901 97133.4249
[128,] 108837.8583 34853.8901
[129,] 30799.6250 108837.8583
[130,] 31362.2250 30799.6250
[131,] 83529.8754 31362.2250
[132,] 105944.9960 83529.8754
[133,] 37654.3602 105944.9960
[134,] -5657.6658 37654.3602
[135,] 66006.0499 -5657.6658
[136,] 34271.6014 66006.0499
[137,] 122236.1478 34271.6014
[138,] -31976.9224 122236.1478
[139,] 21866.1329 -31976.9224
[140,] 121694.4761 21866.1329
[141,] 92736.6837 121694.4761
[142,] -111375.4345 92736.6837
[143,] -248386.0315 -111375.4345
[144,] 130437.9464 -248386.0315
[145,] 74148.3704 130437.9464
[146,] 130974.6342 74148.3704
[147,] 40615.9476 130974.6342
[148,] -17114.3772 40615.9476
[149,] 7813.6856 -17114.3772
[150,] 140288.8159 7813.6856
[151,] 123778.7897 140288.8159
[152,] 45033.2574 123778.7897
[153,] 54005.7392 45033.2574
[154,] -85032.5698 54005.7392
[155,] 79098.9801 -85032.5698
[156,] 115031.2701 79098.9801
[157,] 42902.6339 115031.2701
[158,] 16371.0585 42902.6339
[159,] -251985.7870 16371.0585
[160,] 19402.5665 -251985.7870
[161,] 2171.2951 19402.5665
[162,] 76705.5087 2171.2951
[163,] 126664.9797 76705.5087
[164,] 83514.9431 126664.9797
[165,] 23196.3749 83514.9431
[166,] -9279.2649 23196.3749
[167,] 85063.5039 -9279.2649
[168,] 104005.9545 85063.5039
[169,] -19182.2198 104005.9545
[170,] 84592.6901 -19182.2198
[171,] -83762.6200 84592.6901
[172,] -157828.8037 -83762.6200
[173,] 93454.5965 -157828.8037
[174,] 153982.9185 93454.5965
[175,] -1063.9769 153982.9185
[176,] -88654.7749 -1063.9769
[177,] -14717.3841 -88654.7749
[178,] 27726.1118 -14717.3841
[179,] 150907.0199 27726.1118
[180,] 50410.5136 150907.0199
[181,] 136340.8610 50410.5136
[182,] 107702.1801 136340.8610
[183,] 37168.5810 107702.1801
[184,] 52125.9209 37168.5810
[185,] -3204.6641 52125.9209
[186,] 94774.2108 -3204.6641
[187,] -7863.8557 94774.2108
[188,] 116718.5415 -7863.8557
[189,] -79191.7184 116718.5415
[190,] 109593.7106 -79191.7184
[191,] 122227.8025 109593.7106
[192,] 73201.0440 122227.8025
[193,] 120165.1651 73201.0440
[194,] 16561.5414 120165.1651
[195,] 128971.2808 16561.5414
[196,] -233667.5167 128971.2808
[197,] 21662.5213 -233667.5167
[198,] 48826.6673 21662.5213
[199,] 38481.8023 48826.6673
[200,] -78403.2373 38481.8023
[201,] 117537.1159 -78403.2373
[202,] 136088.8357 117537.1159
[203,] 6946.1219 136088.8357
[204,] 93266.7355 6946.1219
[205,] -668.5685 93266.7355
[206,] 118420.8990 -668.5685
[207,] 139610.0784 118420.8990
[208,] 12404.8292 139610.0784
[209,] 35055.2748 12404.8292
[210,] 18783.7467 35055.2748
[211,] -30962.5639 18783.7467
[212,] 114267.8417 -30962.5639
[213,] 43675.8982 114267.8417
[214,] 28218.3265 43675.8982
[215,] 99945.4123 28218.3265
[216,] 132700.5090 99945.4123
[217,] 125365.3889 132700.5090
[218,] 18329.0421 125365.3889
[219,] 2471.4575 18329.0421
[220,] -10213.5335 2471.4575
[221,] 154648.0856 -10213.5335
[222,] 41763.4498 154648.0856
[223,] 15150.3021 41763.4498
[224,] 138207.4323 15150.3021
[225,] 44422.5117 138207.4323
[226,] 76679.5657 44422.5117
[227,] 21093.2038 76679.5657
[228,] 47729.4636 21093.2038
[229,] 49632.7602 47729.4636
[230,] 25610.6683 49632.7602
[231,] 130428.6334 25610.6683
[232,] 324.4753 130428.6334
[233,] 12274.3581 324.4753
[234,] 412.3646 12274.3581
[235,] 38221.7388 412.3646
[236,] -174553.0343 38221.7388
[237,] 123848.6213 -174553.0343
[238,] 104615.0985 123848.6213
[239,] 1868.7387 104615.0985
[240,] 89744.2163 1868.7387
[241,] -55925.2029 89744.2163
[242,] 73957.3433 -55925.2029
[243,] 128466.9310 73957.3433
[244,] -167325.2007 128466.9310
[245,] 42024.1763 -167325.2007
[246,] 132326.6958 42024.1763
[247,] 133507.6425 132326.6958
[248,] 117088.6655 133507.6425
[249,] 11943.5225 117088.6655
[250,] 47244.5077 11943.5225
[251,] 118025.5208 47244.5077
[252,] 22294.6438 118025.5208
[253,] -10291.5560 22294.6438
[254,] 18965.2072 -10291.5560
[255,] 31857.6418 18965.2072
[256,] 83370.6252 31857.6418
[257,] 38998.8280 83370.6252
[258,] 59227.3013 38998.8280
[259,] 50603.8573 59227.3013
[260,] 2776.1763 50603.8573
[261,] -151092.5860 2776.1763
[262,] 152638.1761 -151092.5860
[263,] 23380.0287 152638.1761
[264,] -5569.0632 23380.0287
[265,] 28679.7829 -5569.0632
[266,] 156884.4521 28679.7829
[267,] 42333.3868 156884.4521
[268,] 34539.9273 42333.3868
[269,] 9483.0331 34539.9273
[270,] 104009.3825 9483.0331
[271,] 45912.0144 104009.3825
[272,] 130341.1689 45912.0144
[273,] 19409.9300 130341.1689
[274,] 124519.3013 19409.9300
[275,] 104079.9597 124519.3013
[276,] 4616.6885 104079.9597
[277,] 38288.1777 4616.6885
[278,] 131453.3728 38288.1777
[279,] 8544.5835 131453.3728
[280,] -245065.2680 8544.5835
[281,] 119708.2855 -245065.2680
[282,] 79892.2718 119708.2855
[283,] 73565.3265 79892.2718
[284,] 3841.6443 73565.3265
[285,] 119692.0345 3841.6443
[286,] 88489.1070 119692.0345
[287,] 106970.1481 88489.1070
[288,] -2234.1475 106970.1481
[289,] 26377.8586 -2234.1475
[290,] 109111.2094 26377.8586
[291,] 128575.7542 109111.2094
[292,] 68524.4566 128575.7542
[293,] 55507.1617 68524.4566
[294,] -22791.3372 55507.1617
[295,] 107569.0796 -22791.3372
[296,] 8708.4356 107569.0796
[297,] -586.0401 8708.4356
[298,] 63994.9277 -586.0401
[299,] 45714.7608 63994.9277
[300,] 107931.1100 45714.7608
[301,] 88592.0011 107931.1100
[302,] 74084.2139 88592.0011
[303,] 35272.3096 74084.2139
[304,] -29392.3792 35272.3096
[305,] -402263.5485 -29392.3792
[306,] -265592.9731 -402263.5485
[307,] 86319.3857 -265592.9731
[308,] -99418.7196 86319.3857
[309,] -135733.4010 -99418.7196
[310,] -991661.4399 -135733.4010
[311,] -114661.6884 -991661.4399
[312,] -353186.5569 -114661.6884
[313,] -524362.4068 -353186.5569
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 2780931.1541 1362967.1736
2 1317495.4894 2780931.1541
3 -2084201.5018 1317495.4894
4 -445713.3330 -2084201.5018
5 73375.2662 -445713.3330
6 -181641.9261 73375.2662
7 118317.2713 -181641.9261
8 -400659.6024 118317.2713
9 615160.3371 -400659.6024
10 -147065.6848 615160.3371
11 256126.7917 -147065.6848
12 300419.3313 256126.7917
13 -992377.4418 300419.3313
14 -284064.4725 -992377.4418
15 425089.7262 -284064.4725
16 -213230.7692 425089.7262
17 -298381.9459 -213230.7692
18 -375587.3613 -298381.9459
19 -631095.5993 -375587.3613
20 146575.4045 -631095.5993
21 -185929.8003 146575.4045
22 -144652.0777 -185929.8003
23 -667955.7109 -144652.0777
24 -167653.6538 -667955.7109
25 -137159.5009 -167653.6538
26 -534929.5535 -137159.5009
27 -195277.5256 -534929.5535
28 246581.8931 -195277.5256
29 -311679.0206 246581.8931
30 201870.8099 -311679.0206
31 -124892.6240 201870.8099
32 182402.4323 -124892.6240
33 79700.1927 182402.4323
34 -59423.9376 79700.1927
35 -598182.5702 -59423.9376
36 -14898.6672 -598182.5702
37 106459.0671 -14898.6672
38 -64003.9842 106459.0671
39 -157210.3445 -64003.9842
40 -334464.9916 -157210.3445
41 -324247.1903 -334464.9916
42 -39455.2658 -324247.1903
43 -222597.4140 -39455.2658
44 -58031.0557 -222597.4140
45 -76332.5003 -58031.0557
46 -243886.2669 -76332.5003
47 -287196.3969 -243886.2669
48 -79173.5418 -287196.3969
49 -18280.0960 -79173.5418
50 -6585.8409 -18280.0960
51 -122531.9144 -6585.8409
52 -231818.6310 -122531.9144
53 23874.2287 -231818.6310
54 -218084.7311 23874.2287
55 -28892.8951 -218084.7311
56 131753.3868 -28892.8951
57 -58190.4800 131753.3868
58 -56790.4135 -58190.4800
59 -231033.4100 -56790.4135
60 -28907.1432 -231033.4100
61 -15528.2888 -28907.1432
62 -26858.7533 -15528.2888
63 -189832.5784 -26858.7533
64 -37078.0150 -189832.5784
65 -528319.2771 -37078.0150
66 88995.1516 -528319.2771
67 -40530.5083 88995.1516
68 -106062.4759 -40530.5083
69 -106911.7828 -106062.4759
70 20023.1890 -106911.7828
71 -40068.0969 20023.1890
72 29806.4332 -40068.0969
73 -2114.3821 29806.4332
74 -55930.2579 -2114.3821
75 52003.2369 -55930.2579
76 -34640.2824 52003.2369
77 83825.5884 -34640.2824
78 62657.2124 83825.5884
79 -18699.6284 62657.2124
80 2156.0473 -18699.6284
81 -83372.4103 2156.0473
82 -70123.4819 -83372.4103
83 144347.9482 -70123.4819
84 -12073.3311 144347.9482
85 -123408.2490 -12073.3311
86 -75403.0055 -123408.2490
87 -21812.5294 -75403.0055
88 10975.8713 -21812.5294
89 -22779.3065 10975.8713
90 -26893.5165 -22779.3065
91 -49481.5892 -26893.5165
92 19588.7263 -49481.5892
93 -163558.3325 19588.7263
94 14364.8869 -163558.3325
95 -11763.4129 14364.8869
96 -19704.5514 -11763.4129
97 7394.5410 -19704.5514
98 -66224.2995 7394.5410
99 -4994.4797 -66224.2995
100 77735.4451 -4994.4797
101 48326.0561 77735.4451
102 -26277.7460 48326.0561
103 57943.2431 -26277.7460
104 4081.3474 57943.2431
105 86456.0662 4081.3474
106 -90424.7736 86456.0662
107 63560.9019 -90424.7736
108 -152800.5120 63560.9019
109 59041.3533 -152800.5120
110 47968.1096 59041.3533
111 -32747.7456 47968.1096
112 -22746.1688 -32747.7456
113 72301.3139 -22746.1688
114 9813.1066 72301.3139
115 2657.6602 9813.1066
116 15653.1987 2657.6602
117 18280.8119 15653.1987
118 31920.6663 18280.8119
119 37182.9258 31920.6663
120 -298720.4283 37182.9258
121 9017.0950 -298720.4283
122 -121840.2138 9017.0950
123 38069.6663 -121840.2138
124 -4045.3824 38069.6663
125 -955.7540 -4045.3824
126 97133.4249 -955.7540
127 34853.8901 97133.4249
128 108837.8583 34853.8901
129 30799.6250 108837.8583
130 31362.2250 30799.6250
131 83529.8754 31362.2250
132 105944.9960 83529.8754
133 37654.3602 105944.9960
134 -5657.6658 37654.3602
135 66006.0499 -5657.6658
136 34271.6014 66006.0499
137 122236.1478 34271.6014
138 -31976.9224 122236.1478
139 21866.1329 -31976.9224
140 121694.4761 21866.1329
141 92736.6837 121694.4761
142 -111375.4345 92736.6837
143 -248386.0315 -111375.4345
144 130437.9464 -248386.0315
145 74148.3704 130437.9464
146 130974.6342 74148.3704
147 40615.9476 130974.6342
148 -17114.3772 40615.9476
149 7813.6856 -17114.3772
150 140288.8159 7813.6856
151 123778.7897 140288.8159
152 45033.2574 123778.7897
153 54005.7392 45033.2574
154 -85032.5698 54005.7392
155 79098.9801 -85032.5698
156 115031.2701 79098.9801
157 42902.6339 115031.2701
158 16371.0585 42902.6339
159 -251985.7870 16371.0585
160 19402.5665 -251985.7870
161 2171.2951 19402.5665
162 76705.5087 2171.2951
163 126664.9797 76705.5087
164 83514.9431 126664.9797
165 23196.3749 83514.9431
166 -9279.2649 23196.3749
167 85063.5039 -9279.2649
168 104005.9545 85063.5039
169 -19182.2198 104005.9545
170 84592.6901 -19182.2198
171 -83762.6200 84592.6901
172 -157828.8037 -83762.6200
173 93454.5965 -157828.8037
174 153982.9185 93454.5965
175 -1063.9769 153982.9185
176 -88654.7749 -1063.9769
177 -14717.3841 -88654.7749
178 27726.1118 -14717.3841
179 150907.0199 27726.1118
180 50410.5136 150907.0199
181 136340.8610 50410.5136
182 107702.1801 136340.8610
183 37168.5810 107702.1801
184 52125.9209 37168.5810
185 -3204.6641 52125.9209
186 94774.2108 -3204.6641
187 -7863.8557 94774.2108
188 116718.5415 -7863.8557
189 -79191.7184 116718.5415
190 109593.7106 -79191.7184
191 122227.8025 109593.7106
192 73201.0440 122227.8025
193 120165.1651 73201.0440
194 16561.5414 120165.1651
195 128971.2808 16561.5414
196 -233667.5167 128971.2808
197 21662.5213 -233667.5167
198 48826.6673 21662.5213
199 38481.8023 48826.6673
200 -78403.2373 38481.8023
201 117537.1159 -78403.2373
202 136088.8357 117537.1159
203 6946.1219 136088.8357
204 93266.7355 6946.1219
205 -668.5685 93266.7355
206 118420.8990 -668.5685
207 139610.0784 118420.8990
208 12404.8292 139610.0784
209 35055.2748 12404.8292
210 18783.7467 35055.2748
211 -30962.5639 18783.7467
212 114267.8417 -30962.5639
213 43675.8982 114267.8417
214 28218.3265 43675.8982
215 99945.4123 28218.3265
216 132700.5090 99945.4123
217 125365.3889 132700.5090
218 18329.0421 125365.3889
219 2471.4575 18329.0421
220 -10213.5335 2471.4575
221 154648.0856 -10213.5335
222 41763.4498 154648.0856
223 15150.3021 41763.4498
224 138207.4323 15150.3021
225 44422.5117 138207.4323
226 76679.5657 44422.5117
227 21093.2038 76679.5657
228 47729.4636 21093.2038
229 49632.7602 47729.4636
230 25610.6683 49632.7602
231 130428.6334 25610.6683
232 324.4753 130428.6334
233 12274.3581 324.4753
234 412.3646 12274.3581
235 38221.7388 412.3646
236 -174553.0343 38221.7388
237 123848.6213 -174553.0343
238 104615.0985 123848.6213
239 1868.7387 104615.0985
240 89744.2163 1868.7387
241 -55925.2029 89744.2163
242 73957.3433 -55925.2029
243 128466.9310 73957.3433
244 -167325.2007 128466.9310
245 42024.1763 -167325.2007
246 132326.6958 42024.1763
247 133507.6425 132326.6958
248 117088.6655 133507.6425
249 11943.5225 117088.6655
250 47244.5077 11943.5225
251 118025.5208 47244.5077
252 22294.6438 118025.5208
253 -10291.5560 22294.6438
254 18965.2072 -10291.5560
255 31857.6418 18965.2072
256 83370.6252 31857.6418
257 38998.8280 83370.6252
258 59227.3013 38998.8280
259 50603.8573 59227.3013
260 2776.1763 50603.8573
261 -151092.5860 2776.1763
262 152638.1761 -151092.5860
263 23380.0287 152638.1761
264 -5569.0632 23380.0287
265 28679.7829 -5569.0632
266 156884.4521 28679.7829
267 42333.3868 156884.4521
268 34539.9273 42333.3868
269 9483.0331 34539.9273
270 104009.3825 9483.0331
271 45912.0144 104009.3825
272 130341.1689 45912.0144
273 19409.9300 130341.1689
274 124519.3013 19409.9300
275 104079.9597 124519.3013
276 4616.6885 104079.9597
277 38288.1777 4616.6885
278 131453.3728 38288.1777
279 8544.5835 131453.3728
280 -245065.2680 8544.5835
281 119708.2855 -245065.2680
282 79892.2718 119708.2855
283 73565.3265 79892.2718
284 3841.6443 73565.3265
285 119692.0345 3841.6443
286 88489.1070 119692.0345
287 106970.1481 88489.1070
288 -2234.1475 106970.1481
289 26377.8586 -2234.1475
290 109111.2094 26377.8586
291 128575.7542 109111.2094
292 68524.4566 128575.7542
293 55507.1617 68524.4566
294 -22791.3372 55507.1617
295 107569.0796 -22791.3372
296 8708.4356 107569.0796
297 -586.0401 8708.4356
298 63994.9277 -586.0401
299 45714.7608 63994.9277
300 107931.1100 45714.7608
301 88592.0011 107931.1100
302 74084.2139 88592.0011
303 35272.3096 74084.2139
304 -29392.3792 35272.3096
305 -402263.5485 -29392.3792
306 -265592.9731 -402263.5485
307 86319.3857 -265592.9731
308 -99418.7196 86319.3857
309 -135733.4010 -99418.7196
310 -991661.4399 -135733.4010
311 -114661.6884 -991661.4399
312 -353186.5569 -114661.6884
313 -524362.4068 -353186.5569
> plot(z,main=paste('Residual Lag plot, lowess, and regression line'), ylab='values of Residuals', xlab='lagged values of Residuals')
> lines(lowess(z))
> abline(lm(z))
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/7plq91291232200.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> acf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Autocorrelation Function')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/8plq91291232200.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> pacf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Partial Autocorrelation Function')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/9plq91291232200.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> opar <- par(mfrow = c(2,2), oma = c(0, 0, 1.1, 0))
> plot(mylm, las = 1, sub='Residual Diagnostics')
> par(opar)
> dev.off()
null device
1
> if (n > n25) {
+ postscript(file="/var/www/html/rcomp/tmp/10icpc1291232200.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
+ plot(kp3:nmkm3,gqarr[,2], main='Goldfeld-Quandt test',ylab='2-sided p-value',xlab='breakpoint')
+ grid()
+ dev.off()
+ }
null device
1
>
> #Note: the /var/www/html/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/www/html/rcomp/createtable")
>
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Estimated Regression Equation', 1, TRUE)
> a<-table.row.end(a)
> myeq <- colnames(x)[1]
> myeq <- paste(myeq, '[t] = ', sep='')
> for (i in 1:k){
+ if (mysum$coefficients[i,1] > 0) myeq <- paste(myeq, '+', '')
+ myeq <- paste(myeq, mysum$coefficients[i,1], sep=' ')
+ if (rownames(mysum$coefficients)[i] != '(Intercept)') {
+ myeq <- paste(myeq, rownames(mysum$coefficients)[i], sep='')
+ if (rownames(mysum$coefficients)[i] != 't') myeq <- paste(myeq, '[t]', sep='')
+ }
+ }
> myeq <- paste(myeq, ' + e[t]')
> a<-table.row.start(a)
> a<-table.element(a, myeq)
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/11luo01291232200.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,hyperlink('http://www.xycoon.com/ols1.htm','Multiple Linear Regression - Ordinary Least Squares',''), 6, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Variable',header=TRUE)
> a<-table.element(a,'Parameter',header=TRUE)
> a<-table.element(a,'S.D.',header=TRUE)
> a<-table.element(a,'T-STAT
H0: parameter = 0',header=TRUE)
> a<-table.element(a,'2-tail p-value',header=TRUE)
> a<-table.element(a,'1-tail p-value',header=TRUE)
> a<-table.row.end(a)
> for (i in 1:k){
+ a<-table.row.start(a)
+ a<-table.element(a,rownames(mysum$coefficients)[i],header=TRUE)
+ a<-table.element(a,mysum$coefficients[i,1])
+ a<-table.element(a, round(mysum$coefficients[i,2],6))
+ a<-table.element(a, round(mysum$coefficients[i,3],4))
+ a<-table.element(a, round(mysum$coefficients[i,4],6))
+ a<-table.element(a, round(mysum$coefficients[i,4]/2,6))
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/12pv5o1291232200.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Regression Statistics', 2, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple R',1,TRUE)
> a<-table.element(a, sqrt(mysum$r.squared))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'R-squared',1,TRUE)
> a<-table.element(a, mysum$r.squared)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Adjusted R-squared',1,TRUE)
> a<-table.element(a, mysum$adj.r.squared)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (value)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[1])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (DF numerator)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[2])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (DF denominator)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[3])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'p-value',1,TRUE)
> a<-table.element(a, 1-pf(mysum$fstatistic[1],mysum$fstatistic[2],mysum$fstatistic[3]))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Residual Statistics', 2, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Residual Standard Deviation',1,TRUE)
> a<-table.element(a, mysum$sigma)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Sum Squared Residuals',1,TRUE)
> a<-table.element(a, sum(myerror*myerror))
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/13lnke1291232200.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Actuals, Interpolation, and Residuals', 4, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Time or Index', 1, TRUE)
> a<-table.element(a, 'Actuals', 1, TRUE)
> a<-table.element(a, 'Interpolation
Forecast', 1, TRUE)
> a<-table.element(a, 'Residuals
Prediction Error', 1, TRUE)
> a<-table.row.end(a)
> for (i in 1:n) {
+ a<-table.row.start(a)
+ a<-table.element(a,i, 1, TRUE)
+ a<-table.element(a,x[i])
+ a<-table.element(a,x[i]-mysum$resid[i])
+ a<-table.element(a,mysum$resid[i])
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/14o5jk1291232200.tab")
> if (n > n25) {
+ a<-table.start()
+ a<-table.row.start(a)
+ a<-table.element(a,'Goldfeld-Quandt test for Heteroskedasticity',4,TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'p-values',header=TRUE)
+ a<-table.element(a,'Alternative Hypothesis',3,header=TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'breakpoint index',header=TRUE)
+ a<-table.element(a,'greater',header=TRUE)
+ a<-table.element(a,'2-sided',header=TRUE)
+ a<-table.element(a,'less',header=TRUE)
+ a<-table.row.end(a)
+ for (mypoint in kp3:nmkm3) {
+ a<-table.row.start(a)
+ a<-table.element(a,mypoint,header=TRUE)
+ a<-table.element(a,gqarr[mypoint-kp3+1,1])
+ a<-table.element(a,gqarr[mypoint-kp3+1,2])
+ a<-table.element(a,gqarr[mypoint-kp3+1,3])
+ a<-table.row.end(a)
+ }
+ a<-table.end(a)
+ table.save(a,file="/var/www/html/rcomp/tmp/15r6hq1291232200.tab")
+ a<-table.start()
+ a<-table.row.start(a)
+ a<-table.element(a,'Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity',4,TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'Description',header=TRUE)
+ a<-table.element(a,'# significant tests',header=TRUE)
+ a<-table.element(a,'% significant tests',header=TRUE)
+ a<-table.element(a,'OK/NOK',header=TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'1% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant1)
+ a<-table.element(a,numsignificant1/numgqtests)
+ if (numsignificant1/numgqtests < 0.01) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'5% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant5)
+ a<-table.element(a,numsignificant5/numgqtests)
+ if (numsignificant5/numgqtests < 0.05) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'10% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant10)
+ a<-table.element(a,numsignificant10/numgqtests)
+ if (numsignificant10/numgqtests < 0.1) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.end(a)
+ table.save(a,file="/var/www/html/rcomp/tmp/16d6gw1291232200.tab")
+ }
>
> try(system("convert tmp/1bbai1291232200.ps tmp/1bbai1291232200.png",intern=TRUE))
character(0)
> try(system("convert tmp/2bbai1291232200.ps tmp/2bbai1291232200.png",intern=TRUE))
character(0)
> try(system("convert tmp/3mkal1291232200.ps tmp/3mkal1291232200.png",intern=TRUE))
character(0)
> try(system("convert tmp/4mkal1291232200.ps tmp/4mkal1291232200.png",intern=TRUE))
character(0)
> try(system("convert tmp/5mkal1291232200.ps tmp/5mkal1291232200.png",intern=TRUE))
character(0)
> try(system("convert tmp/6xcr61291232200.ps tmp/6xcr61291232200.png",intern=TRUE))
character(0)
> try(system("convert tmp/7plq91291232200.ps tmp/7plq91291232200.png",intern=TRUE))
character(0)
> try(system("convert tmp/8plq91291232200.ps tmp/8plq91291232200.png",intern=TRUE))
character(0)
> try(system("convert tmp/9plq91291232200.ps tmp/9plq91291232200.png",intern=TRUE))
character(0)
> try(system("convert tmp/10icpc1291232200.ps tmp/10icpc1291232200.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
8.491 2.111 21.299