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(1
+ ,1595
+ ,17
+ ,60720
+ ,319369
+ ,0
+ ,5565
+ ,47
+ ,94355
+ ,493408
+ ,0
+ ,601
+ ,6
+ ,60720
+ ,319210
+ ,1
+ ,188
+ ,11
+ ,77655
+ ,381180
+ ,1
+ ,7146
+ ,105
+ ,134028
+ ,297978
+ ,0
+ ,1135
+ ,17
+ ,62285
+ ,290476
+ ,1
+ ,450
+ ,8
+ ,59325
+ ,292136
+ ,1
+ ,34
+ ,8
+ ,60630
+ ,314353
+ ,1
+ ,133
+ ,14
+ ,65990
+ ,339445
+ ,1
+ ,119
+ ,6
+ ,59118
+ ,303677
+ ,0
+ ,2053
+ ,53
+ ,100423
+ ,397144
+ ,0
+ ,4036
+ ,63
+ ,100269
+ ,424898
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,4655
+ ,393
+ ,60720
+ ,341570
+ ,1
+ ,131
+ ,11
+ ,61808
+ ,308989
+ ,1
+ ,1766
+ ,73
+ ,60438
+ ,305959
+ ,1
+ ,312
+ ,14
+ ,58598
+ ,318690
+ ,1
+ ,448
+ ,40
+ ,59781
+ ,323361
+ ,1
+ ,115
+ ,13
+ ,60945
+ ,318903
+ ,1
+ ,60
+ ,10
+ ,61124
+ ,314049
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,364
+ ,35
+ ,47705
+ ,298466
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1442
+ ,118
+ ,121173
+ ,438493
+ ,0
+ ,1389
+ ,9
+ ,57530
+ ,378049
+ ,1
+ ,149
+ ,5
+ ,62041
+ ,313332
+ ,1
+ ,2212
+ ,14
+ ,60720
+ ,319864
+ ,0
+ ,7489
+ ,322
+ ,136996
+ ,430866
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,402
+ ,26
+ ,84990
+ ,332743
+ ,0
+ ,7419
+ ,29
+ ,63255
+ ,286963
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,307
+ ,15
+ ,58856
+ ,331955
+ ,0
+ ,1134
+ ,162
+ ,146216
+ ,527021
+ ,0
+ ,7561
+ ,87
+ ,82425
+ ,364304
+ ,1
+ ,5131
+ ,71
+ ,86111
+ ,292154
+ ,1
+ ,9
+ ,6
+ ,60835
+ ,314987
+ ,1
+ ,2264
+ ,24
+ ,55174
+ ,272713
+ ,1
+ ,40949
+ ,481
+ ,129352
+ ,1398893
+ ,0
+ ,4980
+ ,91
+ ,113521
+ ,409642
+ ,1
+ ,272
+ ,8
+ ,112995
+ ,429112
+ ,1
+ ,757
+ ,19
+ ,45689
+ ,382712
+ ,1
+ ,1424
+ ,94
+ ,131116
+ ,374943
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,378
+ ,52
+ ,71873
+ ,297413
+ ,1
+ ,1061
+ ,5
+ ,70415
+ ,304252
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,203
+ ,9
+ ,61938
+ ,309596
+ ,0
+ ,3689
+ ,25
+ ,101792
+ ,377305
+ ,1
+ ,567
+ ,11
+ ,103345
+ ,386212
+ ,0
+ ,27664
+ ,282
+ ,74996
+ ,657954
+ ,1
+ ,1686
+ ,40
+ ,68696
+ ,368186
+ ,0
+ ,1164
+ ,17
+ ,57320
+ ,269753
+ ,1
+ ,15824
+ ,64
+ ,63838
+ ,585715
+ ,1
+ ,6575
+ ,159
+ ,151352
+ ,1071292
+ ,0
+ ,37597
+ ,248
+ ,37884
+ ,992426
+ ,1
+ ,314
+ ,5
+ ,60720
+ ,315371
+ ,0
+ ,1932
+ ,36
+ ,72954
+ ,249898
+ ,0
+ ,5692
+ ,72
+ ,64239
+ ,357312
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,6533
+ ,75
+ ,126910
+ ,364839
+ ,0
+ ,2126
+ ,138
+ ,21001
+ ,230621
+ ,1
+ ,14
+ ,1
+ ,61016
+ ,315877
+ ,0
+ ,3192
+ ,60
+ ,59831
+ ,125390
+ ,1
+ ,47
+ ,3
+ ,60720
+ ,314882
+ ,0
+ ,4891
+ ,217
+ ,210568
+ ,370837
+ ,0
+ ,531
+ ,10
+ ,71641
+ ,330068
+ ,1
+ ,2332
+ ,4
+ ,72680
+ ,324385
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,16275
+ ,109
+ ,47818
+ ,1443586
+ ,1
+ ,9252
+ ,48
+ ,79510
+ ,253537
+ ,0
+ ,29790
+ ,309
+ ,81767
+ ,4321023
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,2300
+ ,29
+ ,64873
+ ,352108
+ ,1
+ ,267
+ ,32
+ ,57640
+ ,330059
+ ,0
+ ,382
+ ,21
+ ,59365
+ ,340968
+ ,1
+ ,397
+ ,28
+ ,60646
+ ,317736
+ ,1
+ ,11920
+ ,182
+ ,31080
+ ,209458
+ ,0
+ ,54660
+ ,302
+ ,126942
+ ,1491348
+ ,1
+ ,5
+ ,2
+ ,60720
+ ,314887
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,2647
+ ,37
+ ,62920
+ ,353058
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,94
+ ,7
+ ,60793
+ ,314533
+ ,1
+ ,3
+ ,2
+ ,60698
+ ,315354
+ ,1
+ ,454
+ ,51
+ ,63261
+ ,302187
+ ,1
+ ,227
+ ,9
+ ,76644
+ ,336639
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,206
+ ,9
+ ,53110
+ ,296702
+ ,1
+ ,25830
+ ,115
+ ,245546
+ ,1073089
+ ,1
+ ,1528
+ ,38
+ ,60326
+ ,146494
+ ,1
+ ,271
+ ,2
+ ,60720
+ ,325249
+ ,1
+ ,138
+ ,6
+ ,69817
+ ,331420
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,278
+ ,26
+ ,61404
+ ,314922
+ ,1
+ ,282
+ ,30
+ ,62452
+ ,320016
+ ,1
+ ,571
+ ,80
+ ,41477
+ ,280398
+ ,1
+ ,2253
+ ,99
+ ,63593
+ ,452469
+ ,1
+ ,290
+ ,2
+ ,58790
+ ,301164
+ ,1
+ ,78
+ ,16
+ ,62700
+ ,317330
+ ,1
+ ,20
+ ,2
+ ,60805
+ ,315576
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,18073
+ ,180
+ ,27284
+ ,-7170
+ ,1
+ ,1866
+ ,29
+ ,56225
+ ,322331
+ ,0
+ ,42634
+ ,156
+ ,157214
+ ,1629616
+ ,1
+ ,249
+ ,62
+ ,54323
+ ,292754
+ ,1
+ ,422
+ ,28
+ ,57935
+ ,318056
+ ,1
+ ,2675
+ ,79
+ ,59017
+ ,355178
+ ,1
+ ,965
+ ,3
+ ,73490
+ ,204325
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,621
+ ,16
+ ,68005
+ ,317046
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,365
+ ,8
+ ,54820
+ ,309560
+ ,0
+ ,3122
+ ,34
+ ,94670
+ ,414462
+ ,0
+ ,12988
+ ,119
+ ,82340
+ ,857217
+ ,0
+ ,5336
+ ,81
+ ,112477
+ ,697458
+ ,0
+ ,3160
+ ,108
+ ,108094
+ ,530670
+ ,1
+ ,2489
+ ,40
+ ,67804
+ ,238125
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,6984
+ ,45
+ ,80570
+ ,741409
+ ,1
+ ,2001
+ ,52
+ ,95551
+ ,393343
+ ,1
+ ,15236
+ ,90
+ ,77440
+ ,372631
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,530
+ ,60
+ ,73433
+ ,317291
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,35624
+ ,223
+ ,157278
+ ,306275
+ ,0
+ ,1383
+ ,11
+ ,73221
+ ,317892
+ ,1
+ ,875
+ ,14
+ ,67000
+ ,334280
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,72
+ ,3
+ ,60398
+ ,314210
+ ,0
+ ,265
+ ,31
+ ,60720
+ ,306948
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,335
+ ,33
+ ,64175
+ ,320398
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,4525
+ ,226
+ ,93811
+ ,501749
+ ,0
+ ,3045
+ ,58
+ ,27330
+ ,202055
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,638
+ ,14
+ ,60370
+ ,333210
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,607
+ ,11
+ ,60436
+ ,322340
+ ,1
+ ,1558
+ ,63
+ ,55637
+ ,369448
+ ,0
+ ,1324
+ ,117
+ ,67440
+ ,291841
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,611
+ ,8
+ ,59190
+ ,296919
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,923
+ ,8
+ ,58620
+ ,309038
+ ,1
+ ,661
+ ,3
+ ,65920
+ ,246541
+ ,0
+ ,2397
+ ,14
+ ,74020
+ ,289513
+ ,1
+ ,366
+ ,53
+ ,61808
+ ,344425
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,135
+ ,2
+ ,62065
+ ,314210
+ ,0
+ ,1659
+ ,29
+ ,107577
+ ,480382
+ ,1
+ ,316
+ ,9
+ ,60505
+ ,315009
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,309
+ ,11
+ ,56535
+ ,312878
+ ,0
+ ,49
+ ,8
+ ,64107
+ ,322031
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,4519
+ ,49
+ ,102129
+ ,597793
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,837
+ ,69
+ ,61262
+ ,315688
+ ,1
+ ,5119
+ ,49
+ ,39039
+ ,378525
+ ,1
+ ,1280
+ ,117
+ ,69465
+ ,312378
+ ,1
+ ,2564
+ ,22
+ ,130140
+ ,403560
+ ,0
+ ,2045
+ ,20
+ ,97890
+ ,510834
+ ,0
+ ,2234
+ ,16
+ ,77200
+ ,214215
+ ,1
+ ,975
+ ,23
+ ,90534
+ ,235133
+ ,0
+ ,1136
+ ,32
+ ,48522
+ ,343613
+ ,1
+ ,453
+ ,21
+ ,81125
+ ,365959
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,61
+ ,14
+ ,60720
+ ,314551
+ ,1
+ ,368
+ ,17
+ ,61686
+ ,303230
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,4901
+ ,84
+ ,121920
+ ,469107
+ ,1
+ ,540
+ ,14
+ ,103960
+ ,354228
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,2
+ ,9
+ ,60735
+ ,315394
+ ,1
+ ,36
+ ,4
+ ,61564
+ ,312412
+ ,1
+ ,776
+ ,9
+ ,64230
+ ,333505
+ ,1
+ ,84738
+ ,588
+ ,-26007
+ ,223193
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,3
+ ,4
+ ,60761
+ ,315656
+ ,1
+ ,529
+ ,9
+ ,63870
+ ,296261
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,405
+ ,7
+ ,60845
+ ,336425
+ ,1
+ ,972
+ ,34
+ ,71642
+ ,359335
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,2099
+ ,64
+ ,106611
+ ,308636
+ ,0
+ ,3437
+ ,47
+ ,48022
+ ,158492
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,22330
+ ,84
+ ,79801
+ ,711969
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,483
+ ,27
+ ,60830
+ ,306268
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,2239
+ ,21
+ ,88590
+ ,442882
+ ,0
+ ,2949
+ ,41
+ ,82903
+ ,378509
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,365
+ ,83
+ ,87192
+ ,346611
+ ,1
+ ,2461
+ ,57
+ ,55792
+ ,314289
+ ,0
+ ,21950
+ ,519
+ ,114337
+ ,856956
+ ,0
+ ,3294
+ ,23
+ ,75832
+ ,217193
+ ,1
+ ,141
+ ,6
+ ,61630
+ ,315366
+ ,1
+ ,572
+ ,16
+ ,58580
+ ,307930
+ ,1
+ ,13326
+ ,102
+ ,165548
+ ,702380
+ ,1
+ ,2284
+ ,33
+ ,76403
+ ,194493
+ ,0
+ ,10
+ ,2
+ ,61656
+ ,316155
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1414
+ ,198
+ ,70184
+ ,330546
+ ,1
+ ,1975
+ ,35
+ ,118881
+ ,394510
+ ,1
+ ,43
+ ,3
+ ,60887
+ ,312846
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,844
+ ,78
+ ,60925
+ ,296139
+ ,1
+ ,304
+ ,15
+ ,62969
+ ,295580
+ ,1
+ ,458
+ ,11
+ ,58625
+ ,297765
+ ,1
+ ,18562
+ ,155
+ ,102313
+ ,377934
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,7123
+ ,109
+ ,288170
+ ,638830
+ ,1
+ ,622
+ ,73
+ ,73007
+ ,304376
+ ,1
+ ,174
+ ,1
+ ,64820
+ ,307424
+ ,1
+ ,2220
+ ,22
+ ,301670
+ ,644190
+ ,1
+ ,121
+ ,24
+ ,56178
+ ,295370
+ ,0
+ ,11819
+ ,85
+ ,106113
+ ,574339
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,125
+ ,12
+ ,60798
+ ,310201
+ ,1
+ ,1182
+ ,13
+ ,70694
+ ,327007
+ ,1
+ ,1503
+ ,39
+ ,56364
+ ,343466
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,30
+ ,4
+ ,62045
+ ,318098
+ ,0
+ ,3310
+ ,33
+ ,75230
+ ,448243
+ ,1
+ ,554
+ ,11
+ ,79285
+ ,325738
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,468
+ ,21
+ ,60720
+ ,312161
+ ,1
+ ,4917
+ ,42
+ ,52811
+ ,243650
+ ,1
+ ,3256
+ ,66
+ ,35250
+ ,407159
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,125
+ ,13
+ ,59734
+ ,317698
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,22
+ ,2
+ ,60722
+ ,312502
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,514
+ ,52
+ ,78780
+ ,322378
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,10327
+ ,91
+ ,88577
+ ,640273
+ ,1
+ ,13718
+ ,61
+ ,96448
+ ,345783
+ ,1
+ ,3748
+ ,9
+ ,50350
+ ,652925
+ ,1
+ ,14416
+ ,361
+ ,49857
+ ,439798
+ ,0
+ ,1526
+ ,25
+ ,69351
+ ,278990
+ ,0
+ ,666
+ ,171
+ ,117869
+ ,339836
+ ,1
+ ,2844
+ ,67
+ ,72683
+ ,240897
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,368
+ ,17
+ ,61167
+ ,297141
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,333
+ ,13
+ ,70811
+ ,331323
+ ,1
+ ,26
+ ,1
+ ,60896
+ ,313880
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,1303
+ ,61
+ ,69863
+ ,309422
+ ,1
+ ,20
+ ,3
+ ,60938
+ ,315245
+ ,1
+ ,2384
+ ,97
+ ,61348
+ ,405972
+ ,0
+ ,203
+ ,27
+ ,50804
+ ,300962
+ ,0
+ ,71
+ ,1
+ ,60745
+ ,316176
+ ,1
+ ,53
+ ,14
+ ,59506
+ ,302409
+ ,1
+ ,562
+ ,13
+ ,58456
+ ,283587
+ ,1
+ ,622
+ ,7
+ ,60950
+ ,263276
+ ,1
+ ,645
+ ,5
+ ,60720
+ ,312075
+ ,0
+ ,1763
+ ,24
+ ,61600
+ ,308336
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,317
+ ,9
+ ,63915
+ ,298700
+ ,1
+ ,1
+ ,6
+ ,60719
+ ,315372
+ ,1
+ ,275
+ ,4
+ ,59500
+ ,318745
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,936
+ ,30
+ ,67939
+ ,408881
+ ,1
+ ,8568
+ ,37
+ ,32168
+ ,786690
+ ,0
+ ,11528
+ ,70
+ ,-14545
+ ,-83265
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,738
+ ,10
+ ,60720
+ ,321376
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,592
+ ,7
+ ,64270
+ ,276898
+ ,1
+ ,126
+ ,10
+ ,60951
+ ,315547
+ ,1
+ ,2
+ ,1
+ ,60743
+ ,315487
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,18014
+ ,239
+ ,-1710
+ ,1405225
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,37704
+ ,452
+ ,60448
+ ,983660
+ ,0
+ ,63
+ ,4
+ ,65688
+ ,318574
+ ,1
+ ,1431
+ ,34
+ ,106885
+ ,310768
+ ,1
+ ,94
+ ,3
+ ,61360
+ ,312887
+ ,1
+ ,192
+ ,4
+ ,65276
+ ,312339
+ ,1
+ ,32
+ ,1
+ ,59988
+ ,314964
+ ,1
+ ,6869
+ ,105
+ ,117520
+ ,379983
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,1
+ ,4
+ ,60722
+ ,315398
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,2328
+ ,94
+ ,82732
+ ,253588
+ ,0
+ ,209
+ ,13
+ ,64016
+ ,316647
+ ,1
+ ,28
+ ,2
+ ,60890
+ ,315688
+ ,1
+ ,176
+ ,28
+ ,68136
+ ,310670
+ ,1
+ ,1920
+ ,15
+ ,79420
+ ,165404
+ ,0
+ ,87550
+ ,458
+ ,153198
+ ,4111912
+ ,1
+ ,520
+ ,10
+ ,58650
+ ,291650
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,1013
+ ,2
+ ,69770
+ ,253468
+ ,1
+ ,15
+ ,8
+ ,60831
+ ,315688
+ ,1
+ ,587
+ ,12
+ ,59595
+ ,325699
+ ,1
+ ,5371
+ ,66
+ ,87720
+ ,446211
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,1012
+ ,47
+ ,114768
+ ,368078
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,876
+ ,59
+ ,138971
+ ,352850
+ ,0
+ ,162556
+ ,1081
+ ,213118
+ ,6282154
+ ,1
+ ,43556
+ ,247
+ ,32648
+ ,227132
+ ,0
+ ,3425
+ ,12
+ ,83620
+ ,283910
+ ,1
+ ,810
+ ,43
+ ,74015
+ ,236761
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,2365
+ ,79
+ ,191778
+ ,550608
+ ,0
+ ,1261
+ ,33
+ ,76114
+ ,307528
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,1585
+ ,107
+ ,93099
+ ,355864
+ ,1
+ ,16189
+ ,295
+ ,116384
+ ,358589
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,10579
+ ,92
+ ,110309
+ ,375195
+ ,1
+ ,474
+ ,20
+ ,61977
+ ,288985
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,3642
+ ,86
+ ,90262
+ ,458343
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,472
+ ,20
+ ,80045
+ ,269587
+ ,1
+ ,98
+ ,2
+ ,61490
+ ,315236
+ ,1
+ ,3999
+ ,30
+ ,51252
+ ,42754
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,621
+ ,48
+ ,60720
+ ,308256
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,30
+ ,4
+ ,60798
+ ,313164
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,746
+ ,23
+ ,71561
+ ,269661
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,4150
+ ,150
+ ,134759
+ ,518365
+ ,0
+ ,4658
+ ,158
+ ,156608
+ ,233773
+ ,1
+ ,814
+ ,18
+ ,39625
+ ,301881
+ ,1
+ ,1002
+ ,15
+ ,56750
+ ,298568
+ ,0
+ ,496
+ ,17
+ ,87390
+ ,325479
+ ,1
+ ,389
+ ,14
+ ,58990
+ ,325506
+ ,1
+ ,12679
+ ,110
+ ,48020
+ ,984885
+ ,1
+ ,400
+ ,152
+ ,60720
+ ,313267
+ ,1
+ ,53
+ ,6
+ ,60349
+ ,315793
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,5109
+ ,97
+ ,67038
+ ,215362
+ ,1
+ ,576
+ ,18
+ ,113761
+ ,314073
+ ,1
+ ,438
+ ,11
+ ,58320
+ ,298096
+ ,1
+ ,165
+ ,24
+ ,62841
+ ,325176
+ ,1
+ ,4069
+ ,145
+ ,79804
+ ,207393
+ ,0
+ ,23
+ ,5
+ ,62555
+ ,314806
+ ,1
+ ,1285
+ ,107
+ ,99489
+ ,341340
+ ,0
+ ,4677
+ ,75
+ ,90131
+ ,426280
+ ,0
+ ,24811
+ ,505
+ ,95350
+ ,929118
+ ,1
+ ,167
+ ,2
+ ,64245
+ ,307322
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,4628
+ ,62
+ ,69159
+ ,387475
+ ,1
+ ,226
+ ,5
+ ,65745
+ ,291787
+ ,0
+ ,1765
+ ,20
+ ,77623
+ ,247060
+ ,1
+ ,460
+ ,15
+ ,63346
+ ,329784
+ ,1
+ ,36
+ ,3
+ ,60894
+ ,315834
+ ,1
+ ,989
+ ,11
+ ,58930
+ ,304555
+ ,0
+ ,3055
+ ,43
+ ,60247
+ ,376641
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,13253
+ ,947
+ ,90829
+ ,357760
+ ,1
+ ,24
+ ,8
+ ,59818
+ ,315637
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,131
+ ,6
+ ,59661
+ ,327071
+ ,1
+ ,514
+ ,15
+ ,102725
+ ,377516
+ ,1
+ ,3366
+ ,80
+ ,108479
+ ,299243
+ ,0
+ ,9327
+ ,101
+ ,87419
+ ,387699
+ ,1
+ ,384
+ ,22
+ ,54683
+ ,309836
+ ,1
+ ,17821
+ ,373
+ ,122844
+ ,444477
+ ,1
+ ,397
+ ,10
+ ,62710
+ ,322327
+ ,1
+ ,897
+ ,5
+ ,60720
+ ,314913
+ ,1
+ ,218
+ ,21
+ ,60720
+ ,315553
+ ,0
+ ,3369
+ ,28
+ ,87161
+ ,688779
+ ,1
+ ,5702
+ ,61
+ ,101481
+ ,321896
+ ,1
+ ,8636
+ ,855
+ ,128294
+ ,301607
+ ,1
+ ,534
+ ,7
+ ,62620
+ ,304485
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,726
+ ,3
+ ,69980
+ ,231861
+ ,1
+ ,1380
+ ,26
+ ,60982
+ ,347385
+ ,1
+ ,180
+ ,11
+ ,59635
+ ,316386
+ ,0
+ ,7285
+ ,115
+ ,188873
+ ,491303
+ ,1
+ ,880
+ ,74
+ ,80791
+ ,261216
+ ,1
+ ,1
+ ,2
+ ,60727
+ ,315388
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,96
+ ,5
+ ,60379
+ ,313729
+ ,1
+ ,1889
+ ,45
+ ,37527
+ ,358649
+ ,0
+ ,45187
+ ,353
+ ,234817
+ ,1926517
+ ,1
+ ,288
+ ,4
+ ,60510
+ ,296656
+ ,1
+ ,1270
+ ,26
+ ,69206
+ ,275311
+ ,1
+ ,6526
+ ,26
+ ,55830
+ ,-42143
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,226
+ ,12
+ ,72835
+ ,343929
+ ,1
+ ,694
+ ,8
+ ,68060
+ ,367655
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,249
+ ,72
+ ,56726
+ ,313491)
+ ,dim=c(5
+ ,431)
+ ,dimnames=list(c('group'
+ ,'costs'
+ ,'trades'
+ ,'dividends'
+ ,'wealth')
+ ,1:431))
> y <- array(NA,dim=c(5,431),dimnames=list(c('group','costs','trades','dividends','wealth'),1:431))
> for (i in 1:dim(x)[1])
+ {
+ for (j in 1:dim(x)[2])
+ {
+ y[i,j] <- as.numeric(x[i,j])
+ }
+ }
> par20 = ''
> par19 = ''
> par18 = ''
> par17 = ''
> par16 = ''
> par15 = ''
> par14 = ''
> par13 = ''
> par12 = ''
> par11 = ''
> par10 = ''
> par9 = ''
> par8 = ''
> par7 = ''
> par6 = ''
> par5 = ''
> par4 = ''
> par3 = 'No Linear Trend'
> par2 = 'Do not include Seasonal Dummies'
> par1 = '5'
> ylab = ''
> xlab = ''
> main = ''
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Prof. Dr. P. Wessa
> #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #Technical description: Write here your technical program description (don't use hard returns!)
> library(lattice)
> library(lmtest)
Loading required package: zoo
> n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test
> par1 <- as.numeric(par1)
> x <- t(y)
> k <- length(x[1,])
> n <- length(x[,1])
> x1 <- cbind(x[,par1], x[,1:k!=par1])
> mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1])
> colnames(x1) <- mycolnames #colnames(x)[par1]
> x <- x1
> if (par3 == 'First Differences'){
+ x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep='')))
+ for (i in 1:n-1) {
+ for (j in 1:k) {
+ x2[i,j] <- x[i+1,j] - x[i,j]
+ }
+ }
+ x <- x2
+ }
> if (par2 == 'Include Monthly Dummies'){
+ x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep ='')))
+ for (i in 1:11){
+ x2[seq(i,n,12),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> if (par2 == 'Include Quarterly Dummies'){
+ x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep ='')))
+ for (i in 1:3){
+ x2[seq(i,n,4),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> k <- length(x[1,])
> if (par3 == 'Linear Trend'){
+ x <- cbind(x, c(1:n))
+ colnames(x)[k+1] <- 't'
+ }
> x
wealth group costs trades dividends
1 319369 1 1595 17 60720
2 493408 0 5565 47 94355
3 319210 0 601 6 60720
4 381180 1 188 11 77655
5 297978 1 7146 105 134028
6 290476 0 1135 17 62285
7 292136 1 450 8 59325
8 314353 1 34 8 60630
9 339445 1 133 14 65990
10 303677 1 119 6 59118
11 397144 0 2053 53 100423
12 424898 0 4036 63 100269
13 315380 1 0 0 60720
14 341570 1 4655 393 60720
15 308989 1 131 11 61808
16 305959 1 1766 73 60438
17 318690 1 312 14 58598
18 323361 1 448 40 59781
19 318903 1 115 13 60945
20 314049 1 60 10 61124
21 315380 1 0 0 60720
22 298466 1 364 35 47705
23 315380 1 0 0 60720
24 438493 0 1442 118 121173
25 378049 0 1389 9 57530
26 313332 1 149 5 62041
27 319864 1 2212 14 60720
28 430866 0 7489 322 136996
29 315380 1 0 0 60720
30 332743 1 402 26 84990
31 286963 0 7419 29 63255
32 315380 1 0 0 60720
33 331955 1 307 15 58856
34 527021 0 1134 162 146216
35 364304 0 7561 87 82425
36 292154 1 5131 71 86111
37 314987 1 9 6 60835
38 272713 1 2264 24 55174
39 1398893 1 40949 481 129352
40 409642 0 4980 91 113521
41 429112 1 272 8 112995
42 382712 1 757 19 45689
43 374943 1 1424 94 131116
44 315380 1 0 0 60720
45 297413 0 378 52 71873
46 304252 1 1061 5 70415
47 315380 1 0 0 60720
48 309596 1 203 9 61938
49 377305 0 3689 25 101792
50 386212 1 567 11 103345
51 657954 0 27664 282 74996
52 368186 1 1686 40 68696
53 269753 0 1164 17 57320
54 585715 1 15824 64 63838
55 1071292 1 6575 159 151352
56 992426 0 37597 248 37884
57 315371 1 314 5 60720
58 249898 0 1932 36 72954
59 357312 0 5692 72 64239
60 315380 1 0 0 60720
61 315380 1 0 0 60720
62 364839 0 6533 75 126910
63 230621 0 2126 138 21001
64 315877 1 14 1 61016
65 125390 0 3192 60 59831
66 314882 1 47 3 60720
67 370837 0 4891 217 210568
68 330068 0 531 10 71641
69 324385 1 2332 4 72680
70 315380 1 0 0 60720
71 1443586 0 16275 109 47818
72 253537 1 9252 48 79510
73 4321023 0 29790 309 81767
74 315380 1 0 0 60720
75 352108 0 2300 29 64873
76 330059 1 267 32 57640
77 340968 0 382 21 59365
78 317736 1 397 28 60646
79 209458 1 11920 182 31080
80 1491348 0 54660 302 126942
81 314887 1 5 2 60720
82 315380 1 0 0 60720
83 315380 1 0 0 60720
84 353058 0 2647 37 62920
85 315380 1 0 0 60720
86 315380 0 0 0 60720
87 314533 1 94 7 60793
88 315354 1 3 2 60698
89 302187 1 454 51 63261
90 336639 1 227 9 76644
91 315380 0 0 0 60720
92 296702 1 206 9 53110
93 1073089 1 25830 115 245546
94 146494 1 1528 38 60326
95 325249 1 271 2 60720
96 331420 1 138 6 69817
97 315380 1 0 0 60720
98 314922 1 278 26 61404
99 320016 1 282 30 62452
100 280398 1 571 80 41477
101 452469 1 2253 99 63593
102 301164 1 290 2 58790
103 317330 1 78 16 62700
104 315576 1 20 2 60805
105 315380 1 0 0 60720
106 -7170 1 18073 180 27284
107 322331 1 1866 29 56225
108 1629616 0 42634 156 157214
109 292754 1 249 62 54323
110 318056 1 422 28 57935
111 355178 1 2675 79 59017
112 204325 1 965 3 73490
113 315380 1 0 0 60720
114 317046 1 621 16 68005
115 315380 0 0 0 60720
116 309560 1 365 8 54820
117 414462 0 3122 34 94670
118 857217 0 12988 119 82340
119 697458 0 5336 81 112477
120 530670 0 3160 108 108094
121 238125 1 2489 40 67804
122 315380 1 0 0 60720
123 741409 0 6984 45 80570
124 393343 1 2001 52 95551
125 372631 1 15236 90 77440
126 315380 0 0 0 60720
127 317291 1 530 60 73433
128 315380 1 0 0 60720
129 306275 1 35624 223 157278
130 317892 0 1383 11 73221
131 334280 1 875 14 67000
132 315380 0 0 0 60720
133 315380 0 0 0 60720
134 315380 0 0 0 60720
135 314210 1 72 3 60398
136 306948 0 265 31 60720
137 315380 1 0 0 60720
138 320398 1 335 33 64175
139 315380 0 0 0 60720
140 315380 1 0 0 60720
141 501749 1 4525 226 93811
142 202055 0 3045 58 27330
143 315380 1 0 0 60720
144 315380 1 0 0 60720
145 333210 0 638 14 60370
146 315380 0 0 0 60720
147 322340 1 607 11 60436
148 369448 1 1558 63 55637
149 291841 0 1324 117 67440
150 315380 1 0 0 60720
151 315380 0 0 0 60720
152 296919 1 611 8 59190
153 315380 1 0 0 60720
154 309038 1 923 8 58620
155 246541 1 661 3 65920
156 289513 0 2397 14 74020
157 344425 1 366 53 61808
158 315380 1 0 0 60720
159 314210 1 135 2 62065
160 480382 0 1659 29 107577
161 315009 1 316 9 60505
162 315380 1 0 0 60720
163 312878 1 309 11 56535
164 322031 0 49 8 64107
165 315380 1 0 0 60720
166 597793 0 4519 49 102129
167 315380 1 0 0 60720
168 315688 1 837 69 61262
169 378525 1 5119 49 39039
170 312378 1 1280 117 69465
171 403560 1 2564 22 130140
172 510834 0 2045 20 97890
173 214215 0 2234 16 77200
174 235133 1 975 23 90534
175 343613 0 1136 32 48522
176 365959 1 453 21 81125
177 315380 1 0 0 60720
178 314551 1 61 14 60720
179 303230 1 368 17 61686
180 315380 0 0 0 60720
181 469107 1 4901 84 121920
182 354228 1 540 14 103960
183 315380 1 0 0 60720
184 315380 1 0 0 60720
185 315380 0 0 0 60720
186 315394 1 2 9 60735
187 312412 1 36 4 61564
188 333505 1 776 9 64230
189 223193 1 84738 588 -26007
190 315380 1 0 0 60720
191 315656 1 3 4 60761
192 296261 1 529 9 63870
193 315380 1 0 0 60720
194 336425 1 405 7 60845
195 359335 1 972 34 71642
196 315380 1 0 0 60720
197 308636 1 2099 64 106611
198 158492 0 3437 47 48022
199 315380 0 0 0 60720
200 315380 1 0 0 60720
201 711969 0 22330 84 79801
202 315380 0 0 0 60720
203 315380 0 0 0 60720
204 306268 1 483 27 60830
205 315380 1 0 0 60720
206 442882 1 2239 21 88590
207 378509 0 2949 41 82903
208 315380 1 0 0 60720
209 346611 1 365 83 87192
210 314289 1 2461 57 55792
211 856956 0 21950 519 114337
212 217193 0 3294 23 75832
213 315366 1 141 6 61630
214 307930 1 572 16 58580
215 702380 1 13326 102 165548
216 194493 1 2284 33 76403
217 316155 0 10 2 61656
218 315380 1 0 0 60720
219 330546 0 1414 198 70184
220 394510 1 1975 35 118881
221 312846 1 43 3 60887
222 315380 1 0 0 60720
223 296139 1 844 78 60925
224 295580 1 304 15 62969
225 297765 1 458 11 58625
226 377934 1 18562 155 102313
227 315380 1 0 0 60720
228 638830 1 7123 109 288170
229 304376 1 622 73 73007
230 307424 1 174 1 64820
231 644190 1 2220 22 301670
232 295370 1 121 24 56178
233 574339 0 11819 85 106113
234 315380 1 0 0 60720
235 310201 1 125 12 60798
236 327007 1 1182 13 70694
237 343466 1 1503 39 56364
238 315380 1 0 0 60720
239 315380 1 0 0 60720
240 318098 1 30 4 62045
241 448243 0 3310 33 75230
242 325738 1 554 11 79285
243 315380 0 0 0 60720
244 312161 1 468 21 60720
245 243650 1 4917 42 52811
246 407159 1 3256 66 35250
247 315380 1 0 0 60720
248 317698 1 125 13 59734
249 315380 1 0 0 60720
250 312502 1 22 2 60722
251 315380 0 0 0 60720
252 322378 1 514 52 78780
253 315380 0 0 0 60720
254 315380 1 0 0 60720
255 315380 0 0 0 60720
256 315380 1 0 0 60720
257 315380 1 0 0 60720
258 315380 1 0 0 60720
259 640273 1 10327 91 88577
260 345783 1 13718 61 96448
261 652925 1 3748 9 50350
262 439798 1 14416 361 49857
263 278990 0 1526 25 69351
264 339836 0 666 171 117869
265 240897 1 2844 67 72683
266 315380 1 0 0 60720
267 297141 1 368 17 61167
268 315380 1 0 0 60720
269 331323 1 333 13 70811
270 313880 1 26 1 60896
271 315380 1 0 0 60720
272 315380 1 0 0 60720
273 309422 1 1303 61 69863
274 315245 1 20 3 60938
275 405972 1 2384 97 61348
276 300962 0 203 27 50804
277 316176 0 71 1 60745
278 302409 1 53 14 59506
279 283587 1 562 13 58456
280 263276 1 622 7 60950
281 312075 1 645 5 60720
282 308336 0 1763 24 61600
283 315380 1 0 0 60720
284 298700 1 317 9 63915
285 315372 1 1 6 60719
286 318745 1 275 4 59500
287 315380 1 0 0 60720
288 408881 1 936 30 67939
289 786690 1 8568 37 32168
290 -83265 0 11528 70 -14545
291 315380 0 0 0 60720
292 321376 1 738 10 60720
293 315380 1 0 0 60720
294 276898 1 592 7 64270
295 315547 1 126 10 60951
296 315487 1 2 1 60743
297 315380 1 0 0 60720
298 1405225 0 18014 239 -1710
299 315380 0 0 0 60720
300 315380 0 0 0 60720
301 983660 0 37704 452 60448
302 318574 0 63 4 65688
303 310768 1 1431 34 106885
304 312887 1 94 3 61360
305 312339 1 192 4 65276
306 314964 1 32 1 59988
307 379983 1 6869 105 117520
308 315380 1 0 0 60720
309 315380 1 0 0 60720
310 315398 1 1 4 60722
311 315380 0 0 0 60720
312 253588 1 2328 94 82732
313 316647 0 209 13 64016
314 315688 1 28 2 60890
315 310670 1 176 28 68136
316 165404 1 1920 15 79420
317 4111912 0 87550 458 153198
318 291650 1 520 10 58650
319 315380 0 0 0 60720
320 253468 1 1013 2 69770
321 315688 1 15 8 60831
322 325699 1 587 12 59595
323 446211 1 5371 66 87720
324 315380 1 0 0 60720
325 368078 1 1012 47 114768
326 315380 1 0 0 60720
327 352850 1 876 59 138971
328 6282154 0 162556 1081 213118
329 227132 1 43556 247 32648
330 283910 0 3425 12 83620
331 236761 1 810 43 74015
332 315380 1 0 0 60720
333 315380 0 0 0 60720
334 550608 0 2365 79 191778
335 307528 0 1261 33 76114
336 315380 1 0 0 60720
337 355864 1 1585 107 93099
338 358589 1 16189 295 116384
339 315380 1 0 0 60720
340 315380 0 0 0 60720
341 375195 1 10579 92 110309
342 288985 1 474 20 61977
343 315380 1 0 0 60720
344 458343 1 3642 86 90262
345 315380 1 0 0 60720
346 269587 0 472 20 80045
347 315236 1 98 2 61490
348 42754 1 3999 30 51252
349 315380 1 0 0 60720
350 308256 1 621 48 60720
351 315380 1 0 0 60720
352 315380 1 0 0 60720
353 313164 0 30 4 60798
354 315380 1 0 0 60720
355 269661 1 746 23 71561
356 315380 0 0 0 60720
357 315380 1 0 0 60720
358 315380 0 0 0 60720
359 315380 1 0 0 60720
360 518365 0 4150 150 134759
361 233773 0 4658 158 156608
362 301881 1 814 18 39625
363 298568 1 1002 15 56750
364 325479 0 496 17 87390
365 325506 1 389 14 58990
366 984885 1 12679 110 48020
367 313267 1 400 152 60720
368 315793 1 53 6 60349
369 315380 1 0 0 60720
370 215362 0 5109 97 67038
371 314073 1 576 18 113761
372 298096 1 438 11 58320
373 325176 1 165 24 62841
374 207393 1 4069 145 79804
375 314806 0 23 5 62555
376 341340 1 1285 107 99489
377 426280 0 4677 75 90131
378 929118 0 24811 505 95350
379 307322 1 167 2 64245
380 315380 1 0 0 60720
381 387475 1 4628 62 69159
382 291787 1 226 5 65745
383 247060 0 1765 20 77623
384 329784 1 460 15 63346
385 315834 1 36 3 60894
386 304555 1 989 11 58930
387 376641 0 3055 43 60247
388 315380 0 0 0 60720
389 315380 1 0 0 60720
390 315380 0 0 0 60720
391 315380 1 0 0 60720
392 315380 1 0 0 60720
393 357760 1 13253 947 90829
394 315637 1 24 8 59818
395 315380 1 0 0 60720
396 315380 0 0 0 60720
397 315380 1 0 0 60720
398 315380 1 0 0 60720
399 327071 1 131 6 59661
400 377516 1 514 15 102725
401 299243 1 3366 80 108479
402 387699 0 9327 101 87419
403 309836 1 384 22 54683
404 444477 1 17821 373 122844
405 322327 1 397 10 62710
406 314913 1 897 5 60720
407 315553 1 218 21 60720
408 688779 0 3369 28 87161
409 321896 1 5702 61 101481
410 301607 1 8636 855 128294
411 304485 1 534 7 62620
412 315380 1 0 0 60720
413 315380 0 0 0 60720
414 231861 1 726 3 69980
415 347385 1 1380 26 60982
416 316386 1 180 11 59635
417 491303 0 7285 115 188873
418 261216 1 880 74 80791
419 315388 1 1 2 60727
420 315380 1 0 0 60720
421 313729 1 96 5 60379
422 358649 1 1889 45 37527
423 1926517 0 45187 353 234817
424 296656 1 288 4 60510
425 275311 1 1270 26 69206
426 -42143 1 6526 26 55830
427 315380 0 0 0 60720
428 343929 1 226 12 72835
429 367655 1 694 8 68060
430 315380 1 0 0 60720
431 313491 1 249 72 56726
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) group costs trades dividends
205978.971 -55485.819 29.811 -407.150 1.926
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-2163938 -14344 35201 47931 3195284
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 2.060e+05 3.854e+04 5.345 1.48e-07 ***
group -5.549e+04 2.784e+04 -1.993 0.0469 *
costs 2.981e+01 1.578e+00 18.889 < 2e-16 ***
trades -4.072e+02 1.663e+02 -2.448 0.0148 *
dividends 1.926e+00 4.105e-01 4.692 3.65e-06 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 253400 on 426 degrees of freedom
Multiple R-squared: 0.6569, Adjusted R-squared: 0.6537
F-statistic: 203.9 on 4 and 426 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,] 8.484368e-04 1.696874e-03 9.991516e-01
[2,] 5.665972e-04 1.133194e-03 9.994334e-01
[3,] 7.190595e-05 1.438119e-04 9.999281e-01
[4,] 9.710950e-05 1.942190e-04 9.999029e-01
[5,] 4.832291e-05 9.664581e-05 9.999517e-01
[6,] 9.861619e-06 1.972324e-05 9.999901e-01
[7,] 1.611924e-05 3.223847e-05 9.999839e-01
[8,] 3.240916e-06 6.481832e-06 9.999968e-01
[9,] 6.401563e-07 1.280313e-06 9.999994e-01
[10,] 1.197364e-07 2.394728e-07 9.999999e-01
[11,] 2.181095e-08 4.362190e-08 1.000000e+00
[12,] 3.780405e-09 7.560810e-09 1.000000e+00
[13,] 6.273737e-10 1.254747e-09 1.000000e+00
[14,] 1.010434e-10 2.020868e-10 1.000000e+00
[15,] 1.615546e-11 3.231092e-11 1.000000e+00
[16,] 2.459365e-12 4.918730e-12 1.000000e+00
[17,] 4.223568e-13 8.447135e-13 1.000000e+00
[18,] 6.406709e-14 1.281342e-13 1.000000e+00
[19,] 9.077889e-15 1.815578e-14 1.000000e+00
[20,] 1.279272e-15 2.558545e-15 1.000000e+00
[21,] 1.855800e-16 3.711600e-16 1.000000e+00
[22,] 2.469671e-17 4.939342e-17 1.000000e+00
[23,] 3.209765e-18 6.419529e-18 1.000000e+00
[24,] 1.583470e-18 3.166941e-18 1.000000e+00
[25,] 2.086646e-19 4.173293e-19 1.000000e+00
[26,] 3.129255e-20 6.258509e-20 1.000000e+00
[27,] 7.759474e-21 1.551895e-20 1.000000e+00
[28,] 1.146319e-21 2.292638e-21 1.000000e+00
[29,] 1.806038e-22 3.612076e-22 1.000000e+00
[30,] 2.268877e-23 4.537755e-23 1.000000e+00
[31,] 3.258163e-24 6.516326e-24 1.000000e+00
[32,] 1.981634e-13 3.963269e-13 1.000000e+00
[33,] 6.034170e-14 1.206834e-13 1.000000e+00
[34,] 2.678091e-14 5.356181e-14 1.000000e+00
[35,] 1.362866e-14 2.725733e-14 1.000000e+00
[36,] 3.854280e-15 7.708560e-15 1.000000e+00
[37,] 1.057449e-15 2.114899e-15 1.000000e+00
[38,] 2.799644e-16 5.599288e-16 1.000000e+00
[39,] 7.883505e-17 1.576701e-16 1.000000e+00
[40,] 2.089755e-17 4.179509e-17 1.000000e+00
[41,] 5.343678e-18 1.068736e-17 1.000000e+00
[42,] 1.465547e-18 2.931094e-18 1.000000e+00
[43,] 4.083408e-19 8.166816e-19 1.000000e+00
[44,] 5.765247e-18 1.153049e-17 1.000000e+00
[45,] 1.622127e-18 3.244253e-18 1.000000e+00
[46,] 4.383994e-19 8.767988e-19 1.000000e+00
[47,] 1.250605e-19 2.501209e-19 1.000000e+00
[48,] 5.163040e-14 1.032608e-13 1.000000e+00
[49,] 2.802594e-14 5.605188e-14 1.000000e+00
[50,] 9.907984e-15 1.981597e-14 1.000000e+00
[51,] 4.019064e-15 8.038128e-15 1.000000e+00
[52,] 1.414009e-15 2.828018e-15 1.000000e+00
[53,] 4.893971e-16 9.787943e-16 1.000000e+00
[54,] 1.676641e-16 3.353283e-16 1.000000e+00
[55,] 1.897401e-16 3.794802e-16 1.000000e+00
[56,] 7.157465e-17 1.431493e-16 1.000000e+00
[57,] 2.439192e-17 4.878383e-17 1.000000e+00
[58,] 3.525940e-17 7.051880e-17 1.000000e+00
[59,] 1.211744e-17 2.423487e-17 1.000000e+00
[60,] 7.198683e-17 1.439737e-16 1.000000e+00
[61,] 3.299063e-17 6.598125e-17 1.000000e+00
[62,] 1.320711e-17 2.641422e-17 1.000000e+00
[63,] 4.637037e-18 9.274073e-18 1.000000e+00
[64,] 4.366882e-09 8.733765e-09 1.000000e+00
[65,] 1.001800e-08 2.003599e-08 1.000000e+00
[66,] 1.000000e+00 9.921018e-16 4.960509e-16
[67,] 1.000000e+00 1.895110e-15 9.475552e-16
[68,] 1.000000e+00 3.642158e-15 1.821079e-15
[69,] 1.000000e+00 6.757272e-15 3.378636e-15
[70,] 1.000000e+00 1.271375e-14 6.356875e-15
[71,] 1.000000e+00 2.362851e-14 1.181425e-14
[72,] 1.000000e+00 9.671149e-15 4.835575e-15
[73,] 1.000000e+00 8.345993e-16 4.172996e-16
[74,] 1.000000e+00 1.571996e-15 7.859979e-16
[75,] 1.000000e+00 2.937336e-15 1.468668e-15
[76,] 1.000000e+00 5.448409e-15 2.724204e-15
[77,] 1.000000e+00 1.008237e-14 5.041186e-15
[78,] 1.000000e+00 1.842538e-14 9.212692e-15
[79,] 1.000000e+00 3.396701e-14 1.698350e-14
[80,] 1.000000e+00 6.133035e-14 3.066518e-14
[81,] 1.000000e+00 1.096820e-13 5.484099e-14
[82,] 1.000000e+00 1.963150e-13 9.815749e-14
[83,] 1.000000e+00 3.462971e-13 1.731485e-13
[84,] 1.000000e+00 6.153513e-13 3.076757e-13
[85,] 1.000000e+00 1.077475e-12 5.387375e-13
[86,] 1.000000e+00 1.313099e-12 6.565497e-13
[87,] 1.000000e+00 1.836230e-12 9.181151e-13
[88,] 1.000000e+00 3.136701e-12 1.568351e-12
[89,] 1.000000e+00 5.324004e-12 2.662002e-12
[90,] 1.000000e+00 8.989909e-12 4.494954e-12
[91,] 1.000000e+00 1.514019e-11 7.570096e-12
[92,] 1.000000e+00 2.527376e-11 1.263688e-11
[93,] 1.000000e+00 4.088254e-11 2.044127e-11
[94,] 1.000000e+00 6.143946e-11 3.071973e-11
[95,] 1.000000e+00 1.015079e-10 5.075395e-11
[96,] 1.000000e+00 1.651878e-10 8.259389e-11
[97,] 1.000000e+00 2.669176e-10 1.334588e-10
[98,] 1.000000e+00 4.285435e-10 2.142717e-10
[99,] 1.000000e+00 1.176614e-11 5.883070e-12
[100,] 1.000000e+00 1.969451e-11 9.847255e-12
[101,] 1.000000e+00 3.119349e-11 1.559675e-11
[102,] 1.000000e+00 5.087096e-11 2.543548e-11
[103,] 1.000000e+00 8.247451e-11 4.123725e-11
[104,] 1.000000e+00 1.334495e-10 6.672476e-11
[105,] 1.000000e+00 2.006704e-10 1.003352e-10
[106,] 1.000000e+00 3.195677e-10 1.597839e-10
[107,] 1.000000e+00 5.123797e-10 2.561899e-10
[108,] 1.000000e+00 8.196487e-10 4.098244e-10
[109,] 1.000000e+00 1.287450e-09 6.437252e-10
[110,] 1.000000e+00 2.010781e-09 1.005391e-09
[111,] 1.000000e+00 2.762543e-09 1.381271e-09
[112,] 1.000000e+00 3.681272e-09 1.840636e-09
[113,] 1.000000e+00 5.603798e-09 2.801899e-09
[114,] 1.000000e+00 8.074612e-09 4.037306e-09
[115,] 1.000000e+00 1.225828e-08 6.129138e-09
[116,] 1.000000e+00 1.486480e-08 7.432400e-09
[117,] 1.000000e+00 2.263544e-08 1.131772e-08
[118,] 1.000000e+00 1.652429e-08 8.262144e-09
[119,] 1.000000e+00 2.516410e-08 1.258205e-08
[120,] 1.000000e+00 3.794774e-08 1.897387e-08
[121,] 1.000000e+00 5.632909e-08 2.816454e-08
[122,] 1.000000e+00 3.982383e-11 1.991191e-11
[123,] 1.000000e+00 6.242953e-11 3.121476e-11
[124,] 1.000000e+00 9.976004e-11 4.988002e-11
[125,] 1.000000e+00 1.591813e-10 7.959063e-11
[126,] 1.000000e+00 2.527101e-10 1.263551e-10
[127,] 1.000000e+00 3.991536e-10 1.995768e-10
[128,] 1.000000e+00 6.204676e-10 3.102338e-10
[129,] 1.000000e+00 9.594216e-10 4.797108e-10
[130,] 1.000000e+00 1.473770e-09 7.368852e-10
[131,] 1.000000e+00 2.260209e-09 1.130105e-09
[132,] 1.000000e+00 3.479437e-09 1.739719e-09
[133,] 1.000000e+00 5.264166e-09 2.632083e-09
[134,] 1.000000e+00 7.481755e-09 3.740877e-09
[135,] 1.000000e+00 9.381567e-09 4.690784e-09
[136,] 1.000000e+00 1.401844e-08 7.009222e-09
[137,] 1.000000e+00 2.084636e-08 1.042318e-08
[138,] 1.000000e+00 3.120861e-08 1.560431e-08
[139,] 1.000000e+00 4.654986e-08 2.327493e-08
[140,] 1.000000e+00 6.851278e-08 3.425639e-08
[141,] 1.000000e+00 9.736479e-08 4.868239e-08
[142,] 9.999999e-01 1.355086e-07 6.775431e-08
[143,] 9.999999e-01 1.958388e-07 9.791941e-08
[144,] 9.999999e-01 2.852404e-07 1.426202e-07
[145,] 9.999998e-01 4.133511e-07 2.066756e-07
[146,] 9.999997e-01 5.887539e-07 2.943769e-07
[147,] 9.999996e-01 8.439698e-07 4.219849e-07
[148,] 9.999994e-01 1.193783e-06 5.968914e-07
[149,] 9.999992e-01 1.576432e-06 7.882159e-07
[150,] 9.999989e-01 2.163169e-06 1.081585e-06
[151,] 9.999985e-01 3.012143e-06 1.506072e-06
[152,] 9.999979e-01 4.189007e-06 2.094504e-06
[153,] 9.999971e-01 5.770129e-06 2.885064e-06
[154,] 9.999960e-01 7.952795e-06 3.976398e-06
[155,] 9.999946e-01 1.087242e-05 5.436209e-06
[156,] 9.999926e-01 1.480999e-05 7.404995e-06
[157,] 9.999899e-01 2.027990e-05 1.013995e-05
[158,] 9.999863e-01 2.734521e-05 1.367261e-05
[159,] 9.999822e-01 3.559695e-05 1.779848e-05
[160,] 9.999762e-01 4.759168e-05 2.379584e-05
[161,] 9.999683e-01 6.331389e-05 3.165695e-05
[162,] 9.999578e-01 8.433054e-05 4.216527e-05
[163,] 9.999445e-01 1.109703e-04 5.548515e-05
[164,] 9.999285e-01 1.429787e-04 7.148937e-05
[165,] 9.999080e-01 1.840707e-04 9.203536e-05
[166,] 9.998970e-01 2.059166e-04 1.029583e-04
[167,] 9.998708e-01 2.584229e-04 1.292114e-04
[168,] 9.998323e-01 3.353787e-04 1.676894e-04
[169,] 9.997847e-01 4.306850e-04 2.153425e-04
[170,] 9.997235e-01 5.529107e-04 2.764554e-04
[171,] 9.996467e-01 7.066665e-04 3.533332e-04
[172,] 9.995475e-01 9.050856e-04 4.525428e-04
[173,] 9.994216e-01 1.156865e-03 5.784324e-04
[174,] 9.992685e-01 1.462979e-03 7.314896e-04
[175,] 9.990778e-01 1.844470e-03 9.222352e-04
[176,] 9.988463e-01 2.307397e-03 1.153699e-03
[177,] 9.985628e-01 2.874349e-03 1.437175e-03
[178,] 9.982033e-01 3.593447e-03 1.796723e-03
[179,] 9.977823e-01 4.435495e-03 2.217748e-03
[180,] 9.972682e-01 5.463658e-03 2.731829e-03
[181,] 9.966456e-01 6.708709e-03 3.354354e-03
[182,] 1.000000e+00 6.486876e-17 3.243438e-17
[183,] 1.000000e+00 1.145452e-16 5.727261e-17
[184,] 1.000000e+00 2.010375e-16 1.005187e-16
[185,] 1.000000e+00 3.590393e-16 1.795196e-16
[186,] 1.000000e+00 6.254204e-16 3.127102e-16
[187,] 1.000000e+00 1.071454e-15 5.357270e-16
[188,] 1.000000e+00 1.827067e-15 9.135334e-16
[189,] 1.000000e+00 3.140341e-15 1.570170e-15
[190,] 1.000000e+00 5.086919e-15 2.543460e-15
[191,] 1.000000e+00 5.017429e-15 2.508715e-15
[192,] 1.000000e+00 8.608374e-15 4.304187e-15
[193,] 1.000000e+00 1.461258e-14 7.306292e-15
[194,] 1.000000e+00 7.335428e-15 3.667714e-15
[195,] 1.000000e+00 1.258696e-14 6.293480e-15
[196,] 1.000000e+00 2.150346e-14 1.075173e-14
[197,] 1.000000e+00 3.663851e-14 1.831926e-14
[198,] 1.000000e+00 6.147306e-14 3.073653e-14
[199,] 1.000000e+00 9.945364e-14 4.972682e-14
[200,] 1.000000e+00 1.639913e-13 8.199563e-14
[201,] 1.000000e+00 2.716200e-13 1.358100e-13
[202,] 1.000000e+00 4.312368e-13 2.156184e-13
[203,] 1.000000e+00 7.215039e-13 3.607519e-13
[204,] 1.000000e+00 1.186368e-12 5.931840e-13
[205,] 1.000000e+00 1.220945e-12 6.104724e-13
[206,] 1.000000e+00 1.993336e-12 9.966680e-13
[207,] 1.000000e+00 3.268640e-12 1.634320e-12
[208,] 1.000000e+00 4.798555e-12 2.399278e-12
[209,] 1.000000e+00 6.016286e-12 3.008143e-12
[210,] 1.000000e+00 9.757872e-12 4.878936e-12
[211,] 1.000000e+00 1.560175e-11 7.800875e-12
[212,] 1.000000e+00 2.280856e-11 1.140428e-11
[213,] 1.000000e+00 3.648137e-11 1.824068e-11
[214,] 1.000000e+00 5.775923e-11 2.887961e-11
[215,] 1.000000e+00 9.083114e-11 4.541557e-11
[216,] 1.000000e+00 1.429135e-10 7.145675e-11
[217,] 1.000000e+00 2.260035e-10 1.130018e-10
[218,] 1.000000e+00 3.556400e-10 1.778200e-10
[219,] 1.000000e+00 5.220834e-11 2.610417e-11
[220,] 1.000000e+00 8.242160e-11 4.121080e-11
[221,] 1.000000e+00 1.002645e-10 5.013225e-11
[222,] 1.000000e+00 1.580578e-10 7.902889e-11
[223,] 1.000000e+00 2.498562e-10 1.249281e-10
[224,] 1.000000e+00 3.737215e-10 1.868607e-10
[225,] 1.000000e+00 5.816809e-10 2.908404e-10
[226,] 1.000000e+00 7.356029e-10 3.678015e-10
[227,] 1.000000e+00 1.132015e-09 5.660074e-10
[228,] 1.000000e+00 1.739411e-09 8.697056e-10
[229,] 1.000000e+00 2.694508e-09 1.347254e-09
[230,] 1.000000e+00 4.081513e-09 2.040756e-09
[231,] 1.000000e+00 6.169855e-09 3.084927e-09
[232,] 1.000000e+00 9.285235e-09 4.642618e-09
[233,] 1.000000e+00 1.389478e-08 6.947392e-09
[234,] 1.000000e+00 2.091167e-08 1.045583e-08
[235,] 1.000000e+00 3.148867e-08 1.574433e-08
[236,] 1.000000e+00 4.686588e-08 2.343294e-08
[237,] 1.000000e+00 6.936867e-08 3.468434e-08
[238,] 1.000000e+00 8.318332e-08 4.159166e-08
[239,] 9.999999e-01 1.143899e-07 5.719495e-08
[240,] 9.999999e-01 1.667603e-07 8.338016e-08
[241,] 9.999999e-01 2.410357e-07 1.205179e-07
[242,] 9.999998e-01 3.482596e-07 1.741298e-07
[243,] 9.999997e-01 5.018126e-07 2.509063e-07
[244,] 9.999996e-01 7.232469e-07 3.616235e-07
[245,] 9.999995e-01 1.035168e-06 5.175841e-07
[246,] 9.999993e-01 1.478584e-06 7.392921e-07
[247,] 9.999990e-01 2.088956e-06 1.044478e-06
[248,] 9.999985e-01 2.956987e-06 1.478493e-06
[249,] 9.999979e-01 4.140453e-06 2.070227e-06
[250,] 9.999971e-01 5.771556e-06 2.885778e-06
[251,] 9.999960e-01 8.009095e-06 4.004547e-06
[252,] 9.999946e-01 1.081374e-05 5.406868e-06
[253,] 9.999977e-01 4.615971e-06 2.307985e-06
[254,] 9.999982e-01 3.631730e-06 1.815865e-06
[255,] 9.999977e-01 4.605737e-06 2.302868e-06
[256,] 9.999969e-01 6.203360e-06 3.101680e-06
[257,] 9.999959e-01 8.165397e-06 4.082698e-06
[258,] 9.999948e-01 1.046507e-05 5.232537e-06
[259,] 9.999928e-01 1.447155e-05 7.235777e-06
[260,] 9.999899e-01 2.011157e-05 1.005578e-05
[261,] 9.999862e-01 2.755454e-05 1.377727e-05
[262,] 9.999812e-01 3.765978e-05 1.882989e-05
[263,] 9.999744e-01 5.118334e-05 2.559167e-05
[264,] 9.999654e-01 6.916072e-05 3.458036e-05
[265,] 9.999535e-01 9.302223e-05 4.651111e-05
[266,] 9.999371e-01 1.257906e-04 6.289529e-05
[267,] 9.999162e-01 1.675937e-04 8.379683e-05
[268,] 9.998942e-01 2.115774e-04 1.057887e-04
[269,] 9.998600e-01 2.800421e-04 1.400211e-04
[270,] 9.998149e-01 3.702998e-04 1.851499e-04
[271,] 9.997569e-01 4.861448e-04 2.430724e-04
[272,] 9.996796e-01 6.407712e-04 3.203856e-04
[273,] 9.995811e-01 8.377625e-04 4.188813e-04
[274,] 9.994546e-01 1.090767e-03 5.453835e-04
[275,] 9.992996e-01 1.400709e-03 7.003545e-04
[276,] 9.991039e-01 1.792185e-03 8.960926e-04
[277,] 9.988486e-01 2.302798e-03 1.151399e-03
[278,] 9.985426e-01 2.914834e-03 1.457417e-03
[279,] 9.981595e-01 3.681073e-03 1.840536e-03
[280,] 9.976886e-01 4.622843e-03 2.311422e-03
[281,] 9.972482e-01 5.503583e-03 2.751792e-03
[282,] 9.978107e-01 4.378580e-03 2.189290e-03
[283,] 9.996868e-01 6.264513e-04 3.132256e-04
[284,] 9.995868e-01 8.264326e-04 4.132163e-04
[285,] 9.994602e-01 1.079672e-03 5.398359e-04
[286,] 9.993029e-01 1.394196e-03 6.970981e-04
[287,] 9.990934e-01 1.813229e-03 9.066147e-04
[288,] 9.988408e-01 2.318448e-03 1.159224e-03
[289,] 9.985249e-01 2.950238e-03 1.475119e-03
[290,] 9.981314e-01 3.737125e-03 1.868562e-03
[291,] 9.998542e-01 2.915175e-04 1.457587e-04
[292,] 9.998048e-01 3.904402e-04 1.952201e-04
[293,] 9.997399e-01 5.201471e-04 2.600735e-04
[294,] 9.998766e-01 2.468940e-04 1.234470e-04
[295,] 9.998331e-01 3.338266e-04 1.669133e-04
[296,] 9.997783e-01 4.433009e-04 2.216504e-04
[297,] 9.997064e-01 5.871813e-04 2.935906e-04
[298,] 9.996113e-01 7.774578e-04 3.887289e-04
[299,] 9.994920e-01 1.016099e-03 5.080495e-04
[300,] 9.994095e-01 1.181021e-03 5.905106e-04
[301,] 9.992334e-01 1.533102e-03 7.665509e-04
[302,] 9.990100e-01 1.980046e-03 9.900230e-04
[303,] 9.987295e-01 2.541049e-03 1.270525e-03
[304,] 9.983548e-01 3.290451e-03 1.645225e-03
[305,] 9.979449e-01 4.110107e-03 2.055053e-03
[306,] 9.973663e-01 5.267438e-03 2.633719e-03
[307,] 9.966800e-01 6.639974e-03 3.319987e-03
[308,] 9.958227e-01 8.354700e-03 4.177350e-03
[309,] 9.954982e-01 9.003552e-03 4.501776e-03
[310,] 9.999903e-01 1.934470e-05 9.672349e-06
[311,] 9.999860e-01 2.800932e-05 1.400466e-05
[312,] 9.999799e-01 4.010688e-05 2.005344e-05
[313,] 9.999718e-01 5.641322e-05 2.820661e-05
[314,] 9.999607e-01 7.858445e-05 3.929222e-05
[315,] 9.999453e-01 1.093790e-04 5.468950e-05
[316,] 9.999227e-01 1.545355e-04 7.726775e-05
[317,] 9.998938e-01 2.123160e-04 1.061580e-04
[318,] 9.998522e-01 2.956329e-04 1.478164e-04
[319,] 9.997992e-01 4.015649e-04 2.007825e-04
[320,] 9.997250e-01 5.499780e-04 2.749890e-04
[321,] 1.000000e+00 9.631258e-13 4.815629e-13
[322,] 1.000000e+00 1.299596e-20 6.497978e-21
[323,] 1.000000e+00 1.726870e-20 8.634350e-21
[324,] 1.000000e+00 3.829884e-20 1.914942e-20
[325,] 1.000000e+00 9.285787e-20 4.642894e-20
[326,] 1.000000e+00 2.312863e-19 1.156431e-19
[327,] 1.000000e+00 2.987356e-19 1.493678e-19
[328,] 1.000000e+00 7.117202e-19 3.558601e-19
[329,] 1.000000e+00 1.697864e-18 8.489321e-19
[330,] 1.000000e+00 3.749497e-18 1.874748e-18
[331,] 1.000000e+00 2.041115e-19 1.020557e-19
[332,] 1.000000e+00 5.034865e-19 2.517433e-19
[333,] 1.000000e+00 1.269140e-18 6.345702e-19
[334,] 1.000000e+00 4.378779e-19 2.189390e-19
[335,] 1.000000e+00 1.123281e-18 5.616406e-19
[336,] 1.000000e+00 2.783901e-18 1.391951e-18
[337,] 1.000000e+00 5.795603e-18 2.897802e-18
[338,] 1.000000e+00 1.417015e-17 7.085075e-18
[339,] 1.000000e+00 3.297251e-17 1.648626e-17
[340,] 1.000000e+00 7.999336e-17 3.999668e-17
[341,] 1.000000e+00 1.736588e-18 8.682939e-19
[342,] 1.000000e+00 4.473357e-18 2.236679e-18
[343,] 1.000000e+00 1.157456e-17 5.787279e-18
[344,] 1.000000e+00 2.935581e-17 1.467791e-17
[345,] 1.000000e+00 7.386488e-17 3.693244e-17
[346,] 1.000000e+00 1.879568e-16 9.397839e-17
[347,] 1.000000e+00 4.655715e-16 2.327857e-16
[348,] 1.000000e+00 1.118305e-15 5.591527e-16
[349,] 1.000000e+00 2.774310e-15 1.387155e-15
[350,] 1.000000e+00 6.722435e-15 3.361217e-15
[351,] 1.000000e+00 1.639046e-14 8.195231e-15
[352,] 1.000000e+00 3.906479e-14 1.953239e-14
[353,] 1.000000e+00 5.036680e-14 2.518340e-14
[354,] 1.000000e+00 4.458427e-14 2.229214e-14
[355,] 1.000000e+00 1.088889e-13 5.444446e-14
[356,] 1.000000e+00 2.622499e-13 1.311249e-13
[357,] 1.000000e+00 6.275484e-13 3.137742e-13
[358,] 1.000000e+00 1.453913e-12 7.269567e-13
[359,] 1.000000e+00 5.679096e-14 2.839548e-14
[360,] 1.000000e+00 1.118570e-13 5.592850e-14
[361,] 1.000000e+00 2.688117e-13 1.344058e-13
[362,] 1.000000e+00 6.425918e-13 3.212959e-13
[363,] 1.000000e+00 3.013306e-13 1.506653e-13
[364,] 1.000000e+00 7.514672e-13 3.757336e-13
[365,] 1.000000e+00 1.890680e-12 9.453402e-13
[366,] 1.000000e+00 4.330187e-12 2.165093e-12
[367,] 1.000000e+00 4.670571e-12 2.335285e-12
[368,] 1.000000e+00 1.171694e-11 5.858472e-12
[369,] 1.000000e+00 2.832921e-11 1.416460e-11
[370,] 1.000000e+00 6.706326e-11 3.353163e-11
[371,] 1.000000e+00 1.551195e-10 7.755977e-11
[372,] 1.000000e+00 3.685980e-10 1.842990e-10
[373,] 1.000000e+00 8.404425e-10 4.202213e-10
[374,] 1.000000e+00 1.972418e-09 9.862089e-10
[375,] 1.000000e+00 4.590103e-09 2.295052e-09
[376,] 1.000000e+00 5.883611e-09 2.941806e-09
[377,] 1.000000e+00 1.291454e-08 6.457270e-09
[378,] 1.000000e+00 2.823164e-08 1.411582e-08
[379,] 1.000000e+00 6.367452e-08 3.183726e-08
[380,] 9.999999e-01 1.386605e-07 6.933026e-08
[381,] 9.999998e-01 3.009328e-07 1.504664e-07
[382,] 9.999997e-01 6.274968e-07 3.137484e-07
[383,] 9.999993e-01 1.321403e-06 6.607016e-07
[384,] 9.999987e-01 2.684638e-06 1.342319e-06
[385,] 9.999973e-01 5.375235e-06 2.687617e-06
[386,] 9.999948e-01 1.038571e-05 5.192855e-06
[387,] 9.999899e-01 2.011283e-05 1.005642e-05
[388,] 9.999807e-01 3.866451e-05 1.933225e-05
[389,] 9.999625e-01 7.508213e-05 3.754107e-05
[390,] 9.999299e-01 1.401252e-04 7.006262e-05
[391,] 9.998714e-01 2.571669e-04 1.285835e-04
[392,] 9.997735e-01 4.530684e-04 2.265342e-04
[393,] 9.996045e-01 7.909579e-04 3.954789e-04
[394,] 9.993499e-01 1.300221e-03 6.501103e-04
[395,] 9.995763e-01 8.473422e-04 4.236711e-04
[396,] 9.992353e-01 1.529455e-03 7.647274e-04
[397,] 9.996841e-01 6.318893e-04 3.159447e-04
[398,] 9.994032e-01 1.193665e-03 5.968325e-04
[399,] 9.988576e-01 2.284762e-03 1.142381e-03
[400,] 9.979459e-01 4.108281e-03 2.054141e-03
[401,] 9.983773e-01 3.245341e-03 1.622670e-03
[402,] 9.978334e-01 4.333206e-03 2.166603e-03
[403,] 9.969624e-01 6.075206e-03 3.037603e-03
[404,] 9.943040e-01 1.139194e-02 5.695970e-03
[405,] 9.902485e-01 1.950303e-02 9.751513e-03
[406,] 9.824527e-01 3.509458e-02 1.754729e-02
[407,] 9.695200e-01 6.096005e-02 3.048003e-02
[408,] 9.491105e-01 1.017791e-01 5.088953e-02
[409,] 9.195370e-01 1.609259e-01 8.046297e-02
[410,] 9.640072e-01 7.198563e-02 3.599282e-02
[411,] 9.844260e-01 3.114798e-02 1.557399e-02
[412,] 9.682975e-01 6.340502e-02 3.170251e-02
[413,] 9.400735e-01 1.198531e-01 5.992655e-02
[414,] 8.891404e-01 2.217192e-01 1.108596e-01
[415,] 9.810401e-01 3.791974e-02 1.895987e-02
[416,] 9.967602e-01 6.479602e-03 3.239801e-03
> postscript(file="/var/www/rcomp/tmp/1kz5g1293276484.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/2c84j1293276484.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/3c84j1293276484.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/4c84j1293276484.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/5n0m41293276484.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> qqnorm(mysum$resid, main='Residual Normal Q-Q Plot')
> qqline(mysum$resid)
> grid()
> dev.off()
null device
1
> (myerror <- as.ts(mysum$resid))
Time Series:
Start = 1
End = 431
Frequency = 1
1 2 3 4 5 6
11292.552 -41076.166 -19198.651 79985.494 -280953.361 -62387.577
7 8 9 10 11 12
17215.934 49320.735 63580.147 38208.833 -41888.465 -68881.815
13 14 15 16 17 18
47930.758 95359.943 40017.496 16128.433 51727.086 60651.036
19 20 21 22 23 24
52885.045 48104.426 47930.758 59484.659 47930.758 4172.197
25 26 27 28 29 30
23515.058 40932.200 -7827.369 -131141.638 47930.758 17147.822
31 32 33 34 35 36
-250216.313 47930.758 65051.344 71559.952 -190418.162 -148255.495
37 38 39 40 41 42
49490.851 -41774.594 -25648.977 -126404.718 56121.615 129383.608
43 44 45 46 47 48
-32278.740 47930.758 -37101.259 -11465.171 47930.758 37413.394
49 50 51 52 53 54
-124535.591 24236.170 -402357.714 51398.186 -74411.743 -133413.771
55 56 57 58 59 60
498000.209 -306359.379 40596.813 -139539.338 -112771.362 47930.758
61 62 63 64 65 66
47930.758 -249808.092 -23000.886 47847.411 -266560.840 47253.085
67 68 69 70 71 72
-298183.103 -25660.802 -33992.012 47930.758 704705.338 -306374.089
73 74 75 76 77 78
3195283.878 47930.758 -35584.633 73611.546 17805.179 49994.478
79 80 81 82 83 84
-282147.420 -465658.107 48103.003 47930.758 47930.758 -37960.116
85 86 87 88 89 90
47930.758 -7555.060 46990.954 48672.001 37073.804 35414.901
91 92 93 94 95 96
-7555.060 41434.050 -273562.970 -150276.044 50535.241 44777.497
97 98 99 100 101 102
47930.758 48453.678 53038.424 65563.601 152629.298 29601.307
103 104 105 106 107 108
50256.107 48181.113 47930.758 -675705.948 19719.600 -86633.872
109 110 111 112 113 114
55446.705 54791.004 43429.078 -115267.526 47930.758 23566.412
115 116 117 118 119 120
-7555.060 45851.206 -53093.255 153902.340 148737.922 66254.339
121 122 123 124 125 126
-100883.026 47930.758 190360.572 20323.606 -344582.470 -7555.060
127 128 129 130 131 132
33983.667 47930.758 -1118357.214 -65872.062 34349.869 -7555.060
133 134 135 136 137 138
-7555.060 -7555.060 46456.029 -11265.356 47930.758 49743.121
139 140 141 142 143 144
-7555.060 47930.758 127681.941 -123725.963 47930.758 47930.758
145 146 147 148 149 150
-2370.308 -7555.060 41821.079 90994.114 -35871.190 47930.758
151 152 153 154 155 156
-7555.060 17459.373 47930.758 21375.207 -49407.953 -124797.099
157 158 159 160 161 162
85548.186 47930.758 40959.878 29543.810 42217.914 47930.758
163 164 165 166 167 168
48756.725 -5631.489 47930.758 80331.658 47930.758 50336.225
169 170 171 172 173 174
20183.871 37562.856 -65081.319 63483.018 -200546.755 -109444.008
175 176 177 178 179 180
23336.527 54252.289 47930.758 50983.381 29871.149 -7555.060
181 182 183 184 185 186
-28126.641 -6906.064 47930.758 47930.758 -7555.060 51520.595
187 188 189 190 191 192
43892.484 39825.868 -2163938.266 47930.758 49666.953 10638.634
193 194 195 196 197 198
47930.758 59511.531 55714.984 47930.758 -83722.349 -223309.564
199 200 201 202 203 204
-7555.060 47930.758 -279201.028 -7555.060 -7555.060 35201.158
205 206 207 208 209 210
47930.758 63553.861 -58373.820 47930.758 51085.001 6174.204
211 212 213 214 215 216
-12297.139 -223683.527 44403.489 34065.162 -122718.005 -157816.795
217 218 219 220 221 222
-8066.752 47930.758 27844.604 -29593.037 45014.662 47930.758
223 224 225 226 227 228
34892.012 20843.505 25176.204 -459875.779 47930.758 -234688.381
229 230 231 232 233 234
24439.537 27297.538 -144589.537 42833.808 -153760.016 47930.758
235 236 237 238 239 240
43760.928 10402.487 55479.808 47930.758 47930.758 48830.870
241 242 243 244 245 246
12120.541 10492.987 -7555.060 39310.301 -138046.328 118575.766
247 248 249 250 251 252
47930.758 53714.507 47930.758 45207.361 -7555.060 25991.297
253 254 255 256 257 258
-7555.060 47930.758 -7555.060 47930.758 47930.758 47930.758
259 260 261 262 263 264
48357.958 -374596.864 297382.198 -89503.531 -95882.735 -43408.390
265 266 267 268 269 270
-107098.631 47930.758 24781.823 47930.758 39802.780 45723.816
271 272 273 274 275 276
47930.758 47930.758 10354.181 48001.084 105736.955 2068.080
277 278 279 280 281 282
-8516.655 41418.222 9037.666 -20308.732 27433.327 -59079.502
283 284 285 286 287 288
47930.758 19310.917 50337.774 47076.205 47930.758 111838.134
289 290 291 292 293 294
333879.073 -576390.302 -7555.060 35997.643 47930.758 -12187.230
295 296 297 298 299 300
47968.115 48340.985 47930.758 762830.868 -7555.060 -7555.060
301 302 303 304 305 306
-278718.282 -14179.696 -74418.783 42624.223 32019.062 48377.897
307 308 309 310 311 312
-158893.721 47930.758 47930.758 49543.695 -7555.060 -87387.956
313 314 315 316 317 318
-13574.240 47890.901 35089.843 -189194.453 1187360.002 16757.609
319 320 321 322 323 324
-7555.060 -60797.317 50834.989 47803.348 -6488.100 47930.758
325 326 327 328 329 330
-14508.840 47930.758 -67415.436 1259827.443 -1184133.917 -180351.327
331 332 333 334 335 336
-62936.025 47930.758 -7555.060 -63103.461 -69214.165 47930.758
337 338 339 340 341 342
22362.228 -378580.851 47930.758 -7555.060 -265684.483 13127.108
343 344 345 346 347 348
47930.758 60434.071 47930.758 -96498.844 44196.429 -313458.634
349 350 351 352 353 354
47930.758 41837.252 47930.758 47930.758 -9187.034 47930.758
355 356 357 358 359 360
-31544.333 -7555.060 47930.758 -7555.060 47930.758 -9824.282
361 362 363 364 365 366
-348387.682 58126.423 15002.087 -56691.365 55492.576 408709.063
367 368 369 370 371 372
95780.130 49921.272 47930.758 -232554.023 -65383.903 26690.904
373 374 375 376 377 378
58494.152 -159079.702 -10313.459 4473.442 -62195.598 5446.994
379 380 381 382 383 384
28918.905 47930.758 -8951.678 9957.267 -152906.495 49670.808
385 386 387 388 389 390
48197.857 15549.014 -18948.548 -7555.060 47930.758 -7555.060
391 392 393 394 395 396
47930.758 47930.758 22800.421 52466.884 47930.758 -7555.060
397 398 399 400 401 402
47930.758 47930.758 60199.198 19942.976 -127969.710 -223588.738
403 404 405 406 407 408
51524.780 -322029.866 43281.192 22758.922 50155.084 225880.987
409 410 411 412 413 414
-169212.146 -5335.757 20306.970 47930.758 -7555.060 -73845.864
415 416 417 418 419 420
48877.645 50139.283 -248826.358 -40997.771 48709.764 47930.758
421 422 423 424 425 426
46110.459 97881.580 64892.532 22654.245 -35757.824 -484134.906
427 428 429 430 431
-7555.060 51292.884 68636.060 47930.758 75626.658
> postscript(file="/var/www/rcomp/tmp/6n0m41293276484.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 11292.552 NA
1 -41076.166 11292.552
2 -19198.651 -41076.166
3 79985.494 -19198.651
4 -280953.361 79985.494
5 -62387.577 -280953.361
6 17215.934 -62387.577
7 49320.735 17215.934
8 63580.147 49320.735
9 38208.833 63580.147
10 -41888.465 38208.833
11 -68881.815 -41888.465
12 47930.758 -68881.815
13 95359.943 47930.758
14 40017.496 95359.943
15 16128.433 40017.496
16 51727.086 16128.433
17 60651.036 51727.086
18 52885.045 60651.036
19 48104.426 52885.045
20 47930.758 48104.426
21 59484.659 47930.758
22 47930.758 59484.659
23 4172.197 47930.758
24 23515.058 4172.197
25 40932.200 23515.058
26 -7827.369 40932.200
27 -131141.638 -7827.369
28 47930.758 -131141.638
29 17147.822 47930.758
30 -250216.313 17147.822
31 47930.758 -250216.313
32 65051.344 47930.758
33 71559.952 65051.344
34 -190418.162 71559.952
35 -148255.495 -190418.162
36 49490.851 -148255.495
37 -41774.594 49490.851
38 -25648.977 -41774.594
39 -126404.718 -25648.977
40 56121.615 -126404.718
41 129383.608 56121.615
42 -32278.740 129383.608
43 47930.758 -32278.740
44 -37101.259 47930.758
45 -11465.171 -37101.259
46 47930.758 -11465.171
47 37413.394 47930.758
48 -124535.591 37413.394
49 24236.170 -124535.591
50 -402357.714 24236.170
51 51398.186 -402357.714
52 -74411.743 51398.186
53 -133413.771 -74411.743
54 498000.209 -133413.771
55 -306359.379 498000.209
56 40596.813 -306359.379
57 -139539.338 40596.813
58 -112771.362 -139539.338
59 47930.758 -112771.362
60 47930.758 47930.758
61 -249808.092 47930.758
62 -23000.886 -249808.092
63 47847.411 -23000.886
64 -266560.840 47847.411
65 47253.085 -266560.840
66 -298183.103 47253.085
67 -25660.802 -298183.103
68 -33992.012 -25660.802
69 47930.758 -33992.012
70 704705.338 47930.758
71 -306374.089 704705.338
72 3195283.878 -306374.089
73 47930.758 3195283.878
74 -35584.633 47930.758
75 73611.546 -35584.633
76 17805.179 73611.546
77 49994.478 17805.179
78 -282147.420 49994.478
79 -465658.107 -282147.420
80 48103.003 -465658.107
81 47930.758 48103.003
82 47930.758 47930.758
83 -37960.116 47930.758
84 47930.758 -37960.116
85 -7555.060 47930.758
86 46990.954 -7555.060
87 48672.001 46990.954
88 37073.804 48672.001
89 35414.901 37073.804
90 -7555.060 35414.901
91 41434.050 -7555.060
92 -273562.970 41434.050
93 -150276.044 -273562.970
94 50535.241 -150276.044
95 44777.497 50535.241
96 47930.758 44777.497
97 48453.678 47930.758
98 53038.424 48453.678
99 65563.601 53038.424
100 152629.298 65563.601
101 29601.307 152629.298
102 50256.107 29601.307
103 48181.113 50256.107
104 47930.758 48181.113
105 -675705.948 47930.758
106 19719.600 -675705.948
107 -86633.872 19719.600
108 55446.705 -86633.872
109 54791.004 55446.705
110 43429.078 54791.004
111 -115267.526 43429.078
112 47930.758 -115267.526
113 23566.412 47930.758
114 -7555.060 23566.412
115 45851.206 -7555.060
116 -53093.255 45851.206
117 153902.340 -53093.255
118 148737.922 153902.340
119 66254.339 148737.922
120 -100883.026 66254.339
121 47930.758 -100883.026
122 190360.572 47930.758
123 20323.606 190360.572
124 -344582.470 20323.606
125 -7555.060 -344582.470
126 33983.667 -7555.060
127 47930.758 33983.667
128 -1118357.214 47930.758
129 -65872.062 -1118357.214
130 34349.869 -65872.062
131 -7555.060 34349.869
132 -7555.060 -7555.060
133 -7555.060 -7555.060
134 46456.029 -7555.060
135 -11265.356 46456.029
136 47930.758 -11265.356
137 49743.121 47930.758
138 -7555.060 49743.121
139 47930.758 -7555.060
140 127681.941 47930.758
141 -123725.963 127681.941
142 47930.758 -123725.963
143 47930.758 47930.758
144 -2370.308 47930.758
145 -7555.060 -2370.308
146 41821.079 -7555.060
147 90994.114 41821.079
148 -35871.190 90994.114
149 47930.758 -35871.190
150 -7555.060 47930.758
151 17459.373 -7555.060
152 47930.758 17459.373
153 21375.207 47930.758
154 -49407.953 21375.207
155 -124797.099 -49407.953
156 85548.186 -124797.099
157 47930.758 85548.186
158 40959.878 47930.758
159 29543.810 40959.878
160 42217.914 29543.810
161 47930.758 42217.914
162 48756.725 47930.758
163 -5631.489 48756.725
164 47930.758 -5631.489
165 80331.658 47930.758
166 47930.758 80331.658
167 50336.225 47930.758
168 20183.871 50336.225
169 37562.856 20183.871
170 -65081.319 37562.856
171 63483.018 -65081.319
172 -200546.755 63483.018
173 -109444.008 -200546.755
174 23336.527 -109444.008
175 54252.289 23336.527
176 47930.758 54252.289
177 50983.381 47930.758
178 29871.149 50983.381
179 -7555.060 29871.149
180 -28126.641 -7555.060
181 -6906.064 -28126.641
182 47930.758 -6906.064
183 47930.758 47930.758
184 -7555.060 47930.758
185 51520.595 -7555.060
186 43892.484 51520.595
187 39825.868 43892.484
188 -2163938.266 39825.868
189 47930.758 -2163938.266
190 49666.953 47930.758
191 10638.634 49666.953
192 47930.758 10638.634
193 59511.531 47930.758
194 55714.984 59511.531
195 47930.758 55714.984
196 -83722.349 47930.758
197 -223309.564 -83722.349
198 -7555.060 -223309.564
199 47930.758 -7555.060
200 -279201.028 47930.758
201 -7555.060 -279201.028
202 -7555.060 -7555.060
203 35201.158 -7555.060
204 47930.758 35201.158
205 63553.861 47930.758
206 -58373.820 63553.861
207 47930.758 -58373.820
208 51085.001 47930.758
209 6174.204 51085.001
210 -12297.139 6174.204
211 -223683.527 -12297.139
212 44403.489 -223683.527
213 34065.162 44403.489
214 -122718.005 34065.162
215 -157816.795 -122718.005
216 -8066.752 -157816.795
217 47930.758 -8066.752
218 27844.604 47930.758
219 -29593.037 27844.604
220 45014.662 -29593.037
221 47930.758 45014.662
222 34892.012 47930.758
223 20843.505 34892.012
224 25176.204 20843.505
225 -459875.779 25176.204
226 47930.758 -459875.779
227 -234688.381 47930.758
228 24439.537 -234688.381
229 27297.538 24439.537
230 -144589.537 27297.538
231 42833.808 -144589.537
232 -153760.016 42833.808
233 47930.758 -153760.016
234 43760.928 47930.758
235 10402.487 43760.928
236 55479.808 10402.487
237 47930.758 55479.808
238 47930.758 47930.758
239 48830.870 47930.758
240 12120.541 48830.870
241 10492.987 12120.541
242 -7555.060 10492.987
243 39310.301 -7555.060
244 -138046.328 39310.301
245 118575.766 -138046.328
246 47930.758 118575.766
247 53714.507 47930.758
248 47930.758 53714.507
249 45207.361 47930.758
250 -7555.060 45207.361
251 25991.297 -7555.060
252 -7555.060 25991.297
253 47930.758 -7555.060
254 -7555.060 47930.758
255 47930.758 -7555.060
256 47930.758 47930.758
257 47930.758 47930.758
258 48357.958 47930.758
259 -374596.864 48357.958
260 297382.198 -374596.864
261 -89503.531 297382.198
262 -95882.735 -89503.531
263 -43408.390 -95882.735
264 -107098.631 -43408.390
265 47930.758 -107098.631
266 24781.823 47930.758
267 47930.758 24781.823
268 39802.780 47930.758
269 45723.816 39802.780
270 47930.758 45723.816
271 47930.758 47930.758
272 10354.181 47930.758
273 48001.084 10354.181
274 105736.955 48001.084
275 2068.080 105736.955
276 -8516.655 2068.080
277 41418.222 -8516.655
278 9037.666 41418.222
279 -20308.732 9037.666
280 27433.327 -20308.732
281 -59079.502 27433.327
282 47930.758 -59079.502
283 19310.917 47930.758
284 50337.774 19310.917
285 47076.205 50337.774
286 47930.758 47076.205
287 111838.134 47930.758
288 333879.073 111838.134
289 -576390.302 333879.073
290 -7555.060 -576390.302
291 35997.643 -7555.060
292 47930.758 35997.643
293 -12187.230 47930.758
294 47968.115 -12187.230
295 48340.985 47968.115
296 47930.758 48340.985
297 762830.868 47930.758
298 -7555.060 762830.868
299 -7555.060 -7555.060
300 -278718.282 -7555.060
301 -14179.696 -278718.282
302 -74418.783 -14179.696
303 42624.223 -74418.783
304 32019.062 42624.223
305 48377.897 32019.062
306 -158893.721 48377.897
307 47930.758 -158893.721
308 47930.758 47930.758
309 49543.695 47930.758
310 -7555.060 49543.695
311 -87387.956 -7555.060
312 -13574.240 -87387.956
313 47890.901 -13574.240
314 35089.843 47890.901
315 -189194.453 35089.843
316 1187360.002 -189194.453
317 16757.609 1187360.002
318 -7555.060 16757.609
319 -60797.317 -7555.060
320 50834.989 -60797.317
321 47803.348 50834.989
322 -6488.100 47803.348
323 47930.758 -6488.100
324 -14508.840 47930.758
325 47930.758 -14508.840
326 -67415.436 47930.758
327 1259827.443 -67415.436
328 -1184133.917 1259827.443
329 -180351.327 -1184133.917
330 -62936.025 -180351.327
331 47930.758 -62936.025
332 -7555.060 47930.758
333 -63103.461 -7555.060
334 -69214.165 -63103.461
335 47930.758 -69214.165
336 22362.228 47930.758
337 -378580.851 22362.228
338 47930.758 -378580.851
339 -7555.060 47930.758
340 -265684.483 -7555.060
341 13127.108 -265684.483
342 47930.758 13127.108
343 60434.071 47930.758
344 47930.758 60434.071
345 -96498.844 47930.758
346 44196.429 -96498.844
347 -313458.634 44196.429
348 47930.758 -313458.634
349 41837.252 47930.758
350 47930.758 41837.252
351 47930.758 47930.758
352 -9187.034 47930.758
353 47930.758 -9187.034
354 -31544.333 47930.758
355 -7555.060 -31544.333
356 47930.758 -7555.060
357 -7555.060 47930.758
358 47930.758 -7555.060
359 -9824.282 47930.758
360 -348387.682 -9824.282
361 58126.423 -348387.682
362 15002.087 58126.423
363 -56691.365 15002.087
364 55492.576 -56691.365
365 408709.063 55492.576
366 95780.130 408709.063
367 49921.272 95780.130
368 47930.758 49921.272
369 -232554.023 47930.758
370 -65383.903 -232554.023
371 26690.904 -65383.903
372 58494.152 26690.904
373 -159079.702 58494.152
374 -10313.459 -159079.702
375 4473.442 -10313.459
376 -62195.598 4473.442
377 5446.994 -62195.598
378 28918.905 5446.994
379 47930.758 28918.905
380 -8951.678 47930.758
381 9957.267 -8951.678
382 -152906.495 9957.267
383 49670.808 -152906.495
384 48197.857 49670.808
385 15549.014 48197.857
386 -18948.548 15549.014
387 -7555.060 -18948.548
388 47930.758 -7555.060
389 -7555.060 47930.758
390 47930.758 -7555.060
391 47930.758 47930.758
392 22800.421 47930.758
393 52466.884 22800.421
394 47930.758 52466.884
395 -7555.060 47930.758
396 47930.758 -7555.060
397 47930.758 47930.758
398 60199.198 47930.758
399 19942.976 60199.198
400 -127969.710 19942.976
401 -223588.738 -127969.710
402 51524.780 -223588.738
403 -322029.866 51524.780
404 43281.192 -322029.866
405 22758.922 43281.192
406 50155.084 22758.922
407 225880.987 50155.084
408 -169212.146 225880.987
409 -5335.757 -169212.146
410 20306.970 -5335.757
411 47930.758 20306.970
412 -7555.060 47930.758
413 -73845.864 -7555.060
414 48877.645 -73845.864
415 50139.283 48877.645
416 -248826.358 50139.283
417 -40997.771 -248826.358
418 48709.764 -40997.771
419 47930.758 48709.764
420 46110.459 47930.758
421 97881.580 46110.459
422 64892.532 97881.580
423 22654.245 64892.532
424 -35757.824 22654.245
425 -484134.906 -35757.824
426 -7555.060 -484134.906
427 51292.884 -7555.060
428 68636.060 51292.884
429 47930.758 68636.060
430 75626.658 47930.758
431 NA 75626.658
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] -41076.166 11292.552
[2,] -19198.651 -41076.166
[3,] 79985.494 -19198.651
[4,] -280953.361 79985.494
[5,] -62387.577 -280953.361
[6,] 17215.934 -62387.577
[7,] 49320.735 17215.934
[8,] 63580.147 49320.735
[9,] 38208.833 63580.147
[10,] -41888.465 38208.833
[11,] -68881.815 -41888.465
[12,] 47930.758 -68881.815
[13,] 95359.943 47930.758
[14,] 40017.496 95359.943
[15,] 16128.433 40017.496
[16,] 51727.086 16128.433
[17,] 60651.036 51727.086
[18,] 52885.045 60651.036
[19,] 48104.426 52885.045
[20,] 47930.758 48104.426
[21,] 59484.659 47930.758
[22,] 47930.758 59484.659
[23,] 4172.197 47930.758
[24,] 23515.058 4172.197
[25,] 40932.200 23515.058
[26,] -7827.369 40932.200
[27,] -131141.638 -7827.369
[28,] 47930.758 -131141.638
[29,] 17147.822 47930.758
[30,] -250216.313 17147.822
[31,] 47930.758 -250216.313
[32,] 65051.344 47930.758
[33,] 71559.952 65051.344
[34,] -190418.162 71559.952
[35,] -148255.495 -190418.162
[36,] 49490.851 -148255.495
[37,] -41774.594 49490.851
[38,] -25648.977 -41774.594
[39,] -126404.718 -25648.977
[40,] 56121.615 -126404.718
[41,] 129383.608 56121.615
[42,] -32278.740 129383.608
[43,] 47930.758 -32278.740
[44,] -37101.259 47930.758
[45,] -11465.171 -37101.259
[46,] 47930.758 -11465.171
[47,] 37413.394 47930.758
[48,] -124535.591 37413.394
[49,] 24236.170 -124535.591
[50,] -402357.714 24236.170
[51,] 51398.186 -402357.714
[52,] -74411.743 51398.186
[53,] -133413.771 -74411.743
[54,] 498000.209 -133413.771
[55,] -306359.379 498000.209
[56,] 40596.813 -306359.379
[57,] -139539.338 40596.813
[58,] -112771.362 -139539.338
[59,] 47930.758 -112771.362
[60,] 47930.758 47930.758
[61,] -249808.092 47930.758
[62,] -23000.886 -249808.092
[63,] 47847.411 -23000.886
[64,] -266560.840 47847.411
[65,] 47253.085 -266560.840
[66,] -298183.103 47253.085
[67,] -25660.802 -298183.103
[68,] -33992.012 -25660.802
[69,] 47930.758 -33992.012
[70,] 704705.338 47930.758
[71,] -306374.089 704705.338
[72,] 3195283.878 -306374.089
[73,] 47930.758 3195283.878
[74,] -35584.633 47930.758
[75,] 73611.546 -35584.633
[76,] 17805.179 73611.546
[77,] 49994.478 17805.179
[78,] -282147.420 49994.478
[79,] -465658.107 -282147.420
[80,] 48103.003 -465658.107
[81,] 47930.758 48103.003
[82,] 47930.758 47930.758
[83,] -37960.116 47930.758
[84,] 47930.758 -37960.116
[85,] -7555.060 47930.758
[86,] 46990.954 -7555.060
[87,] 48672.001 46990.954
[88,] 37073.804 48672.001
[89,] 35414.901 37073.804
[90,] -7555.060 35414.901
[91,] 41434.050 -7555.060
[92,] -273562.970 41434.050
[93,] -150276.044 -273562.970
[94,] 50535.241 -150276.044
[95,] 44777.497 50535.241
[96,] 47930.758 44777.497
[97,] 48453.678 47930.758
[98,] 53038.424 48453.678
[99,] 65563.601 53038.424
[100,] 152629.298 65563.601
[101,] 29601.307 152629.298
[102,] 50256.107 29601.307
[103,] 48181.113 50256.107
[104,] 47930.758 48181.113
[105,] -675705.948 47930.758
[106,] 19719.600 -675705.948
[107,] -86633.872 19719.600
[108,] 55446.705 -86633.872
[109,] 54791.004 55446.705
[110,] 43429.078 54791.004
[111,] -115267.526 43429.078
[112,] 47930.758 -115267.526
[113,] 23566.412 47930.758
[114,] -7555.060 23566.412
[115,] 45851.206 -7555.060
[116,] -53093.255 45851.206
[117,] 153902.340 -53093.255
[118,] 148737.922 153902.340
[119,] 66254.339 148737.922
[120,] -100883.026 66254.339
[121,] 47930.758 -100883.026
[122,] 190360.572 47930.758
[123,] 20323.606 190360.572
[124,] -344582.470 20323.606
[125,] -7555.060 -344582.470
[126,] 33983.667 -7555.060
[127,] 47930.758 33983.667
[128,] -1118357.214 47930.758
[129,] -65872.062 -1118357.214
[130,] 34349.869 -65872.062
[131,] -7555.060 34349.869
[132,] -7555.060 -7555.060
[133,] -7555.060 -7555.060
[134,] 46456.029 -7555.060
[135,] -11265.356 46456.029
[136,] 47930.758 -11265.356
[137,] 49743.121 47930.758
[138,] -7555.060 49743.121
[139,] 47930.758 -7555.060
[140,] 127681.941 47930.758
[141,] -123725.963 127681.941
[142,] 47930.758 -123725.963
[143,] 47930.758 47930.758
[144,] -2370.308 47930.758
[145,] -7555.060 -2370.308
[146,] 41821.079 -7555.060
[147,] 90994.114 41821.079
[148,] -35871.190 90994.114
[149,] 47930.758 -35871.190
[150,] -7555.060 47930.758
[151,] 17459.373 -7555.060
[152,] 47930.758 17459.373
[153,] 21375.207 47930.758
[154,] -49407.953 21375.207
[155,] -124797.099 -49407.953
[156,] 85548.186 -124797.099
[157,] 47930.758 85548.186
[158,] 40959.878 47930.758
[159,] 29543.810 40959.878
[160,] 42217.914 29543.810
[161,] 47930.758 42217.914
[162,] 48756.725 47930.758
[163,] -5631.489 48756.725
[164,] 47930.758 -5631.489
[165,] 80331.658 47930.758
[166,] 47930.758 80331.658
[167,] 50336.225 47930.758
[168,] 20183.871 50336.225
[169,] 37562.856 20183.871
[170,] -65081.319 37562.856
[171,] 63483.018 -65081.319
[172,] -200546.755 63483.018
[173,] -109444.008 -200546.755
[174,] 23336.527 -109444.008
[175,] 54252.289 23336.527
[176,] 47930.758 54252.289
[177,] 50983.381 47930.758
[178,] 29871.149 50983.381
[179,] -7555.060 29871.149
[180,] -28126.641 -7555.060
[181,] -6906.064 -28126.641
[182,] 47930.758 -6906.064
[183,] 47930.758 47930.758
[184,] -7555.060 47930.758
[185,] 51520.595 -7555.060
[186,] 43892.484 51520.595
[187,] 39825.868 43892.484
[188,] -2163938.266 39825.868
[189,] 47930.758 -2163938.266
[190,] 49666.953 47930.758
[191,] 10638.634 49666.953
[192,] 47930.758 10638.634
[193,] 59511.531 47930.758
[194,] 55714.984 59511.531
[195,] 47930.758 55714.984
[196,] -83722.349 47930.758
[197,] -223309.564 -83722.349
[198,] -7555.060 -223309.564
[199,] 47930.758 -7555.060
[200,] -279201.028 47930.758
[201,] -7555.060 -279201.028
[202,] -7555.060 -7555.060
[203,] 35201.158 -7555.060
[204,] 47930.758 35201.158
[205,] 63553.861 47930.758
[206,] -58373.820 63553.861
[207,] 47930.758 -58373.820
[208,] 51085.001 47930.758
[209,] 6174.204 51085.001
[210,] -12297.139 6174.204
[211,] -223683.527 -12297.139
[212,] 44403.489 -223683.527
[213,] 34065.162 44403.489
[214,] -122718.005 34065.162
[215,] -157816.795 -122718.005
[216,] -8066.752 -157816.795
[217,] 47930.758 -8066.752
[218,] 27844.604 47930.758
[219,] -29593.037 27844.604
[220,] 45014.662 -29593.037
[221,] 47930.758 45014.662
[222,] 34892.012 47930.758
[223,] 20843.505 34892.012
[224,] 25176.204 20843.505
[225,] -459875.779 25176.204
[226,] 47930.758 -459875.779
[227,] -234688.381 47930.758
[228,] 24439.537 -234688.381
[229,] 27297.538 24439.537
[230,] -144589.537 27297.538
[231,] 42833.808 -144589.537
[232,] -153760.016 42833.808
[233,] 47930.758 -153760.016
[234,] 43760.928 47930.758
[235,] 10402.487 43760.928
[236,] 55479.808 10402.487
[237,] 47930.758 55479.808
[238,] 47930.758 47930.758
[239,] 48830.870 47930.758
[240,] 12120.541 48830.870
[241,] 10492.987 12120.541
[242,] -7555.060 10492.987
[243,] 39310.301 -7555.060
[244,] -138046.328 39310.301
[245,] 118575.766 -138046.328
[246,] 47930.758 118575.766
[247,] 53714.507 47930.758
[248,] 47930.758 53714.507
[249,] 45207.361 47930.758
[250,] -7555.060 45207.361
[251,] 25991.297 -7555.060
[252,] -7555.060 25991.297
[253,] 47930.758 -7555.060
[254,] -7555.060 47930.758
[255,] 47930.758 -7555.060
[256,] 47930.758 47930.758
[257,] 47930.758 47930.758
[258,] 48357.958 47930.758
[259,] -374596.864 48357.958
[260,] 297382.198 -374596.864
[261,] -89503.531 297382.198
[262,] -95882.735 -89503.531
[263,] -43408.390 -95882.735
[264,] -107098.631 -43408.390
[265,] 47930.758 -107098.631
[266,] 24781.823 47930.758
[267,] 47930.758 24781.823
[268,] 39802.780 47930.758
[269,] 45723.816 39802.780
[270,] 47930.758 45723.816
[271,] 47930.758 47930.758
[272,] 10354.181 47930.758
[273,] 48001.084 10354.181
[274,] 105736.955 48001.084
[275,] 2068.080 105736.955
[276,] -8516.655 2068.080
[277,] 41418.222 -8516.655
[278,] 9037.666 41418.222
[279,] -20308.732 9037.666
[280,] 27433.327 -20308.732
[281,] -59079.502 27433.327
[282,] 47930.758 -59079.502
[283,] 19310.917 47930.758
[284,] 50337.774 19310.917
[285,] 47076.205 50337.774
[286,] 47930.758 47076.205
[287,] 111838.134 47930.758
[288,] 333879.073 111838.134
[289,] -576390.302 333879.073
[290,] -7555.060 -576390.302
[291,] 35997.643 -7555.060
[292,] 47930.758 35997.643
[293,] -12187.230 47930.758
[294,] 47968.115 -12187.230
[295,] 48340.985 47968.115
[296,] 47930.758 48340.985
[297,] 762830.868 47930.758
[298,] -7555.060 762830.868
[299,] -7555.060 -7555.060
[300,] -278718.282 -7555.060
[301,] -14179.696 -278718.282
[302,] -74418.783 -14179.696
[303,] 42624.223 -74418.783
[304,] 32019.062 42624.223
[305,] 48377.897 32019.062
[306,] -158893.721 48377.897
[307,] 47930.758 -158893.721
[308,] 47930.758 47930.758
[309,] 49543.695 47930.758
[310,] -7555.060 49543.695
[311,] -87387.956 -7555.060
[312,] -13574.240 -87387.956
[313,] 47890.901 -13574.240
[314,] 35089.843 47890.901
[315,] -189194.453 35089.843
[316,] 1187360.002 -189194.453
[317,] 16757.609 1187360.002
[318,] -7555.060 16757.609
[319,] -60797.317 -7555.060
[320,] 50834.989 -60797.317
[321,] 47803.348 50834.989
[322,] -6488.100 47803.348
[323,] 47930.758 -6488.100
[324,] -14508.840 47930.758
[325,] 47930.758 -14508.840
[326,] -67415.436 47930.758
[327,] 1259827.443 -67415.436
[328,] -1184133.917 1259827.443
[329,] -180351.327 -1184133.917
[330,] -62936.025 -180351.327
[331,] 47930.758 -62936.025
[332,] -7555.060 47930.758
[333,] -63103.461 -7555.060
[334,] -69214.165 -63103.461
[335,] 47930.758 -69214.165
[336,] 22362.228 47930.758
[337,] -378580.851 22362.228
[338,] 47930.758 -378580.851
[339,] -7555.060 47930.758
[340,] -265684.483 -7555.060
[341,] 13127.108 -265684.483
[342,] 47930.758 13127.108
[343,] 60434.071 47930.758
[344,] 47930.758 60434.071
[345,] -96498.844 47930.758
[346,] 44196.429 -96498.844
[347,] -313458.634 44196.429
[348,] 47930.758 -313458.634
[349,] 41837.252 47930.758
[350,] 47930.758 41837.252
[351,] 47930.758 47930.758
[352,] -9187.034 47930.758
[353,] 47930.758 -9187.034
[354,] -31544.333 47930.758
[355,] -7555.060 -31544.333
[356,] 47930.758 -7555.060
[357,] -7555.060 47930.758
[358,] 47930.758 -7555.060
[359,] -9824.282 47930.758
[360,] -348387.682 -9824.282
[361,] 58126.423 -348387.682
[362,] 15002.087 58126.423
[363,] -56691.365 15002.087
[364,] 55492.576 -56691.365
[365,] 408709.063 55492.576
[366,] 95780.130 408709.063
[367,] 49921.272 95780.130
[368,] 47930.758 49921.272
[369,] -232554.023 47930.758
[370,] -65383.903 -232554.023
[371,] 26690.904 -65383.903
[372,] 58494.152 26690.904
[373,] -159079.702 58494.152
[374,] -10313.459 -159079.702
[375,] 4473.442 -10313.459
[376,] -62195.598 4473.442
[377,] 5446.994 -62195.598
[378,] 28918.905 5446.994
[379,] 47930.758 28918.905
[380,] -8951.678 47930.758
[381,] 9957.267 -8951.678
[382,] -152906.495 9957.267
[383,] 49670.808 -152906.495
[384,] 48197.857 49670.808
[385,] 15549.014 48197.857
[386,] -18948.548 15549.014
[387,] -7555.060 -18948.548
[388,] 47930.758 -7555.060
[389,] -7555.060 47930.758
[390,] 47930.758 -7555.060
[391,] 47930.758 47930.758
[392,] 22800.421 47930.758
[393,] 52466.884 22800.421
[394,] 47930.758 52466.884
[395,] -7555.060 47930.758
[396,] 47930.758 -7555.060
[397,] 47930.758 47930.758
[398,] 60199.198 47930.758
[399,] 19942.976 60199.198
[400,] -127969.710 19942.976
[401,] -223588.738 -127969.710
[402,] 51524.780 -223588.738
[403,] -322029.866 51524.780
[404,] 43281.192 -322029.866
[405,] 22758.922 43281.192
[406,] 50155.084 22758.922
[407,] 225880.987 50155.084
[408,] -169212.146 225880.987
[409,] -5335.757 -169212.146
[410,] 20306.970 -5335.757
[411,] 47930.758 20306.970
[412,] -7555.060 47930.758
[413,] -73845.864 -7555.060
[414,] 48877.645 -73845.864
[415,] 50139.283 48877.645
[416,] -248826.358 50139.283
[417,] -40997.771 -248826.358
[418,] 48709.764 -40997.771
[419,] 47930.758 48709.764
[420,] 46110.459 47930.758
[421,] 97881.580 46110.459
[422,] 64892.532 97881.580
[423,] 22654.245 64892.532
[424,] -35757.824 22654.245
[425,] -484134.906 -35757.824
[426,] -7555.060 -484134.906
[427,] 51292.884 -7555.060
[428,] 68636.060 51292.884
[429,] 47930.758 68636.060
[430,] 75626.658 47930.758
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 -41076.166 11292.552
2 -19198.651 -41076.166
3 79985.494 -19198.651
4 -280953.361 79985.494
5 -62387.577 -280953.361
6 17215.934 -62387.577
7 49320.735 17215.934
8 63580.147 49320.735
9 38208.833 63580.147
10 -41888.465 38208.833
11 -68881.815 -41888.465
12 47930.758 -68881.815
13 95359.943 47930.758
14 40017.496 95359.943
15 16128.433 40017.496
16 51727.086 16128.433
17 60651.036 51727.086
18 52885.045 60651.036
19 48104.426 52885.045
20 47930.758 48104.426
21 59484.659 47930.758
22 47930.758 59484.659
23 4172.197 47930.758
24 23515.058 4172.197
25 40932.200 23515.058
26 -7827.369 40932.200
27 -131141.638 -7827.369
28 47930.758 -131141.638
29 17147.822 47930.758
30 -250216.313 17147.822
31 47930.758 -250216.313
32 65051.344 47930.758
33 71559.952 65051.344
34 -190418.162 71559.952
35 -148255.495 -190418.162
36 49490.851 -148255.495
37 -41774.594 49490.851
38 -25648.977 -41774.594
39 -126404.718 -25648.977
40 56121.615 -126404.718
41 129383.608 56121.615
42 -32278.740 129383.608
43 47930.758 -32278.740
44 -37101.259 47930.758
45 -11465.171 -37101.259
46 47930.758 -11465.171
47 37413.394 47930.758
48 -124535.591 37413.394
49 24236.170 -124535.591
50 -402357.714 24236.170
51 51398.186 -402357.714
52 -74411.743 51398.186
53 -133413.771 -74411.743
54 498000.209 -133413.771
55 -306359.379 498000.209
56 40596.813 -306359.379
57 -139539.338 40596.813
58 -112771.362 -139539.338
59 47930.758 -112771.362
60 47930.758 47930.758
61 -249808.092 47930.758
62 -23000.886 -249808.092
63 47847.411 -23000.886
64 -266560.840 47847.411
65 47253.085 -266560.840
66 -298183.103 47253.085
67 -25660.802 -298183.103
68 -33992.012 -25660.802
69 47930.758 -33992.012
70 704705.338 47930.758
71 -306374.089 704705.338
72 3195283.878 -306374.089
73 47930.758 3195283.878
74 -35584.633 47930.758
75 73611.546 -35584.633
76 17805.179 73611.546
77 49994.478 17805.179
78 -282147.420 49994.478
79 -465658.107 -282147.420
80 48103.003 -465658.107
81 47930.758 48103.003
82 47930.758 47930.758
83 -37960.116 47930.758
84 47930.758 -37960.116
85 -7555.060 47930.758
86 46990.954 -7555.060
87 48672.001 46990.954
88 37073.804 48672.001
89 35414.901 37073.804
90 -7555.060 35414.901
91 41434.050 -7555.060
92 -273562.970 41434.050
93 -150276.044 -273562.970
94 50535.241 -150276.044
95 44777.497 50535.241
96 47930.758 44777.497
97 48453.678 47930.758
98 53038.424 48453.678
99 65563.601 53038.424
100 152629.298 65563.601
101 29601.307 152629.298
102 50256.107 29601.307
103 48181.113 50256.107
104 47930.758 48181.113
105 -675705.948 47930.758
106 19719.600 -675705.948
107 -86633.872 19719.600
108 55446.705 -86633.872
109 54791.004 55446.705
110 43429.078 54791.004
111 -115267.526 43429.078
112 47930.758 -115267.526
113 23566.412 47930.758
114 -7555.060 23566.412
115 45851.206 -7555.060
116 -53093.255 45851.206
117 153902.340 -53093.255
118 148737.922 153902.340
119 66254.339 148737.922
120 -100883.026 66254.339
121 47930.758 -100883.026
122 190360.572 47930.758
123 20323.606 190360.572
124 -344582.470 20323.606
125 -7555.060 -344582.470
126 33983.667 -7555.060
127 47930.758 33983.667
128 -1118357.214 47930.758
129 -65872.062 -1118357.214
130 34349.869 -65872.062
131 -7555.060 34349.869
132 -7555.060 -7555.060
133 -7555.060 -7555.060
134 46456.029 -7555.060
135 -11265.356 46456.029
136 47930.758 -11265.356
137 49743.121 47930.758
138 -7555.060 49743.121
139 47930.758 -7555.060
140 127681.941 47930.758
141 -123725.963 127681.941
142 47930.758 -123725.963
143 47930.758 47930.758
144 -2370.308 47930.758
145 -7555.060 -2370.308
146 41821.079 -7555.060
147 90994.114 41821.079
148 -35871.190 90994.114
149 47930.758 -35871.190
150 -7555.060 47930.758
151 17459.373 -7555.060
152 47930.758 17459.373
153 21375.207 47930.758
154 -49407.953 21375.207
155 -124797.099 -49407.953
156 85548.186 -124797.099
157 47930.758 85548.186
158 40959.878 47930.758
159 29543.810 40959.878
160 42217.914 29543.810
161 47930.758 42217.914
162 48756.725 47930.758
163 -5631.489 48756.725
164 47930.758 -5631.489
165 80331.658 47930.758
166 47930.758 80331.658
167 50336.225 47930.758
168 20183.871 50336.225
169 37562.856 20183.871
170 -65081.319 37562.856
171 63483.018 -65081.319
172 -200546.755 63483.018
173 -109444.008 -200546.755
174 23336.527 -109444.008
175 54252.289 23336.527
176 47930.758 54252.289
177 50983.381 47930.758
178 29871.149 50983.381
179 -7555.060 29871.149
180 -28126.641 -7555.060
181 -6906.064 -28126.641
182 47930.758 -6906.064
183 47930.758 47930.758
184 -7555.060 47930.758
185 51520.595 -7555.060
186 43892.484 51520.595
187 39825.868 43892.484
188 -2163938.266 39825.868
189 47930.758 -2163938.266
190 49666.953 47930.758
191 10638.634 49666.953
192 47930.758 10638.634
193 59511.531 47930.758
194 55714.984 59511.531
195 47930.758 55714.984
196 -83722.349 47930.758
197 -223309.564 -83722.349
198 -7555.060 -223309.564
199 47930.758 -7555.060
200 -279201.028 47930.758
201 -7555.060 -279201.028
202 -7555.060 -7555.060
203 35201.158 -7555.060
204 47930.758 35201.158
205 63553.861 47930.758
206 -58373.820 63553.861
207 47930.758 -58373.820
208 51085.001 47930.758
209 6174.204 51085.001
210 -12297.139 6174.204
211 -223683.527 -12297.139
212 44403.489 -223683.527
213 34065.162 44403.489
214 -122718.005 34065.162
215 -157816.795 -122718.005
216 -8066.752 -157816.795
217 47930.758 -8066.752
218 27844.604 47930.758
219 -29593.037 27844.604
220 45014.662 -29593.037
221 47930.758 45014.662
222 34892.012 47930.758
223 20843.505 34892.012
224 25176.204 20843.505
225 -459875.779 25176.204
226 47930.758 -459875.779
227 -234688.381 47930.758
228 24439.537 -234688.381
229 27297.538 24439.537
230 -144589.537 27297.538
231 42833.808 -144589.537
232 -153760.016 42833.808
233 47930.758 -153760.016
234 43760.928 47930.758
235 10402.487 43760.928
236 55479.808 10402.487
237 47930.758 55479.808
238 47930.758 47930.758
239 48830.870 47930.758
240 12120.541 48830.870
241 10492.987 12120.541
242 -7555.060 10492.987
243 39310.301 -7555.060
244 -138046.328 39310.301
245 118575.766 -138046.328
246 47930.758 118575.766
247 53714.507 47930.758
248 47930.758 53714.507
249 45207.361 47930.758
250 -7555.060 45207.361
251 25991.297 -7555.060
252 -7555.060 25991.297
253 47930.758 -7555.060
254 -7555.060 47930.758
255 47930.758 -7555.060
256 47930.758 47930.758
257 47930.758 47930.758
258 48357.958 47930.758
259 -374596.864 48357.958
260 297382.198 -374596.864
261 -89503.531 297382.198
262 -95882.735 -89503.531
263 -43408.390 -95882.735
264 -107098.631 -43408.390
265 47930.758 -107098.631
266 24781.823 47930.758
267 47930.758 24781.823
268 39802.780 47930.758
269 45723.816 39802.780
270 47930.758 45723.816
271 47930.758 47930.758
272 10354.181 47930.758
273 48001.084 10354.181
274 105736.955 48001.084
275 2068.080 105736.955
276 -8516.655 2068.080
277 41418.222 -8516.655
278 9037.666 41418.222
279 -20308.732 9037.666
280 27433.327 -20308.732
281 -59079.502 27433.327
282 47930.758 -59079.502
283 19310.917 47930.758
284 50337.774 19310.917
285 47076.205 50337.774
286 47930.758 47076.205
287 111838.134 47930.758
288 333879.073 111838.134
289 -576390.302 333879.073
290 -7555.060 -576390.302
291 35997.643 -7555.060
292 47930.758 35997.643
293 -12187.230 47930.758
294 47968.115 -12187.230
295 48340.985 47968.115
296 47930.758 48340.985
297 762830.868 47930.758
298 -7555.060 762830.868
299 -7555.060 -7555.060
300 -278718.282 -7555.060
301 -14179.696 -278718.282
302 -74418.783 -14179.696
303 42624.223 -74418.783
304 32019.062 42624.223
305 48377.897 32019.062
306 -158893.721 48377.897
307 47930.758 -158893.721
308 47930.758 47930.758
309 49543.695 47930.758
310 -7555.060 49543.695
311 -87387.956 -7555.060
312 -13574.240 -87387.956
313 47890.901 -13574.240
314 35089.843 47890.901
315 -189194.453 35089.843
316 1187360.002 -189194.453
317 16757.609 1187360.002
318 -7555.060 16757.609
319 -60797.317 -7555.060
320 50834.989 -60797.317
321 47803.348 50834.989
322 -6488.100 47803.348
323 47930.758 -6488.100
324 -14508.840 47930.758
325 47930.758 -14508.840
326 -67415.436 47930.758
327 1259827.443 -67415.436
328 -1184133.917 1259827.443
329 -180351.327 -1184133.917
330 -62936.025 -180351.327
331 47930.758 -62936.025
332 -7555.060 47930.758
333 -63103.461 -7555.060
334 -69214.165 -63103.461
335 47930.758 -69214.165
336 22362.228 47930.758
337 -378580.851 22362.228
338 47930.758 -378580.851
339 -7555.060 47930.758
340 -265684.483 -7555.060
341 13127.108 -265684.483
342 47930.758 13127.108
343 60434.071 47930.758
344 47930.758 60434.071
345 -96498.844 47930.758
346 44196.429 -96498.844
347 -313458.634 44196.429
348 47930.758 -313458.634
349 41837.252 47930.758
350 47930.758 41837.252
351 47930.758 47930.758
352 -9187.034 47930.758
353 47930.758 -9187.034
354 -31544.333 47930.758
355 -7555.060 -31544.333
356 47930.758 -7555.060
357 -7555.060 47930.758
358 47930.758 -7555.060
359 -9824.282 47930.758
360 -348387.682 -9824.282
361 58126.423 -348387.682
362 15002.087 58126.423
363 -56691.365 15002.087
364 55492.576 -56691.365
365 408709.063 55492.576
366 95780.130 408709.063
367 49921.272 95780.130
368 47930.758 49921.272
369 -232554.023 47930.758
370 -65383.903 -232554.023
371 26690.904 -65383.903
372 58494.152 26690.904
373 -159079.702 58494.152
374 -10313.459 -159079.702
375 4473.442 -10313.459
376 -62195.598 4473.442
377 5446.994 -62195.598
378 28918.905 5446.994
379 47930.758 28918.905
380 -8951.678 47930.758
381 9957.267 -8951.678
382 -152906.495 9957.267
383 49670.808 -152906.495
384 48197.857 49670.808
385 15549.014 48197.857
386 -18948.548 15549.014
387 -7555.060 -18948.548
388 47930.758 -7555.060
389 -7555.060 47930.758
390 47930.758 -7555.060
391 47930.758 47930.758
392 22800.421 47930.758
393 52466.884 22800.421
394 47930.758 52466.884
395 -7555.060 47930.758
396 47930.758 -7555.060
397 47930.758 47930.758
398 60199.198 47930.758
399 19942.976 60199.198
400 -127969.710 19942.976
401 -223588.738 -127969.710
402 51524.780 -223588.738
403 -322029.866 51524.780
404 43281.192 -322029.866
405 22758.922 43281.192
406 50155.084 22758.922
407 225880.987 50155.084
408 -169212.146 225880.987
409 -5335.757 -169212.146
410 20306.970 -5335.757
411 47930.758 20306.970
412 -7555.060 47930.758
413 -73845.864 -7555.060
414 48877.645 -73845.864
415 50139.283 48877.645
416 -248826.358 50139.283
417 -40997.771 -248826.358
418 48709.764 -40997.771
419 47930.758 48709.764
420 46110.459 47930.758
421 97881.580 46110.459
422 64892.532 97881.580
423 22654.245 64892.532
424 -35757.824 22654.245
425 -484134.906 -35757.824
426 -7555.060 -484134.906
427 51292.884 -7555.060
428 68636.060 51292.884
429 47930.758 68636.060
430 75626.658 47930.758
> 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/7grlp1293276484.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/8grlp1293276484.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/9r0ks1293276484.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/10r0ks1293276484.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/11cj1y1293276484.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/12f1hm1293276484.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/13cbfc1293276484.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/14ftdi1293276484.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/150uc61293276484.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/164dtu1293276484.tab")
+ }
>
> try(system("convert tmp/1kz5g1293276484.ps tmp/1kz5g1293276484.png",intern=TRUE))
character(0)
> try(system("convert tmp/2c84j1293276484.ps tmp/2c84j1293276484.png",intern=TRUE))
character(0)
> try(system("convert tmp/3c84j1293276484.ps tmp/3c84j1293276484.png",intern=TRUE))
character(0)
> try(system("convert tmp/4c84j1293276484.ps tmp/4c84j1293276484.png",intern=TRUE))
character(0)
> try(system("convert tmp/5n0m41293276484.ps tmp/5n0m41293276484.png",intern=TRUE))
character(0)
> try(system("convert tmp/6n0m41293276484.ps tmp/6n0m41293276484.png",intern=TRUE))
character(0)
> try(system("convert tmp/7grlp1293276484.ps tmp/7grlp1293276484.png",intern=TRUE))
character(0)
> try(system("convert tmp/8grlp1293276484.ps tmp/8grlp1293276484.png",intern=TRUE))
character(0)
> try(system("convert tmp/9r0ks1293276484.ps tmp/9r0ks1293276484.png",intern=TRUE))
character(0)
> try(system("convert tmp/10r0ks1293276484.ps tmp/10r0ks1293276484.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
12.230 1.810 14.224