R version 2.15.2 (2012-10-26) -- "Trick or Treat"
Copyright (C) 2012 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i686-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(210907
+ ,56
+ ,396
+ ,3
+ ,30
+ ,112285
+ ,1
+ ,120982
+ ,56
+ ,297
+ ,4
+ ,28
+ ,84786
+ ,1
+ ,176508
+ ,54
+ ,559
+ ,12
+ ,38
+ ,83123
+ ,1
+ ,179321
+ ,89
+ ,967
+ ,2
+ ,30
+ ,101193
+ ,1
+ ,123185
+ ,40
+ ,270
+ ,1
+ ,22
+ ,38361
+ ,1
+ ,52746
+ ,25
+ ,143
+ ,3
+ ,26
+ ,68504
+ ,1
+ ,385534
+ ,92
+ ,1562
+ ,0
+ ,25
+ ,119182
+ ,1
+ ,33170
+ ,18
+ ,109
+ ,0
+ ,18
+ ,22807
+ ,1
+ ,101645
+ ,63
+ ,371
+ ,0
+ ,11
+ ,17140
+ ,0
+ ,149061
+ ,44
+ ,656
+ ,5
+ ,26
+ ,116174
+ ,1
+ ,165446
+ ,33
+ ,511
+ ,0
+ ,25
+ ,57635
+ ,1
+ ,237213
+ ,84
+ ,655
+ ,0
+ ,38
+ ,66198
+ ,1
+ ,173326
+ ,88
+ ,465
+ ,7
+ ,44
+ ,71701
+ ,1
+ ,133131
+ ,55
+ ,525
+ ,7
+ ,30
+ ,57793
+ ,1
+ ,258873
+ ,60
+ ,885
+ ,3
+ ,40
+ ,80444
+ ,1
+ ,180083
+ ,66
+ ,497
+ ,9
+ ,34
+ ,53855
+ ,1
+ ,324799
+ ,154
+ ,1436
+ ,0
+ ,47
+ ,97668
+ ,1
+ ,230964
+ ,53
+ ,612
+ ,4
+ ,30
+ ,133824
+ ,1
+ ,236785
+ ,119
+ ,865
+ ,3
+ ,31
+ ,101481
+ ,1
+ ,135473
+ ,41
+ ,385
+ ,0
+ ,23
+ ,99645
+ ,1
+ ,202925
+ ,61
+ ,567
+ ,7
+ ,36
+ ,114789
+ ,1
+ ,215147
+ ,58
+ ,639
+ ,0
+ ,36
+ ,99052
+ ,1
+ ,344297
+ ,75
+ ,963
+ ,1
+ ,30
+ ,67654
+ ,1
+ ,153935
+ ,33
+ ,398
+ ,5
+ ,25
+ ,65553
+ ,1
+ ,132943
+ ,40
+ ,410
+ ,7
+ ,39
+ ,97500
+ ,1
+ ,174724
+ ,92
+ ,966
+ ,0
+ ,34
+ ,69112
+ ,1
+ ,174415
+ ,100
+ ,801
+ ,0
+ ,31
+ ,82753
+ ,1
+ ,225548
+ ,112
+ ,892
+ ,5
+ ,31
+ ,85323
+ ,1
+ ,223632
+ ,73
+ ,513
+ ,0
+ ,33
+ ,72654
+ ,1
+ ,124817
+ ,40
+ ,469
+ ,0
+ ,25
+ ,30727
+ ,1
+ ,221698
+ ,45
+ ,683
+ ,0
+ ,33
+ ,77873
+ ,1
+ ,210767
+ ,60
+ ,643
+ ,3
+ ,35
+ ,117478
+ ,1
+ ,170266
+ ,62
+ ,535
+ ,4
+ ,42
+ ,74007
+ ,1
+ ,260561
+ ,75
+ ,625
+ ,1
+ ,43
+ ,90183
+ ,1
+ ,84853
+ ,31
+ ,264
+ ,4
+ ,30
+ ,61542
+ ,1
+ ,294424
+ ,77
+ ,992
+ ,2
+ ,33
+ ,101494
+ ,1
+ ,101011
+ ,34
+ ,238
+ ,0
+ ,13
+ ,27570
+ ,1
+ ,215641
+ ,46
+ ,818
+ ,0
+ ,32
+ ,55813
+ ,1
+ ,325107
+ ,99
+ ,937
+ ,0
+ ,36
+ ,79215
+ ,1
+ ,7176
+ ,17
+ ,70
+ ,0
+ ,0
+ ,1423
+ ,1
+ ,167542
+ ,66
+ ,507
+ ,2
+ ,28
+ ,55461
+ ,1
+ ,106408
+ ,30
+ ,260
+ ,1
+ ,14
+ ,31081
+ ,1
+ ,96560
+ ,76
+ ,503
+ ,0
+ ,17
+ ,22996
+ ,1
+ ,265769
+ ,146
+ ,927
+ ,2
+ ,32
+ ,83122
+ ,1
+ ,269651
+ ,67
+ ,1269
+ ,10
+ ,30
+ ,70106
+ ,1
+ ,149112
+ ,56
+ ,537
+ ,6
+ ,35
+ ,60578
+ ,1
+ ,175824
+ ,107
+ ,910
+ ,0
+ ,20
+ ,39992
+ ,1
+ ,152871
+ ,58
+ ,532
+ ,5
+ ,28
+ ,79892
+ ,1
+ ,111665
+ ,34
+ ,345
+ ,4
+ ,28
+ ,49810
+ ,1
+ ,116408
+ ,61
+ ,918
+ ,1
+ ,39
+ ,71570
+ ,1
+ ,362301
+ ,119
+ ,1635
+ ,2
+ ,34
+ ,100708
+ ,1
+ ,78800
+ ,42
+ ,330
+ ,2
+ ,26
+ ,33032
+ ,1
+ ,183167
+ ,66
+ ,557
+ ,0
+ ,39
+ ,82875
+ ,1
+ ,277965
+ ,89
+ ,1178
+ ,8
+ ,39
+ ,139077
+ ,1
+ ,150629
+ ,44
+ ,740
+ ,3
+ ,33
+ ,71595
+ ,1
+ ,168809
+ ,66
+ ,452
+ ,0
+ ,28
+ ,72260
+ ,1
+ ,24188
+ ,24
+ ,218
+ ,0
+ ,4
+ ,5950
+ ,1
+ ,329267
+ ,259
+ ,764
+ ,8
+ ,39
+ ,115762
+ ,1
+ ,65029
+ ,17
+ ,255
+ ,5
+ ,18
+ ,32551
+ ,1
+ ,101097
+ ,64
+ ,454
+ ,3
+ ,14
+ ,31701
+ ,1
+ ,218946
+ ,41
+ ,866
+ ,1
+ ,29
+ ,80670
+ ,1
+ ,244052
+ ,68
+ ,574
+ ,5
+ ,44
+ ,143558
+ ,1
+ ,341570
+ ,168
+ ,1276
+ ,1
+ ,21
+ ,117105
+ ,1
+ ,103597
+ ,43
+ ,379
+ ,1
+ ,16
+ ,23789
+ ,1
+ ,233328
+ ,132
+ ,825
+ ,5
+ ,28
+ ,120733
+ ,1
+ ,256462
+ ,105
+ ,798
+ ,0
+ ,35
+ ,105195
+ ,1
+ ,206161
+ ,71
+ ,663
+ ,12
+ ,28
+ ,73107
+ ,1
+ ,311473
+ ,112
+ ,1069
+ ,8
+ ,38
+ ,132068
+ ,1
+ ,235800
+ ,94
+ ,921
+ ,8
+ ,23
+ ,149193
+ ,1
+ ,177939
+ ,82
+ ,858
+ ,8
+ ,36
+ ,46821
+ ,1
+ ,207176
+ ,70
+ ,711
+ ,8
+ ,32
+ ,87011
+ ,1
+ ,196553
+ ,57
+ ,503
+ ,2
+ ,29
+ ,95260
+ ,1
+ ,174184
+ ,53
+ ,382
+ ,0
+ ,25
+ ,55183
+ ,1
+ ,143246
+ ,103
+ ,464
+ ,5
+ ,27
+ ,106671
+ ,1
+ ,187559
+ ,121
+ ,717
+ ,8
+ ,36
+ ,73511
+ ,1
+ ,187681
+ ,62
+ ,690
+ ,2
+ ,28
+ ,92945
+ ,1
+ ,119016
+ ,52
+ ,462
+ ,5
+ ,23
+ ,78664
+ ,1
+ ,182192
+ ,52
+ ,657
+ ,12
+ ,40
+ ,70054
+ ,1
+ ,73566
+ ,32
+ ,385
+ ,6
+ ,23
+ ,22618
+ ,1
+ ,194979
+ ,62
+ ,577
+ ,7
+ ,40
+ ,74011
+ ,1
+ ,167488
+ ,45
+ ,619
+ ,2
+ ,28
+ ,83737
+ ,1
+ ,143756
+ ,46
+ ,479
+ ,0
+ ,34
+ ,69094
+ ,1
+ ,275541
+ ,63
+ ,817
+ ,4
+ ,33
+ ,93133
+ ,1
+ ,243199
+ ,75
+ ,752
+ ,3
+ ,28
+ ,95536
+ ,1
+ ,182999
+ ,88
+ ,430
+ ,6
+ ,34
+ ,225920
+ ,1
+ ,135649
+ ,46
+ ,451
+ ,2
+ ,30
+ ,62133
+ ,1
+ ,152299
+ ,53
+ ,537
+ ,0
+ ,33
+ ,61370
+ ,1
+ ,120221
+ ,37
+ ,519
+ ,1
+ ,22
+ ,43836
+ ,1
+ ,346485
+ ,90
+ ,1000
+ ,0
+ ,38
+ ,106117
+ ,1
+ ,145790
+ ,63
+ ,637
+ ,5
+ ,26
+ ,38692
+ ,1
+ ,193339
+ ,78
+ ,465
+ ,2
+ ,35
+ ,84651
+ ,1
+ ,80953
+ ,25
+ ,437
+ ,0
+ ,8
+ ,56622
+ ,1
+ ,122774
+ ,45
+ ,711
+ ,0
+ ,24
+ ,15986
+ ,1
+ ,130585
+ ,46
+ ,299
+ ,5
+ ,29
+ ,95364
+ ,1
+ ,112611
+ ,41
+ ,248
+ ,0
+ ,20
+ ,26706
+ ,0
+ ,286468
+ ,144
+ ,1162
+ ,1
+ ,29
+ ,89691
+ ,1
+ ,241066
+ ,82
+ ,714
+ ,0
+ ,45
+ ,67267
+ ,1
+ ,148446
+ ,91
+ ,905
+ ,1
+ ,37
+ ,126846
+ ,1
+ ,204713
+ ,71
+ ,649
+ ,1
+ ,33
+ ,41140
+ ,1
+ ,182079
+ ,63
+ ,512
+ ,2
+ ,33
+ ,102860
+ ,1
+ ,140344
+ ,53
+ ,472
+ ,6
+ ,25
+ ,51715
+ ,1
+ ,220516
+ ,62
+ ,905
+ ,1
+ ,32
+ ,55801
+ ,1
+ ,243060
+ ,63
+ ,786
+ ,4
+ ,29
+ ,111813
+ ,1
+ ,162765
+ ,32
+ ,489
+ ,2
+ ,28
+ ,120293
+ ,1
+ ,182613
+ ,39
+ ,479
+ ,3
+ ,28
+ ,138599
+ ,1
+ ,232138
+ ,62
+ ,617
+ ,0
+ ,31
+ ,161647
+ ,1
+ ,265318
+ ,117
+ ,925
+ ,10
+ ,52
+ ,115929
+ ,1
+ ,85574
+ ,34
+ ,351
+ ,0
+ ,21
+ ,24266
+ ,0
+ ,310839
+ ,92
+ ,1144
+ ,9
+ ,24
+ ,162901
+ ,1
+ ,225060
+ ,93
+ ,669
+ ,7
+ ,41
+ ,109825
+ ,1
+ ,232317
+ ,54
+ ,707
+ ,0
+ ,33
+ ,129838
+ ,1
+ ,144966
+ ,144
+ ,458
+ ,0
+ ,32
+ ,37510
+ ,1
+ ,43287
+ ,14
+ ,214
+ ,4
+ ,19
+ ,43750
+ ,1
+ ,155754
+ ,61
+ ,599
+ ,4
+ ,20
+ ,40652
+ ,1
+ ,164709
+ ,109
+ ,572
+ ,0
+ ,31
+ ,87771
+ ,1
+ ,201940
+ ,38
+ ,897
+ ,0
+ ,31
+ ,85872
+ ,1
+ ,235454
+ ,73
+ ,819
+ ,0
+ ,32
+ ,89275
+ ,1
+ ,220801
+ ,75
+ ,720
+ ,1
+ ,18
+ ,44418
+ ,0
+ ,99466
+ ,50
+ ,273
+ ,0
+ ,23
+ ,192565
+ ,1
+ ,92661
+ ,61
+ ,508
+ ,1
+ ,17
+ ,35232
+ ,0
+ ,133328
+ ,55
+ ,506
+ ,0
+ ,20
+ ,40909
+ ,0
+ ,61361
+ ,77
+ ,451
+ ,0
+ ,12
+ ,13294
+ ,0
+ ,125930
+ ,75
+ ,699
+ ,4
+ ,17
+ ,32387
+ ,0
+ ,100750
+ ,72
+ ,407
+ ,0
+ ,30
+ ,140867
+ ,1
+ ,224549
+ ,50
+ ,465
+ ,4
+ ,31
+ ,120662
+ ,1
+ ,82316
+ ,32
+ ,245
+ ,4
+ ,10
+ ,21233
+ ,0
+ ,102010
+ ,53
+ ,370
+ ,3
+ ,13
+ ,44332
+ ,0
+ ,101523
+ ,42
+ ,316
+ ,0
+ ,22
+ ,61056
+ ,0
+ ,243511
+ ,71
+ ,603
+ ,0
+ ,42
+ ,101338
+ ,1
+ ,22938
+ ,10
+ ,154
+ ,0
+ ,1
+ ,1168
+ ,1
+ ,41566
+ ,35
+ ,229
+ ,5
+ ,9
+ ,13497
+ ,0
+ ,152474
+ ,65
+ ,577
+ ,0
+ ,32
+ ,65567
+ ,1
+ ,61857
+ ,25
+ ,192
+ ,4
+ ,11
+ ,25162
+ ,1
+ ,99923
+ ,66
+ ,617
+ ,0
+ ,25
+ ,32334
+ ,0
+ ,132487
+ ,41
+ ,411
+ ,0
+ ,36
+ ,40735
+ ,1
+ ,317394
+ ,86
+ ,975
+ ,1
+ ,31
+ ,91413
+ ,1
+ ,21054
+ ,16
+ ,146
+ ,0
+ ,0
+ ,855
+ ,1
+ ,209641
+ ,42
+ ,705
+ ,5
+ ,24
+ ,97068
+ ,1
+ ,22648
+ ,19
+ ,184
+ ,0
+ ,13
+ ,44339
+ ,0
+ ,31414
+ ,19
+ ,200
+ ,0
+ ,8
+ ,14116
+ ,1
+ ,46698
+ ,45
+ ,274
+ ,0
+ ,13
+ ,10288
+ ,0
+ ,131698
+ ,65
+ ,502
+ ,0
+ ,19
+ ,65622
+ ,0
+ ,91735
+ ,35
+ ,382
+ ,0
+ ,18
+ ,16563
+ ,0
+ ,244749
+ ,95
+ ,964
+ ,2
+ ,33
+ ,76643
+ ,1
+ ,184510
+ ,49
+ ,537
+ ,7
+ ,40
+ ,110681
+ ,1
+ ,79863
+ ,37
+ ,438
+ ,1
+ ,22
+ ,29011
+ ,0
+ ,128423
+ ,64
+ ,369
+ ,8
+ ,38
+ ,92696
+ ,1
+ ,97839
+ ,38
+ ,417
+ ,2
+ ,24
+ ,94785
+ ,1
+ ,38214
+ ,34
+ ,276
+ ,0
+ ,8
+ ,8773
+ ,1
+ ,151101
+ ,32
+ ,514
+ ,2
+ ,35
+ ,83209
+ ,1
+ ,272458
+ ,65
+ ,822
+ ,0
+ ,43
+ ,93815
+ ,1
+ ,172494
+ ,52
+ ,389
+ ,0
+ ,43
+ ,86687
+ ,1
+ ,108043
+ ,62
+ ,466
+ ,1
+ ,14
+ ,34553
+ ,0
+ ,328107
+ ,65
+ ,1255
+ ,3
+ ,41
+ ,105547
+ ,1
+ ,250579
+ ,83
+ ,694
+ ,0
+ ,38
+ ,103487
+ ,1
+ ,351067
+ ,95
+ ,1024
+ ,3
+ ,45
+ ,213688
+ ,1
+ ,158015
+ ,29
+ ,400
+ ,0
+ ,31
+ ,71220
+ ,1
+ ,98866
+ ,18
+ ,397
+ ,0
+ ,13
+ ,23517
+ ,0
+ ,85439
+ ,33
+ ,350
+ ,0
+ ,28
+ ,56926
+ ,1
+ ,229242
+ ,247
+ ,719
+ ,4
+ ,31
+ ,91721
+ ,1
+ ,351619
+ ,139
+ ,1277
+ ,4
+ ,40
+ ,115168
+ ,1
+ ,84207
+ ,29
+ ,356
+ ,11
+ ,30
+ ,111194
+ ,1
+ ,120445
+ ,118
+ ,457
+ ,0
+ ,16
+ ,51009
+ ,0
+ ,324598
+ ,110
+ ,1402
+ ,0
+ ,37
+ ,135777
+ ,1
+ ,131069
+ ,67
+ ,600
+ ,4
+ ,30
+ ,51513
+ ,1
+ ,204271
+ ,42
+ ,480
+ ,0
+ ,35
+ ,74163
+ ,1
+ ,165543
+ ,65
+ ,595
+ ,1
+ ,32
+ ,51633
+ ,1
+ ,141722
+ ,94
+ ,436
+ ,0
+ ,27
+ ,75345
+ ,1
+ ,116048
+ ,64
+ ,230
+ ,0
+ ,20
+ ,33416
+ ,0
+ ,250047
+ ,81
+ ,651
+ ,0
+ ,18
+ ,83305
+ ,0
+ ,299775
+ ,95
+ ,1367
+ ,9
+ ,31
+ ,98952
+ ,1
+ ,195838
+ ,67
+ ,564
+ ,1
+ ,31
+ ,102372
+ ,1
+ ,173260
+ ,63
+ ,716
+ ,3
+ ,21
+ ,37238
+ ,1
+ ,254488
+ ,83
+ ,747
+ ,10
+ ,39
+ ,103772
+ ,1
+ ,104389
+ ,45
+ ,467
+ ,5
+ ,41
+ ,123969
+ ,1
+ ,136084
+ ,30
+ ,671
+ ,0
+ ,13
+ ,27142
+ ,0
+ ,199476
+ ,70
+ ,861
+ ,2
+ ,32
+ ,135400
+ ,1
+ ,92499
+ ,32
+ ,319
+ ,0
+ ,18
+ ,21399
+ ,0
+ ,224330
+ ,83
+ ,612
+ ,1
+ ,39
+ ,130115
+ ,1
+ ,135781
+ ,31
+ ,433
+ ,2
+ ,14
+ ,24874
+ ,0
+ ,74408
+ ,67
+ ,434
+ ,4
+ ,7
+ ,34988
+ ,0
+ ,81240
+ ,66
+ ,503
+ ,0
+ ,17
+ ,45549
+ ,0
+ ,14688
+ ,10
+ ,85
+ ,0
+ ,0
+ ,6023
+ ,1
+ ,181633
+ ,70
+ ,564
+ ,2
+ ,30
+ ,64466
+ ,1
+ ,271856
+ ,103
+ ,824
+ ,1
+ ,37
+ ,54990
+ ,1
+ ,7199
+ ,5
+ ,74
+ ,0
+ ,0
+ ,1644
+ ,1
+ ,46660
+ ,20
+ ,259
+ ,0
+ ,5
+ ,6179
+ ,1
+ ,17547
+ ,5
+ ,69
+ ,0
+ ,1
+ ,3926
+ ,1
+ ,133368
+ ,36
+ ,535
+ ,1
+ ,16
+ ,32755
+ ,0
+ ,95227
+ ,34
+ ,239
+ ,0
+ ,32
+ ,34777
+ ,1
+ ,152601
+ ,48
+ ,438
+ ,2
+ ,24
+ ,73224
+ ,1
+ ,98146
+ ,40
+ ,459
+ ,0
+ ,17
+ ,27114
+ ,0
+ ,79619
+ ,43
+ ,426
+ ,3
+ ,11
+ ,20760
+ ,0
+ ,59194
+ ,31
+ ,288
+ ,6
+ ,24
+ ,37636
+ ,0
+ ,139942
+ ,42
+ ,498
+ ,0
+ ,22
+ ,65461
+ ,0
+ ,118612
+ ,46
+ ,454
+ ,2
+ ,12
+ ,30080
+ ,0
+ ,72880
+ ,33
+ ,376
+ ,0
+ ,19
+ ,24094
+ ,0
+ ,65475
+ ,18
+ ,225
+ ,2
+ ,13
+ ,69008
+ ,0
+ ,99643
+ ,55
+ ,555
+ ,1
+ ,17
+ ,54968
+ ,0
+ ,71965
+ ,35
+ ,252
+ ,1
+ ,15
+ ,46090
+ ,0
+ ,77272
+ ,59
+ ,208
+ ,2
+ ,16
+ ,27507
+ ,0
+ ,49289
+ ,19
+ ,130
+ ,1
+ ,24
+ ,10672
+ ,0
+ ,135131
+ ,66
+ ,481
+ ,0
+ ,15
+ ,34029
+ ,0
+ ,108446
+ ,60
+ ,389
+ ,1
+ ,17
+ ,46300
+ ,0
+ ,89746
+ ,36
+ ,565
+ ,3
+ ,18
+ ,24760
+ ,0
+ ,44296
+ ,25
+ ,173
+ ,0
+ ,20
+ ,18779
+ ,0
+ ,77648
+ ,47
+ ,278
+ ,0
+ ,16
+ ,21280
+ ,0
+ ,181528
+ ,54
+ ,609
+ ,0
+ ,16
+ ,40662
+ ,0
+ ,134019
+ ,53
+ ,422
+ ,0
+ ,18
+ ,28987
+ ,0
+ ,124064
+ ,40
+ ,445
+ ,1
+ ,22
+ ,22827
+ ,0
+ ,92630
+ ,40
+ ,387
+ ,4
+ ,8
+ ,18513
+ ,0
+ ,121848
+ ,39
+ ,339
+ ,0
+ ,17
+ ,30594
+ ,0
+ ,52915
+ ,14
+ ,181
+ ,0
+ ,18
+ ,24006
+ ,0
+ ,81872
+ ,45
+ ,245
+ ,0
+ ,16
+ ,27913
+ ,0
+ ,58981
+ ,36
+ ,384
+ ,7
+ ,23
+ ,42744
+ ,0
+ ,53515
+ ,28
+ ,212
+ ,2
+ ,22
+ ,12934
+ ,0
+ ,60812
+ ,44
+ ,399
+ ,0
+ ,13
+ ,22574
+ ,0
+ ,56375
+ ,30
+ ,229
+ ,7
+ ,13
+ ,41385
+ ,0
+ ,65490
+ ,22
+ ,224
+ ,3
+ ,16
+ ,18653
+ ,0
+ ,80949
+ ,17
+ ,203
+ ,0
+ ,16
+ ,18472
+ ,0
+ ,76302
+ ,31
+ ,333
+ ,0
+ ,20
+ ,30976
+ ,0
+ ,104011
+ ,55
+ ,384
+ ,6
+ ,22
+ ,63339
+ ,0
+ ,98104
+ ,54
+ ,636
+ ,2
+ ,17
+ ,25568
+ ,0
+ ,67989
+ ,21
+ ,185
+ ,0
+ ,18
+ ,33747
+ ,0
+ ,30989
+ ,14
+ ,93
+ ,0
+ ,17
+ ,4154
+ ,0
+ ,135458
+ ,81
+ ,581
+ ,3
+ ,12
+ ,19474
+ ,0
+ ,73504
+ ,35
+ ,248
+ ,0
+ ,7
+ ,35130
+ ,0
+ ,63123
+ ,43
+ ,304
+ ,1
+ ,17
+ ,39067
+ ,0
+ ,61254
+ ,46
+ ,344
+ ,1
+ ,14
+ ,13310
+ ,0
+ ,74914
+ ,30
+ ,407
+ ,0
+ ,23
+ ,65892
+ ,0
+ ,31774
+ ,23
+ ,170
+ ,1
+ ,17
+ ,4143
+ ,0
+ ,81437
+ ,38
+ ,312
+ ,0
+ ,14
+ ,28579
+ ,0
+ ,87186
+ ,54
+ ,507
+ ,0
+ ,15
+ ,51776
+ ,0
+ ,50090
+ ,20
+ ,224
+ ,0
+ ,17
+ ,21152
+ ,0
+ ,65745
+ ,53
+ ,340
+ ,0
+ ,21
+ ,38084
+ ,0
+ ,56653
+ ,45
+ ,168
+ ,0
+ ,18
+ ,27717
+ ,0
+ ,158399
+ ,39
+ ,443
+ ,0
+ ,18
+ ,32928
+ ,0
+ ,46455
+ ,20
+ ,204
+ ,0
+ ,17
+ ,11342
+ ,0
+ ,73624
+ ,24
+ ,367
+ ,0
+ ,17
+ ,19499
+ ,0
+ ,38395
+ ,31
+ ,210
+ ,0
+ ,16
+ ,16380
+ ,0
+ ,91899
+ ,35
+ ,335
+ ,0
+ ,15
+ ,36874
+ ,0
+ ,139526
+ ,151
+ ,364
+ ,0
+ ,21
+ ,48259
+ ,0
+ ,52164
+ ,52
+ ,178
+ ,0
+ ,16
+ ,16734
+ ,0
+ ,51567
+ ,30
+ ,206
+ ,2
+ ,14
+ ,28207
+ ,0
+ ,70551
+ ,31
+ ,279
+ ,0
+ ,15
+ ,30143
+ ,0
+ ,84856
+ ,29
+ ,387
+ ,1
+ ,17
+ ,41369
+ ,0
+ ,102538
+ ,57
+ ,490
+ ,1
+ ,15
+ ,45833
+ ,0
+ ,86678
+ ,40
+ ,238
+ ,0
+ ,15
+ ,29156
+ ,0
+ ,85709
+ ,44
+ ,343
+ ,0
+ ,10
+ ,35944
+ ,0
+ ,34662
+ ,25
+ ,232
+ ,0
+ ,6
+ ,36278
+ ,0
+ ,150580
+ ,77
+ ,530
+ ,0
+ ,22
+ ,45588
+ ,0
+ ,99611
+ ,35
+ ,291
+ ,0
+ ,21
+ ,45097
+ ,0
+ ,19349
+ ,11
+ ,67
+ ,0
+ ,1
+ ,3895
+ ,0
+ ,99373
+ ,63
+ ,397
+ ,1
+ ,18
+ ,28394
+ ,0
+ ,86230
+ ,44
+ ,467
+ ,0
+ ,17
+ ,18632
+ ,0
+ ,30837
+ ,19
+ ,178
+ ,0
+ ,4
+ ,2325
+ ,0
+ ,31706
+ ,13
+ ,175
+ ,0
+ ,10
+ ,25139
+ ,0
+ ,89806
+ ,42
+ ,299
+ ,0
+ ,16
+ ,27975
+ ,0
+ ,62088
+ ,38
+ ,154
+ ,1
+ ,16
+ ,14483
+ ,0
+ ,40151
+ ,29
+ ,106
+ ,0
+ ,9
+ ,13127
+ ,0
+ ,27634
+ ,20
+ ,189
+ ,0
+ ,16
+ ,5839
+ ,0
+ ,76990
+ ,27
+ ,194
+ ,0
+ ,17
+ ,24069
+ ,0
+ ,37460
+ ,20
+ ,135
+ ,0
+ ,7
+ ,3738
+ ,0
+ ,54157
+ ,19
+ ,201
+ ,0
+ ,15
+ ,18625
+ ,0
+ ,49862
+ ,37
+ ,207
+ ,0
+ ,14
+ ,36341
+ ,0
+ ,84337
+ ,26
+ ,280
+ ,0
+ ,14
+ ,24548
+ ,0
+ ,64175
+ ,42
+ ,260
+ ,0
+ ,18
+ ,21792
+ ,0
+ ,59382
+ ,49
+ ,227
+ ,0
+ ,12
+ ,26263
+ ,0
+ ,119308
+ ,30
+ ,239
+ ,0
+ ,16
+ ,23686
+ ,0
+ ,76702
+ ,49
+ ,333
+ ,0
+ ,21
+ ,49303
+ ,0
+ ,103425
+ ,67
+ ,428
+ ,1
+ ,19
+ ,25659
+ ,0
+ ,70344
+ ,28
+ ,230
+ ,0
+ ,16
+ ,28904
+ ,0
+ ,43410
+ ,19
+ ,292
+ ,0
+ ,1
+ ,2781
+ ,0
+ ,104838
+ ,49
+ ,350
+ ,1
+ ,16
+ ,29236
+ ,0
+ ,62215
+ ,27
+ ,186
+ ,0
+ ,10
+ ,19546
+ ,0
+ ,69304
+ ,30
+ ,326
+ ,6
+ ,19
+ ,22818
+ ,0
+ ,53117
+ ,22
+ ,155
+ ,3
+ ,12
+ ,32689
+ ,0
+ ,19764
+ ,12
+ ,75
+ ,1
+ ,2
+ ,5752
+ ,0
+ ,86680
+ ,31
+ ,361
+ ,2
+ ,14
+ ,22197
+ ,0
+ ,84105
+ ,20
+ ,261
+ ,0
+ ,17
+ ,20055
+ ,0
+ ,77945
+ ,20
+ ,299
+ ,0
+ ,19
+ ,25272
+ ,0
+ ,89113
+ ,39
+ ,300
+ ,0
+ ,14
+ ,82206
+ ,0
+ ,91005
+ ,29
+ ,450
+ ,3
+ ,11
+ ,32073
+ ,0
+ ,40248
+ ,16
+ ,183
+ ,1
+ ,4
+ ,5444
+ ,0
+ ,64187
+ ,27
+ ,238
+ ,0
+ ,16
+ ,20154
+ ,0
+ ,50857
+ ,21
+ ,165
+ ,0
+ ,20
+ ,36944
+ ,0
+ ,56613
+ ,19
+ ,234
+ ,1
+ ,12
+ ,8019
+ ,0
+ ,62792
+ ,35
+ ,176
+ ,0
+ ,15
+ ,30884
+ ,0
+ ,72535
+ ,14
+ ,329
+ ,0
+ ,16
+ ,19540
+ ,0)
+ ,dim=c(7
+ ,289)
+ ,dimnames=list(c('Time'
+ ,'Logins'
+ ,'Views'
+ ,'Shared'
+ ,'Reviewed'
+ ,'TotalSize'
+ ,'Pop')
+ ,1:289))
> y <- array(NA,dim=c(7,289),dimnames=list(c('Time','Logins','Views','Shared','Reviewed','TotalSize','Pop'),1:289))
> for (i in 1:dim(x)[1])
+ {
+ for (j in 1:dim(x)[2])
+ {
+ y[i,j] <- as.numeric(x[i,j])
+ }
+ }
> par3 = 'No Linear Trend'
> par2 = 'Do not include Seasonal Dummies'
> par1 = '1'
> library(lattice)
> library(lmtest)
Loading required package: zoo
Attaching package: 'zoo'
The following object(s) are masked from 'package:base':
as.Date, as.Date.numeric
> n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test
> par1 <- as.numeric(par1)
> x <- t(y)
> k <- length(x[1,])
> n <- length(x[,1])
> x1 <- cbind(x[,par1], x[,1:k!=par1])
> mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1])
> colnames(x1) <- mycolnames #colnames(x)[par1]
> x <- x1
> if (par3 == 'First Differences'){
+ x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep='')))
+ for (i in 1:n-1) {
+ for (j in 1:k) {
+ x2[i,j] <- x[i+1,j] - x[i,j]
+ }
+ }
+ x <- x2
+ }
> if (par2 == 'Include Monthly Dummies'){
+ x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep ='')))
+ for (i in 1:11){
+ x2[seq(i,n,12),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> if (par2 == 'Include Quarterly Dummies'){
+ x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep ='')))
+ for (i in 1:3){
+ x2[seq(i,n,4),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> k <- length(x[1,])
> if (par3 == 'Linear Trend'){
+ x <- cbind(x, c(1:n))
+ colnames(x)[k+1] <- 't'
+ }
> x
Time Logins Views Shared Reviewed TotalSize Pop
1 210907 56 396 3 30 112285 1
2 120982 56 297 4 28 84786 1
3 176508 54 559 12 38 83123 1
4 179321 89 967 2 30 101193 1
5 123185 40 270 1 22 38361 1
6 52746 25 143 3 26 68504 1
7 385534 92 1562 0 25 119182 1
8 33170 18 109 0 18 22807 1
9 101645 63 371 0 11 17140 0
10 149061 44 656 5 26 116174 1
11 165446 33 511 0 25 57635 1
12 237213 84 655 0 38 66198 1
13 173326 88 465 7 44 71701 1
14 133131 55 525 7 30 57793 1
15 258873 60 885 3 40 80444 1
16 180083 66 497 9 34 53855 1
17 324799 154 1436 0 47 97668 1
18 230964 53 612 4 30 133824 1
19 236785 119 865 3 31 101481 1
20 135473 41 385 0 23 99645 1
21 202925 61 567 7 36 114789 1
22 215147 58 639 0 36 99052 1
23 344297 75 963 1 30 67654 1
24 153935 33 398 5 25 65553 1
25 132943 40 410 7 39 97500 1
26 174724 92 966 0 34 69112 1
27 174415 100 801 0 31 82753 1
28 225548 112 892 5 31 85323 1
29 223632 73 513 0 33 72654 1
30 124817 40 469 0 25 30727 1
31 221698 45 683 0 33 77873 1
32 210767 60 643 3 35 117478 1
33 170266 62 535 4 42 74007 1
34 260561 75 625 1 43 90183 1
35 84853 31 264 4 30 61542 1
36 294424 77 992 2 33 101494 1
37 101011 34 238 0 13 27570 1
38 215641 46 818 0 32 55813 1
39 325107 99 937 0 36 79215 1
40 7176 17 70 0 0 1423 1
41 167542 66 507 2 28 55461 1
42 106408 30 260 1 14 31081 1
43 96560 76 503 0 17 22996 1
44 265769 146 927 2 32 83122 1
45 269651 67 1269 10 30 70106 1
46 149112 56 537 6 35 60578 1
47 175824 107 910 0 20 39992 1
48 152871 58 532 5 28 79892 1
49 111665 34 345 4 28 49810 1
50 116408 61 918 1 39 71570 1
51 362301 119 1635 2 34 100708 1
52 78800 42 330 2 26 33032 1
53 183167 66 557 0 39 82875 1
54 277965 89 1178 8 39 139077 1
55 150629 44 740 3 33 71595 1
56 168809 66 452 0 28 72260 1
57 24188 24 218 0 4 5950 1
58 329267 259 764 8 39 115762 1
59 65029 17 255 5 18 32551 1
60 101097 64 454 3 14 31701 1
61 218946 41 866 1 29 80670 1
62 244052 68 574 5 44 143558 1
63 341570 168 1276 1 21 117105 1
64 103597 43 379 1 16 23789 1
65 233328 132 825 5 28 120733 1
66 256462 105 798 0 35 105195 1
67 206161 71 663 12 28 73107 1
68 311473 112 1069 8 38 132068 1
69 235800 94 921 8 23 149193 1
70 177939 82 858 8 36 46821 1
71 207176 70 711 8 32 87011 1
72 196553 57 503 2 29 95260 1
73 174184 53 382 0 25 55183 1
74 143246 103 464 5 27 106671 1
75 187559 121 717 8 36 73511 1
76 187681 62 690 2 28 92945 1
77 119016 52 462 5 23 78664 1
78 182192 52 657 12 40 70054 1
79 73566 32 385 6 23 22618 1
80 194979 62 577 7 40 74011 1
81 167488 45 619 2 28 83737 1
82 143756 46 479 0 34 69094 1
83 275541 63 817 4 33 93133 1
84 243199 75 752 3 28 95536 1
85 182999 88 430 6 34 225920 1
86 135649 46 451 2 30 62133 1
87 152299 53 537 0 33 61370 1
88 120221 37 519 1 22 43836 1
89 346485 90 1000 0 38 106117 1
90 145790 63 637 5 26 38692 1
91 193339 78 465 2 35 84651 1
92 80953 25 437 0 8 56622 1
93 122774 45 711 0 24 15986 1
94 130585 46 299 5 29 95364 1
95 112611 41 248 0 20 26706 0
96 286468 144 1162 1 29 89691 1
97 241066 82 714 0 45 67267 1
98 148446 91 905 1 37 126846 1
99 204713 71 649 1 33 41140 1
100 182079 63 512 2 33 102860 1
101 140344 53 472 6 25 51715 1
102 220516 62 905 1 32 55801 1
103 243060 63 786 4 29 111813 1
104 162765 32 489 2 28 120293 1
105 182613 39 479 3 28 138599 1
106 232138 62 617 0 31 161647 1
107 265318 117 925 10 52 115929 1
108 85574 34 351 0 21 24266 0
109 310839 92 1144 9 24 162901 1
110 225060 93 669 7 41 109825 1
111 232317 54 707 0 33 129838 1
112 144966 144 458 0 32 37510 1
113 43287 14 214 4 19 43750 1
114 155754 61 599 4 20 40652 1
115 164709 109 572 0 31 87771 1
116 201940 38 897 0 31 85872 1
117 235454 73 819 0 32 89275 1
118 220801 75 720 1 18 44418 0
119 99466 50 273 0 23 192565 1
120 92661 61 508 1 17 35232 0
121 133328 55 506 0 20 40909 0
122 61361 77 451 0 12 13294 0
123 125930 75 699 4 17 32387 0
124 100750 72 407 0 30 140867 1
125 224549 50 465 4 31 120662 1
126 82316 32 245 4 10 21233 0
127 102010 53 370 3 13 44332 0
128 101523 42 316 0 22 61056 0
129 243511 71 603 0 42 101338 1
130 22938 10 154 0 1 1168 1
131 41566 35 229 5 9 13497 0
132 152474 65 577 0 32 65567 1
133 61857 25 192 4 11 25162 1
134 99923 66 617 0 25 32334 0
135 132487 41 411 0 36 40735 1
136 317394 86 975 1 31 91413 1
137 21054 16 146 0 0 855 1
138 209641 42 705 5 24 97068 1
139 22648 19 184 0 13 44339 0
140 31414 19 200 0 8 14116 1
141 46698 45 274 0 13 10288 0
142 131698 65 502 0 19 65622 0
143 91735 35 382 0 18 16563 0
144 244749 95 964 2 33 76643 1
145 184510 49 537 7 40 110681 1
146 79863 37 438 1 22 29011 0
147 128423 64 369 8 38 92696 1
148 97839 38 417 2 24 94785 1
149 38214 34 276 0 8 8773 1
150 151101 32 514 2 35 83209 1
151 272458 65 822 0 43 93815 1
152 172494 52 389 0 43 86687 1
153 108043 62 466 1 14 34553 0
154 328107 65 1255 3 41 105547 1
155 250579 83 694 0 38 103487 1
156 351067 95 1024 3 45 213688 1
157 158015 29 400 0 31 71220 1
158 98866 18 397 0 13 23517 0
159 85439 33 350 0 28 56926 1
160 229242 247 719 4 31 91721 1
161 351619 139 1277 4 40 115168 1
162 84207 29 356 11 30 111194 1
163 120445 118 457 0 16 51009 0
164 324598 110 1402 0 37 135777 1
165 131069 67 600 4 30 51513 1
166 204271 42 480 0 35 74163 1
167 165543 65 595 1 32 51633 1
168 141722 94 436 0 27 75345 1
169 116048 64 230 0 20 33416 0
170 250047 81 651 0 18 83305 0
171 299775 95 1367 9 31 98952 1
172 195838 67 564 1 31 102372 1
173 173260 63 716 3 21 37238 1
174 254488 83 747 10 39 103772 1
175 104389 45 467 5 41 123969 1
176 136084 30 671 0 13 27142 0
177 199476 70 861 2 32 135400 1
178 92499 32 319 0 18 21399 0
179 224330 83 612 1 39 130115 1
180 135781 31 433 2 14 24874 0
181 74408 67 434 4 7 34988 0
182 81240 66 503 0 17 45549 0
183 14688 10 85 0 0 6023 1
184 181633 70 564 2 30 64466 1
185 271856 103 824 1 37 54990 1
186 7199 5 74 0 0 1644 1
187 46660 20 259 0 5 6179 1
188 17547 5 69 0 1 3926 1
189 133368 36 535 1 16 32755 0
190 95227 34 239 0 32 34777 1
191 152601 48 438 2 24 73224 1
192 98146 40 459 0 17 27114 0
193 79619 43 426 3 11 20760 0
194 59194 31 288 6 24 37636 0
195 139942 42 498 0 22 65461 0
196 118612 46 454 2 12 30080 0
197 72880 33 376 0 19 24094 0
198 65475 18 225 2 13 69008 0
199 99643 55 555 1 17 54968 0
200 71965 35 252 1 15 46090 0
201 77272 59 208 2 16 27507 0
202 49289 19 130 1 24 10672 0
203 135131 66 481 0 15 34029 0
204 108446 60 389 1 17 46300 0
205 89746 36 565 3 18 24760 0
206 44296 25 173 0 20 18779 0
207 77648 47 278 0 16 21280 0
208 181528 54 609 0 16 40662 0
209 134019 53 422 0 18 28987 0
210 124064 40 445 1 22 22827 0
211 92630 40 387 4 8 18513 0
212 121848 39 339 0 17 30594 0
213 52915 14 181 0 18 24006 0
214 81872 45 245 0 16 27913 0
215 58981 36 384 7 23 42744 0
216 53515 28 212 2 22 12934 0
217 60812 44 399 0 13 22574 0
218 56375 30 229 7 13 41385 0
219 65490 22 224 3 16 18653 0
220 80949 17 203 0 16 18472 0
221 76302 31 333 0 20 30976 0
222 104011 55 384 6 22 63339 0
223 98104 54 636 2 17 25568 0
224 67989 21 185 0 18 33747 0
225 30989 14 93 0 17 4154 0
226 135458 81 581 3 12 19474 0
227 73504 35 248 0 7 35130 0
228 63123 43 304 1 17 39067 0
229 61254 46 344 1 14 13310 0
230 74914 30 407 0 23 65892 0
231 31774 23 170 1 17 4143 0
232 81437 38 312 0 14 28579 0
233 87186 54 507 0 15 51776 0
234 50090 20 224 0 17 21152 0
235 65745 53 340 0 21 38084 0
236 56653 45 168 0 18 27717 0
237 158399 39 443 0 18 32928 0
238 46455 20 204 0 17 11342 0
239 73624 24 367 0 17 19499 0
240 38395 31 210 0 16 16380 0
241 91899 35 335 0 15 36874 0
242 139526 151 364 0 21 48259 0
243 52164 52 178 0 16 16734 0
244 51567 30 206 2 14 28207 0
245 70551 31 279 0 15 30143 0
246 84856 29 387 1 17 41369 0
247 102538 57 490 1 15 45833 0
248 86678 40 238 0 15 29156 0
249 85709 44 343 0 10 35944 0
250 34662 25 232 0 6 36278 0
251 150580 77 530 0 22 45588 0
252 99611 35 291 0 21 45097 0
253 19349 11 67 0 1 3895 0
254 99373 63 397 1 18 28394 0
255 86230 44 467 0 17 18632 0
256 30837 19 178 0 4 2325 0
257 31706 13 175 0 10 25139 0
258 89806 42 299 0 16 27975 0
259 62088 38 154 1 16 14483 0
260 40151 29 106 0 9 13127 0
261 27634 20 189 0 16 5839 0
262 76990 27 194 0 17 24069 0
263 37460 20 135 0 7 3738 0
264 54157 19 201 0 15 18625 0
265 49862 37 207 0 14 36341 0
266 84337 26 280 0 14 24548 0
267 64175 42 260 0 18 21792 0
268 59382 49 227 0 12 26263 0
269 119308 30 239 0 16 23686 0
270 76702 49 333 0 21 49303 0
271 103425 67 428 1 19 25659 0
272 70344 28 230 0 16 28904 0
273 43410 19 292 0 1 2781 0
274 104838 49 350 1 16 29236 0
275 62215 27 186 0 10 19546 0
276 69304 30 326 6 19 22818 0
277 53117 22 155 3 12 32689 0
278 19764 12 75 1 2 5752 0
279 86680 31 361 2 14 22197 0
280 84105 20 261 0 17 20055 0
281 77945 20 299 0 19 25272 0
282 89113 39 300 0 14 82206 0
283 91005 29 450 3 11 32073 0
284 40248 16 183 1 4 5444 0
285 64187 27 238 0 16 20154 0
286 50857 21 165 0 20 36944 0
287 56613 19 234 1 12 8019 0
288 62792 35 176 0 15 30884 0
289 72535 14 329 0 16 19540 0
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) Logins Views Shared Reviewed TotalSize
-1.516e+04 2.831e+02 1.639e+02 -1.655e+03 1.499e+03 3.693e-01
Pop
6.459e+03
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-125849 -15485 625 15770 105643
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -1.516e+04 4.286e+03 -3.538 0.000472 ***
Logins 2.831e+02 7.085e+01 3.995 8.26e-05 ***
Views 1.639e+02 9.291e+00 17.639 < 2e-16 ***
Shared -1.655e+03 6.833e+02 -2.422 0.016087 *
Reviewed 1.499e+03 2.672e+02 5.610 4.82e-08 ***
TotalSize 3.693e-01 6.740e-02 5.479 9.50e-08 ***
Pop 6.459e+03 4.682e+03 1.380 0.168783
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 28650 on 282 degrees of freedom
Multiple R-squared: 0.8814, Adjusted R-squared: 0.8789
F-statistic: 349.4 on 6 and 282 DF, p-value: < 2.2e-16
> if (n > n25) {
+ kp3 <- k + 3
+ nmkm3 <- n - k - 3
+ gqarr <- array(NA, dim=c(nmkm3-kp3+1,3))
+ numgqtests <- 0
+ numsignificant1 <- 0
+ numsignificant5 <- 0
+ numsignificant10 <- 0
+ for (mypoint in kp3:nmkm3) {
+ j <- 0
+ numgqtests <- numgqtests + 1
+ for (myalt in c('greater', 'two.sided', 'less')) {
+ j <- j + 1
+ gqarr[mypoint-kp3+1,j] <- gqtest(mylm, point=mypoint, alternative=myalt)$p.value
+ }
+ if (gqarr[mypoint-kp3+1,2] < 0.01) numsignificant1 <- numsignificant1 + 1
+ if (gqarr[mypoint-kp3+1,2] < 0.05) numsignificant5 <- numsignificant5 + 1
+ if (gqarr[mypoint-kp3+1,2] < 0.10) numsignificant10 <- numsignificant10 + 1
+ }
+ gqarr
+ }
[,1] [,2] [,3]
[1,] 0.9987200 2.559993e-03 1.279996e-03
[2,] 0.9966509 6.698286e-03 3.349143e-03
[3,] 0.9926542 1.469167e-02 7.345835e-03
[4,] 0.9860749 2.785019e-02 1.392510e-02
[5,] 0.9744625 5.107507e-02 2.553753e-02
[6,] 0.9608103 7.837947e-02 3.918973e-02
[7,] 0.9660898 6.782035e-02 3.391017e-02
[8,] 0.9793712 4.125769e-02 2.062884e-02
[9,] 0.9760960 4.780807e-02 2.390403e-02
[10,] 0.9671504 6.569916e-02 3.284958e-02
[11,] 0.9525469 9.490622e-02 4.745311e-02
[12,] 0.9321930 1.356140e-01 6.780701e-02
[13,] 0.9058463 1.883073e-01 9.415366e-02
[14,] 0.9967806 6.438794e-03 3.219397e-03
[15,] 0.9957175 8.564955e-03 4.282477e-03
[16,] 0.9957841 8.431722e-03 4.215861e-03
[17,] 0.9995608 8.783659e-04 4.391829e-04
[18,] 0.9996493 7.014486e-04 3.507243e-04
[19,] 0.9994300 1.140077e-03 5.700387e-04
[20,] 0.9997594 4.812001e-04 2.406001e-04
[21,] 0.9996388 7.223503e-04 3.611752e-04
[22,] 0.9994646 1.070869e-03 5.354346e-04
[23,] 0.9991516 1.696886e-03 8.484430e-04
[24,] 0.9988258 2.348395e-03 1.174198e-03
[25,] 0.9991889 1.622177e-03 8.110887e-04
[26,] 0.9991647 1.670667e-03 8.353334e-04
[27,] 0.9990961 1.807867e-03 9.039336e-04
[28,] 0.9990364 1.927238e-03 9.636188e-04
[29,] 0.9986219 2.756294e-03 1.378147e-03
[30,] 0.9997004 5.992457e-04 2.996229e-04
[31,] 0.9995379 9.241288e-04 4.620644e-04
[32,] 0.9993460 1.307967e-03 6.539833e-04
[33,] 0.9992664 1.467168e-03 7.335838e-04
[34,] 0.9992002 1.599604e-03 7.998021e-04
[35,] 0.9990570 1.885931e-03 9.429654e-04
[36,] 0.9986275 2.744955e-03 1.372477e-03
[37,] 0.9981775 3.645002e-03 1.822501e-03
[38,] 0.9981722 3.655576e-03 1.827788e-03
[39,] 0.9974453 5.109464e-03 2.554732e-03
[40,] 0.9964984 7.003108e-03 3.501554e-03
[41,] 0.9999959 8.253026e-06 4.126513e-06
[42,] 0.9999938 1.242476e-05 6.212381e-06
[43,] 0.9999928 1.439800e-05 7.199002e-06
[44,] 0.9999890 2.199257e-05 1.099629e-05
[45,] 0.9999882 2.354674e-05 1.177337e-05
[46,] 0.9999929 1.418736e-05 7.093680e-06
[47,] 0.9999895 2.090773e-05 1.045386e-05
[48,] 0.9999861 2.778799e-05 1.389399e-05
[49,] 0.9999914 1.722657e-05 8.613283e-06
[50,] 0.9999865 2.702709e-05 1.351354e-05
[51,] 0.9999799 4.019701e-05 2.009851e-05
[52,] 0.9999700 6.004433e-05 3.002216e-05
[53,] 0.9999627 7.462685e-05 3.731343e-05
[54,] 0.9999479 1.041143e-04 5.205717e-05
[55,] 0.9999247 1.506990e-04 7.534948e-05
[56,] 0.9999066 1.868516e-04 9.342581e-05
[57,] 0.9998674 2.651745e-04 1.325873e-04
[58,] 0.9998947 2.105412e-04 1.052706e-04
[59,] 0.9998636 2.728702e-04 1.364351e-04
[60,] 0.9998356 3.288187e-04 1.644094e-04
[61,] 0.9998299 3.402187e-04 1.701093e-04
[62,] 0.9997699 4.601703e-04 2.300852e-04
[63,] 0.9997509 4.982729e-04 2.491365e-04
[64,] 0.9998402 3.196801e-04 1.598400e-04
[65,] 0.9998644 2.712980e-04 1.356490e-04
[66,] 0.9998387 3.225698e-04 1.612849e-04
[67,] 0.9997793 4.413119e-04 2.206559e-04
[68,] 0.9997357 5.286151e-04 2.643075e-04
[69,] 0.9996372 7.256995e-04 3.628497e-04
[70,] 0.9995558 8.884492e-04 4.442246e-04
[71,] 0.9994421 1.115735e-03 5.578674e-04
[72,] 0.9992543 1.491439e-03 7.457197e-04
[73,] 0.9990841 1.831808e-03 9.159040e-04
[74,] 0.9995757 8.486682e-04 4.243341e-04
[75,] 0.9996069 7.862637e-04 3.931319e-04
[76,] 0.9997600 4.800976e-04 2.400488e-04
[77,] 0.9996721 6.558333e-04 3.279167e-04
[78,] 0.9995804 8.391973e-04 4.195987e-04
[79,] 0.9994594 1.081143e-03 5.405713e-04
[80,] 0.9998708 2.583371e-04 1.291685e-04
[81,] 0.9998258 3.484854e-04 1.742427e-04
[82,] 0.9997987 4.025035e-04 2.012517e-04
[83,] 0.9997640 4.719564e-04 2.359782e-04
[84,] 0.9998155 3.690959e-04 1.845480e-04
[85,] 0.9997507 4.985642e-04 2.492821e-04
[86,] 0.9997374 5.251677e-04 2.625838e-04
[87,] 0.9996536 6.927747e-04 3.463873e-04
[88,] 0.9995644 8.712606e-04 4.356303e-04
[89,] 0.9999988 2.359031e-06 1.179516e-06
[90,] 0.9999986 2.752758e-06 1.376379e-06
[91,] 0.9999979 4.110388e-06 2.055194e-06
[92,] 0.9999971 5.853512e-06 2.926756e-06
[93,] 0.9999957 8.546289e-06 4.273144e-06
[94,] 0.9999955 9.072055e-06 4.536028e-06
[95,] 0.9999933 1.346397e-05 6.731987e-06
[96,] 0.9999911 1.787574e-05 8.937869e-06
[97,] 0.9999886 2.278848e-05 1.139424e-05
[98,] 0.9999848 3.042648e-05 1.521324e-05
[99,] 0.9999808 3.832182e-05 1.916091e-05
[100,] 0.9999797 4.064590e-05 2.032295e-05
[101,] 0.9999724 5.517753e-05 2.758877e-05
[102,] 0.9999633 7.347224e-05 3.673612e-05
[103,] 0.9999592 8.165430e-05 4.082715e-05
[104,] 0.9999542 9.166666e-05 4.583333e-05
[105,] 0.9999382 1.236688e-04 6.183441e-05
[106,] 0.9999400 1.200711e-04 6.003557e-05
[107,] 0.9999388 1.224715e-04 6.123577e-05
[108,] 0.9999152 1.696650e-04 8.483248e-05
[109,] 0.9999510 9.797229e-05 4.898615e-05
[110,] 0.9999809 3.814179e-05 1.907089e-05
[111,] 0.9999861 2.786399e-05 1.393200e-05
[112,] 0.9999801 3.989023e-05 1.994512e-05
[113,] 0.9999892 2.160975e-05 1.080488e-05
[114,] 0.9999889 2.224236e-05 1.112118e-05
[115,] 0.9999989 2.181898e-06 1.090949e-06
[116,] 0.9999998 4.565769e-07 2.282885e-07
[117,] 0.9999998 3.787828e-07 1.893914e-07
[118,] 0.9999997 5.511937e-07 2.755969e-07
[119,] 0.9999996 8.456919e-07 4.228460e-07
[120,] 0.9999996 7.054938e-07 3.527469e-07
[121,] 0.9999995 1.081331e-06 5.406656e-07
[122,] 0.9999992 1.640255e-06 8.201274e-07
[123,] 0.9999991 1.784986e-06 8.924931e-07
[124,] 0.9999988 2.397076e-06 1.198538e-06
[125,] 0.9999997 5.738172e-07 2.869086e-07
[126,] 0.9999996 8.610743e-07 4.305372e-07
[127,] 0.9999999 1.484519e-07 7.422594e-08
[128,] 0.9999999 2.359066e-07 1.179533e-07
[129,] 0.9999999 2.023076e-07 1.011538e-07
[130,] 0.9999999 1.677164e-07 8.385821e-08
[131,] 0.9999999 2.296261e-07 1.148131e-07
[132,] 0.9999999 2.846984e-07 1.423492e-07
[133,] 0.9999998 4.377504e-07 2.188752e-07
[134,] 0.9999997 6.800348e-07 3.400174e-07
[135,] 0.9999995 1.024610e-06 5.123049e-07
[136,] 0.9999993 1.406416e-06 7.032081e-07
[137,] 0.9999994 1.299024e-06 6.495120e-07
[138,] 0.9999991 1.797254e-06 8.986268e-07
[139,] 0.9999994 1.248047e-06 6.240236e-07
[140,] 0.9999994 1.259733e-06 6.298666e-07
[141,] 0.9999991 1.785923e-06 8.929615e-07
[142,] 0.9999991 1.789120e-06 8.945602e-07
[143,] 0.9999987 2.640374e-06 1.320187e-06
[144,] 0.9999980 3.988206e-06 1.994103e-06
[145,] 0.9999974 5.127955e-06 2.563978e-06
[146,] 0.9999975 4.937430e-06 2.468715e-06
[147,] 0.9999984 3.145984e-06 1.572992e-06
[148,] 0.9999983 3.368484e-06 1.684242e-06
[149,] 0.9999978 4.493998e-06 2.246999e-06
[150,] 0.9999982 3.537131e-06 1.768565e-06
[151,] 0.9999983 3.479915e-06 1.739958e-06
[152,] 0.9999979 4.225194e-06 2.112597e-06
[153,] 0.9999980 4.004088e-06 2.002044e-06
[154,] 0.9999977 4.646530e-06 2.323265e-06
[155,] 0.9999979 4.255864e-06 2.127932e-06
[156,] 0.9999985 3.005035e-06 1.502517e-06
[157,] 0.9999994 1.143740e-06 5.718702e-07
[158,] 0.9999991 1.722931e-06 8.614657e-07
[159,] 0.9999990 2.016586e-06 1.008293e-06
[160,] 0.9999992 1.697484e-06 8.487420e-07
[161,] 1.0000000 2.077579e-08 1.038790e-08
[162,] 1.0000000 3.381893e-08 1.690946e-08
[163,] 1.0000000 4.266477e-08 2.133238e-08
[164,] 1.0000000 6.958600e-08 3.479300e-08
[165,] 1.0000000 1.147686e-08 5.738430e-09
[166,] 1.0000000 1.319759e-09 6.598794e-10
[167,] 1.0000000 2.311449e-09 1.155724e-09
[168,] 1.0000000 1.057180e-09 5.285900e-10
[169,] 1.0000000 1.699891e-09 8.499454e-10
[170,] 1.0000000 2.760897e-09 1.380448e-09
[171,] 1.0000000 6.575299e-10 3.287649e-10
[172,] 1.0000000 8.719129e-10 4.359564e-10
[173,] 1.0000000 1.434277e-10 7.171387e-11
[174,] 1.0000000 2.254112e-10 1.127056e-10
[175,] 1.0000000 3.779802e-10 1.889901e-10
[176,] 1.0000000 7.167197e-11 3.583598e-11
[177,] 1.0000000 1.085619e-10 5.428096e-11
[178,] 1.0000000 1.386072e-10 6.930359e-11
[179,] 1.0000000 1.704848e-10 8.524239e-11
[180,] 1.0000000 2.011115e-10 1.005557e-10
[181,] 1.0000000 2.889156e-10 1.444578e-10
[182,] 1.0000000 5.234176e-10 2.617088e-10
[183,] 1.0000000 9.348223e-10 4.674111e-10
[184,] 1.0000000 1.624339e-09 8.121695e-10
[185,] 1.0000000 2.632402e-09 1.316201e-09
[186,] 1.0000000 3.477874e-09 1.738937e-09
[187,] 1.0000000 3.987585e-09 1.993793e-09
[188,] 1.0000000 5.321290e-09 2.660645e-09
[189,] 1.0000000 9.474837e-09 4.737418e-09
[190,] 1.0000000 7.070885e-09 3.535443e-09
[191,] 1.0000000 1.295911e-08 6.479553e-09
[192,] 1.0000000 2.161659e-08 1.080830e-08
[193,] 1.0000000 3.849009e-08 1.924505e-08
[194,] 1.0000000 5.095810e-08 2.547905e-08
[195,] 1.0000000 8.845449e-08 4.422724e-08
[196,] 1.0000000 8.736445e-08 4.368223e-08
[197,] 0.9999999 1.344234e-07 6.721169e-08
[198,] 0.9999999 2.372196e-07 1.186098e-07
[199,] 1.0000000 3.573318e-08 1.786659e-08
[200,] 1.0000000 2.144854e-08 1.072427e-08
[201,] 1.0000000 2.039669e-08 1.019835e-08
[202,] 1.0000000 2.444289e-08 1.222145e-08
[203,] 1.0000000 7.845522e-09 3.922761e-09
[204,] 1.0000000 1.504052e-08 7.520259e-09
[205,] 1.0000000 2.573093e-08 1.286547e-08
[206,] 1.0000000 1.807702e-08 9.038512e-09
[207,] 1.0000000 3.116808e-08 1.558404e-08
[208,] 1.0000000 2.308020e-08 1.154010e-08
[209,] 1.0000000 4.400686e-08 2.200343e-08
[210,] 1.0000000 7.511171e-08 3.755585e-08
[211,] 1.0000000 5.667909e-08 2.833954e-08
[212,] 0.9999999 1.019756e-07 5.098781e-08
[213,] 0.9999999 1.877732e-07 9.388661e-08
[214,] 0.9999999 1.089234e-07 5.446168e-08
[215,] 0.9999999 1.750331e-07 8.751657e-08
[216,] 0.9999998 3.264709e-07 1.632354e-07
[217,] 0.9999997 5.445549e-07 2.722775e-07
[218,] 0.9999996 8.436940e-07 4.218470e-07
[219,] 0.9999995 1.052287e-06 5.261435e-07
[220,] 0.9999994 1.171083e-06 5.855415e-07
[221,] 0.9999997 6.780617e-07 3.390308e-07
[222,] 0.9999996 8.456853e-07 4.228427e-07
[223,] 0.9999992 1.582200e-06 7.910998e-07
[224,] 0.9999996 7.410761e-07 3.705381e-07
[225,] 0.9999994 1.197626e-06 5.988128e-07
[226,] 0.9999998 4.735410e-07 2.367705e-07
[227,] 0.9999996 8.613557e-07 4.306779e-07
[228,] 1.0000000 1.214814e-08 6.074068e-09
[229,] 1.0000000 2.165686e-08 1.082843e-08
[230,] 1.0000000 4.038017e-08 2.019009e-08
[231,] 1.0000000 2.592709e-08 1.296355e-08
[232,] 1.0000000 5.159209e-08 2.579604e-08
[233,] 0.9999999 1.123592e-07 5.617959e-08
[234,] 0.9999999 1.608400e-07 8.042002e-08
[235,] 0.9999998 3.083562e-07 1.541781e-07
[236,] 0.9999997 6.602484e-07 3.301242e-07
[237,] 0.9999993 1.336971e-06 6.684857e-07
[238,] 0.9999988 2.356361e-06 1.178181e-06
[239,] 0.9999983 3.360170e-06 1.680085e-06
[240,] 0.9999966 6.835874e-06 3.417937e-06
[241,] 0.9999967 6.597967e-06 3.298983e-06
[242,] 0.9999958 8.346788e-06 4.173394e-06
[243,] 0.9999932 1.354031e-05 6.770153e-06
[244,] 0.9999860 2.798781e-05 1.399390e-05
[245,] 0.9999710 5.805939e-05 2.902969e-05
[246,] 0.9999647 7.066222e-05 3.533111e-05
[247,] 0.9999357 1.285493e-04 6.427464e-05
[248,] 0.9999208 1.584448e-04 7.922238e-05
[249,] 0.9998579 2.841229e-04 1.420614e-04
[250,] 0.9997462 5.075406e-04 2.537703e-04
[251,] 0.9995041 9.917396e-04 4.958698e-04
[252,] 0.9997655 4.690350e-04 2.345175e-04
[253,] 0.9996334 7.331755e-04 3.665878e-04
[254,] 0.9992609 1.478150e-03 7.390750e-04
[255,] 0.9986414 2.717262e-03 1.358631e-03
[256,] 0.9981870 3.626054e-03 1.813027e-03
[257,] 0.9969997 6.000669e-03 3.000335e-03
[258,] 0.9957145 8.571022e-03 4.285511e-03
[259,] 0.9931683 1.366334e-02 6.831671e-03
[260,] 0.9999505 9.903193e-05 4.951597e-05
[261,] 0.9999839 3.229061e-05 1.614531e-05
[262,] 0.9999861 2.778975e-05 1.389488e-05
[263,] 0.9999459 1.082352e-04 5.411761e-05
[264,] 0.9999552 8.967855e-05 4.483928e-05
[265,] 0.9998405 3.190098e-04 1.595049e-04
[266,] 0.9994868 1.026442e-03 5.132211e-04
[267,] 0.9985044 2.991207e-03 1.495603e-03
[268,] 0.9968047 6.390590e-03 3.195295e-03
[269,] 0.9876807 2.463865e-02 1.231932e-02
[270,] 0.9553589 8.928219e-02 4.464109e-02
> postscript(file="/var/wessaorg/rcomp/tmp/16j6x1355218833.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/wessaorg/rcomp/tmp/2y2ko1355218833.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/wessaorg/rcomp/tmp/3joe61355218833.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/wessaorg/rcomp/tmp/4g21p1355218833.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/wessaorg/rcomp/tmp/5son71355218833.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 = 289
Frequency = 1
1 2 3 4 5
57381.11794 -1511.90811 10501.17082 -74681.58651 30819.96881
6 7 8 9 10
-28378.61573 30717.59077 -16496.24686 15351.12634 -35806.27187
11 12 13 14 15
22297.04832 33374.39258 47.29319 -14512.12552 20840.28937
16 17 18 19 20
32680.77750 -51961.21467 36590.20525 -8944.82804 -1805.12382
21 22 23 24 25
16657.91439 12157.59221 105642.71844 34654.50653 -19764.26410
26 27 28 29 30
-77424.81227 -53496.76887 -13349.23944 51291.41395 -3494.46935
31 32 33 34 35
27496.39450 6215.55758 -9941.44218 49487.62786 -19571.80669
36 37 38 39 40
35111.90998 31413.06822 8677.49503 68999.86930 -930.54728
41 42 43 44 45
15321.23850 33195.82855 -32664.16493 5860.64248 -2899.41017
46 47 48 49 50
-10960.83946 -39649.27269 -5239.37136 447.46891 -125849.14401
51 52 53 54 55
-15485.14617 -26339.15190 -7173.40119 -28171.72815 -45348.42941
56 57 58 59 60
16089.14399 -17823.72420 51464.29413 -3605.32157 -10452.14157
61 62 63 64 65
2505.82322 28728.51718 20529.52422 6896.81673 -8827.80248
66 67 68 69 70
13342.09226 36990.65941 20776.23197 -9379.63009 -35209.79356
71 72 73 74 75
12671.81419 31339.97577 47419.76163 -24847.57652 -23377.73801
76 77 78 79 80
-7239.02624 -17973.81838 2517.94397 -22793.42358 15850.00152
81 82 83 84 85
-7584.09708 -15554.51361 55267.99624 35136.35832 -28151.02707
86 87 88 89 90
-7195.42067 -14146.61998 -14123.02164 69668.44704 -12731.14567
91 92 93 94 95
23330.06097 -21940.41106 -39669.57886 6843.37884 35675.30501
96 97 98 99 100
-10968.58808 17234.03830 -117585.07544 23942.12876 4888.69166
101 102 103 104 105
10039.32580 -3575.20141 26967.18332 -821.09625 13579.32886
106 107 108 109 110
16005.04785 -14917.78955 -6858.40724 24773.46998 7354.72167
111 112 113 114 115
12446.17238 -23982.64205 -25068.37781 10644.38922 -30074.45043
116 117 118 119 120
-25305.67103 8327.61408 55001.97162 -56316.34712 -29540.73584
121 122 123 124 125
4903.28961 -42085.63943 -25521.43149 -74625.98305 58475.78260
126 127 128 129 130
32053.34571 10635.28023 -2522.46536 32900.92429 1641.66583
131 132 133 134 135
-913.08066 -23974.12311 12852.10903 -54136.94656 -6791.24183
136 137 138 139 140
63386.79098 985.25829 27363.47235 -33587.04807 -15245.41479
141 142 143 144 145
-19072.70354 -6527.94620 1281.75086 -5893.30000 2075.56041
146 147 148 149 150
-29273.11738 -19430.81573 -40229.48323 -23173.55634 -13384.37261
151 152 153 154 155
28934.80243 6243.54112 -2809.87691 17253.88781 26862.87247
156 157 158 159 160
23649.04916 20176.49159 15695.89678 -35561.38334 -23535.62236
161 162 163 164 165
15814.68234 -41480.04988 -15514.94985 -33213.15976 -34907.34146
166 167 168 169 170
42557.73770 -7055.04988 -15941.52180 33073.81360 77843.37210
171 172 173 174 175
-10567.50642 10517.82044 6516.62322 37028.73077 -75156.95534
176 177 178 179 180
3274.91471 -47406.71555 11433.78388 4374.93135 44340.19551
181 182 183 184 185
-17316.73907 -47022.02775 4406.14997 12614.75872 42235.77627
186 187 188 189 190
1752.18093 -2522.53629 10577.55326 16232.29357 -5684.49964
191 192 193 194 195
16221.82643 -8738.12690 -6399.47350 -21571.45769 4443.50282
196 197 198 199 200
20560.04069 -20304.07951 -2995.87050 -35850.38990 -1933.77753
201 202 203 204 205
10807.51736 -503.27315 17726.98844 1942.31810 -29043.29067
206 207 208 209 210
-12892.32868 2098.58134 42595.32910 27327.55140 15214.97874
211 212 213 214 215
20834.74317 33627.72752 -1401.91244 9847.53970 -37664.50139
216 217 218 219 220
-8444.96432 -29697.52180 2325.32435 11801.18346 27219.72710
221 222 223 224 225
-13309.64206 -5773.03200 -37869.86553 7438.13087 -76.45526
226 227 228 229 230
12256.62081 14648.21807 -21967.69165 -17232.05350 -43932.95034
231 232 233 234 235
-12799.20274 3166.83219 -37635.26986 -10419.00216 -35365.45895
236 237 238 239 240
-5678.98516 50773.80086 -7153.96261 -10841.98093 -19671.96340
241 242 243 244 245
6146.39423 2984.16950 -6733.97382 -3620.33496 -2406.48171
246 247 248 249 250
-10723.89172 -16497.28395 18256.47603 3938.06610 -17665.47305
251 252 253 254 255
7268.01690 9036.43523 17479.59684 -4178.43422 -19965.43252
256 257 258 259 260
4593.71816 -9767.78007 9758.25289 13572.10719 11390.80740
261 262 263 264 265
-19985.32693 18338.82340 12960.62004 1632.24302 -13783.69448
266 267 268 269 270
16196.37801 -10197.03482 -4218.03514 54073.75254 -26271.69019
271 272 273 274 275
-6828.52716 5224.04983 2814.00833 15639.87683 17040.77017
276 277 278 279 280
-4437.28245 11550.69999 15770.09184 8027.17270 17937.45199
281 282 283 284 285
624.75388 -7275.78206 838.97164 14537.53539 1270.06348
286 287 288 289
-10595.53443 8749.38291 5308.42024 -1388.50917
> postscript(file="/var/wessaorg/rcomp/tmp/66d1e1355218833.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 = 289
Frequency = 1
lag(myerror, k = 1) myerror
0 57381.11794 NA
1 -1511.90811 57381.11794
2 10501.17082 -1511.90811
3 -74681.58651 10501.17082
4 30819.96881 -74681.58651
5 -28378.61573 30819.96881
6 30717.59077 -28378.61573
7 -16496.24686 30717.59077
8 15351.12634 -16496.24686
9 -35806.27187 15351.12634
10 22297.04832 -35806.27187
11 33374.39258 22297.04832
12 47.29319 33374.39258
13 -14512.12552 47.29319
14 20840.28937 -14512.12552
15 32680.77750 20840.28937
16 -51961.21467 32680.77750
17 36590.20525 -51961.21467
18 -8944.82804 36590.20525
19 -1805.12382 -8944.82804
20 16657.91439 -1805.12382
21 12157.59221 16657.91439
22 105642.71844 12157.59221
23 34654.50653 105642.71844
24 -19764.26410 34654.50653
25 -77424.81227 -19764.26410
26 -53496.76887 -77424.81227
27 -13349.23944 -53496.76887
28 51291.41395 -13349.23944
29 -3494.46935 51291.41395
30 27496.39450 -3494.46935
31 6215.55758 27496.39450
32 -9941.44218 6215.55758
33 49487.62786 -9941.44218
34 -19571.80669 49487.62786
35 35111.90998 -19571.80669
36 31413.06822 35111.90998
37 8677.49503 31413.06822
38 68999.86930 8677.49503
39 -930.54728 68999.86930
40 15321.23850 -930.54728
41 33195.82855 15321.23850
42 -32664.16493 33195.82855
43 5860.64248 -32664.16493
44 -2899.41017 5860.64248
45 -10960.83946 -2899.41017
46 -39649.27269 -10960.83946
47 -5239.37136 -39649.27269
48 447.46891 -5239.37136
49 -125849.14401 447.46891
50 -15485.14617 -125849.14401
51 -26339.15190 -15485.14617
52 -7173.40119 -26339.15190
53 -28171.72815 -7173.40119
54 -45348.42941 -28171.72815
55 16089.14399 -45348.42941
56 -17823.72420 16089.14399
57 51464.29413 -17823.72420
58 -3605.32157 51464.29413
59 -10452.14157 -3605.32157
60 2505.82322 -10452.14157
61 28728.51718 2505.82322
62 20529.52422 28728.51718
63 6896.81673 20529.52422
64 -8827.80248 6896.81673
65 13342.09226 -8827.80248
66 36990.65941 13342.09226
67 20776.23197 36990.65941
68 -9379.63009 20776.23197
69 -35209.79356 -9379.63009
70 12671.81419 -35209.79356
71 31339.97577 12671.81419
72 47419.76163 31339.97577
73 -24847.57652 47419.76163
74 -23377.73801 -24847.57652
75 -7239.02624 -23377.73801
76 -17973.81838 -7239.02624
77 2517.94397 -17973.81838
78 -22793.42358 2517.94397
79 15850.00152 -22793.42358
80 -7584.09708 15850.00152
81 -15554.51361 -7584.09708
82 55267.99624 -15554.51361
83 35136.35832 55267.99624
84 -28151.02707 35136.35832
85 -7195.42067 -28151.02707
86 -14146.61998 -7195.42067
87 -14123.02164 -14146.61998
88 69668.44704 -14123.02164
89 -12731.14567 69668.44704
90 23330.06097 -12731.14567
91 -21940.41106 23330.06097
92 -39669.57886 -21940.41106
93 6843.37884 -39669.57886
94 35675.30501 6843.37884
95 -10968.58808 35675.30501
96 17234.03830 -10968.58808
97 -117585.07544 17234.03830
98 23942.12876 -117585.07544
99 4888.69166 23942.12876
100 10039.32580 4888.69166
101 -3575.20141 10039.32580
102 26967.18332 -3575.20141
103 -821.09625 26967.18332
104 13579.32886 -821.09625
105 16005.04785 13579.32886
106 -14917.78955 16005.04785
107 -6858.40724 -14917.78955
108 24773.46998 -6858.40724
109 7354.72167 24773.46998
110 12446.17238 7354.72167
111 -23982.64205 12446.17238
112 -25068.37781 -23982.64205
113 10644.38922 -25068.37781
114 -30074.45043 10644.38922
115 -25305.67103 -30074.45043
116 8327.61408 -25305.67103
117 55001.97162 8327.61408
118 -56316.34712 55001.97162
119 -29540.73584 -56316.34712
120 4903.28961 -29540.73584
121 -42085.63943 4903.28961
122 -25521.43149 -42085.63943
123 -74625.98305 -25521.43149
124 58475.78260 -74625.98305
125 32053.34571 58475.78260
126 10635.28023 32053.34571
127 -2522.46536 10635.28023
128 32900.92429 -2522.46536
129 1641.66583 32900.92429
130 -913.08066 1641.66583
131 -23974.12311 -913.08066
132 12852.10903 -23974.12311
133 -54136.94656 12852.10903
134 -6791.24183 -54136.94656
135 63386.79098 -6791.24183
136 985.25829 63386.79098
137 27363.47235 985.25829
138 -33587.04807 27363.47235
139 -15245.41479 -33587.04807
140 -19072.70354 -15245.41479
141 -6527.94620 -19072.70354
142 1281.75086 -6527.94620
143 -5893.30000 1281.75086
144 2075.56041 -5893.30000
145 -29273.11738 2075.56041
146 -19430.81573 -29273.11738
147 -40229.48323 -19430.81573
148 -23173.55634 -40229.48323
149 -13384.37261 -23173.55634
150 28934.80243 -13384.37261
151 6243.54112 28934.80243
152 -2809.87691 6243.54112
153 17253.88781 -2809.87691
154 26862.87247 17253.88781
155 23649.04916 26862.87247
156 20176.49159 23649.04916
157 15695.89678 20176.49159
158 -35561.38334 15695.89678
159 -23535.62236 -35561.38334
160 15814.68234 -23535.62236
161 -41480.04988 15814.68234
162 -15514.94985 -41480.04988
163 -33213.15976 -15514.94985
164 -34907.34146 -33213.15976
165 42557.73770 -34907.34146
166 -7055.04988 42557.73770
167 -15941.52180 -7055.04988
168 33073.81360 -15941.52180
169 77843.37210 33073.81360
170 -10567.50642 77843.37210
171 10517.82044 -10567.50642
172 6516.62322 10517.82044
173 37028.73077 6516.62322
174 -75156.95534 37028.73077
175 3274.91471 -75156.95534
176 -47406.71555 3274.91471
177 11433.78388 -47406.71555
178 4374.93135 11433.78388
179 44340.19551 4374.93135
180 -17316.73907 44340.19551
181 -47022.02775 -17316.73907
182 4406.14997 -47022.02775
183 12614.75872 4406.14997
184 42235.77627 12614.75872
185 1752.18093 42235.77627
186 -2522.53629 1752.18093
187 10577.55326 -2522.53629
188 16232.29357 10577.55326
189 -5684.49964 16232.29357
190 16221.82643 -5684.49964
191 -8738.12690 16221.82643
192 -6399.47350 -8738.12690
193 -21571.45769 -6399.47350
194 4443.50282 -21571.45769
195 20560.04069 4443.50282
196 -20304.07951 20560.04069
197 -2995.87050 -20304.07951
198 -35850.38990 -2995.87050
199 -1933.77753 -35850.38990
200 10807.51736 -1933.77753
201 -503.27315 10807.51736
202 17726.98844 -503.27315
203 1942.31810 17726.98844
204 -29043.29067 1942.31810
205 -12892.32868 -29043.29067
206 2098.58134 -12892.32868
207 42595.32910 2098.58134
208 27327.55140 42595.32910
209 15214.97874 27327.55140
210 20834.74317 15214.97874
211 33627.72752 20834.74317
212 -1401.91244 33627.72752
213 9847.53970 -1401.91244
214 -37664.50139 9847.53970
215 -8444.96432 -37664.50139
216 -29697.52180 -8444.96432
217 2325.32435 -29697.52180
218 11801.18346 2325.32435
219 27219.72710 11801.18346
220 -13309.64206 27219.72710
221 -5773.03200 -13309.64206
222 -37869.86553 -5773.03200
223 7438.13087 -37869.86553
224 -76.45526 7438.13087
225 12256.62081 -76.45526
226 14648.21807 12256.62081
227 -21967.69165 14648.21807
228 -17232.05350 -21967.69165
229 -43932.95034 -17232.05350
230 -12799.20274 -43932.95034
231 3166.83219 -12799.20274
232 -37635.26986 3166.83219
233 -10419.00216 -37635.26986
234 -35365.45895 -10419.00216
235 -5678.98516 -35365.45895
236 50773.80086 -5678.98516
237 -7153.96261 50773.80086
238 -10841.98093 -7153.96261
239 -19671.96340 -10841.98093
240 6146.39423 -19671.96340
241 2984.16950 6146.39423
242 -6733.97382 2984.16950
243 -3620.33496 -6733.97382
244 -2406.48171 -3620.33496
245 -10723.89172 -2406.48171
246 -16497.28395 -10723.89172
247 18256.47603 -16497.28395
248 3938.06610 18256.47603
249 -17665.47305 3938.06610
250 7268.01690 -17665.47305
251 9036.43523 7268.01690
252 17479.59684 9036.43523
253 -4178.43422 17479.59684
254 -19965.43252 -4178.43422
255 4593.71816 -19965.43252
256 -9767.78007 4593.71816
257 9758.25289 -9767.78007
258 13572.10719 9758.25289
259 11390.80740 13572.10719
260 -19985.32693 11390.80740
261 18338.82340 -19985.32693
262 12960.62004 18338.82340
263 1632.24302 12960.62004
264 -13783.69448 1632.24302
265 16196.37801 -13783.69448
266 -10197.03482 16196.37801
267 -4218.03514 -10197.03482
268 54073.75254 -4218.03514
269 -26271.69019 54073.75254
270 -6828.52716 -26271.69019
271 5224.04983 -6828.52716
272 2814.00833 5224.04983
273 15639.87683 2814.00833
274 17040.77017 15639.87683
275 -4437.28245 17040.77017
276 11550.69999 -4437.28245
277 15770.09184 11550.69999
278 8027.17270 15770.09184
279 17937.45199 8027.17270
280 624.75388 17937.45199
281 -7275.78206 624.75388
282 838.97164 -7275.78206
283 14537.53539 838.97164
284 1270.06348 14537.53539
285 -10595.53443 1270.06348
286 8749.38291 -10595.53443
287 5308.42024 8749.38291
288 -1388.50917 5308.42024
289 NA -1388.50917
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] -1511.90811 57381.11794
[2,] 10501.17082 -1511.90811
[3,] -74681.58651 10501.17082
[4,] 30819.96881 -74681.58651
[5,] -28378.61573 30819.96881
[6,] 30717.59077 -28378.61573
[7,] -16496.24686 30717.59077
[8,] 15351.12634 -16496.24686
[9,] -35806.27187 15351.12634
[10,] 22297.04832 -35806.27187
[11,] 33374.39258 22297.04832
[12,] 47.29319 33374.39258
[13,] -14512.12552 47.29319
[14,] 20840.28937 -14512.12552
[15,] 32680.77750 20840.28937
[16,] -51961.21467 32680.77750
[17,] 36590.20525 -51961.21467
[18,] -8944.82804 36590.20525
[19,] -1805.12382 -8944.82804
[20,] 16657.91439 -1805.12382
[21,] 12157.59221 16657.91439
[22,] 105642.71844 12157.59221
[23,] 34654.50653 105642.71844
[24,] -19764.26410 34654.50653
[25,] -77424.81227 -19764.26410
[26,] -53496.76887 -77424.81227
[27,] -13349.23944 -53496.76887
[28,] 51291.41395 -13349.23944
[29,] -3494.46935 51291.41395
[30,] 27496.39450 -3494.46935
[31,] 6215.55758 27496.39450
[32,] -9941.44218 6215.55758
[33,] 49487.62786 -9941.44218
[34,] -19571.80669 49487.62786
[35,] 35111.90998 -19571.80669
[36,] 31413.06822 35111.90998
[37,] 8677.49503 31413.06822
[38,] 68999.86930 8677.49503
[39,] -930.54728 68999.86930
[40,] 15321.23850 -930.54728
[41,] 33195.82855 15321.23850
[42,] -32664.16493 33195.82855
[43,] 5860.64248 -32664.16493
[44,] -2899.41017 5860.64248
[45,] -10960.83946 -2899.41017
[46,] -39649.27269 -10960.83946
[47,] -5239.37136 -39649.27269
[48,] 447.46891 -5239.37136
[49,] -125849.14401 447.46891
[50,] -15485.14617 -125849.14401
[51,] -26339.15190 -15485.14617
[52,] -7173.40119 -26339.15190
[53,] -28171.72815 -7173.40119
[54,] -45348.42941 -28171.72815
[55,] 16089.14399 -45348.42941
[56,] -17823.72420 16089.14399
[57,] 51464.29413 -17823.72420
[58,] -3605.32157 51464.29413
[59,] -10452.14157 -3605.32157
[60,] 2505.82322 -10452.14157
[61,] 28728.51718 2505.82322
[62,] 20529.52422 28728.51718
[63,] 6896.81673 20529.52422
[64,] -8827.80248 6896.81673
[65,] 13342.09226 -8827.80248
[66,] 36990.65941 13342.09226
[67,] 20776.23197 36990.65941
[68,] -9379.63009 20776.23197
[69,] -35209.79356 -9379.63009
[70,] 12671.81419 -35209.79356
[71,] 31339.97577 12671.81419
[72,] 47419.76163 31339.97577
[73,] -24847.57652 47419.76163
[74,] -23377.73801 -24847.57652
[75,] -7239.02624 -23377.73801
[76,] -17973.81838 -7239.02624
[77,] 2517.94397 -17973.81838
[78,] -22793.42358 2517.94397
[79,] 15850.00152 -22793.42358
[80,] -7584.09708 15850.00152
[81,] -15554.51361 -7584.09708
[82,] 55267.99624 -15554.51361
[83,] 35136.35832 55267.99624
[84,] -28151.02707 35136.35832
[85,] -7195.42067 -28151.02707
[86,] -14146.61998 -7195.42067
[87,] -14123.02164 -14146.61998
[88,] 69668.44704 -14123.02164
[89,] -12731.14567 69668.44704
[90,] 23330.06097 -12731.14567
[91,] -21940.41106 23330.06097
[92,] -39669.57886 -21940.41106
[93,] 6843.37884 -39669.57886
[94,] 35675.30501 6843.37884
[95,] -10968.58808 35675.30501
[96,] 17234.03830 -10968.58808
[97,] -117585.07544 17234.03830
[98,] 23942.12876 -117585.07544
[99,] 4888.69166 23942.12876
[100,] 10039.32580 4888.69166
[101,] -3575.20141 10039.32580
[102,] 26967.18332 -3575.20141
[103,] -821.09625 26967.18332
[104,] 13579.32886 -821.09625
[105,] 16005.04785 13579.32886
[106,] -14917.78955 16005.04785
[107,] -6858.40724 -14917.78955
[108,] 24773.46998 -6858.40724
[109,] 7354.72167 24773.46998
[110,] 12446.17238 7354.72167
[111,] -23982.64205 12446.17238
[112,] -25068.37781 -23982.64205
[113,] 10644.38922 -25068.37781
[114,] -30074.45043 10644.38922
[115,] -25305.67103 -30074.45043
[116,] 8327.61408 -25305.67103
[117,] 55001.97162 8327.61408
[118,] -56316.34712 55001.97162
[119,] -29540.73584 -56316.34712
[120,] 4903.28961 -29540.73584
[121,] -42085.63943 4903.28961
[122,] -25521.43149 -42085.63943
[123,] -74625.98305 -25521.43149
[124,] 58475.78260 -74625.98305
[125,] 32053.34571 58475.78260
[126,] 10635.28023 32053.34571
[127,] -2522.46536 10635.28023
[128,] 32900.92429 -2522.46536
[129,] 1641.66583 32900.92429
[130,] -913.08066 1641.66583
[131,] -23974.12311 -913.08066
[132,] 12852.10903 -23974.12311
[133,] -54136.94656 12852.10903
[134,] -6791.24183 -54136.94656
[135,] 63386.79098 -6791.24183
[136,] 985.25829 63386.79098
[137,] 27363.47235 985.25829
[138,] -33587.04807 27363.47235
[139,] -15245.41479 -33587.04807
[140,] -19072.70354 -15245.41479
[141,] -6527.94620 -19072.70354
[142,] 1281.75086 -6527.94620
[143,] -5893.30000 1281.75086
[144,] 2075.56041 -5893.30000
[145,] -29273.11738 2075.56041
[146,] -19430.81573 -29273.11738
[147,] -40229.48323 -19430.81573
[148,] -23173.55634 -40229.48323
[149,] -13384.37261 -23173.55634
[150,] 28934.80243 -13384.37261
[151,] 6243.54112 28934.80243
[152,] -2809.87691 6243.54112
[153,] 17253.88781 -2809.87691
[154,] 26862.87247 17253.88781
[155,] 23649.04916 26862.87247
[156,] 20176.49159 23649.04916
[157,] 15695.89678 20176.49159
[158,] -35561.38334 15695.89678
[159,] -23535.62236 -35561.38334
[160,] 15814.68234 -23535.62236
[161,] -41480.04988 15814.68234
[162,] -15514.94985 -41480.04988
[163,] -33213.15976 -15514.94985
[164,] -34907.34146 -33213.15976
[165,] 42557.73770 -34907.34146
[166,] -7055.04988 42557.73770
[167,] -15941.52180 -7055.04988
[168,] 33073.81360 -15941.52180
[169,] 77843.37210 33073.81360
[170,] -10567.50642 77843.37210
[171,] 10517.82044 -10567.50642
[172,] 6516.62322 10517.82044
[173,] 37028.73077 6516.62322
[174,] -75156.95534 37028.73077
[175,] 3274.91471 -75156.95534
[176,] -47406.71555 3274.91471
[177,] 11433.78388 -47406.71555
[178,] 4374.93135 11433.78388
[179,] 44340.19551 4374.93135
[180,] -17316.73907 44340.19551
[181,] -47022.02775 -17316.73907
[182,] 4406.14997 -47022.02775
[183,] 12614.75872 4406.14997
[184,] 42235.77627 12614.75872
[185,] 1752.18093 42235.77627
[186,] -2522.53629 1752.18093
[187,] 10577.55326 -2522.53629
[188,] 16232.29357 10577.55326
[189,] -5684.49964 16232.29357
[190,] 16221.82643 -5684.49964
[191,] -8738.12690 16221.82643
[192,] -6399.47350 -8738.12690
[193,] -21571.45769 -6399.47350
[194,] 4443.50282 -21571.45769
[195,] 20560.04069 4443.50282
[196,] -20304.07951 20560.04069
[197,] -2995.87050 -20304.07951
[198,] -35850.38990 -2995.87050
[199,] -1933.77753 -35850.38990
[200,] 10807.51736 -1933.77753
[201,] -503.27315 10807.51736
[202,] 17726.98844 -503.27315
[203,] 1942.31810 17726.98844
[204,] -29043.29067 1942.31810
[205,] -12892.32868 -29043.29067
[206,] 2098.58134 -12892.32868
[207,] 42595.32910 2098.58134
[208,] 27327.55140 42595.32910
[209,] 15214.97874 27327.55140
[210,] 20834.74317 15214.97874
[211,] 33627.72752 20834.74317
[212,] -1401.91244 33627.72752
[213,] 9847.53970 -1401.91244
[214,] -37664.50139 9847.53970
[215,] -8444.96432 -37664.50139
[216,] -29697.52180 -8444.96432
[217,] 2325.32435 -29697.52180
[218,] 11801.18346 2325.32435
[219,] 27219.72710 11801.18346
[220,] -13309.64206 27219.72710
[221,] -5773.03200 -13309.64206
[222,] -37869.86553 -5773.03200
[223,] 7438.13087 -37869.86553
[224,] -76.45526 7438.13087
[225,] 12256.62081 -76.45526
[226,] 14648.21807 12256.62081
[227,] -21967.69165 14648.21807
[228,] -17232.05350 -21967.69165
[229,] -43932.95034 -17232.05350
[230,] -12799.20274 -43932.95034
[231,] 3166.83219 -12799.20274
[232,] -37635.26986 3166.83219
[233,] -10419.00216 -37635.26986
[234,] -35365.45895 -10419.00216
[235,] -5678.98516 -35365.45895
[236,] 50773.80086 -5678.98516
[237,] -7153.96261 50773.80086
[238,] -10841.98093 -7153.96261
[239,] -19671.96340 -10841.98093
[240,] 6146.39423 -19671.96340
[241,] 2984.16950 6146.39423
[242,] -6733.97382 2984.16950
[243,] -3620.33496 -6733.97382
[244,] -2406.48171 -3620.33496
[245,] -10723.89172 -2406.48171
[246,] -16497.28395 -10723.89172
[247,] 18256.47603 -16497.28395
[248,] 3938.06610 18256.47603
[249,] -17665.47305 3938.06610
[250,] 7268.01690 -17665.47305
[251,] 9036.43523 7268.01690
[252,] 17479.59684 9036.43523
[253,] -4178.43422 17479.59684
[254,] -19965.43252 -4178.43422
[255,] 4593.71816 -19965.43252
[256,] -9767.78007 4593.71816
[257,] 9758.25289 -9767.78007
[258,] 13572.10719 9758.25289
[259,] 11390.80740 13572.10719
[260,] -19985.32693 11390.80740
[261,] 18338.82340 -19985.32693
[262,] 12960.62004 18338.82340
[263,] 1632.24302 12960.62004
[264,] -13783.69448 1632.24302
[265,] 16196.37801 -13783.69448
[266,] -10197.03482 16196.37801
[267,] -4218.03514 -10197.03482
[268,] 54073.75254 -4218.03514
[269,] -26271.69019 54073.75254
[270,] -6828.52716 -26271.69019
[271,] 5224.04983 -6828.52716
[272,] 2814.00833 5224.04983
[273,] 15639.87683 2814.00833
[274,] 17040.77017 15639.87683
[275,] -4437.28245 17040.77017
[276,] 11550.69999 -4437.28245
[277,] 15770.09184 11550.69999
[278,] 8027.17270 15770.09184
[279,] 17937.45199 8027.17270
[280,] 624.75388 17937.45199
[281,] -7275.78206 624.75388
[282,] 838.97164 -7275.78206
[283,] 14537.53539 838.97164
[284,] 1270.06348 14537.53539
[285,] -10595.53443 1270.06348
[286,] 8749.38291 -10595.53443
[287,] 5308.42024 8749.38291
[288,] -1388.50917 5308.42024
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 -1511.90811 57381.11794
2 10501.17082 -1511.90811
3 -74681.58651 10501.17082
4 30819.96881 -74681.58651
5 -28378.61573 30819.96881
6 30717.59077 -28378.61573
7 -16496.24686 30717.59077
8 15351.12634 -16496.24686
9 -35806.27187 15351.12634
10 22297.04832 -35806.27187
11 33374.39258 22297.04832
12 47.29319 33374.39258
13 -14512.12552 47.29319
14 20840.28937 -14512.12552
15 32680.77750 20840.28937
16 -51961.21467 32680.77750
17 36590.20525 -51961.21467
18 -8944.82804 36590.20525
19 -1805.12382 -8944.82804
20 16657.91439 -1805.12382
21 12157.59221 16657.91439
22 105642.71844 12157.59221
23 34654.50653 105642.71844
24 -19764.26410 34654.50653
25 -77424.81227 -19764.26410
26 -53496.76887 -77424.81227
27 -13349.23944 -53496.76887
28 51291.41395 -13349.23944
29 -3494.46935 51291.41395
30 27496.39450 -3494.46935
31 6215.55758 27496.39450
32 -9941.44218 6215.55758
33 49487.62786 -9941.44218
34 -19571.80669 49487.62786
35 35111.90998 -19571.80669
36 31413.06822 35111.90998
37 8677.49503 31413.06822
38 68999.86930 8677.49503
39 -930.54728 68999.86930
40 15321.23850 -930.54728
41 33195.82855 15321.23850
42 -32664.16493 33195.82855
43 5860.64248 -32664.16493
44 -2899.41017 5860.64248
45 -10960.83946 -2899.41017
46 -39649.27269 -10960.83946
47 -5239.37136 -39649.27269
48 447.46891 -5239.37136
49 -125849.14401 447.46891
50 -15485.14617 -125849.14401
51 -26339.15190 -15485.14617
52 -7173.40119 -26339.15190
53 -28171.72815 -7173.40119
54 -45348.42941 -28171.72815
55 16089.14399 -45348.42941
56 -17823.72420 16089.14399
57 51464.29413 -17823.72420
58 -3605.32157 51464.29413
59 -10452.14157 -3605.32157
60 2505.82322 -10452.14157
61 28728.51718 2505.82322
62 20529.52422 28728.51718
63 6896.81673 20529.52422
64 -8827.80248 6896.81673
65 13342.09226 -8827.80248
66 36990.65941 13342.09226
67 20776.23197 36990.65941
68 -9379.63009 20776.23197
69 -35209.79356 -9379.63009
70 12671.81419 -35209.79356
71 31339.97577 12671.81419
72 47419.76163 31339.97577
73 -24847.57652 47419.76163
74 -23377.73801 -24847.57652
75 -7239.02624 -23377.73801
76 -17973.81838 -7239.02624
77 2517.94397 -17973.81838
78 -22793.42358 2517.94397
79 15850.00152 -22793.42358
80 -7584.09708 15850.00152
81 -15554.51361 -7584.09708
82 55267.99624 -15554.51361
83 35136.35832 55267.99624
84 -28151.02707 35136.35832
85 -7195.42067 -28151.02707
86 -14146.61998 -7195.42067
87 -14123.02164 -14146.61998
88 69668.44704 -14123.02164
89 -12731.14567 69668.44704
90 23330.06097 -12731.14567
91 -21940.41106 23330.06097
92 -39669.57886 -21940.41106
93 6843.37884 -39669.57886
94 35675.30501 6843.37884
95 -10968.58808 35675.30501
96 17234.03830 -10968.58808
97 -117585.07544 17234.03830
98 23942.12876 -117585.07544
99 4888.69166 23942.12876
100 10039.32580 4888.69166
101 -3575.20141 10039.32580
102 26967.18332 -3575.20141
103 -821.09625 26967.18332
104 13579.32886 -821.09625
105 16005.04785 13579.32886
106 -14917.78955 16005.04785
107 -6858.40724 -14917.78955
108 24773.46998 -6858.40724
109 7354.72167 24773.46998
110 12446.17238 7354.72167
111 -23982.64205 12446.17238
112 -25068.37781 -23982.64205
113 10644.38922 -25068.37781
114 -30074.45043 10644.38922
115 -25305.67103 -30074.45043
116 8327.61408 -25305.67103
117 55001.97162 8327.61408
118 -56316.34712 55001.97162
119 -29540.73584 -56316.34712
120 4903.28961 -29540.73584
121 -42085.63943 4903.28961
122 -25521.43149 -42085.63943
123 -74625.98305 -25521.43149
124 58475.78260 -74625.98305
125 32053.34571 58475.78260
126 10635.28023 32053.34571
127 -2522.46536 10635.28023
128 32900.92429 -2522.46536
129 1641.66583 32900.92429
130 -913.08066 1641.66583
131 -23974.12311 -913.08066
132 12852.10903 -23974.12311
133 -54136.94656 12852.10903
134 -6791.24183 -54136.94656
135 63386.79098 -6791.24183
136 985.25829 63386.79098
137 27363.47235 985.25829
138 -33587.04807 27363.47235
139 -15245.41479 -33587.04807
140 -19072.70354 -15245.41479
141 -6527.94620 -19072.70354
142 1281.75086 -6527.94620
143 -5893.30000 1281.75086
144 2075.56041 -5893.30000
145 -29273.11738 2075.56041
146 -19430.81573 -29273.11738
147 -40229.48323 -19430.81573
148 -23173.55634 -40229.48323
149 -13384.37261 -23173.55634
150 28934.80243 -13384.37261
151 6243.54112 28934.80243
152 -2809.87691 6243.54112
153 17253.88781 -2809.87691
154 26862.87247 17253.88781
155 23649.04916 26862.87247
156 20176.49159 23649.04916
157 15695.89678 20176.49159
158 -35561.38334 15695.89678
159 -23535.62236 -35561.38334
160 15814.68234 -23535.62236
161 -41480.04988 15814.68234
162 -15514.94985 -41480.04988
163 -33213.15976 -15514.94985
164 -34907.34146 -33213.15976
165 42557.73770 -34907.34146
166 -7055.04988 42557.73770
167 -15941.52180 -7055.04988
168 33073.81360 -15941.52180
169 77843.37210 33073.81360
170 -10567.50642 77843.37210
171 10517.82044 -10567.50642
172 6516.62322 10517.82044
173 37028.73077 6516.62322
174 -75156.95534 37028.73077
175 3274.91471 -75156.95534
176 -47406.71555 3274.91471
177 11433.78388 -47406.71555
178 4374.93135 11433.78388
179 44340.19551 4374.93135
180 -17316.73907 44340.19551
181 -47022.02775 -17316.73907
182 4406.14997 -47022.02775
183 12614.75872 4406.14997
184 42235.77627 12614.75872
185 1752.18093 42235.77627
186 -2522.53629 1752.18093
187 10577.55326 -2522.53629
188 16232.29357 10577.55326
189 -5684.49964 16232.29357
190 16221.82643 -5684.49964
191 -8738.12690 16221.82643
192 -6399.47350 -8738.12690
193 -21571.45769 -6399.47350
194 4443.50282 -21571.45769
195 20560.04069 4443.50282
196 -20304.07951 20560.04069
197 -2995.87050 -20304.07951
198 -35850.38990 -2995.87050
199 -1933.77753 -35850.38990
200 10807.51736 -1933.77753
201 -503.27315 10807.51736
202 17726.98844 -503.27315
203 1942.31810 17726.98844
204 -29043.29067 1942.31810
205 -12892.32868 -29043.29067
206 2098.58134 -12892.32868
207 42595.32910 2098.58134
208 27327.55140 42595.32910
209 15214.97874 27327.55140
210 20834.74317 15214.97874
211 33627.72752 20834.74317
212 -1401.91244 33627.72752
213 9847.53970 -1401.91244
214 -37664.50139 9847.53970
215 -8444.96432 -37664.50139
216 -29697.52180 -8444.96432
217 2325.32435 -29697.52180
218 11801.18346 2325.32435
219 27219.72710 11801.18346
220 -13309.64206 27219.72710
221 -5773.03200 -13309.64206
222 -37869.86553 -5773.03200
223 7438.13087 -37869.86553
224 -76.45526 7438.13087
225 12256.62081 -76.45526
226 14648.21807 12256.62081
227 -21967.69165 14648.21807
228 -17232.05350 -21967.69165
229 -43932.95034 -17232.05350
230 -12799.20274 -43932.95034
231 3166.83219 -12799.20274
232 -37635.26986 3166.83219
233 -10419.00216 -37635.26986
234 -35365.45895 -10419.00216
235 -5678.98516 -35365.45895
236 50773.80086 -5678.98516
237 -7153.96261 50773.80086
238 -10841.98093 -7153.96261
239 -19671.96340 -10841.98093
240 6146.39423 -19671.96340
241 2984.16950 6146.39423
242 -6733.97382 2984.16950
243 -3620.33496 -6733.97382
244 -2406.48171 -3620.33496
245 -10723.89172 -2406.48171
246 -16497.28395 -10723.89172
247 18256.47603 -16497.28395
248 3938.06610 18256.47603
249 -17665.47305 3938.06610
250 7268.01690 -17665.47305
251 9036.43523 7268.01690
252 17479.59684 9036.43523
253 -4178.43422 17479.59684
254 -19965.43252 -4178.43422
255 4593.71816 -19965.43252
256 -9767.78007 4593.71816
257 9758.25289 -9767.78007
258 13572.10719 9758.25289
259 11390.80740 13572.10719
260 -19985.32693 11390.80740
261 18338.82340 -19985.32693
262 12960.62004 18338.82340
263 1632.24302 12960.62004
264 -13783.69448 1632.24302
265 16196.37801 -13783.69448
266 -10197.03482 16196.37801
267 -4218.03514 -10197.03482
268 54073.75254 -4218.03514
269 -26271.69019 54073.75254
270 -6828.52716 -26271.69019
271 5224.04983 -6828.52716
272 2814.00833 5224.04983
273 15639.87683 2814.00833
274 17040.77017 15639.87683
275 -4437.28245 17040.77017
276 11550.69999 -4437.28245
277 15770.09184 11550.69999
278 8027.17270 15770.09184
279 17937.45199 8027.17270
280 624.75388 17937.45199
281 -7275.78206 624.75388
282 838.97164 -7275.78206
283 14537.53539 838.97164
284 1270.06348 14537.53539
285 -10595.53443 1270.06348
286 8749.38291 -10595.53443
287 5308.42024 8749.38291
288 -1388.50917 5308.42024
> 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/wessaorg/rcomp/tmp/7hjkl1355218833.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/wessaorg/rcomp/tmp/82jju1355218833.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/wessaorg/rcomp/tmp/9fx901355218833.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/wessaorg/rcomp/tmp/10e5q21355218833.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/wessaorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/wessaorg/rcomp/createtable")
>
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Estimated Regression Equation', 1, TRUE)
> a<-table.row.end(a)
> myeq <- colnames(x)[1]
> myeq <- paste(myeq, '[t] = ', sep='')
> for (i in 1:k){
+ if (mysum$coefficients[i,1] > 0) myeq <- paste(myeq, '+', '')
+ myeq <- paste(myeq, mysum$coefficients[i,1], sep=' ')
+ if (rownames(mysum$coefficients)[i] != '(Intercept)') {
+ myeq <- paste(myeq, rownames(mysum$coefficients)[i], sep='')
+ if (rownames(mysum$coefficients)[i] != 't') myeq <- paste(myeq, '[t]', sep='')
+ }
+ }
> myeq <- paste(myeq, ' + e[t]')
> a<-table.row.start(a)
> a<-table.element(a, myeq)
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/wessaorg/rcomp/tmp/111hde1355218833.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/wessaorg/rcomp/tmp/12d6k11355218833.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/wessaorg/rcomp/tmp/13dztz1355218833.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/wessaorg/rcomp/tmp/14ycou1355218833.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/wessaorg/rcomp/tmp/15stqa1355218833.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/wessaorg/rcomp/tmp/160cpb1355218833.tab")
+ }
>
> try(system("convert tmp/16j6x1355218833.ps tmp/16j6x1355218833.png",intern=TRUE))
character(0)
> try(system("convert tmp/2y2ko1355218833.ps tmp/2y2ko1355218833.png",intern=TRUE))
character(0)
> try(system("convert tmp/3joe61355218833.ps tmp/3joe61355218833.png",intern=TRUE))
character(0)
> try(system("convert tmp/4g21p1355218833.ps tmp/4g21p1355218833.png",intern=TRUE))
character(0)
> try(system("convert tmp/5son71355218833.ps tmp/5son71355218833.png",intern=TRUE))
character(0)
> try(system("convert tmp/66d1e1355218833.ps tmp/66d1e1355218833.png",intern=TRUE))
character(0)
> try(system("convert tmp/7hjkl1355218833.ps tmp/7hjkl1355218833.png",intern=TRUE))
character(0)
> try(system("convert tmp/82jju1355218833.ps tmp/82jju1355218833.png",intern=TRUE))
character(0)
> try(system("convert tmp/9fx901355218833.ps tmp/9fx901355218833.png",intern=TRUE))
character(0)
> try(system("convert tmp/10e5q21355218833.ps tmp/10e5q21355218833.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
12.129 1.078 13.220