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
+ ,112285
+ ,24188
+ ,146283
+ ,144
+ ,145
+ ,120982
+ ,56
+ ,84786
+ ,18273
+ ,98364
+ ,103
+ ,101
+ ,176508
+ ,54
+ ,83123
+ ,14130
+ ,86146
+ ,98
+ ,98
+ ,179321
+ ,89
+ ,101193
+ ,32287
+ ,96933
+ ,135
+ ,132
+ ,123185
+ ,40
+ ,38361
+ ,8654
+ ,79234
+ ,61
+ ,60
+ ,52746
+ ,25
+ ,68504
+ ,9245
+ ,42551
+ ,39
+ ,38
+ ,385534
+ ,92
+ ,119182
+ ,33251
+ ,195663
+ ,150
+ ,144
+ ,33170
+ ,18
+ ,22807
+ ,1271
+ ,6853
+ ,5
+ ,5
+ ,101645
+ ,63
+ ,17140
+ ,5279
+ ,21529
+ ,28
+ ,28
+ ,149061
+ ,44
+ ,116174
+ ,27101
+ ,95757
+ ,84
+ ,84
+ ,165446
+ ,33
+ ,57635
+ ,16373
+ ,85584
+ ,80
+ ,79
+ ,237213
+ ,84
+ ,66198
+ ,19716
+ ,143983
+ ,130
+ ,127
+ ,173326
+ ,88
+ ,71701
+ ,17753
+ ,75851
+ ,82
+ ,78
+ ,133131
+ ,55
+ ,57793
+ ,9028
+ ,59238
+ ,60
+ ,60
+ ,258873
+ ,60
+ ,80444
+ ,18653
+ ,93163
+ ,131
+ ,131
+ ,180083
+ ,66
+ ,53855
+ ,8828
+ ,96037
+ ,84
+ ,84
+ ,324799
+ ,154
+ ,97668
+ ,29498
+ ,151511
+ ,140
+ ,133
+ ,230964
+ ,53
+ ,133824
+ ,27563
+ ,136368
+ ,151
+ ,150
+ ,236785
+ ,119
+ ,101481
+ ,18293
+ ,112642
+ ,91
+ ,91
+ ,135473
+ ,41
+ ,99645
+ ,22530
+ ,94728
+ ,138
+ ,132
+ ,202925
+ ,61
+ ,114789
+ ,15977
+ ,105499
+ ,150
+ ,136
+ ,215147
+ ,58
+ ,99052
+ ,35082
+ ,121527
+ ,124
+ ,124
+ ,344297
+ ,75
+ ,67654
+ ,16116
+ ,127766
+ ,119
+ ,118
+ ,153935
+ ,33
+ ,65553
+ ,15849
+ ,98958
+ ,73
+ ,70
+ ,132943
+ ,40
+ ,97500
+ ,16026
+ ,77900
+ ,110
+ ,107
+ ,174724
+ ,92
+ ,69112
+ ,26569
+ ,85646
+ ,123
+ ,119
+ ,174415
+ ,100
+ ,82753
+ ,24785
+ ,98579
+ ,90
+ ,89
+ ,225548
+ ,112
+ ,85323
+ ,17569
+ ,130767
+ ,116
+ ,112
+ ,223632
+ ,73
+ ,72654
+ ,23825
+ ,131741
+ ,113
+ ,108
+ ,124817
+ ,40
+ ,30727
+ ,7869
+ ,53907
+ ,56
+ ,52
+ ,221698
+ ,45
+ ,77873
+ ,14975
+ ,178812
+ ,115
+ ,112
+ ,210767
+ ,60
+ ,117478
+ ,37791
+ ,146761
+ ,119
+ ,116
+ ,170266
+ ,62
+ ,74007
+ ,9605
+ ,82036
+ ,129
+ ,123
+ ,260561
+ ,75
+ ,90183
+ ,27295
+ ,163253
+ ,127
+ ,125
+ ,84853
+ ,31
+ ,61542
+ ,2746
+ ,27032
+ ,27
+ ,27
+ ,294424
+ ,77
+ ,101494
+ ,34461
+ ,171975
+ ,175
+ ,162
+ ,101011
+ ,34
+ ,27570
+ ,8098
+ ,65990
+ ,35
+ ,32
+ ,215641
+ ,46
+ ,55813
+ ,4787
+ ,86572
+ ,64
+ ,64
+ ,325107
+ ,99
+ ,79215
+ ,24919
+ ,159676
+ ,96
+ ,92
+ ,7176
+ ,17
+ ,1423
+ ,603
+ ,1929
+ ,0
+ ,0
+ ,167542
+ ,66
+ ,55461
+ ,16329
+ ,85371
+ ,84
+ ,83
+ ,106408
+ ,30
+ ,31081
+ ,12558
+ ,58391
+ ,41
+ ,41
+ ,96560
+ ,76
+ ,22996
+ ,7784
+ ,31580
+ ,47
+ ,47
+ ,265769
+ ,146
+ ,83122
+ ,28522
+ ,136815
+ ,126
+ ,120
+ ,269651
+ ,67
+ ,70106
+ ,22265
+ ,120642
+ ,105
+ ,105
+ ,149112
+ ,56
+ ,60578
+ ,14459
+ ,69107
+ ,80
+ ,79
+ ,175824
+ ,107
+ ,39992
+ ,14526
+ ,50495
+ ,70
+ ,65
+ ,152871
+ ,58
+ ,79892
+ ,22240
+ ,108016
+ ,73
+ ,70
+ ,111665
+ ,34
+ ,49810
+ ,11802
+ ,46341
+ ,57
+ ,55
+ ,116408
+ ,61
+ ,71570
+ ,7623
+ ,78348
+ ,40
+ ,39
+ ,362301
+ ,119
+ ,100708
+ ,11912
+ ,79336
+ ,68
+ ,67
+ ,78800
+ ,42
+ ,33032
+ ,7935
+ ,56968
+ ,21
+ ,21
+ ,183167
+ ,66
+ ,82875
+ ,18220
+ ,93176
+ ,127
+ ,127
+ ,277965
+ ,89
+ ,139077
+ ,19199
+ ,161632
+ ,154
+ ,152
+ ,150629
+ ,44
+ ,71595
+ ,19918
+ ,87850
+ ,116
+ ,113
+ ,168809
+ ,66
+ ,72260
+ ,21884
+ ,127969
+ ,102
+ ,99
+ ,24188
+ ,24
+ ,5950
+ ,2694
+ ,15049
+ ,7
+ ,7
+ ,329267
+ ,259
+ ,115762
+ ,15808
+ ,155135
+ ,148
+ ,141
+ ,65029
+ ,17
+ ,32551
+ ,3597
+ ,25109
+ ,21
+ ,21
+ ,101097
+ ,64
+ ,31701
+ ,5296
+ ,45824
+ ,35
+ ,35
+ ,218946
+ ,41
+ ,80670
+ ,25239
+ ,102996
+ ,112
+ ,109
+ ,244052
+ ,68
+ ,143558
+ ,29801
+ ,160604
+ ,137
+ ,133
+ ,341570
+ ,168
+ ,117105
+ ,18450
+ ,158051
+ ,135
+ ,123
+ ,103597
+ ,43
+ ,23789
+ ,7132
+ ,44547
+ ,26
+ ,26
+ ,233328
+ ,132
+ ,120733
+ ,34861
+ ,162647
+ ,230
+ ,230
+ ,256462
+ ,105
+ ,105195
+ ,35940
+ ,174141
+ ,181
+ ,166
+ ,206161
+ ,71
+ ,73107
+ ,16688
+ ,60622
+ ,71
+ ,68
+ ,311473
+ ,112
+ ,132068
+ ,24683
+ ,179566
+ ,147
+ ,147
+ ,235800
+ ,94
+ ,149193
+ ,46230
+ ,184301
+ ,190
+ ,179
+ ,177939
+ ,82
+ ,46821
+ ,10387
+ ,75661
+ ,64
+ ,61
+ ,207176
+ ,70
+ ,87011
+ ,21436
+ ,96144
+ ,105
+ ,101
+ ,196553
+ ,57
+ ,95260
+ ,30546
+ ,129847
+ ,107
+ ,108
+ ,174184
+ ,53
+ ,55183
+ ,19746
+ ,117286
+ ,94
+ ,90
+ ,143246
+ ,103
+ ,106671
+ ,15977
+ ,71180
+ ,116
+ ,114
+ ,187559
+ ,121
+ ,73511
+ ,22583
+ ,109377
+ ,106
+ ,103
+ ,187681
+ ,62
+ ,92945
+ ,17274
+ ,85298
+ ,143
+ ,142
+ ,119016
+ ,52
+ ,78664
+ ,16469
+ ,73631
+ ,81
+ ,79
+ ,182192
+ ,52
+ ,70054
+ ,14251
+ ,86767
+ ,89
+ ,88
+ ,73566
+ ,32
+ ,22618
+ ,3007
+ ,23824
+ ,26
+ ,25
+ ,194979
+ ,62
+ ,74011
+ ,16851
+ ,93487
+ ,84
+ ,83
+ ,167488
+ ,45
+ ,83737
+ ,21113
+ ,82981
+ ,113
+ ,113
+ ,143756
+ ,46
+ ,69094
+ ,17401
+ ,73815
+ ,120
+ ,118
+ ,275541
+ ,63
+ ,93133
+ ,23958
+ ,94552
+ ,110
+ ,110
+ ,243199
+ ,75
+ ,95536
+ ,23567
+ ,132190
+ ,134
+ ,129
+ ,182999
+ ,88
+ ,225920
+ ,13065
+ ,128754
+ ,54
+ ,51
+ ,135649
+ ,46
+ ,62133
+ ,15358
+ ,66363
+ ,96
+ ,93
+ ,152299
+ ,53
+ ,61370
+ ,14587
+ ,67808
+ ,78
+ ,76
+ ,120221
+ ,37
+ ,43836
+ ,12770
+ ,61724
+ ,51
+ ,49
+ ,346485
+ ,90
+ ,106117
+ ,24021
+ ,131722
+ ,121
+ ,118
+ ,145790
+ ,63
+ ,38692
+ ,9648
+ ,68580
+ ,38
+ ,38
+ ,193339
+ ,78
+ ,84651
+ ,20537
+ ,106175
+ ,145
+ ,141
+ ,80953
+ ,25
+ ,56622
+ ,7905
+ ,55792
+ ,59
+ ,58
+ ,122774
+ ,45
+ ,15986
+ ,4527
+ ,25157
+ ,27
+ ,27
+ ,130585
+ ,46
+ ,95364
+ ,30495
+ ,76669
+ ,91
+ ,91
+ ,112611
+ ,41
+ ,26706
+ ,7117
+ ,57283
+ ,48
+ ,48
+ ,286468
+ ,144
+ ,89691
+ ,17719
+ ,105805
+ ,68
+ ,63
+ ,241066
+ ,82
+ ,67267
+ ,27056
+ ,129484
+ ,58
+ ,56
+ ,148446
+ ,91
+ ,126846
+ ,33473
+ ,72413
+ ,150
+ ,144
+ ,204713
+ ,71
+ ,41140
+ ,9758
+ ,87831
+ ,74
+ ,73
+ ,182079
+ ,63
+ ,102860
+ ,21115
+ ,96971
+ ,181
+ ,168
+ ,140344
+ ,53
+ ,51715
+ ,7236
+ ,71299
+ ,65
+ ,64
+ ,220516
+ ,62
+ ,55801
+ ,13790
+ ,77494
+ ,97
+ ,97
+ ,243060
+ ,63
+ ,111813
+ ,32902
+ ,120336
+ ,121
+ ,117
+ ,162765
+ ,32
+ ,120293
+ ,25131
+ ,93913
+ ,99
+ ,100
+ ,182613
+ ,39
+ ,138599
+ ,30910
+ ,136048
+ ,152
+ ,149
+ ,232138
+ ,62
+ ,161647
+ ,35947
+ ,181248
+ ,188
+ ,187
+ ,265318
+ ,117
+ ,115929
+ ,29848
+ ,146123
+ ,138
+ ,127
+ ,85574
+ ,34
+ ,24266
+ ,6943
+ ,32036
+ ,40
+ ,37
+ ,310839
+ ,92
+ ,162901
+ ,42705
+ ,186646
+ ,254
+ ,245
+ ,225060
+ ,93
+ ,109825
+ ,31808
+ ,102255
+ ,87
+ ,87
+ ,232317
+ ,54
+ ,129838
+ ,26675
+ ,168237
+ ,178
+ ,177
+ ,144966
+ ,144
+ ,37510
+ ,8435
+ ,64219
+ ,51
+ ,49
+ ,43287
+ ,14
+ ,43750
+ ,7409
+ ,19630
+ ,49
+ ,49
+ ,155754
+ ,61
+ ,40652
+ ,14993
+ ,76825
+ ,73
+ ,73
+ ,164709
+ ,109
+ ,87771
+ ,36867
+ ,115338
+ ,176
+ ,177
+ ,201940
+ ,38
+ ,85872
+ ,33835
+ ,109427
+ ,94
+ ,94
+ ,235454
+ ,73
+ ,89275
+ ,24164
+ ,118168
+ ,120
+ ,117
+ ,220801
+ ,75
+ ,44418
+ ,12607
+ ,84845
+ ,66
+ ,60
+ ,99466
+ ,50
+ ,192565
+ ,22609
+ ,153197
+ ,56
+ ,55
+ ,92661
+ ,61
+ ,35232
+ ,5892
+ ,29877
+ ,39
+ ,39
+ ,133328
+ ,55
+ ,40909
+ ,17014
+ ,63506
+ ,66
+ ,64
+ ,61361
+ ,77
+ ,13294
+ ,5394
+ ,22445
+ ,27
+ ,26
+ ,125930
+ ,75
+ ,32387
+ ,9178
+ ,47695
+ ,65
+ ,64
+ ,100750
+ ,72
+ ,140867
+ ,6440
+ ,68370
+ ,58
+ ,58
+ ,224549
+ ,50
+ ,120662
+ ,21916
+ ,146304
+ ,98
+ ,95
+ ,82316
+ ,32
+ ,21233
+ ,4011
+ ,38233
+ ,25
+ ,25
+ ,102010
+ ,53
+ ,44332
+ ,5818
+ ,42071
+ ,26
+ ,26
+ ,101523
+ ,42
+ ,61056
+ ,18647
+ ,50517
+ ,77
+ ,76
+ ,243511
+ ,71
+ ,101338
+ ,20556
+ ,103950
+ ,130
+ ,129
+ ,22938
+ ,10
+ ,1168
+ ,238
+ ,5841
+ ,11
+ ,11
+ ,41566
+ ,35
+ ,13497
+ ,70
+ ,2341
+ ,2
+ ,2
+ ,152474
+ ,65
+ ,65567
+ ,22392
+ ,84396
+ ,101
+ ,101
+ ,61857
+ ,25
+ ,25162
+ ,3913
+ ,24610
+ ,31
+ ,28
+ ,99923
+ ,66
+ ,32334
+ ,12237
+ ,35753
+ ,36
+ ,36
+ ,132487
+ ,41
+ ,40735
+ ,8388
+ ,55515
+ ,120
+ ,89
+ ,317394
+ ,86
+ ,91413
+ ,22120
+ ,209056
+ ,195
+ ,193
+ ,21054
+ ,16
+ ,855
+ ,338
+ ,6622
+ ,4
+ ,4
+ ,209641
+ ,42
+ ,97068
+ ,11727
+ ,115814
+ ,89
+ ,84
+ ,22648
+ ,19
+ ,44339
+ ,3704
+ ,11609
+ ,24
+ ,23
+ ,31414
+ ,19
+ ,14116
+ ,3988
+ ,13155
+ ,39
+ ,39
+ ,46698
+ ,45
+ ,10288
+ ,3030
+ ,18274
+ ,14
+ ,14
+ ,131698
+ ,65
+ ,65622
+ ,13520
+ ,72875
+ ,78
+ ,78
+ ,91735
+ ,35
+ ,16563
+ ,1421
+ ,10112
+ ,15
+ ,14
+ ,244749
+ ,95
+ ,76643
+ ,20923
+ ,142775
+ ,106
+ ,101
+ ,184510
+ ,49
+ ,110681
+ ,20237
+ ,68847
+ ,83
+ ,82
+ ,79863
+ ,37
+ ,29011
+ ,3219
+ ,17659
+ ,24
+ ,24
+ ,128423
+ ,64
+ ,92696
+ ,3769
+ ,20112
+ ,37
+ ,36
+ ,97839
+ ,38
+ ,94785
+ ,12252
+ ,61023
+ ,77
+ ,75
+ ,38214
+ ,34
+ ,8773
+ ,1888
+ ,13983
+ ,16
+ ,16
+ ,151101
+ ,32
+ ,83209
+ ,14497
+ ,65176
+ ,56
+ ,55
+ ,272458
+ ,65
+ ,93815
+ ,28864
+ ,132432
+ ,132
+ ,131
+ ,172494
+ ,52
+ ,86687
+ ,21721
+ ,112494
+ ,144
+ ,131
+ ,108043
+ ,62
+ ,34553
+ ,4821
+ ,45109
+ ,40
+ ,39
+ ,328107
+ ,65
+ ,105547
+ ,33644
+ ,170875
+ ,153
+ ,144
+ ,250579
+ ,83
+ ,103487
+ ,15923
+ ,180759
+ ,143
+ ,139
+ ,351067
+ ,95
+ ,213688
+ ,42935
+ ,214921
+ ,220
+ ,211
+ ,158015
+ ,29
+ ,71220
+ ,18864
+ ,100226
+ ,79
+ ,78
+ ,98866
+ ,18
+ ,23517
+ ,4977
+ ,32043
+ ,50
+ ,50
+ ,85439
+ ,33
+ ,56926
+ ,7785
+ ,54454
+ ,39
+ ,39
+ ,229242
+ ,247
+ ,91721
+ ,17939
+ ,78876
+ ,95
+ ,90
+ ,351619
+ ,139
+ ,115168
+ ,23436
+ ,170745
+ ,169
+ ,166
+ ,84207
+ ,29
+ ,111194
+ ,325
+ ,6940
+ ,12
+ ,12
+ ,120445
+ ,118
+ ,51009
+ ,13539
+ ,49025
+ ,63
+ ,57
+ ,324598
+ ,110
+ ,135777
+ ,34538
+ ,122037
+ ,134
+ ,133
+ ,131069
+ ,67
+ ,51513
+ ,12198
+ ,53782
+ ,69
+ ,69
+ ,204271
+ ,42
+ ,74163
+ ,26924
+ ,127748
+ ,119
+ ,119
+ ,165543
+ ,65
+ ,51633
+ ,12716
+ ,86839
+ ,119
+ ,119
+ ,141722
+ ,94
+ ,75345
+ ,8172
+ ,44830
+ ,75
+ ,65
+ ,116048
+ ,64
+ ,33416
+ ,10855
+ ,77395
+ ,63
+ ,61
+ ,250047
+ ,81
+ ,83305
+ ,11932
+ ,89324
+ ,55
+ ,49
+ ,299775
+ ,95
+ ,98952
+ ,14300
+ ,103300
+ ,103
+ ,101
+ ,195838
+ ,67
+ ,102372
+ ,25515
+ ,112283
+ ,197
+ ,196
+ ,173260
+ ,63
+ ,37238
+ ,2805
+ ,10901
+ ,16
+ ,15
+ ,254488
+ ,83
+ ,103772
+ ,29402
+ ,120691
+ ,140
+ ,136
+ ,104389
+ ,45
+ ,123969
+ ,16440
+ ,58106
+ ,89
+ ,89
+ ,136084
+ ,30
+ ,27142
+ ,11221
+ ,57140
+ ,40
+ ,40
+ ,199476
+ ,70
+ ,135400
+ ,28732
+ ,122422
+ ,125
+ ,123
+ ,92499
+ ,32
+ ,21399
+ ,5250
+ ,25899
+ ,21
+ ,21
+ ,224330
+ ,83
+ ,130115
+ ,28608
+ ,139296
+ ,167
+ ,163
+ ,135781
+ ,31
+ ,24874
+ ,8092
+ ,52678
+ ,32
+ ,29
+ ,74408
+ ,67
+ ,34988
+ ,4473
+ ,23853
+ ,36
+ ,35
+ ,81240
+ ,66
+ ,45549
+ ,1572
+ ,17306
+ ,13
+ ,13
+ ,14688
+ ,10
+ ,6023
+ ,2065
+ ,7953
+ ,5
+ ,5
+ ,181633
+ ,70
+ ,64466
+ ,14817
+ ,89455
+ ,96
+ ,96
+ ,271856
+ ,103
+ ,54990
+ ,16714
+ ,147866
+ ,151
+ ,151
+ ,7199
+ ,5
+ ,1644
+ ,556
+ ,4245
+ ,6
+ ,6
+ ,46660
+ ,20
+ ,6179
+ ,2089
+ ,21509
+ ,13
+ ,13
+ ,17547
+ ,5
+ ,3926
+ ,2658
+ ,7670
+ ,3
+ ,3
+ ,133368
+ ,36
+ ,32755
+ ,10695
+ ,66675
+ ,57
+ ,56
+ ,95227
+ ,34
+ ,34777
+ ,1669
+ ,14336
+ ,23
+ ,23
+ ,152601
+ ,48
+ ,73224
+ ,16267
+ ,53608
+ ,61
+ ,57
+ ,98146
+ ,40
+ ,27114
+ ,7768
+ ,30059
+ ,21
+ ,14
+ ,79619
+ ,43
+ ,20760
+ ,7252
+ ,29668
+ ,43
+ ,43
+ ,59194
+ ,31
+ ,37636
+ ,6387
+ ,22097
+ ,20
+ ,20
+ ,139942
+ ,42
+ ,65461
+ ,18715
+ ,96841
+ ,82
+ ,72
+ ,118612
+ ,46
+ ,30080
+ ,7936
+ ,41907
+ ,90
+ ,87
+ ,72880
+ ,33
+ ,24094
+ ,8643
+ ,27080
+ ,25
+ ,21
+ ,65475
+ ,18
+ ,69008
+ ,7294
+ ,35885
+ ,60
+ ,56
+ ,99643
+ ,55
+ ,54968
+ ,4570
+ ,41247
+ ,61
+ ,59
+ ,71965
+ ,35
+ ,46090
+ ,7185
+ ,28313
+ ,85
+ ,82
+ ,77272
+ ,59
+ ,27507
+ ,10058
+ ,36845
+ ,43
+ ,43
+ ,49289
+ ,19
+ ,10672
+ ,2342
+ ,16548
+ ,25
+ ,25
+ ,135131
+ ,66
+ ,34029
+ ,8509
+ ,36134
+ ,41
+ ,38
+ ,108446
+ ,60
+ ,46300
+ ,13275
+ ,55764
+ ,26
+ ,25
+ ,89746
+ ,36
+ ,24760
+ ,6816
+ ,28910
+ ,38
+ ,38
+ ,44296
+ ,25
+ ,18779
+ ,1930
+ ,13339
+ ,12
+ ,12
+ ,77648
+ ,47
+ ,21280
+ ,8086
+ ,25319
+ ,29
+ ,29
+ ,181528
+ ,54
+ ,40662
+ ,10737
+ ,66956
+ ,49
+ ,47
+ ,134019
+ ,53
+ ,28987
+ ,8033
+ ,47487
+ ,46
+ ,45
+ ,124064
+ ,40
+ ,22827
+ ,7058
+ ,52785
+ ,41
+ ,40
+ ,92630
+ ,40
+ ,18513
+ ,6782
+ ,44683
+ ,31
+ ,30
+ ,121848
+ ,39
+ ,30594
+ ,5401
+ ,35619
+ ,41
+ ,41
+ ,52915
+ ,14
+ ,24006
+ ,6521
+ ,21920
+ ,26
+ ,25
+ ,81872
+ ,45
+ ,27913
+ ,10856
+ ,45608
+ ,23
+ ,23
+ ,58981
+ ,36
+ ,42744
+ ,2154
+ ,7721
+ ,14
+ ,14
+ ,53515
+ ,28
+ ,12934
+ ,6117
+ ,20634
+ ,16
+ ,16
+ ,60812
+ ,44
+ ,22574
+ ,5238
+ ,29788
+ ,25
+ ,26
+ ,56375
+ ,30
+ ,41385
+ ,4820
+ ,31931
+ ,21
+ ,21
+ ,65490
+ ,22
+ ,18653
+ ,5615
+ ,37754
+ ,32
+ ,27
+ ,80949
+ ,17
+ ,18472
+ ,4272
+ ,32505
+ ,9
+ ,9
+ ,76302
+ ,31
+ ,30976
+ ,8702
+ ,40557
+ ,35
+ ,33
+ ,104011
+ ,55
+ ,63339
+ ,15340
+ ,94238
+ ,42
+ ,42
+ ,98104
+ ,54
+ ,25568
+ ,8030
+ ,44197
+ ,68
+ ,68
+ ,67989
+ ,21
+ ,33747
+ ,9526
+ ,43228
+ ,32
+ ,32
+ ,30989
+ ,14
+ ,4154
+ ,1278
+ ,4103
+ ,6
+ ,6
+ ,135458
+ ,81
+ ,19474
+ ,4236
+ ,44144
+ ,68
+ ,67
+ ,73504
+ ,35
+ ,35130
+ ,3023
+ ,32868
+ ,33
+ ,33
+ ,63123
+ ,43
+ ,39067
+ ,7196
+ ,27640
+ ,84
+ ,77
+ ,61254
+ ,46
+ ,13310
+ ,3394
+ ,14063
+ ,46
+ ,46
+ ,74914
+ ,30
+ ,65892
+ ,6371
+ ,28990
+ ,30
+ ,30
+ ,31774
+ ,23
+ ,4143
+ ,1574
+ ,4694
+ ,0
+ ,0
+ ,81437
+ ,38
+ ,28579
+ ,9620
+ ,42648
+ ,36
+ ,36
+ ,87186
+ ,54
+ ,51776
+ ,6978
+ ,64329
+ ,47
+ ,46
+ ,50090
+ ,20
+ ,21152
+ ,4911
+ ,21928
+ ,20
+ ,18
+ ,65745
+ ,53
+ ,38084
+ ,8645
+ ,25836
+ ,50
+ ,48
+ ,56653
+ ,45
+ ,27717
+ ,8987
+ ,22779
+ ,30
+ ,29
+ ,158399
+ ,39
+ ,32928
+ ,5544
+ ,40820
+ ,30
+ ,28
+ ,46455
+ ,20
+ ,11342
+ ,3083
+ ,27530
+ ,34
+ ,34
+ ,73624
+ ,24
+ ,19499
+ ,6909
+ ,32378
+ ,33
+ ,33
+ ,38395
+ ,31
+ ,16380
+ ,3189
+ ,10824
+ ,34
+ ,34
+ ,91899
+ ,35
+ ,36874
+ ,6745
+ ,39613
+ ,37
+ ,33
+ ,139526
+ ,151
+ ,48259
+ ,16724
+ ,60865
+ ,83
+ ,80
+ ,52164
+ ,52
+ ,16734
+ ,4850
+ ,19787
+ ,32
+ ,32
+ ,51567
+ ,30
+ ,28207
+ ,7025
+ ,20107
+ ,30
+ ,30
+ ,70551
+ ,31
+ ,30143
+ ,6047
+ ,36605
+ ,43
+ ,41
+ ,84856
+ ,29
+ ,41369
+ ,7377
+ ,40961
+ ,41
+ ,41
+ ,102538
+ ,57
+ ,45833
+ ,9078
+ ,48231
+ ,51
+ ,51
+ ,86678
+ ,40
+ ,29156
+ ,4605
+ ,39725
+ ,19
+ ,18
+ ,85709
+ ,44
+ ,35944
+ ,3238
+ ,21455
+ ,37
+ ,34
+ ,34662
+ ,25
+ ,36278
+ ,8100
+ ,23430
+ ,33
+ ,31
+ ,150580
+ ,77
+ ,45588
+ ,9653
+ ,62991
+ ,41
+ ,39
+ ,99611
+ ,35
+ ,45097
+ ,8914
+ ,49363
+ ,54
+ ,54
+ ,19349
+ ,11
+ ,3895
+ ,786
+ ,9604
+ ,14
+ ,14
+ ,99373
+ ,63
+ ,28394
+ ,6700
+ ,24552
+ ,25
+ ,24
+ ,86230
+ ,44
+ ,18632
+ ,5788
+ ,31493
+ ,25
+ ,24
+ ,30837
+ ,19
+ ,2325
+ ,593
+ ,3439
+ ,8
+ ,8
+ ,31706
+ ,13
+ ,25139
+ ,4506
+ ,19555
+ ,26
+ ,26
+ ,89806
+ ,42
+ ,27975
+ ,6382
+ ,21228
+ ,20
+ ,19
+ ,62088
+ ,38
+ ,14483
+ ,5621
+ ,23177
+ ,11
+ ,11
+ ,40151
+ ,29
+ ,13127
+ ,3997
+ ,22094
+ ,14
+ ,14
+ ,27634
+ ,20
+ ,5839
+ ,520
+ ,2342
+ ,3
+ ,1
+ ,76990
+ ,27
+ ,24069
+ ,8891
+ ,38798
+ ,40
+ ,39
+ ,37460
+ ,20
+ ,3738
+ ,999
+ ,3255
+ ,5
+ ,5
+ ,54157
+ ,19
+ ,18625
+ ,7067
+ ,24261
+ ,38
+ ,37
+ ,49862
+ ,37
+ ,36341
+ ,4639
+ ,18511
+ ,32
+ ,32
+ ,84337
+ ,26
+ ,24548
+ ,5654
+ ,40798
+ ,41
+ ,38
+ ,64175
+ ,42
+ ,21792
+ ,6928
+ ,28893
+ ,46
+ ,47
+ ,59382
+ ,49
+ ,26263
+ ,1514
+ ,21425
+ ,47
+ ,47
+ ,119308
+ ,30
+ ,23686
+ ,9238
+ ,50276
+ ,37
+ ,37
+ ,76702
+ ,49
+ ,49303
+ ,8204
+ ,37643
+ ,51
+ ,51
+ ,103425
+ ,67
+ ,25659
+ ,5926
+ ,30377
+ ,49
+ ,45
+ ,70344
+ ,28
+ ,28904
+ ,5785
+ ,27126
+ ,21
+ ,21
+ ,43410
+ ,19
+ ,2781
+ ,4
+ ,13
+ ,1
+ ,1
+ ,104838
+ ,49
+ ,29236
+ ,5930
+ ,42097
+ ,44
+ ,42
+ ,62215
+ ,27
+ ,19546
+ ,3710
+ ,24451
+ ,26
+ ,26
+ ,69304
+ ,30
+ ,22818
+ ,705
+ ,14335
+ ,21
+ ,21
+ ,53117
+ ,22
+ ,32689
+ ,443
+ ,5084
+ ,4
+ ,4
+ ,19764
+ ,12
+ ,5752
+ ,2416
+ ,9927
+ ,10
+ ,10
+ ,86680
+ ,31
+ ,22197
+ ,7747
+ ,43527
+ ,43
+ ,43
+ ,84105
+ ,20
+ ,20055
+ ,5432
+ ,27184
+ ,34
+ ,34
+ ,77945
+ ,20
+ ,25272
+ ,4913
+ ,21610
+ ,32
+ ,31
+ ,89113
+ ,39
+ ,82206
+ ,2650
+ ,20484
+ ,20
+ ,19
+ ,91005
+ ,29
+ ,32073
+ ,2370
+ ,20156
+ ,34
+ ,34
+ ,40248
+ ,16
+ ,5444
+ ,775
+ ,6012
+ ,6
+ ,6
+ ,64187
+ ,27
+ ,20154
+ ,5576
+ ,18475
+ ,12
+ ,11
+ ,50857
+ ,21
+ ,36944
+ ,1352
+ ,12645
+ ,24
+ ,24
+ ,56613
+ ,19
+ ,8019
+ ,3080
+ ,11017
+ ,16
+ ,16
+ ,62792
+ ,35
+ ,30884
+ ,10205
+ ,37623
+ ,72
+ ,72
+ ,72535
+ ,14
+ ,19540
+ ,6095
+ ,35873
+ ,27
+ ,21)
+ ,dim=c(7
+ ,289)
+ ,dimnames=list(c('TimeRFCSEC'
+ ,'#Logins'
+ ,'CW#characters'
+ ,'CW#revisions'
+ ,'CW#seconds'
+ ,'CWIncludedHyperlinks'
+ ,'CWIncludedBlogs')
+ ,1:289))
> y <- array(NA,dim=c(7,289),dimnames=list(c('TimeRFCSEC','#Logins','CW#characters','CW#revisions','CW#seconds','CWIncludedHyperlinks','CWIncludedBlogs'),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 = '3'
> 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
CW#characters TimeRFCSEC #Logins CW#revisions CW#seconds
1 112285 210907 56 24188 146283
2 84786 120982 56 18273 98364
3 83123 176508 54 14130 86146
4 101193 179321 89 32287 96933
5 38361 123185 40 8654 79234
6 68504 52746 25 9245 42551
7 119182 385534 92 33251 195663
8 22807 33170 18 1271 6853
9 17140 101645 63 5279 21529
10 116174 149061 44 27101 95757
11 57635 165446 33 16373 85584
12 66198 237213 84 19716 143983
13 71701 173326 88 17753 75851
14 57793 133131 55 9028 59238
15 80444 258873 60 18653 93163
16 53855 180083 66 8828 96037
17 97668 324799 154 29498 151511
18 133824 230964 53 27563 136368
19 101481 236785 119 18293 112642
20 99645 135473 41 22530 94728
21 114789 202925 61 15977 105499
22 99052 215147 58 35082 121527
23 67654 344297 75 16116 127766
24 65553 153935 33 15849 98958
25 97500 132943 40 16026 77900
26 69112 174724 92 26569 85646
27 82753 174415 100 24785 98579
28 85323 225548 112 17569 130767
29 72654 223632 73 23825 131741
30 30727 124817 40 7869 53907
31 77873 221698 45 14975 178812
32 117478 210767 60 37791 146761
33 74007 170266 62 9605 82036
34 90183 260561 75 27295 163253
35 61542 84853 31 2746 27032
36 101494 294424 77 34461 171975
37 27570 101011 34 8098 65990
38 55813 215641 46 4787 86572
39 79215 325107 99 24919 159676
40 1423 7176 17 603 1929
41 55461 167542 66 16329 85371
42 31081 106408 30 12558 58391
43 22996 96560 76 7784 31580
44 83122 265769 146 28522 136815
45 70106 269651 67 22265 120642
46 60578 149112 56 14459 69107
47 39992 175824 107 14526 50495
48 79892 152871 58 22240 108016
49 49810 111665 34 11802 46341
50 71570 116408 61 7623 78348
51 100708 362301 119 11912 79336
52 33032 78800 42 7935 56968
53 82875 183167 66 18220 93176
54 139077 277965 89 19199 161632
55 71595 150629 44 19918 87850
56 72260 168809 66 21884 127969
57 5950 24188 24 2694 15049
58 115762 329267 259 15808 155135
59 32551 65029 17 3597 25109
60 31701 101097 64 5296 45824
61 80670 218946 41 25239 102996
62 143558 244052 68 29801 160604
63 117105 341570 168 18450 158051
64 23789 103597 43 7132 44547
65 120733 233328 132 34861 162647
66 105195 256462 105 35940 174141
67 73107 206161 71 16688 60622
68 132068 311473 112 24683 179566
69 149193 235800 94 46230 184301
70 46821 177939 82 10387 75661
71 87011 207176 70 21436 96144
72 95260 196553 57 30546 129847
73 55183 174184 53 19746 117286
74 106671 143246 103 15977 71180
75 73511 187559 121 22583 109377
76 92945 187681 62 17274 85298
77 78664 119016 52 16469 73631
78 70054 182192 52 14251 86767
79 22618 73566 32 3007 23824
80 74011 194979 62 16851 93487
81 83737 167488 45 21113 82981
82 69094 143756 46 17401 73815
83 93133 275541 63 23958 94552
84 95536 243199 75 23567 132190
85 225920 182999 88 13065 128754
86 62133 135649 46 15358 66363
87 61370 152299 53 14587 67808
88 43836 120221 37 12770 61724
89 106117 346485 90 24021 131722
90 38692 145790 63 9648 68580
91 84651 193339 78 20537 106175
92 56622 80953 25 7905 55792
93 15986 122774 45 4527 25157
94 95364 130585 46 30495 76669
95 26706 112611 41 7117 57283
96 89691 286468 144 17719 105805
97 67267 241066 82 27056 129484
98 126846 148446 91 33473 72413
99 41140 204713 71 9758 87831
100 102860 182079 63 21115 96971
101 51715 140344 53 7236 71299
102 55801 220516 62 13790 77494
103 111813 243060 63 32902 120336
104 120293 162765 32 25131 93913
105 138599 182613 39 30910 136048
106 161647 232138 62 35947 181248
107 115929 265318 117 29848 146123
108 24266 85574 34 6943 32036
109 162901 310839 92 42705 186646
110 109825 225060 93 31808 102255
111 129838 232317 54 26675 168237
112 37510 144966 144 8435 64219
113 43750 43287 14 7409 19630
114 40652 155754 61 14993 76825
115 87771 164709 109 36867 115338
116 85872 201940 38 33835 109427
117 89275 235454 73 24164 118168
118 44418 220801 75 12607 84845
119 192565 99466 50 22609 153197
120 35232 92661 61 5892 29877
121 40909 133328 55 17014 63506
122 13294 61361 77 5394 22445
123 32387 125930 75 9178 47695
124 140867 100750 72 6440 68370
125 120662 224549 50 21916 146304
126 21233 82316 32 4011 38233
127 44332 102010 53 5818 42071
128 61056 101523 42 18647 50517
129 101338 243511 71 20556 103950
130 1168 22938 10 238 5841
131 13497 41566 35 70 2341
132 65567 152474 65 22392 84396
133 25162 61857 25 3913 24610
134 32334 99923 66 12237 35753
135 40735 132487 41 8388 55515
136 91413 317394 86 22120 209056
137 855 21054 16 338 6622
138 97068 209641 42 11727 115814
139 44339 22648 19 3704 11609
140 14116 31414 19 3988 13155
141 10288 46698 45 3030 18274
142 65622 131698 65 13520 72875
143 16563 91735 35 1421 10112
144 76643 244749 95 20923 142775
145 110681 184510 49 20237 68847
146 29011 79863 37 3219 17659
147 92696 128423 64 3769 20112
148 94785 97839 38 12252 61023
149 8773 38214 34 1888 13983
150 83209 151101 32 14497 65176
151 93815 272458 65 28864 132432
152 86687 172494 52 21721 112494
153 34553 108043 62 4821 45109
154 105547 328107 65 33644 170875
155 103487 250579 83 15923 180759
156 213688 351067 95 42935 214921
157 71220 158015 29 18864 100226
158 23517 98866 18 4977 32043
159 56926 85439 33 7785 54454
160 91721 229242 247 17939 78876
161 115168 351619 139 23436 170745
162 111194 84207 29 325 6940
163 51009 120445 118 13539 49025
164 135777 324598 110 34538 122037
165 51513 131069 67 12198 53782
166 74163 204271 42 26924 127748
167 51633 165543 65 12716 86839
168 75345 141722 94 8172 44830
169 33416 116048 64 10855 77395
170 83305 250047 81 11932 89324
171 98952 299775 95 14300 103300
172 102372 195838 67 25515 112283
173 37238 173260 63 2805 10901
174 103772 254488 83 29402 120691
175 123969 104389 45 16440 58106
176 27142 136084 30 11221 57140
177 135400 199476 70 28732 122422
178 21399 92499 32 5250 25899
179 130115 224330 83 28608 139296
180 24874 135781 31 8092 52678
181 34988 74408 67 4473 23853
182 45549 81240 66 1572 17306
183 6023 14688 10 2065 7953
184 64466 181633 70 14817 89455
185 54990 271856 103 16714 147866
186 1644 7199 5 556 4245
187 6179 46660 20 2089 21509
188 3926 17547 5 2658 7670
189 32755 133368 36 10695 66675
190 34777 95227 34 1669 14336
191 73224 152601 48 16267 53608
192 27114 98146 40 7768 30059
193 20760 79619 43 7252 29668
194 37636 59194 31 6387 22097
195 65461 139942 42 18715 96841
196 30080 118612 46 7936 41907
197 24094 72880 33 8643 27080
198 69008 65475 18 7294 35885
199 54968 99643 55 4570 41247
200 46090 71965 35 7185 28313
201 27507 77272 59 10058 36845
202 10672 49289 19 2342 16548
203 34029 135131 66 8509 36134
204 46300 108446 60 13275 55764
205 24760 89746 36 6816 28910
206 18779 44296 25 1930 13339
207 21280 77648 47 8086 25319
208 40662 181528 54 10737 66956
209 28987 134019 53 8033 47487
210 22827 124064 40 7058 52785
211 18513 92630 40 6782 44683
212 30594 121848 39 5401 35619
213 24006 52915 14 6521 21920
214 27913 81872 45 10856 45608
215 42744 58981 36 2154 7721
216 12934 53515 28 6117 20634
217 22574 60812 44 5238 29788
218 41385 56375 30 4820 31931
219 18653 65490 22 5615 37754
220 18472 80949 17 4272 32505
221 30976 76302 31 8702 40557
222 63339 104011 55 15340 94238
223 25568 98104 54 8030 44197
224 33747 67989 21 9526 43228
225 4154 30989 14 1278 4103
226 19474 135458 81 4236 44144
227 35130 73504 35 3023 32868
228 39067 63123 43 7196 27640
229 13310 61254 46 3394 14063
230 65892 74914 30 6371 28990
231 4143 31774 23 1574 4694
232 28579 81437 38 9620 42648
233 51776 87186 54 6978 64329
234 21152 50090 20 4911 21928
235 38084 65745 53 8645 25836
236 27717 56653 45 8987 22779
237 32928 158399 39 5544 40820
238 11342 46455 20 3083 27530
239 19499 73624 24 6909 32378
240 16380 38395 31 3189 10824
241 36874 91899 35 6745 39613
242 48259 139526 151 16724 60865
243 16734 52164 52 4850 19787
244 28207 51567 30 7025 20107
245 30143 70551 31 6047 36605
246 41369 84856 29 7377 40961
247 45833 102538 57 9078 48231
248 29156 86678 40 4605 39725
249 35944 85709 44 3238 21455
250 36278 34662 25 8100 23430
251 45588 150580 77 9653 62991
252 45097 99611 35 8914 49363
253 3895 19349 11 786 9604
254 28394 99373 63 6700 24552
255 18632 86230 44 5788 31493
256 2325 30837 19 593 3439
257 25139 31706 13 4506 19555
258 27975 89806 42 6382 21228
259 14483 62088 38 5621 23177
260 13127 40151 29 3997 22094
261 5839 27634 20 520 2342
262 24069 76990 27 8891 38798
263 3738 37460 20 999 3255
264 18625 54157 19 7067 24261
265 36341 49862 37 4639 18511
266 24548 84337 26 5654 40798
267 21792 64175 42 6928 28893
268 26263 59382 49 1514 21425
269 23686 119308 30 9238 50276
270 49303 76702 49 8204 37643
271 25659 103425 67 5926 30377
272 28904 70344 28 5785 27126
273 2781 43410 19 4 13
274 29236 104838 49 5930 42097
275 19546 62215 27 3710 24451
276 22818 69304 30 705 14335
277 32689 53117 22 443 5084
278 5752 19764 12 2416 9927
279 22197 86680 31 7747 43527
280 20055 84105 20 5432 27184
281 25272 77945 20 4913 21610
282 82206 89113 39 2650 20484
283 32073 91005 29 2370 20156
284 5444 40248 16 775 6012
285 20154 64187 27 5576 18475
286 36944 50857 21 1352 12645
287 8019 56613 19 3080 11017
288 30884 62792 35 10205 37623
289 19540 72535 14 6095 35873
CWIncludedHyperlinks CWIncludedBlogs
1 144 145
2 103 101
3 98 98
4 135 132
5 61 60
6 39 38
7 150 144
8 5 5
9 28 28
10 84 84
11 80 79
12 130 127
13 82 78
14 60 60
15 131 131
16 84 84
17 140 133
18 151 150
19 91 91
20 138 132
21 150 136
22 124 124
23 119 118
24 73 70
25 110 107
26 123 119
27 90 89
28 116 112
29 113 108
30 56 52
31 115 112
32 119 116
33 129 123
34 127 125
35 27 27
36 175 162
37 35 32
38 64 64
39 96 92
40 0 0
41 84 83
42 41 41
43 47 47
44 126 120
45 105 105
46 80 79
47 70 65
48 73 70
49 57 55
50 40 39
51 68 67
52 21 21
53 127 127
54 154 152
55 116 113
56 102 99
57 7 7
58 148 141
59 21 21
60 35 35
61 112 109
62 137 133
63 135 123
64 26 26
65 230 230
66 181 166
67 71 68
68 147 147
69 190 179
70 64 61
71 105 101
72 107 108
73 94 90
74 116 114
75 106 103
76 143 142
77 81 79
78 89 88
79 26 25
80 84 83
81 113 113
82 120 118
83 110 110
84 134 129
85 54 51
86 96 93
87 78 76
88 51 49
89 121 118
90 38 38
91 145 141
92 59 58
93 27 27
94 91 91
95 48 48
96 68 63
97 58 56
98 150 144
99 74 73
100 181 168
101 65 64
102 97 97
103 121 117
104 99 100
105 152 149
106 188 187
107 138 127
108 40 37
109 254 245
110 87 87
111 178 177
112 51 49
113 49 49
114 73 73
115 176 177
116 94 94
117 120 117
118 66 60
119 56 55
120 39 39
121 66 64
122 27 26
123 65 64
124 58 58
125 98 95
126 25 25
127 26 26
128 77 76
129 130 129
130 11 11
131 2 2
132 101 101
133 31 28
134 36 36
135 120 89
136 195 193
137 4 4
138 89 84
139 24 23
140 39 39
141 14 14
142 78 78
143 15 14
144 106 101
145 83 82
146 24 24
147 37 36
148 77 75
149 16 16
150 56 55
151 132 131
152 144 131
153 40 39
154 153 144
155 143 139
156 220 211
157 79 78
158 50 50
159 39 39
160 95 90
161 169 166
162 12 12
163 63 57
164 134 133
165 69 69
166 119 119
167 119 119
168 75 65
169 63 61
170 55 49
171 103 101
172 197 196
173 16 15
174 140 136
175 89 89
176 40 40
177 125 123
178 21 21
179 167 163
180 32 29
181 36 35
182 13 13
183 5 5
184 96 96
185 151 151
186 6 6
187 13 13
188 3 3
189 57 56
190 23 23
191 61 57
192 21 14
193 43 43
194 20 20
195 82 72
196 90 87
197 25 21
198 60 56
199 61 59
200 85 82
201 43 43
202 25 25
203 41 38
204 26 25
205 38 38
206 12 12
207 29 29
208 49 47
209 46 45
210 41 40
211 31 30
212 41 41
213 26 25
214 23 23
215 14 14
216 16 16
217 25 26
218 21 21
219 32 27
220 9 9
221 35 33
222 42 42
223 68 68
224 32 32
225 6 6
226 68 67
227 33 33
228 84 77
229 46 46
230 30 30
231 0 0
232 36 36
233 47 46
234 20 18
235 50 48
236 30 29
237 30 28
238 34 34
239 33 33
240 34 34
241 37 33
242 83 80
243 32 32
244 30 30
245 43 41
246 41 41
247 51 51
248 19 18
249 37 34
250 33 31
251 41 39
252 54 54
253 14 14
254 25 24
255 25 24
256 8 8
257 26 26
258 20 19
259 11 11
260 14 14
261 3 1
262 40 39
263 5 5
264 38 37
265 32 32
266 41 38
267 46 47
268 47 47
269 37 37
270 51 51
271 49 45
272 21 21
273 1 1
274 44 42
275 26 26
276 21 21
277 4 4
278 10 10
279 43 43
280 34 34
281 32 31
282 20 19
283 34 34
284 6 6
285 12 11
286 24 24
287 16 16
288 72 72
289 27 21
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) TimeRFCSEC `#Logins`
8317.09924 -0.00668 60.22712
`CW#revisions` `CW#seconds` CWIncludedHyperlinks
1.22430 0.32161 -63.12720
CWIncludedBlogs
205.58252
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-47244 -12105 -4223 7904 149045
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 8317.09924 2652.01414 3.136 0.00189 **
TimeRFCSEC -0.00668 0.04002 -0.167 0.86756
`#Logins` 60.22712 54.97897 1.095 0.27425
`CW#revisions` 1.22430 0.29898 4.095 5.52e-05 ***
`CW#seconds` 0.32161 0.07976 4.032 7.12e-05 ***
CWIncludedHyperlinks -63.12720 435.95899 -0.145 0.88497
CWIncludedBlogs 205.58252 448.33698 0.459 0.64691
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 21710 on 282 degrees of freedom
Multiple R-squared: 0.7175, Adjusted R-squared: 0.7115
F-statistic: 119.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.2327052375 4.654105e-01 7.672948e-01
[2,] 0.3524972724 7.049945e-01 6.475027e-01
[3,] 0.2625788711 5.251577e-01 7.374211e-01
[4,] 0.2319467886 4.638936e-01 7.680532e-01
[5,] 0.1765841870 3.531684e-01 8.234158e-01
[6,] 0.1103319897 2.206640e-01 8.896680e-01
[7,] 0.0673111790 1.346224e-01 9.326888e-01
[8,] 0.0457700319 9.154006e-02 9.542300e-01
[9,] 0.0484823473 9.696469e-02 9.515177e-01
[10,] 0.0488922471 9.778449e-02 9.511078e-01
[11,] 0.0532880602 1.065761e-01 9.467119e-01
[12,] 0.1649852733 3.299705e-01 8.350147e-01
[13,] 0.1658138208 3.316276e-01 8.341862e-01
[14,] 0.1238027781 2.476056e-01 8.761972e-01
[15,] 0.0925239633 1.850479e-01 9.074760e-01
[16,] 0.0778877803 1.557756e-01 9.221122e-01
[17,] 0.1135013875 2.270028e-01 8.864986e-01
[18,] 0.0825840854 1.651682e-01 9.174159e-01
[19,] 0.0591100467 1.182201e-01 9.408900e-01
[20,] 0.0652947066 1.305894e-01 9.347053e-01
[21,] 0.0610310078 1.220620e-01 9.389690e-01
[22,] 0.0503739713 1.007479e-01 9.496260e-01
[23,] 0.0355736881 7.114738e-02 9.644263e-01
[24,] 0.0249407676 4.988154e-02 9.750592e-01
[25,] 0.0202369435 4.047389e-02 9.797631e-01
[26,] 0.0346207077 6.924142e-02 9.653793e-01
[27,] 0.0356945787 7.138916e-02 9.643054e-01
[28,] 0.0321269398 6.425388e-02 9.678731e-01
[29,] 0.0249717261 4.994345e-02 9.750283e-01
[30,] 0.0194162855 3.883257e-02 9.805837e-01
[31,] 0.0209530938 4.190619e-02 9.790469e-01
[32,] 0.0181301912 3.626038e-02 9.818698e-01
[33,] 0.0177556363 3.551127e-02 9.822444e-01
[34,] 0.0176030883 3.520618e-02 9.823969e-01
[35,] 0.0139451861 2.789037e-02 9.860548e-01
[36,] 0.0118460561 2.369211e-02 9.881539e-01
[37,] 0.0083929559 1.678591e-02 9.916070e-01
[38,] 0.0063254038 1.265081e-02 9.936746e-01
[39,] 0.0047327577 9.465515e-03 9.952672e-01
[40,] 0.0032227907 6.445581e-03 9.967772e-01
[41,] 0.0054852039 1.097041e-02 9.945148e-01
[42,] 0.0352248832 7.044977e-02 9.647751e-01
[43,] 0.0270670551 5.413411e-02 9.729329e-01
[44,] 0.0207622798 4.152456e-02 9.792377e-01
[45,] 0.0372823671 7.456473e-02 9.627176e-01
[46,] 0.0316495575 6.329912e-02 9.683504e-01
[47,] 0.0269393633 5.387873e-02 9.730606e-01
[48,] 0.0238262590 4.765252e-02 9.761737e-01
[49,] 0.0218108868 4.362177e-02 9.781891e-01
[50,] 0.0168492279 3.369846e-02 9.831508e-01
[51,] 0.0130091713 2.601834e-02 9.869908e-01
[52,] 0.0099034597 1.980692e-02 9.900965e-01
[53,] 0.0192334734 3.846695e-02 9.807665e-01
[54,] 0.0196521449 3.930429e-02 9.803479e-01
[55,] 0.0161971546 3.239431e-02 9.838028e-01
[56,] 0.0202015914 4.040318e-02 9.797984e-01
[57,] 0.0212438809 4.248776e-02 9.787561e-01
[58,] 0.0177963490 3.559270e-02 9.822037e-01
[59,] 0.0163012382 3.260248e-02 9.836988e-01
[60,] 0.0145531415 2.910628e-02 9.854469e-01
[61,] 0.0118228822 2.364576e-02 9.881771e-01
[62,] 0.0091378713 1.827574e-02 9.908621e-01
[63,] 0.0070477712 1.409554e-02 9.929522e-01
[64,] 0.0082651096 1.653022e-02 9.917349e-01
[65,] 0.0116838238 2.336765e-02 9.883162e-01
[66,] 0.0102860114 2.057202e-02 9.897140e-01
[67,] 0.0081238862 1.624777e-02 9.918761e-01
[68,] 0.0069878133 1.397563e-02 9.930122e-01
[69,] 0.0052286468 1.045729e-02 9.947714e-01
[70,] 0.0039623450 7.924690e-03 9.960377e-01
[71,] 0.0029310737 5.862147e-03 9.970689e-01
[72,] 0.0021525895 4.305179e-03 9.978474e-01
[73,] 0.0017023663 3.404733e-03 9.982976e-01
[74,] 0.0012566820 2.513364e-03 9.987433e-01
[75,] 0.0009118334 1.823667e-03 9.990882e-01
[76,] 0.9605242534 7.895149e-02 3.947575e-02
[77,] 0.9519984056 9.600319e-02 4.800159e-02
[78,] 0.9420199429 1.159601e-01 5.798006e-02
[79,] 0.9324359619 1.351281e-01 6.756404e-02
[80,] 0.9216561659 1.566877e-01 7.834383e-02
[81,] 0.9128727533 1.742545e-01 8.712725e-02
[82,] 0.8987673155 2.024654e-01 1.012327e-01
[83,] 0.8849495618 2.301009e-01 1.150504e-01
[84,] 0.8737400162 2.525200e-01 1.262600e-01
[85,] 0.8636425159 2.727150e-01 1.363575e-01
[86,] 0.8590987461 2.818025e-01 1.409013e-01
[87,] 0.8431129302 3.137741e-01 1.568871e-01
[88,] 0.8516165460 2.967669e-01 1.483835e-01
[89,] 0.8712319297 2.575361e-01 1.287681e-01
[90,] 0.8704564483 2.590871e-01 1.295436e-01
[91,] 0.8559066772 2.881866e-01 1.440933e-01
[92,] 0.8349023539 3.301953e-01 1.650976e-01
[93,] 0.8177544255 3.644911e-01 1.822456e-01
[94,] 0.7990919653 4.018161e-01 2.009080e-01
[95,] 0.8384317673 3.231365e-01 1.615682e-01
[96,] 0.8482597080 3.034806e-01 1.517403e-01
[97,] 0.8479262985 3.041474e-01 1.520737e-01
[98,] 0.8277435875 3.445128e-01 1.722564e-01
[99,] 0.8094637544 3.810725e-01 1.905362e-01
[100,] 0.7858713874 4.282572e-01 2.141286e-01
[101,] 0.7700165638 4.599669e-01 2.299834e-01
[102,] 0.7465048520 5.069903e-01 2.534951e-01
[103,] 0.7373237166 5.253526e-01 2.626763e-01
[104,] 0.7174671118 5.650658e-01 2.825329e-01
[105,] 0.7247148531 5.505703e-01 2.752851e-01
[106,] 0.7645139889 4.709720e-01 2.354860e-01
[107,] 0.7510499862 4.979000e-01 2.489500e-01
[108,] 0.7257535706 5.484929e-01 2.742464e-01
[109,] 0.7179476438 5.641047e-01 2.820524e-01
[110,] 0.9871645924 2.567082e-02 1.283541e-02
[111,] 0.9840211758 3.195765e-02 1.597882e-02
[112,] 0.9837402325 3.251953e-02 1.625977e-02
[113,] 0.9822494422 3.550112e-02 1.775056e-02
[114,] 0.9805050796 3.898984e-02 1.949492e-02
[115,] 0.9998626692 2.746615e-04 1.373308e-04
[116,] 0.9998996129 2.007742e-04 1.003871e-04
[117,] 0.9998685911 2.628178e-04 1.314089e-04
[118,] 0.9998319129 3.361742e-04 1.680871e-04
[119,] 0.9997698142 4.603716e-04 2.301858e-04
[120,] 0.9997112934 5.774133e-04 2.887066e-04
[121,] 0.9996411294 7.177411e-04 3.588706e-04
[122,] 0.9995131936 9.736128e-04 4.868064e-04
[123,] 0.9994472406 1.105519e-03 5.527594e-04
[124,] 0.9992562141 1.487572e-03 7.437859e-04
[125,] 0.9991082191 1.783562e-03 8.917809e-04
[126,] 0.9990896589 1.820682e-03 9.103411e-04
[127,] 0.9994904833 1.019033e-03 5.095167e-04
[128,] 0.9993668998 1.266200e-03 6.331002e-04
[129,] 0.9994990512 1.001898e-03 5.009488e-04
[130,] 0.9995454717 9.090567e-04 4.545283e-04
[131,] 0.9994302962 1.139408e-03 5.697038e-04
[132,] 0.9992917421 1.416516e-03 7.082579e-04
[133,] 0.9990765201 1.846960e-03 9.234799e-04
[134,] 0.9987930740 2.413852e-03 1.206926e-03
[135,] 0.9986736503 2.652699e-03 1.326350e-03
[136,] 0.9992886619 1.422676e-03 7.113381e-04
[137,] 0.9990576538 1.884692e-03 9.423462e-04
[138,] 0.9999118825 1.762350e-04 8.811751e-05
[139,] 0.9999686078 6.278449e-05 3.139224e-05
[140,] 0.9999585784 8.284324e-05 4.142162e-05
[141,] 0.9999669554 6.608918e-05 3.304459e-05
[142,] 0.9999611503 7.769938e-05 3.884969e-05
[143,] 0.9999475680 1.048640e-04 5.243198e-05
[144,] 0.9999262129 1.475741e-04 7.378707e-05
[145,] 0.9999478116 1.043769e-04 5.218843e-05
[146,] 0.9999425651 1.148699e-04 5.743494e-05
[147,] 0.9999922153 1.556944e-05 7.784722e-06
[148,] 0.9999890059 2.198827e-05 1.099413e-05
[149,] 0.9999853778 2.924432e-05 1.462216e-05
[150,] 0.9999872798 2.544041e-05 1.272021e-05
[151,] 0.9999824918 3.501640e-05 1.750820e-05
[152,] 0.9999748055 5.038900e-05 2.519450e-05
[153,] 0.9999999990 2.010129e-09 1.005065e-09
[154,] 0.9999999985 3.089514e-09 1.544757e-09
[155,] 0.9999999980 4.068744e-09 2.034372e-09
[156,] 0.9999999965 6.931384e-09 3.465692e-09
[157,] 0.9999999969 6.163541e-09 3.081771e-09
[158,] 0.9999999962 7.607927e-09 3.803963e-09
[159,] 0.9999999972 5.538104e-09 2.769052e-09
[160,] 0.9999999968 6.339797e-09 3.169899e-09
[161,] 0.9999999975 4.908341e-09 2.454171e-09
[162,] 0.9999999983 3.458822e-09 1.729411e-09
[163,] 0.9999999975 5.051249e-09 2.525625e-09
[164,] 0.9999999967 6.627309e-09 3.313654e-09
[165,] 0.9999999951 9.889809e-09 4.944905e-09
[166,] 1.0000000000 4.942995e-11 2.471497e-11
[167,] 1.0000000000 5.089820e-11 2.544910e-11
[168,] 1.0000000000 4.688743e-12 2.344372e-12
[169,] 1.0000000000 8.620920e-12 4.310460e-12
[170,] 1.0000000000 1.352291e-12 6.761454e-13
[171,] 1.0000000000 1.870084e-12 9.350419e-13
[172,] 1.0000000000 3.418548e-12 1.709274e-12
[173,] 1.0000000000 1.768967e-12 8.844836e-13
[174,] 1.0000000000 3.013880e-12 1.506940e-12
[175,] 1.0000000000 4.954204e-12 2.477102e-12
[176,] 1.0000000000 1.607975e-12 8.039873e-13
[177,] 1.0000000000 2.657570e-12 1.328785e-12
[178,] 1.0000000000 3.641104e-12 1.820552e-12
[179,] 1.0000000000 5.810302e-12 2.905151e-12
[180,] 1.0000000000 7.318079e-12 3.659039e-12
[181,] 1.0000000000 9.557559e-12 4.778780e-12
[182,] 1.0000000000 4.564369e-12 2.282184e-12
[183,] 1.0000000000 8.936866e-12 4.468433e-12
[184,] 1.0000000000 1.412313e-11 7.061563e-12
[185,] 1.0000000000 1.731946e-11 8.659728e-12
[186,] 1.0000000000 3.379641e-11 1.689820e-11
[187,] 1.0000000000 4.147458e-11 2.073729e-11
[188,] 1.0000000000 7.774215e-11 3.887107e-11
[189,] 1.0000000000 8.562368e-12 4.281184e-12
[190,] 1.0000000000 7.118840e-12 3.559420e-12
[191,] 1.0000000000 9.160790e-12 4.580395e-12
[192,] 1.0000000000 1.611678e-11 8.058388e-12
[193,] 1.0000000000 2.668339e-11 1.334169e-11
[194,] 1.0000000000 5.395642e-11 2.697821e-11
[195,] 1.0000000000 9.747201e-11 4.873600e-11
[196,] 0.9999999999 1.863878e-10 9.319388e-11
[197,] 0.9999999998 3.718996e-10 1.859498e-10
[198,] 0.9999999997 6.461805e-10 3.230902e-10
[199,] 0.9999999994 1.167642e-09 5.838211e-10
[200,] 0.9999999991 1.819853e-09 9.099263e-10
[201,] 0.9999999990 1.995846e-09 9.979228e-10
[202,] 0.9999999989 2.213117e-09 1.106558e-09
[203,] 0.9999999979 4.209185e-09 2.104592e-09
[204,] 0.9999999960 8.073891e-09 4.036946e-09
[205,] 0.9999999932 1.352571e-08 6.762853e-09
[206,] 0.9999999970 5.966044e-09 2.983022e-09
[207,] 0.9999999953 9.350049e-09 4.675025e-09
[208,] 0.9999999914 1.726017e-08 8.630086e-09
[209,] 0.9999999906 1.874377e-08 9.371886e-09
[210,] 0.9999999858 2.832017e-08 1.416009e-08
[211,] 0.9999999761 4.787535e-08 2.393767e-08
[212,] 0.9999999551 8.987382e-08 4.493691e-08
[213,] 0.9999999327 1.346651e-07 6.733256e-08
[214,] 0.9999999162 1.676182e-07 8.380912e-08
[215,] 0.9999998456 3.088424e-07 1.544212e-07
[216,] 0.9999997493 5.014704e-07 2.507352e-07
[217,] 0.9999998977 2.045846e-07 1.022923e-07
[218,] 0.9999998187 3.625840e-07 1.812920e-07
[219,] 0.9999996720 6.560583e-07 3.280292e-07
[220,] 0.9999996143 7.714589e-07 3.857295e-07
[221,] 0.9999999837 3.257289e-08 1.628645e-08
[222,] 0.9999999704 5.913276e-08 2.956638e-08
[223,] 0.9999999415 1.170382e-07 5.851908e-08
[224,] 0.9999999389 1.222663e-07 6.113316e-08
[225,] 0.9999998792 2.416428e-07 1.208214e-07
[226,] 0.9999997900 4.200922e-07 2.100461e-07
[227,] 0.9999996119 7.761050e-07 3.880525e-07
[228,] 0.9999994574 1.085148e-06 5.425739e-07
[229,] 0.9999991969 1.606236e-06 8.031180e-07
[230,] 0.9999986488 2.702394e-06 1.351197e-06
[231,] 0.9999975093 4.981371e-06 2.490685e-06
[232,] 0.9999956866 8.626722e-06 4.313361e-06
[233,] 0.9999932773 1.344544e-05 6.722720e-06
[234,] 0.9999897422 2.051567e-05 1.025784e-05
[235,] 0.9999821455 3.570893e-05 1.785446e-05
[236,] 0.9999672394 6.552125e-05 3.276062e-05
[237,] 0.9999535515 9.289698e-05 4.644849e-05
[238,] 0.9999222954 1.554092e-04 7.770458e-05
[239,] 0.9998651931 2.696139e-04 1.348069e-04
[240,] 0.9997676914 4.646171e-04 2.323086e-04
[241,] 0.9998423048 3.153903e-04 1.576952e-04
[242,] 0.9997188937 5.622126e-04 2.811063e-04
[243,] 0.9995820752 8.358496e-04 4.179248e-04
[244,] 0.9993197410 1.360518e-03 6.802590e-04
[245,] 0.9988766672 2.246666e-03 1.123333e-03
[246,] 0.9984080149 3.183970e-03 1.591985e-03
[247,] 0.9978311596 4.337681e-03 2.168840e-03
[248,] 0.9973418053 5.316389e-03 2.658195e-03
[249,] 0.9955891394 8.821721e-03 4.410861e-03
[250,] 0.9934435319 1.311294e-02 6.556468e-03
[251,] 0.9895914805 2.081704e-02 1.040852e-02
[252,] 0.9840829356 3.183413e-02 1.591706e-02
[253,] 0.9754589917 4.908202e-02 2.454101e-02
[254,] 0.9695816839 6.083663e-02 3.041832e-02
[255,] 0.9543943234 9.121135e-02 4.560568e-02
[256,] 0.9424541727 1.150917e-01 5.754583e-02
[257,] 0.9159013203 1.681974e-01 8.409868e-02
[258,] 0.8852381170 2.295238e-01 1.147619e-01
[259,] 0.8518757489 2.962485e-01 1.481243e-01
[260,] 0.8231380548 3.537239e-01 1.768619e-01
[261,] 0.8055356035 3.889288e-01 1.944644e-01
[262,] 0.8620025597 2.759949e-01 1.379974e-01
[263,] 0.8129962791 3.740074e-01 1.870037e-01
[264,] 0.8201835444 3.596329e-01 1.798165e-01
[265,] 0.9212235199 1.575530e-01 7.877648e-02
[266,] 0.8751726428 2.496547e-01 1.248274e-01
[267,] 0.9431028399 1.137943e-01 5.689716e-02
[268,] 0.8864660184 2.270680e-01 1.135340e-01
[269,] 0.8259290336 3.481419e-01 1.740710e-01
[270,] 0.7238855672 5.522289e-01 2.761144e-01
> postscript(file="/var/wessaorg/rcomp/tmp/1w2t91353339095.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/2ogyx1353339095.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/3strj1353339095.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/4yfjx1353339095.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/5wtzv1353339095.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 6
4625.14380 5635.92984 13767.13349 -604.92763 -16104.21531 28679.82767
7 8 9 10 11 12
-15872.21420 9155.01849 -11668.23176 30260.05356 -10325.56574 -33941.02040
13 14 15 16 17 18
2253.14136 8400.76515 -1218.31091 -10895.17995 -21101.24635 24949.56662
19 20 21 22 23 24
15992.01269 13289.07841 32173.04387 -11020.95753 -20448.75420 -4735.81973
25 26 27 28 29 30
27934.33201 -20351.69750 -5085.56771 -7501.34764 -25173.79613 -13291.73719
31 32 33 34 35 36
-23281.01854 -2847.83656 7806.77379 -24512.65172 36022.63346 -29250.76938
37 38 39 40 41 42
-17626.82621 3345.44351 -27608.31829 -9228.66658 -14920.50815 -18326.77924
43 44 45 46 47 48
-15635.19511 -27849.62532 -21461.81694 -1234.38629 -16562.82744 -2647.25621
49 50 51 52 53 54
3129.30786 20333.50610 38063.37819 -8316.23709 1441.37335 30241.37970
55 56 57 58 59 60
-6913.17620 -20767.09543 -12786.36544 5154.21917 8173.70479 -6002.69443
61 62 63 64 65 66
-8016.94705 25944.13209 10767.48650 -13188.23825 -21729.71406 -30440.62937
67 68 69 70 71 72
12465.47829 10175.16542 -3888.18463 -10796.78614 4561.30496 -9783.18207
73 74 75 76 77 78
-29626.67899 34540.75917 -18149.52387 13400.58101 13038.79612 1996.37577
79 80 81 82 83 84
-1976.76158 804.38879 5194.72829 -2760.55805 7451.29328 -5101.64026
85 86 87 88 89 90
149045.16614 -1253.30656 511.14181 -8245.97212 6301.33852 -11727.04098
91 92 93 94 95 96
-6196.60669 11519.23145 -11700.68538 10192.66568 -17302.27954 10234.21204
97 98 99 100 101 102
-26997.91242 29635.21227 -20616.09690 11814.91185 299.68127 -10401.40368
103 104 105 106 107 108
5927.07008 35855.82849 26518.66477 22269.04063 1402.40997 -9412.12861
109 110 111 112 113 114
4474.69022 13187.62407 7903.79268 -16346.07118 12514.49371 -23761.52858
115 116 117 118 119 120
-33518.63587 -13392.76826 -5931.96874 -17831.67157 97178.85390 1481.88590
121 122 123 124 125 126
-20075.36766 -16713.86090 -15235.84584 90751.02751 23604.75963 -9229.74759
127 128 129 130 131 132
9146.92882 1047.68360 13459.53913 -11335.07485 2226.10522 -14591.60152
133 134 135 136 137 138
-752.50006 -10899.33625 -8011.70449 -41647.54451 -11398.44447 24366.58557
139 140 141 142 143 144
23547.13143 -9804.63758 -12008.52007 3168.35231 -172.38010 -21367.30493
145 146 147 148 149 150
42108.92355 5959.68447 65234.31622 39649.08768 -10424.41851 27492.01397
151 152 153 154 155 156
-13125.27327 -4223.08502 -2679.04094 -20584.33065 -5332.70462 50818.03114
157 158 159 160 161 162
-4165.61552 -8745.28764 14592.13791 10223.61335 -6236.41950 97353.46398
163 164 165 166 167 168
-3694.33435 22586.38705 -2024.13781 -26319.63044 -19941.89321 29262.08157
169 170 171 172 173 174
-24724.90277 21842.23541 21924.05988 -3879.80124 17270.21023 -1777.91575
175 176 177 178 179 180
62145.40836 -19885.86603 32254.79094 -5976.04070 15505.82840 -15193.86536
181 182 183 184 185 186
5062.21370 24457.25208 -8596.48978 -7439.66290 -47243.71062 -9826.83191
187 188 189 190 191 192
-14357.99297 -10723.33760 -19291.12376 15117.83161 18011.27731 -3686.74712
193 194 195 196 197 198
-14160.79509 10071.90947 -8134.34000 -15613.29101 -7753.67920 31848.10582
199 200 201 202 203 204
16864.82388 6751.33027 -14136.69830 -10210.89591 -4623.02995 -2591.51207
205 206 207 208 209 210
-8181.70460 889.77237 -11522.86686 -10943.07281 -13081.44780 -18322.93053
211 212 213 214 215 216
-18478.74542 -3166.62835 -3332.45294 -13802.96372 25538.03500 -13116.44412
217 218 219 220 221 222
-7746.89756 12475.58772 -13098.86912 -7294.53272 -6970.68896 -2667.76182
223 224 225 226 227 228
-19078.40790 -5504.59552 -8538.23321 -21681.40358 6223.12133 355.21231
229 230 231 232 233 234
-12599.41796 34871.30668 -8783.76660 -12104.99425 5067.06339 -3537.82604
235 236 237 238 239 240
1409.20607 -5328.72791 -458.07596 -15341.30441 -13344.60880 -5776.56115
241 242 243 244 245 246
1616.39692 -19477.50168 -11226.61404 -913.45781 -4460.22574 3826.27158
247 248 249 250 251 252
876.82175 -1906.19055 10030.87165 4944.71388 -3867.05394 855.61960
253 254 255 256 257 258
-11000.78507 -2508.42139 -12329.65378 -9902.10033 741.08901 444.08379
259 260 261 262 263 264
-11610.78019 -10662.07525 -4904.09918 -14215.71971 -8515.60811 -12137.11705
265 266 267 268 269 270
9937.11971 -10038.87896 -13158.75994 -48.09113 -18391.23347 9131.32358
271 272 273 274 275 276
-9185.22901 572.25528 -6541.97025 -7987.78570 -6091.36506 4692.03015
277 278 279 280 281 282
20654.46157 -10730.90448 -17017.15667 -9141.40030 -1046.97755 59659.49390
283 284 285 286 287 288
8389.73878 -7304.97409 -3632.81826 18560.88732 -10657.56696 -13972.36400
289
-10747.84240
> postscript(file="/var/wessaorg/rcomp/tmp/62oma1353339095.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 4625.14380 NA
1 5635.92984 4625.14380
2 13767.13349 5635.92984
3 -604.92763 13767.13349
4 -16104.21531 -604.92763
5 28679.82767 -16104.21531
6 -15872.21420 28679.82767
7 9155.01849 -15872.21420
8 -11668.23176 9155.01849
9 30260.05356 -11668.23176
10 -10325.56574 30260.05356
11 -33941.02040 -10325.56574
12 2253.14136 -33941.02040
13 8400.76515 2253.14136
14 -1218.31091 8400.76515
15 -10895.17995 -1218.31091
16 -21101.24635 -10895.17995
17 24949.56662 -21101.24635
18 15992.01269 24949.56662
19 13289.07841 15992.01269
20 32173.04387 13289.07841
21 -11020.95753 32173.04387
22 -20448.75420 -11020.95753
23 -4735.81973 -20448.75420
24 27934.33201 -4735.81973
25 -20351.69750 27934.33201
26 -5085.56771 -20351.69750
27 -7501.34764 -5085.56771
28 -25173.79613 -7501.34764
29 -13291.73719 -25173.79613
30 -23281.01854 -13291.73719
31 -2847.83656 -23281.01854
32 7806.77379 -2847.83656
33 -24512.65172 7806.77379
34 36022.63346 -24512.65172
35 -29250.76938 36022.63346
36 -17626.82621 -29250.76938
37 3345.44351 -17626.82621
38 -27608.31829 3345.44351
39 -9228.66658 -27608.31829
40 -14920.50815 -9228.66658
41 -18326.77924 -14920.50815
42 -15635.19511 -18326.77924
43 -27849.62532 -15635.19511
44 -21461.81694 -27849.62532
45 -1234.38629 -21461.81694
46 -16562.82744 -1234.38629
47 -2647.25621 -16562.82744
48 3129.30786 -2647.25621
49 20333.50610 3129.30786
50 38063.37819 20333.50610
51 -8316.23709 38063.37819
52 1441.37335 -8316.23709
53 30241.37970 1441.37335
54 -6913.17620 30241.37970
55 -20767.09543 -6913.17620
56 -12786.36544 -20767.09543
57 5154.21917 -12786.36544
58 8173.70479 5154.21917
59 -6002.69443 8173.70479
60 -8016.94705 -6002.69443
61 25944.13209 -8016.94705
62 10767.48650 25944.13209
63 -13188.23825 10767.48650
64 -21729.71406 -13188.23825
65 -30440.62937 -21729.71406
66 12465.47829 -30440.62937
67 10175.16542 12465.47829
68 -3888.18463 10175.16542
69 -10796.78614 -3888.18463
70 4561.30496 -10796.78614
71 -9783.18207 4561.30496
72 -29626.67899 -9783.18207
73 34540.75917 -29626.67899
74 -18149.52387 34540.75917
75 13400.58101 -18149.52387
76 13038.79612 13400.58101
77 1996.37577 13038.79612
78 -1976.76158 1996.37577
79 804.38879 -1976.76158
80 5194.72829 804.38879
81 -2760.55805 5194.72829
82 7451.29328 -2760.55805
83 -5101.64026 7451.29328
84 149045.16614 -5101.64026
85 -1253.30656 149045.16614
86 511.14181 -1253.30656
87 -8245.97212 511.14181
88 6301.33852 -8245.97212
89 -11727.04098 6301.33852
90 -6196.60669 -11727.04098
91 11519.23145 -6196.60669
92 -11700.68538 11519.23145
93 10192.66568 -11700.68538
94 -17302.27954 10192.66568
95 10234.21204 -17302.27954
96 -26997.91242 10234.21204
97 29635.21227 -26997.91242
98 -20616.09690 29635.21227
99 11814.91185 -20616.09690
100 299.68127 11814.91185
101 -10401.40368 299.68127
102 5927.07008 -10401.40368
103 35855.82849 5927.07008
104 26518.66477 35855.82849
105 22269.04063 26518.66477
106 1402.40997 22269.04063
107 -9412.12861 1402.40997
108 4474.69022 -9412.12861
109 13187.62407 4474.69022
110 7903.79268 13187.62407
111 -16346.07118 7903.79268
112 12514.49371 -16346.07118
113 -23761.52858 12514.49371
114 -33518.63587 -23761.52858
115 -13392.76826 -33518.63587
116 -5931.96874 -13392.76826
117 -17831.67157 -5931.96874
118 97178.85390 -17831.67157
119 1481.88590 97178.85390
120 -20075.36766 1481.88590
121 -16713.86090 -20075.36766
122 -15235.84584 -16713.86090
123 90751.02751 -15235.84584
124 23604.75963 90751.02751
125 -9229.74759 23604.75963
126 9146.92882 -9229.74759
127 1047.68360 9146.92882
128 13459.53913 1047.68360
129 -11335.07485 13459.53913
130 2226.10522 -11335.07485
131 -14591.60152 2226.10522
132 -752.50006 -14591.60152
133 -10899.33625 -752.50006
134 -8011.70449 -10899.33625
135 -41647.54451 -8011.70449
136 -11398.44447 -41647.54451
137 24366.58557 -11398.44447
138 23547.13143 24366.58557
139 -9804.63758 23547.13143
140 -12008.52007 -9804.63758
141 3168.35231 -12008.52007
142 -172.38010 3168.35231
143 -21367.30493 -172.38010
144 42108.92355 -21367.30493
145 5959.68447 42108.92355
146 65234.31622 5959.68447
147 39649.08768 65234.31622
148 -10424.41851 39649.08768
149 27492.01397 -10424.41851
150 -13125.27327 27492.01397
151 -4223.08502 -13125.27327
152 -2679.04094 -4223.08502
153 -20584.33065 -2679.04094
154 -5332.70462 -20584.33065
155 50818.03114 -5332.70462
156 -4165.61552 50818.03114
157 -8745.28764 -4165.61552
158 14592.13791 -8745.28764
159 10223.61335 14592.13791
160 -6236.41950 10223.61335
161 97353.46398 -6236.41950
162 -3694.33435 97353.46398
163 22586.38705 -3694.33435
164 -2024.13781 22586.38705
165 -26319.63044 -2024.13781
166 -19941.89321 -26319.63044
167 29262.08157 -19941.89321
168 -24724.90277 29262.08157
169 21842.23541 -24724.90277
170 21924.05988 21842.23541
171 -3879.80124 21924.05988
172 17270.21023 -3879.80124
173 -1777.91575 17270.21023
174 62145.40836 -1777.91575
175 -19885.86603 62145.40836
176 32254.79094 -19885.86603
177 -5976.04070 32254.79094
178 15505.82840 -5976.04070
179 -15193.86536 15505.82840
180 5062.21370 -15193.86536
181 24457.25208 5062.21370
182 -8596.48978 24457.25208
183 -7439.66290 -8596.48978
184 -47243.71062 -7439.66290
185 -9826.83191 -47243.71062
186 -14357.99297 -9826.83191
187 -10723.33760 -14357.99297
188 -19291.12376 -10723.33760
189 15117.83161 -19291.12376
190 18011.27731 15117.83161
191 -3686.74712 18011.27731
192 -14160.79509 -3686.74712
193 10071.90947 -14160.79509
194 -8134.34000 10071.90947
195 -15613.29101 -8134.34000
196 -7753.67920 -15613.29101
197 31848.10582 -7753.67920
198 16864.82388 31848.10582
199 6751.33027 16864.82388
200 -14136.69830 6751.33027
201 -10210.89591 -14136.69830
202 -4623.02995 -10210.89591
203 -2591.51207 -4623.02995
204 -8181.70460 -2591.51207
205 889.77237 -8181.70460
206 -11522.86686 889.77237
207 -10943.07281 -11522.86686
208 -13081.44780 -10943.07281
209 -18322.93053 -13081.44780
210 -18478.74542 -18322.93053
211 -3166.62835 -18478.74542
212 -3332.45294 -3166.62835
213 -13802.96372 -3332.45294
214 25538.03500 -13802.96372
215 -13116.44412 25538.03500
216 -7746.89756 -13116.44412
217 12475.58772 -7746.89756
218 -13098.86912 12475.58772
219 -7294.53272 -13098.86912
220 -6970.68896 -7294.53272
221 -2667.76182 -6970.68896
222 -19078.40790 -2667.76182
223 -5504.59552 -19078.40790
224 -8538.23321 -5504.59552
225 -21681.40358 -8538.23321
226 6223.12133 -21681.40358
227 355.21231 6223.12133
228 -12599.41796 355.21231
229 34871.30668 -12599.41796
230 -8783.76660 34871.30668
231 -12104.99425 -8783.76660
232 5067.06339 -12104.99425
233 -3537.82604 5067.06339
234 1409.20607 -3537.82604
235 -5328.72791 1409.20607
236 -458.07596 -5328.72791
237 -15341.30441 -458.07596
238 -13344.60880 -15341.30441
239 -5776.56115 -13344.60880
240 1616.39692 -5776.56115
241 -19477.50168 1616.39692
242 -11226.61404 -19477.50168
243 -913.45781 -11226.61404
244 -4460.22574 -913.45781
245 3826.27158 -4460.22574
246 876.82175 3826.27158
247 -1906.19055 876.82175
248 10030.87165 -1906.19055
249 4944.71388 10030.87165
250 -3867.05394 4944.71388
251 855.61960 -3867.05394
252 -11000.78507 855.61960
253 -2508.42139 -11000.78507
254 -12329.65378 -2508.42139
255 -9902.10033 -12329.65378
256 741.08901 -9902.10033
257 444.08379 741.08901
258 -11610.78019 444.08379
259 -10662.07525 -11610.78019
260 -4904.09918 -10662.07525
261 -14215.71971 -4904.09918
262 -8515.60811 -14215.71971
263 -12137.11705 -8515.60811
264 9937.11971 -12137.11705
265 -10038.87896 9937.11971
266 -13158.75994 -10038.87896
267 -48.09113 -13158.75994
268 -18391.23347 -48.09113
269 9131.32358 -18391.23347
270 -9185.22901 9131.32358
271 572.25528 -9185.22901
272 -6541.97025 572.25528
273 -7987.78570 -6541.97025
274 -6091.36506 -7987.78570
275 4692.03015 -6091.36506
276 20654.46157 4692.03015
277 -10730.90448 20654.46157
278 -17017.15667 -10730.90448
279 -9141.40030 -17017.15667
280 -1046.97755 -9141.40030
281 59659.49390 -1046.97755
282 8389.73878 59659.49390
283 -7304.97409 8389.73878
284 -3632.81826 -7304.97409
285 18560.88732 -3632.81826
286 -10657.56696 18560.88732
287 -13972.36400 -10657.56696
288 -10747.84240 -13972.36400
289 NA -10747.84240
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] 5635.92984 4625.14380
[2,] 13767.13349 5635.92984
[3,] -604.92763 13767.13349
[4,] -16104.21531 -604.92763
[5,] 28679.82767 -16104.21531
[6,] -15872.21420 28679.82767
[7,] 9155.01849 -15872.21420
[8,] -11668.23176 9155.01849
[9,] 30260.05356 -11668.23176
[10,] -10325.56574 30260.05356
[11,] -33941.02040 -10325.56574
[12,] 2253.14136 -33941.02040
[13,] 8400.76515 2253.14136
[14,] -1218.31091 8400.76515
[15,] -10895.17995 -1218.31091
[16,] -21101.24635 -10895.17995
[17,] 24949.56662 -21101.24635
[18,] 15992.01269 24949.56662
[19,] 13289.07841 15992.01269
[20,] 32173.04387 13289.07841
[21,] -11020.95753 32173.04387
[22,] -20448.75420 -11020.95753
[23,] -4735.81973 -20448.75420
[24,] 27934.33201 -4735.81973
[25,] -20351.69750 27934.33201
[26,] -5085.56771 -20351.69750
[27,] -7501.34764 -5085.56771
[28,] -25173.79613 -7501.34764
[29,] -13291.73719 -25173.79613
[30,] -23281.01854 -13291.73719
[31,] -2847.83656 -23281.01854
[32,] 7806.77379 -2847.83656
[33,] -24512.65172 7806.77379
[34,] 36022.63346 -24512.65172
[35,] -29250.76938 36022.63346
[36,] -17626.82621 -29250.76938
[37,] 3345.44351 -17626.82621
[38,] -27608.31829 3345.44351
[39,] -9228.66658 -27608.31829
[40,] -14920.50815 -9228.66658
[41,] -18326.77924 -14920.50815
[42,] -15635.19511 -18326.77924
[43,] -27849.62532 -15635.19511
[44,] -21461.81694 -27849.62532
[45,] -1234.38629 -21461.81694
[46,] -16562.82744 -1234.38629
[47,] -2647.25621 -16562.82744
[48,] 3129.30786 -2647.25621
[49,] 20333.50610 3129.30786
[50,] 38063.37819 20333.50610
[51,] -8316.23709 38063.37819
[52,] 1441.37335 -8316.23709
[53,] 30241.37970 1441.37335
[54,] -6913.17620 30241.37970
[55,] -20767.09543 -6913.17620
[56,] -12786.36544 -20767.09543
[57,] 5154.21917 -12786.36544
[58,] 8173.70479 5154.21917
[59,] -6002.69443 8173.70479
[60,] -8016.94705 -6002.69443
[61,] 25944.13209 -8016.94705
[62,] 10767.48650 25944.13209
[63,] -13188.23825 10767.48650
[64,] -21729.71406 -13188.23825
[65,] -30440.62937 -21729.71406
[66,] 12465.47829 -30440.62937
[67,] 10175.16542 12465.47829
[68,] -3888.18463 10175.16542
[69,] -10796.78614 -3888.18463
[70,] 4561.30496 -10796.78614
[71,] -9783.18207 4561.30496
[72,] -29626.67899 -9783.18207
[73,] 34540.75917 -29626.67899
[74,] -18149.52387 34540.75917
[75,] 13400.58101 -18149.52387
[76,] 13038.79612 13400.58101
[77,] 1996.37577 13038.79612
[78,] -1976.76158 1996.37577
[79,] 804.38879 -1976.76158
[80,] 5194.72829 804.38879
[81,] -2760.55805 5194.72829
[82,] 7451.29328 -2760.55805
[83,] -5101.64026 7451.29328
[84,] 149045.16614 -5101.64026
[85,] -1253.30656 149045.16614
[86,] 511.14181 -1253.30656
[87,] -8245.97212 511.14181
[88,] 6301.33852 -8245.97212
[89,] -11727.04098 6301.33852
[90,] -6196.60669 -11727.04098
[91,] 11519.23145 -6196.60669
[92,] -11700.68538 11519.23145
[93,] 10192.66568 -11700.68538
[94,] -17302.27954 10192.66568
[95,] 10234.21204 -17302.27954
[96,] -26997.91242 10234.21204
[97,] 29635.21227 -26997.91242
[98,] -20616.09690 29635.21227
[99,] 11814.91185 -20616.09690
[100,] 299.68127 11814.91185
[101,] -10401.40368 299.68127
[102,] 5927.07008 -10401.40368
[103,] 35855.82849 5927.07008
[104,] 26518.66477 35855.82849
[105,] 22269.04063 26518.66477
[106,] 1402.40997 22269.04063
[107,] -9412.12861 1402.40997
[108,] 4474.69022 -9412.12861
[109,] 13187.62407 4474.69022
[110,] 7903.79268 13187.62407
[111,] -16346.07118 7903.79268
[112,] 12514.49371 -16346.07118
[113,] -23761.52858 12514.49371
[114,] -33518.63587 -23761.52858
[115,] -13392.76826 -33518.63587
[116,] -5931.96874 -13392.76826
[117,] -17831.67157 -5931.96874
[118,] 97178.85390 -17831.67157
[119,] 1481.88590 97178.85390
[120,] -20075.36766 1481.88590
[121,] -16713.86090 -20075.36766
[122,] -15235.84584 -16713.86090
[123,] 90751.02751 -15235.84584
[124,] 23604.75963 90751.02751
[125,] -9229.74759 23604.75963
[126,] 9146.92882 -9229.74759
[127,] 1047.68360 9146.92882
[128,] 13459.53913 1047.68360
[129,] -11335.07485 13459.53913
[130,] 2226.10522 -11335.07485
[131,] -14591.60152 2226.10522
[132,] -752.50006 -14591.60152
[133,] -10899.33625 -752.50006
[134,] -8011.70449 -10899.33625
[135,] -41647.54451 -8011.70449
[136,] -11398.44447 -41647.54451
[137,] 24366.58557 -11398.44447
[138,] 23547.13143 24366.58557
[139,] -9804.63758 23547.13143
[140,] -12008.52007 -9804.63758
[141,] 3168.35231 -12008.52007
[142,] -172.38010 3168.35231
[143,] -21367.30493 -172.38010
[144,] 42108.92355 -21367.30493
[145,] 5959.68447 42108.92355
[146,] 65234.31622 5959.68447
[147,] 39649.08768 65234.31622
[148,] -10424.41851 39649.08768
[149,] 27492.01397 -10424.41851
[150,] -13125.27327 27492.01397
[151,] -4223.08502 -13125.27327
[152,] -2679.04094 -4223.08502
[153,] -20584.33065 -2679.04094
[154,] -5332.70462 -20584.33065
[155,] 50818.03114 -5332.70462
[156,] -4165.61552 50818.03114
[157,] -8745.28764 -4165.61552
[158,] 14592.13791 -8745.28764
[159,] 10223.61335 14592.13791
[160,] -6236.41950 10223.61335
[161,] 97353.46398 -6236.41950
[162,] -3694.33435 97353.46398
[163,] 22586.38705 -3694.33435
[164,] -2024.13781 22586.38705
[165,] -26319.63044 -2024.13781
[166,] -19941.89321 -26319.63044
[167,] 29262.08157 -19941.89321
[168,] -24724.90277 29262.08157
[169,] 21842.23541 -24724.90277
[170,] 21924.05988 21842.23541
[171,] -3879.80124 21924.05988
[172,] 17270.21023 -3879.80124
[173,] -1777.91575 17270.21023
[174,] 62145.40836 -1777.91575
[175,] -19885.86603 62145.40836
[176,] 32254.79094 -19885.86603
[177,] -5976.04070 32254.79094
[178,] 15505.82840 -5976.04070
[179,] -15193.86536 15505.82840
[180,] 5062.21370 -15193.86536
[181,] 24457.25208 5062.21370
[182,] -8596.48978 24457.25208
[183,] -7439.66290 -8596.48978
[184,] -47243.71062 -7439.66290
[185,] -9826.83191 -47243.71062
[186,] -14357.99297 -9826.83191
[187,] -10723.33760 -14357.99297
[188,] -19291.12376 -10723.33760
[189,] 15117.83161 -19291.12376
[190,] 18011.27731 15117.83161
[191,] -3686.74712 18011.27731
[192,] -14160.79509 -3686.74712
[193,] 10071.90947 -14160.79509
[194,] -8134.34000 10071.90947
[195,] -15613.29101 -8134.34000
[196,] -7753.67920 -15613.29101
[197,] 31848.10582 -7753.67920
[198,] 16864.82388 31848.10582
[199,] 6751.33027 16864.82388
[200,] -14136.69830 6751.33027
[201,] -10210.89591 -14136.69830
[202,] -4623.02995 -10210.89591
[203,] -2591.51207 -4623.02995
[204,] -8181.70460 -2591.51207
[205,] 889.77237 -8181.70460
[206,] -11522.86686 889.77237
[207,] -10943.07281 -11522.86686
[208,] -13081.44780 -10943.07281
[209,] -18322.93053 -13081.44780
[210,] -18478.74542 -18322.93053
[211,] -3166.62835 -18478.74542
[212,] -3332.45294 -3166.62835
[213,] -13802.96372 -3332.45294
[214,] 25538.03500 -13802.96372
[215,] -13116.44412 25538.03500
[216,] -7746.89756 -13116.44412
[217,] 12475.58772 -7746.89756
[218,] -13098.86912 12475.58772
[219,] -7294.53272 -13098.86912
[220,] -6970.68896 -7294.53272
[221,] -2667.76182 -6970.68896
[222,] -19078.40790 -2667.76182
[223,] -5504.59552 -19078.40790
[224,] -8538.23321 -5504.59552
[225,] -21681.40358 -8538.23321
[226,] 6223.12133 -21681.40358
[227,] 355.21231 6223.12133
[228,] -12599.41796 355.21231
[229,] 34871.30668 -12599.41796
[230,] -8783.76660 34871.30668
[231,] -12104.99425 -8783.76660
[232,] 5067.06339 -12104.99425
[233,] -3537.82604 5067.06339
[234,] 1409.20607 -3537.82604
[235,] -5328.72791 1409.20607
[236,] -458.07596 -5328.72791
[237,] -15341.30441 -458.07596
[238,] -13344.60880 -15341.30441
[239,] -5776.56115 -13344.60880
[240,] 1616.39692 -5776.56115
[241,] -19477.50168 1616.39692
[242,] -11226.61404 -19477.50168
[243,] -913.45781 -11226.61404
[244,] -4460.22574 -913.45781
[245,] 3826.27158 -4460.22574
[246,] 876.82175 3826.27158
[247,] -1906.19055 876.82175
[248,] 10030.87165 -1906.19055
[249,] 4944.71388 10030.87165
[250,] -3867.05394 4944.71388
[251,] 855.61960 -3867.05394
[252,] -11000.78507 855.61960
[253,] -2508.42139 -11000.78507
[254,] -12329.65378 -2508.42139
[255,] -9902.10033 -12329.65378
[256,] 741.08901 -9902.10033
[257,] 444.08379 741.08901
[258,] -11610.78019 444.08379
[259,] -10662.07525 -11610.78019
[260,] -4904.09918 -10662.07525
[261,] -14215.71971 -4904.09918
[262,] -8515.60811 -14215.71971
[263,] -12137.11705 -8515.60811
[264,] 9937.11971 -12137.11705
[265,] -10038.87896 9937.11971
[266,] -13158.75994 -10038.87896
[267,] -48.09113 -13158.75994
[268,] -18391.23347 -48.09113
[269,] 9131.32358 -18391.23347
[270,] -9185.22901 9131.32358
[271,] 572.25528 -9185.22901
[272,] -6541.97025 572.25528
[273,] -7987.78570 -6541.97025
[274,] -6091.36506 -7987.78570
[275,] 4692.03015 -6091.36506
[276,] 20654.46157 4692.03015
[277,] -10730.90448 20654.46157
[278,] -17017.15667 -10730.90448
[279,] -9141.40030 -17017.15667
[280,] -1046.97755 -9141.40030
[281,] 59659.49390 -1046.97755
[282,] 8389.73878 59659.49390
[283,] -7304.97409 8389.73878
[284,] -3632.81826 -7304.97409
[285,] 18560.88732 -3632.81826
[286,] -10657.56696 18560.88732
[287,] -13972.36400 -10657.56696
[288,] -10747.84240 -13972.36400
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 5635.92984 4625.14380
2 13767.13349 5635.92984
3 -604.92763 13767.13349
4 -16104.21531 -604.92763
5 28679.82767 -16104.21531
6 -15872.21420 28679.82767
7 9155.01849 -15872.21420
8 -11668.23176 9155.01849
9 30260.05356 -11668.23176
10 -10325.56574 30260.05356
11 -33941.02040 -10325.56574
12 2253.14136 -33941.02040
13 8400.76515 2253.14136
14 -1218.31091 8400.76515
15 -10895.17995 -1218.31091
16 -21101.24635 -10895.17995
17 24949.56662 -21101.24635
18 15992.01269 24949.56662
19 13289.07841 15992.01269
20 32173.04387 13289.07841
21 -11020.95753 32173.04387
22 -20448.75420 -11020.95753
23 -4735.81973 -20448.75420
24 27934.33201 -4735.81973
25 -20351.69750 27934.33201
26 -5085.56771 -20351.69750
27 -7501.34764 -5085.56771
28 -25173.79613 -7501.34764
29 -13291.73719 -25173.79613
30 -23281.01854 -13291.73719
31 -2847.83656 -23281.01854
32 7806.77379 -2847.83656
33 -24512.65172 7806.77379
34 36022.63346 -24512.65172
35 -29250.76938 36022.63346
36 -17626.82621 -29250.76938
37 3345.44351 -17626.82621
38 -27608.31829 3345.44351
39 -9228.66658 -27608.31829
40 -14920.50815 -9228.66658
41 -18326.77924 -14920.50815
42 -15635.19511 -18326.77924
43 -27849.62532 -15635.19511
44 -21461.81694 -27849.62532
45 -1234.38629 -21461.81694
46 -16562.82744 -1234.38629
47 -2647.25621 -16562.82744
48 3129.30786 -2647.25621
49 20333.50610 3129.30786
50 38063.37819 20333.50610
51 -8316.23709 38063.37819
52 1441.37335 -8316.23709
53 30241.37970 1441.37335
54 -6913.17620 30241.37970
55 -20767.09543 -6913.17620
56 -12786.36544 -20767.09543
57 5154.21917 -12786.36544
58 8173.70479 5154.21917
59 -6002.69443 8173.70479
60 -8016.94705 -6002.69443
61 25944.13209 -8016.94705
62 10767.48650 25944.13209
63 -13188.23825 10767.48650
64 -21729.71406 -13188.23825
65 -30440.62937 -21729.71406
66 12465.47829 -30440.62937
67 10175.16542 12465.47829
68 -3888.18463 10175.16542
69 -10796.78614 -3888.18463
70 4561.30496 -10796.78614
71 -9783.18207 4561.30496
72 -29626.67899 -9783.18207
73 34540.75917 -29626.67899
74 -18149.52387 34540.75917
75 13400.58101 -18149.52387
76 13038.79612 13400.58101
77 1996.37577 13038.79612
78 -1976.76158 1996.37577
79 804.38879 -1976.76158
80 5194.72829 804.38879
81 -2760.55805 5194.72829
82 7451.29328 -2760.55805
83 -5101.64026 7451.29328
84 149045.16614 -5101.64026
85 -1253.30656 149045.16614
86 511.14181 -1253.30656
87 -8245.97212 511.14181
88 6301.33852 -8245.97212
89 -11727.04098 6301.33852
90 -6196.60669 -11727.04098
91 11519.23145 -6196.60669
92 -11700.68538 11519.23145
93 10192.66568 -11700.68538
94 -17302.27954 10192.66568
95 10234.21204 -17302.27954
96 -26997.91242 10234.21204
97 29635.21227 -26997.91242
98 -20616.09690 29635.21227
99 11814.91185 -20616.09690
100 299.68127 11814.91185
101 -10401.40368 299.68127
102 5927.07008 -10401.40368
103 35855.82849 5927.07008
104 26518.66477 35855.82849
105 22269.04063 26518.66477
106 1402.40997 22269.04063
107 -9412.12861 1402.40997
108 4474.69022 -9412.12861
109 13187.62407 4474.69022
110 7903.79268 13187.62407
111 -16346.07118 7903.79268
112 12514.49371 -16346.07118
113 -23761.52858 12514.49371
114 -33518.63587 -23761.52858
115 -13392.76826 -33518.63587
116 -5931.96874 -13392.76826
117 -17831.67157 -5931.96874
118 97178.85390 -17831.67157
119 1481.88590 97178.85390
120 -20075.36766 1481.88590
121 -16713.86090 -20075.36766
122 -15235.84584 -16713.86090
123 90751.02751 -15235.84584
124 23604.75963 90751.02751
125 -9229.74759 23604.75963
126 9146.92882 -9229.74759
127 1047.68360 9146.92882
128 13459.53913 1047.68360
129 -11335.07485 13459.53913
130 2226.10522 -11335.07485
131 -14591.60152 2226.10522
132 -752.50006 -14591.60152
133 -10899.33625 -752.50006
134 -8011.70449 -10899.33625
135 -41647.54451 -8011.70449
136 -11398.44447 -41647.54451
137 24366.58557 -11398.44447
138 23547.13143 24366.58557
139 -9804.63758 23547.13143
140 -12008.52007 -9804.63758
141 3168.35231 -12008.52007
142 -172.38010 3168.35231
143 -21367.30493 -172.38010
144 42108.92355 -21367.30493
145 5959.68447 42108.92355
146 65234.31622 5959.68447
147 39649.08768 65234.31622
148 -10424.41851 39649.08768
149 27492.01397 -10424.41851
150 -13125.27327 27492.01397
151 -4223.08502 -13125.27327
152 -2679.04094 -4223.08502
153 -20584.33065 -2679.04094
154 -5332.70462 -20584.33065
155 50818.03114 -5332.70462
156 -4165.61552 50818.03114
157 -8745.28764 -4165.61552
158 14592.13791 -8745.28764
159 10223.61335 14592.13791
160 -6236.41950 10223.61335
161 97353.46398 -6236.41950
162 -3694.33435 97353.46398
163 22586.38705 -3694.33435
164 -2024.13781 22586.38705
165 -26319.63044 -2024.13781
166 -19941.89321 -26319.63044
167 29262.08157 -19941.89321
168 -24724.90277 29262.08157
169 21842.23541 -24724.90277
170 21924.05988 21842.23541
171 -3879.80124 21924.05988
172 17270.21023 -3879.80124
173 -1777.91575 17270.21023
174 62145.40836 -1777.91575
175 -19885.86603 62145.40836
176 32254.79094 -19885.86603
177 -5976.04070 32254.79094
178 15505.82840 -5976.04070
179 -15193.86536 15505.82840
180 5062.21370 -15193.86536
181 24457.25208 5062.21370
182 -8596.48978 24457.25208
183 -7439.66290 -8596.48978
184 -47243.71062 -7439.66290
185 -9826.83191 -47243.71062
186 -14357.99297 -9826.83191
187 -10723.33760 -14357.99297
188 -19291.12376 -10723.33760
189 15117.83161 -19291.12376
190 18011.27731 15117.83161
191 -3686.74712 18011.27731
192 -14160.79509 -3686.74712
193 10071.90947 -14160.79509
194 -8134.34000 10071.90947
195 -15613.29101 -8134.34000
196 -7753.67920 -15613.29101
197 31848.10582 -7753.67920
198 16864.82388 31848.10582
199 6751.33027 16864.82388
200 -14136.69830 6751.33027
201 -10210.89591 -14136.69830
202 -4623.02995 -10210.89591
203 -2591.51207 -4623.02995
204 -8181.70460 -2591.51207
205 889.77237 -8181.70460
206 -11522.86686 889.77237
207 -10943.07281 -11522.86686
208 -13081.44780 -10943.07281
209 -18322.93053 -13081.44780
210 -18478.74542 -18322.93053
211 -3166.62835 -18478.74542
212 -3332.45294 -3166.62835
213 -13802.96372 -3332.45294
214 25538.03500 -13802.96372
215 -13116.44412 25538.03500
216 -7746.89756 -13116.44412
217 12475.58772 -7746.89756
218 -13098.86912 12475.58772
219 -7294.53272 -13098.86912
220 -6970.68896 -7294.53272
221 -2667.76182 -6970.68896
222 -19078.40790 -2667.76182
223 -5504.59552 -19078.40790
224 -8538.23321 -5504.59552
225 -21681.40358 -8538.23321
226 6223.12133 -21681.40358
227 355.21231 6223.12133
228 -12599.41796 355.21231
229 34871.30668 -12599.41796
230 -8783.76660 34871.30668
231 -12104.99425 -8783.76660
232 5067.06339 -12104.99425
233 -3537.82604 5067.06339
234 1409.20607 -3537.82604
235 -5328.72791 1409.20607
236 -458.07596 -5328.72791
237 -15341.30441 -458.07596
238 -13344.60880 -15341.30441
239 -5776.56115 -13344.60880
240 1616.39692 -5776.56115
241 -19477.50168 1616.39692
242 -11226.61404 -19477.50168
243 -913.45781 -11226.61404
244 -4460.22574 -913.45781
245 3826.27158 -4460.22574
246 876.82175 3826.27158
247 -1906.19055 876.82175
248 10030.87165 -1906.19055
249 4944.71388 10030.87165
250 -3867.05394 4944.71388
251 855.61960 -3867.05394
252 -11000.78507 855.61960
253 -2508.42139 -11000.78507
254 -12329.65378 -2508.42139
255 -9902.10033 -12329.65378
256 741.08901 -9902.10033
257 444.08379 741.08901
258 -11610.78019 444.08379
259 -10662.07525 -11610.78019
260 -4904.09918 -10662.07525
261 -14215.71971 -4904.09918
262 -8515.60811 -14215.71971
263 -12137.11705 -8515.60811
264 9937.11971 -12137.11705
265 -10038.87896 9937.11971
266 -13158.75994 -10038.87896
267 -48.09113 -13158.75994
268 -18391.23347 -48.09113
269 9131.32358 -18391.23347
270 -9185.22901 9131.32358
271 572.25528 -9185.22901
272 -6541.97025 572.25528
273 -7987.78570 -6541.97025
274 -6091.36506 -7987.78570
275 4692.03015 -6091.36506
276 20654.46157 4692.03015
277 -10730.90448 20654.46157
278 -17017.15667 -10730.90448
279 -9141.40030 -17017.15667
280 -1046.97755 -9141.40030
281 59659.49390 -1046.97755
282 8389.73878 59659.49390
283 -7304.97409 8389.73878
284 -3632.81826 -7304.97409
285 18560.88732 -3632.81826
286 -10657.56696 18560.88732
287 -13972.36400 -10657.56696
288 -10747.84240 -13972.36400
> 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/7nk3c1353339095.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/8yev21353339095.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/9e8bs1353339095.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/10rnaa1353339095.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/112fwm1353339095.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/128zxz1353339095.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/13psz51353339096.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/14vahi1353339096.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/15p4eg1353339096.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/1624u51353339096.tab")
+ }
>
> try(system("convert tmp/1w2t91353339095.ps tmp/1w2t91353339095.png",intern=TRUE))
character(0)
> try(system("convert tmp/2ogyx1353339095.ps tmp/2ogyx1353339095.png",intern=TRUE))
character(0)
> try(system("convert tmp/3strj1353339095.ps tmp/3strj1353339095.png",intern=TRUE))
character(0)
> try(system("convert tmp/4yfjx1353339095.ps tmp/4yfjx1353339095.png",intern=TRUE))
character(0)
> try(system("convert tmp/5wtzv1353339095.ps tmp/5wtzv1353339095.png",intern=TRUE))
character(0)
> try(system("convert tmp/62oma1353339095.ps tmp/62oma1353339095.png",intern=TRUE))
character(0)
> try(system("convert tmp/7nk3c1353339095.ps tmp/7nk3c1353339095.png",intern=TRUE))
character(0)
> try(system("convert tmp/8yev21353339095.ps tmp/8yev21353339095.png",intern=TRUE))
character(0)
> try(system("convert tmp/9e8bs1353339095.ps tmp/9e8bs1353339095.png",intern=TRUE))
character(0)
> try(system("convert tmp/10rnaa1353339095.ps tmp/10rnaa1353339095.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
12.655 0.931 13.614