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(1418
+ ,210907
+ ,56
+ ,3
+ ,79
+ ,30
+ ,115
+ ,112285
+ ,24188
+ ,144
+ ,145
+ ,869
+ ,120982
+ ,56
+ ,4
+ ,58
+ ,28
+ ,109
+ ,84786
+ ,18273
+ ,103
+ ,101
+ ,1530
+ ,176508
+ ,54
+ ,12
+ ,60
+ ,38
+ ,146
+ ,83123
+ ,14130
+ ,98
+ ,98
+ ,2172
+ ,179321
+ ,89
+ ,2
+ ,108
+ ,30
+ ,116
+ ,101193
+ ,32287
+ ,135
+ ,132
+ ,901
+ ,123185
+ ,40
+ ,1
+ ,49
+ ,22
+ ,68
+ ,38361
+ ,8654
+ ,61
+ ,60
+ ,463
+ ,52746
+ ,25
+ ,3
+ ,0
+ ,26
+ ,101
+ ,68504
+ ,9245
+ ,39
+ ,38
+ ,3201
+ ,385534
+ ,92
+ ,0
+ ,121
+ ,25
+ ,96
+ ,119182
+ ,33251
+ ,150
+ ,144
+ ,371
+ ,33170
+ ,18
+ ,0
+ ,1
+ ,18
+ ,67
+ ,22807
+ ,1271
+ ,5
+ ,5
+ ,1192
+ ,101645
+ ,63
+ ,0
+ ,20
+ ,11
+ ,44
+ ,17140
+ ,5279
+ ,28
+ ,28
+ ,1583
+ ,149061
+ ,44
+ ,5
+ ,43
+ ,26
+ ,100
+ ,116174
+ ,27101
+ ,84
+ ,84
+ ,1439
+ ,165446
+ ,33
+ ,0
+ ,69
+ ,25
+ ,93
+ ,57635
+ ,16373
+ ,80
+ ,79
+ ,1764
+ ,237213
+ ,84
+ ,0
+ ,78
+ ,38
+ ,140
+ ,66198
+ ,19716
+ ,130
+ ,127
+ ,1495
+ ,173326
+ ,88
+ ,7
+ ,86
+ ,44
+ ,166
+ ,71701
+ ,17753
+ ,82
+ ,78
+ ,1373
+ ,133131
+ ,55
+ ,7
+ ,44
+ ,30
+ ,99
+ ,57793
+ ,9028
+ ,60
+ ,60
+ ,2187
+ ,258873
+ ,60
+ ,3
+ ,104
+ ,40
+ ,139
+ ,80444
+ ,18653
+ ,131
+ ,131
+ ,1491
+ ,180083
+ ,66
+ ,9
+ ,63
+ ,34
+ ,130
+ ,53855
+ ,8828
+ ,84
+ ,84
+ ,4041
+ ,324799
+ ,154
+ ,0
+ ,158
+ ,47
+ ,181
+ ,97668
+ ,29498
+ ,140
+ ,133
+ ,1706
+ ,230964
+ ,53
+ ,4
+ ,102
+ ,30
+ ,116
+ ,133824
+ ,27563
+ ,151
+ ,150
+ ,2152
+ ,236785
+ ,119
+ ,3
+ ,77
+ ,31
+ ,116
+ ,101481
+ ,18293
+ ,91
+ ,91
+ ,1036
+ ,135473
+ ,41
+ ,0
+ ,82
+ ,23
+ ,88
+ ,99645
+ ,22530
+ ,138
+ ,132
+ ,1882
+ ,202925
+ ,61
+ ,7
+ ,115
+ ,36
+ ,139
+ ,114789
+ ,15977
+ ,150
+ ,136
+ ,1929
+ ,215147
+ ,58
+ ,0
+ ,101
+ ,36
+ ,135
+ ,99052
+ ,35082
+ ,124
+ ,124
+ ,2242
+ ,344297
+ ,75
+ ,1
+ ,80
+ ,30
+ ,108
+ ,67654
+ ,16116
+ ,119
+ ,118
+ ,1220
+ ,153935
+ ,33
+ ,5
+ ,50
+ ,25
+ ,89
+ ,65553
+ ,15849
+ ,73
+ ,70
+ ,1289
+ ,132943
+ ,40
+ ,7
+ ,83
+ ,39
+ ,156
+ ,97500
+ ,16026
+ ,110
+ ,107
+ ,2515
+ ,174724
+ ,92
+ ,0
+ ,123
+ ,34
+ ,129
+ ,69112
+ ,26569
+ ,123
+ ,119
+ ,2147
+ ,174415
+ ,100
+ ,0
+ ,73
+ ,31
+ ,118
+ ,82753
+ ,24785
+ ,90
+ ,89
+ ,2352
+ ,225548
+ ,112
+ ,5
+ ,81
+ ,31
+ ,118
+ ,85323
+ ,17569
+ ,116
+ ,112
+ ,1638
+ ,223632
+ ,73
+ ,0
+ ,105
+ ,33
+ ,125
+ ,72654
+ ,23825
+ ,113
+ ,108
+ ,1222
+ ,124817
+ ,40
+ ,0
+ ,47
+ ,25
+ ,95
+ ,30727
+ ,7869
+ ,56
+ ,52
+ ,1812
+ ,221698
+ ,45
+ ,0
+ ,105
+ ,33
+ ,126
+ ,77873
+ ,14975
+ ,115
+ ,112
+ ,1677
+ ,210767
+ ,60
+ ,3
+ ,94
+ ,35
+ ,135
+ ,117478
+ ,37791
+ ,119
+ ,116
+ ,1579
+ ,170266
+ ,62
+ ,4
+ ,44
+ ,42
+ ,154
+ ,74007
+ ,9605
+ ,129
+ ,123
+ ,1731
+ ,260561
+ ,75
+ ,1
+ ,114
+ ,43
+ ,165
+ ,90183
+ ,27295
+ ,127
+ ,125
+ ,807
+ ,84853
+ ,31
+ ,4
+ ,38
+ ,30
+ ,113
+ ,61542
+ ,2746
+ ,27
+ ,27
+ ,2452
+ ,294424
+ ,77
+ ,2
+ ,107
+ ,33
+ ,127
+ ,101494
+ ,34461
+ ,175
+ ,162
+ ,829
+ ,101011
+ ,34
+ ,0
+ ,30
+ ,13
+ ,52
+ ,27570
+ ,8098
+ ,35
+ ,32
+ ,1940
+ ,215641
+ ,46
+ ,0
+ ,71
+ ,32
+ ,121
+ ,55813
+ ,4787
+ ,64
+ ,64
+ ,2662
+ ,325107
+ ,99
+ ,0
+ ,84
+ ,36
+ ,136
+ ,79215
+ ,24919
+ ,96
+ ,92
+ ,186
+ ,7176
+ ,17
+ ,0
+ ,0
+ ,0
+ ,0
+ ,1423
+ ,603
+ ,0
+ ,0
+ ,1499
+ ,167542
+ ,66
+ ,2
+ ,59
+ ,28
+ ,108
+ ,55461
+ ,16329
+ ,84
+ ,83
+ ,865
+ ,106408
+ ,30
+ ,1
+ ,33
+ ,14
+ ,46
+ ,31081
+ ,12558
+ ,41
+ ,41
+ ,1793
+ ,96560
+ ,76
+ ,0
+ ,42
+ ,17
+ ,54
+ ,22996
+ ,7784
+ ,47
+ ,47
+ ,2527
+ ,265769
+ ,146
+ ,2
+ ,96
+ ,32
+ ,124
+ ,83122
+ ,28522
+ ,126
+ ,120
+ ,2747
+ ,269651
+ ,67
+ ,10
+ ,106
+ ,30
+ ,115
+ ,70106
+ ,22265
+ ,105
+ ,105
+ ,1324
+ ,149112
+ ,56
+ ,6
+ ,56
+ ,35
+ ,128
+ ,60578
+ ,14459
+ ,80
+ ,79
+ ,2702
+ ,175824
+ ,107
+ ,0
+ ,57
+ ,20
+ ,80
+ ,39992
+ ,14526
+ ,70
+ ,65
+ ,1383
+ ,152871
+ ,58
+ ,5
+ ,59
+ ,28
+ ,97
+ ,79892
+ ,22240
+ ,73
+ ,70
+ ,1179
+ ,111665
+ ,34
+ ,4
+ ,39
+ ,28
+ ,104
+ ,49810
+ ,11802
+ ,57
+ ,55
+ ,2099
+ ,116408
+ ,61
+ ,1
+ ,34
+ ,39
+ ,59
+ ,71570
+ ,7623
+ ,40
+ ,39
+ ,4308
+ ,362301
+ ,119
+ ,2
+ ,76
+ ,34
+ ,125
+ ,100708
+ ,11912
+ ,68
+ ,67
+ ,918
+ ,78800
+ ,42
+ ,2
+ ,20
+ ,26
+ ,82
+ ,33032
+ ,7935
+ ,21
+ ,21
+ ,1831
+ ,183167
+ ,66
+ ,0
+ ,91
+ ,39
+ ,149
+ ,82875
+ ,18220
+ ,127
+ ,127
+ ,3373
+ ,277965
+ ,89
+ ,8
+ ,115
+ ,39
+ ,149
+ ,139077
+ ,19199
+ ,154
+ ,152
+ ,1713
+ ,150629
+ ,44
+ ,3
+ ,85
+ ,33
+ ,122
+ ,71595
+ ,19918
+ ,116
+ ,113
+ ,1438
+ ,168809
+ ,66
+ ,0
+ ,76
+ ,28
+ ,118
+ ,72260
+ ,21884
+ ,102
+ ,99
+ ,496
+ ,24188
+ ,24
+ ,0
+ ,8
+ ,4
+ ,12
+ ,5950
+ ,2694
+ ,7
+ ,7
+ ,2253
+ ,329267
+ ,259
+ ,8
+ ,79
+ ,39
+ ,144
+ ,115762
+ ,15808
+ ,148
+ ,141
+ ,744
+ ,65029
+ ,17
+ ,5
+ ,21
+ ,18
+ ,67
+ ,32551
+ ,3597
+ ,21
+ ,21
+ ,1161
+ ,101097
+ ,64
+ ,3
+ ,30
+ ,14
+ ,52
+ ,31701
+ ,5296
+ ,35
+ ,35
+ ,2352
+ ,218946
+ ,41
+ ,1
+ ,76
+ ,29
+ ,108
+ ,80670
+ ,25239
+ ,112
+ ,109
+ ,2144
+ ,244052
+ ,68
+ ,5
+ ,101
+ ,44
+ ,166
+ ,143558
+ ,29801
+ ,137
+ ,133
+ ,4691
+ ,341570
+ ,168
+ ,1
+ ,94
+ ,21
+ ,80
+ ,117105
+ ,18450
+ ,135
+ ,123
+ ,1112
+ ,103597
+ ,43
+ ,1
+ ,27
+ ,16
+ ,60
+ ,23789
+ ,7132
+ ,26
+ ,26
+ ,2694
+ ,233328
+ ,132
+ ,5
+ ,92
+ ,28
+ ,107
+ ,120733
+ ,34861
+ ,230
+ ,230
+ ,1973
+ ,256462
+ ,105
+ ,0
+ ,123
+ ,35
+ ,127
+ ,105195
+ ,35940
+ ,181
+ ,166
+ ,1769
+ ,206161
+ ,71
+ ,12
+ ,75
+ ,28
+ ,107
+ ,73107
+ ,16688
+ ,71
+ ,68
+ ,3148
+ ,311473
+ ,112
+ ,8
+ ,128
+ ,38
+ ,146
+ ,132068
+ ,24683
+ ,147
+ ,147
+ ,2474
+ ,235800
+ ,94
+ ,8
+ ,105
+ ,23
+ ,84
+ ,149193
+ ,46230
+ ,190
+ ,179
+ ,2084
+ ,177939
+ ,82
+ ,8
+ ,55
+ ,36
+ ,141
+ ,46821
+ ,10387
+ ,64
+ ,61
+ ,1954
+ ,207176
+ ,70
+ ,8
+ ,56
+ ,32
+ ,123
+ ,87011
+ ,21436
+ ,105
+ ,101
+ ,1226
+ ,196553
+ ,57
+ ,2
+ ,41
+ ,29
+ ,111
+ ,95260
+ ,30546
+ ,107
+ ,108
+ ,1389
+ ,174184
+ ,53
+ ,0
+ ,72
+ ,25
+ ,98
+ ,55183
+ ,19746
+ ,94
+ ,90
+ ,1496
+ ,143246
+ ,103
+ ,5
+ ,67
+ ,27
+ ,105
+ ,106671
+ ,15977
+ ,116
+ ,114
+ ,2269
+ ,187559
+ ,121
+ ,8
+ ,75
+ ,36
+ ,135
+ ,73511
+ ,22583
+ ,106
+ ,103
+ ,1833
+ ,187681
+ ,62
+ ,2
+ ,114
+ ,28
+ ,107
+ ,92945
+ ,17274
+ ,143
+ ,142
+ ,1268
+ ,119016
+ ,52
+ ,5
+ ,118
+ ,23
+ ,85
+ ,78664
+ ,16469
+ ,81
+ ,79
+ ,1943
+ ,182192
+ ,52
+ ,12
+ ,77
+ ,40
+ ,155
+ ,70054
+ ,14251
+ ,89
+ ,88
+ ,893
+ ,73566
+ ,32
+ ,6
+ ,22
+ ,23
+ ,88
+ ,22618
+ ,3007
+ ,26
+ ,25
+ ,1762
+ ,194979
+ ,62
+ ,7
+ ,66
+ ,40
+ ,155
+ ,74011
+ ,16851
+ ,84
+ ,83
+ ,1403
+ ,167488
+ ,45
+ ,2
+ ,69
+ ,28
+ ,104
+ ,83737
+ ,21113
+ ,113
+ ,113
+ ,1425
+ ,143756
+ ,46
+ ,0
+ ,105
+ ,34
+ ,132
+ ,69094
+ ,17401
+ ,120
+ ,118
+ ,1857
+ ,275541
+ ,63
+ ,4
+ ,116
+ ,33
+ ,127
+ ,93133
+ ,23958
+ ,110
+ ,110
+ ,1840
+ ,243199
+ ,75
+ ,3
+ ,88
+ ,28
+ ,108
+ ,95536
+ ,23567
+ ,134
+ ,129
+ ,1502
+ ,182999
+ ,88
+ ,6
+ ,73
+ ,34
+ ,129
+ ,225920
+ ,13065
+ ,54
+ ,51
+ ,1441
+ ,135649
+ ,46
+ ,2
+ ,99
+ ,30
+ ,116
+ ,62133
+ ,15358
+ ,96
+ ,93
+ ,1420
+ ,152299
+ ,53
+ ,0
+ ,62
+ ,33
+ ,122
+ ,61370
+ ,14587
+ ,78
+ ,76
+ ,1416
+ ,120221
+ ,37
+ ,1
+ ,53
+ ,22
+ ,85
+ ,43836
+ ,12770
+ ,51
+ ,49
+ ,2970
+ ,346485
+ ,90
+ ,0
+ ,118
+ ,38
+ ,147
+ ,106117
+ ,24021
+ ,121
+ ,118
+ ,1317
+ ,145790
+ ,63
+ ,5
+ ,30
+ ,26
+ ,99
+ ,38692
+ ,9648
+ ,38
+ ,38
+ ,1644
+ ,193339
+ ,78
+ ,2
+ ,100
+ ,35
+ ,87
+ ,84651
+ ,20537
+ ,145
+ ,141
+ ,870
+ ,80953
+ ,25
+ ,0
+ ,49
+ ,8
+ ,28
+ ,56622
+ ,7905
+ ,59
+ ,58
+ ,1654
+ ,122774
+ ,45
+ ,0
+ ,24
+ ,24
+ ,90
+ ,15986
+ ,4527
+ ,27
+ ,27
+ ,1054
+ ,130585
+ ,46
+ ,5
+ ,67
+ ,29
+ ,109
+ ,95364
+ ,30495
+ ,91
+ ,91
+ ,937
+ ,112611
+ ,41
+ ,0
+ ,46
+ ,20
+ ,78
+ ,26706
+ ,7117
+ ,48
+ ,48
+ ,3004
+ ,286468
+ ,144
+ ,1
+ ,57
+ ,29
+ ,111
+ ,89691
+ ,17719
+ ,68
+ ,63
+ ,2008
+ ,241066
+ ,82
+ ,0
+ ,75
+ ,45
+ ,158
+ ,67267
+ ,27056
+ ,58
+ ,56
+ ,2547
+ ,148446
+ ,91
+ ,1
+ ,135
+ ,37
+ ,141
+ ,126846
+ ,33473
+ ,150
+ ,144
+ ,1885
+ ,204713
+ ,71
+ ,1
+ ,68
+ ,33
+ ,122
+ ,41140
+ ,9758
+ ,74
+ ,73
+ ,1626
+ ,182079
+ ,63
+ ,2
+ ,124
+ ,33
+ ,124
+ ,102860
+ ,21115
+ ,181
+ ,168
+ ,1468
+ ,140344
+ ,53
+ ,6
+ ,33
+ ,25
+ ,93
+ ,51715
+ ,7236
+ ,65
+ ,64
+ ,2445
+ ,220516
+ ,62
+ ,1
+ ,98
+ ,32
+ ,124
+ ,55801
+ ,13790
+ ,97
+ ,97
+ ,1964
+ ,243060
+ ,63
+ ,4
+ ,58
+ ,29
+ ,112
+ ,111813
+ ,32902
+ ,121
+ ,117
+ ,1381
+ ,162765
+ ,32
+ ,2
+ ,68
+ ,28
+ ,108
+ ,120293
+ ,25131
+ ,99
+ ,100
+ ,1369
+ ,182613
+ ,39
+ ,3
+ ,81
+ ,28
+ ,99
+ ,138599
+ ,30910
+ ,152
+ ,149
+ ,1659
+ ,232138
+ ,62
+ ,0
+ ,131
+ ,31
+ ,117
+ ,161647
+ ,35947
+ ,188
+ ,187
+ ,2888
+ ,265318
+ ,117
+ ,10
+ ,110
+ ,52
+ ,199
+ ,115929
+ ,29848
+ ,138
+ ,127
+ ,1290
+ ,85574
+ ,34
+ ,0
+ ,37
+ ,21
+ ,78
+ ,24266
+ ,6943
+ ,40
+ ,37
+ ,2845
+ ,310839
+ ,92
+ ,9
+ ,130
+ ,24
+ ,91
+ ,162901
+ ,42705
+ ,254
+ ,245
+ ,1982
+ ,225060
+ ,93
+ ,7
+ ,93
+ ,41
+ ,158
+ ,109825
+ ,31808
+ ,87
+ ,87
+ ,1904
+ ,232317
+ ,54
+ ,0
+ ,118
+ ,33
+ ,126
+ ,129838
+ ,26675
+ ,178
+ ,177
+ ,1391
+ ,144966
+ ,144
+ ,0
+ ,39
+ ,32
+ ,122
+ ,37510
+ ,8435
+ ,51
+ ,49
+ ,602
+ ,43287
+ ,14
+ ,4
+ ,13
+ ,19
+ ,71
+ ,43750
+ ,7409
+ ,49
+ ,49
+ ,1743
+ ,155754
+ ,61
+ ,4
+ ,74
+ ,20
+ ,75
+ ,40652
+ ,14993
+ ,73
+ ,73
+ ,1559
+ ,164709
+ ,109
+ ,0
+ ,81
+ ,31
+ ,115
+ ,87771
+ ,36867
+ ,176
+ ,177
+ ,2014
+ ,201940
+ ,38
+ ,0
+ ,109
+ ,31
+ ,119
+ ,85872
+ ,33835
+ ,94
+ ,94
+ ,2143
+ ,235454
+ ,73
+ ,0
+ ,151
+ ,32
+ ,124
+ ,89275
+ ,24164
+ ,120
+ ,117
+ ,2146
+ ,220801
+ ,75
+ ,1
+ ,51
+ ,18
+ ,72
+ ,44418
+ ,12607
+ ,66
+ ,60
+ ,874
+ ,99466
+ ,50
+ ,0
+ ,28
+ ,23
+ ,91
+ ,192565
+ ,22609
+ ,56
+ ,55
+ ,1590
+ ,92661
+ ,61
+ ,1
+ ,40
+ ,17
+ ,45
+ ,35232
+ ,5892
+ ,39
+ ,39
+ ,1590
+ ,133328
+ ,55
+ ,0
+ ,56
+ ,20
+ ,78
+ ,40909
+ ,17014
+ ,66
+ ,64
+ ,1210
+ ,61361
+ ,77
+ ,0
+ ,27
+ ,12
+ ,39
+ ,13294
+ ,5394
+ ,27
+ ,26
+ ,2072
+ ,125930
+ ,75
+ ,4
+ ,37
+ ,17
+ ,68
+ ,32387
+ ,9178
+ ,65
+ ,64
+ ,1281
+ ,100750
+ ,72
+ ,0
+ ,83
+ ,30
+ ,119
+ ,140867
+ ,6440
+ ,58
+ ,58
+ ,1401
+ ,224549
+ ,50
+ ,4
+ ,54
+ ,31
+ ,117
+ ,120662
+ ,21916
+ ,98
+ ,95
+ ,834
+ ,82316
+ ,32
+ ,4
+ ,27
+ ,10
+ ,39
+ ,21233
+ ,4011
+ ,25
+ ,25
+ ,1105
+ ,102010
+ ,53
+ ,3
+ ,28
+ ,13
+ ,50
+ ,44332
+ ,5818
+ ,26
+ ,26
+ ,1272
+ ,101523
+ ,42
+ ,0
+ ,59
+ ,22
+ ,88
+ ,61056
+ ,18647
+ ,77
+ ,76
+ ,1944
+ ,243511
+ ,71
+ ,0
+ ,133
+ ,42
+ ,155
+ ,101338
+ ,20556
+ ,130
+ ,129
+ ,391
+ ,22938
+ ,10
+ ,0
+ ,12
+ ,1
+ ,0
+ ,1168
+ ,238
+ ,11
+ ,11
+ ,761
+ ,41566
+ ,35
+ ,5
+ ,0
+ ,9
+ ,36
+ ,13497
+ ,70
+ ,2
+ ,2
+ ,1605
+ ,152474
+ ,65
+ ,0
+ ,106
+ ,32
+ ,123
+ ,65567
+ ,22392
+ ,101
+ ,101
+ ,530
+ ,61857
+ ,25
+ ,4
+ ,23
+ ,11
+ ,32
+ ,25162
+ ,3913
+ ,31
+ ,28
+ ,1988
+ ,99923
+ ,66
+ ,0
+ ,44
+ ,25
+ ,99
+ ,32334
+ ,12237
+ ,36
+ ,36
+ ,1386
+ ,132487
+ ,41
+ ,0
+ ,71
+ ,36
+ ,136
+ ,40735
+ ,8388
+ ,120
+ ,89
+ ,2395
+ ,317394
+ ,86
+ ,1
+ ,116
+ ,31
+ ,117
+ ,91413
+ ,22120
+ ,195
+ ,193
+ ,387
+ ,21054
+ ,16
+ ,0
+ ,4
+ ,0
+ ,0
+ ,855
+ ,338
+ ,4
+ ,4
+ ,1742
+ ,209641
+ ,42
+ ,5
+ ,62
+ ,24
+ ,88
+ ,97068
+ ,11727
+ ,89
+ ,84
+ ,620
+ ,22648
+ ,19
+ ,0
+ ,12
+ ,13
+ ,39
+ ,44339
+ ,3704
+ ,24
+ ,23
+ ,449
+ ,31414
+ ,19
+ ,0
+ ,18
+ ,8
+ ,25
+ ,14116
+ ,3988
+ ,39
+ ,39
+ ,800
+ ,46698
+ ,45
+ ,0
+ ,14
+ ,13
+ ,52
+ ,10288
+ ,3030
+ ,14
+ ,14
+ ,1684
+ ,131698
+ ,65
+ ,0
+ ,60
+ ,19
+ ,75
+ ,65622
+ ,13520
+ ,78
+ ,78
+ ,1050
+ ,91735
+ ,35
+ ,0
+ ,7
+ ,18
+ ,71
+ ,16563
+ ,1421
+ ,15
+ ,14
+ ,2699
+ ,244749
+ ,95
+ ,2
+ ,98
+ ,33
+ ,124
+ ,76643
+ ,20923
+ ,106
+ ,101
+ ,1606
+ ,184510
+ ,49
+ ,7
+ ,64
+ ,40
+ ,151
+ ,110681
+ ,20237
+ ,83
+ ,82
+ ,1502
+ ,79863
+ ,37
+ ,1
+ ,29
+ ,22
+ ,71
+ ,29011
+ ,3219
+ ,24
+ ,24
+ ,1204
+ ,128423
+ ,64
+ ,8
+ ,32
+ ,38
+ ,145
+ ,92696
+ ,3769
+ ,37
+ ,36
+ ,1138
+ ,97839
+ ,38
+ ,2
+ ,25
+ ,24
+ ,87
+ ,94785
+ ,12252
+ ,77
+ ,75
+ ,568
+ ,38214
+ ,34
+ ,0
+ ,16
+ ,8
+ ,27
+ ,8773
+ ,1888
+ ,16
+ ,16
+ ,1459
+ ,151101
+ ,32
+ ,2
+ ,48
+ ,35
+ ,131
+ ,83209
+ ,14497
+ ,56
+ ,55
+ ,2158
+ ,272458
+ ,65
+ ,0
+ ,100
+ ,43
+ ,162
+ ,93815
+ ,28864
+ ,132
+ ,131
+ ,1111
+ ,172494
+ ,52
+ ,0
+ ,46
+ ,43
+ ,165
+ ,86687
+ ,21721
+ ,144
+ ,131
+ ,1421
+ ,108043
+ ,62
+ ,1
+ ,45
+ ,14
+ ,54
+ ,34553
+ ,4821
+ ,40
+ ,39
+ ,2833
+ ,328107
+ ,65
+ ,3
+ ,129
+ ,41
+ ,159
+ ,105547
+ ,33644
+ ,153
+ ,144
+ ,1955
+ ,250579
+ ,83
+ ,0
+ ,130
+ ,38
+ ,147
+ ,103487
+ ,15923
+ ,143
+ ,139
+ ,2922
+ ,351067
+ ,95
+ ,3
+ ,136
+ ,45
+ ,170
+ ,213688
+ ,42935
+ ,220
+ ,211
+ ,1002
+ ,158015
+ ,29
+ ,0
+ ,59
+ ,31
+ ,119
+ ,71220
+ ,18864
+ ,79
+ ,78
+ ,1060
+ ,98866
+ ,18
+ ,0
+ ,25
+ ,13
+ ,49
+ ,23517
+ ,4977
+ ,50
+ ,50
+ ,956
+ ,85439
+ ,33
+ ,0
+ ,32
+ ,28
+ ,104
+ ,56926
+ ,7785
+ ,39
+ ,39
+ ,2186
+ ,229242
+ ,247
+ ,4
+ ,63
+ ,31
+ ,120
+ ,91721
+ ,17939
+ ,95
+ ,90
+ ,3604
+ ,351619
+ ,139
+ ,4
+ ,95
+ ,40
+ ,150
+ ,115168
+ ,23436
+ ,169
+ ,166
+ ,1035
+ ,84207
+ ,29
+ ,11
+ ,14
+ ,30
+ ,112
+ ,111194
+ ,325
+ ,12
+ ,12
+ ,1417
+ ,120445
+ ,118
+ ,0
+ ,36
+ ,16
+ ,59
+ ,51009
+ ,13539
+ ,63
+ ,57
+ ,3261
+ ,324598
+ ,110
+ ,0
+ ,113
+ ,37
+ ,136
+ ,135777
+ ,34538
+ ,134
+ ,133
+ ,1587
+ ,131069
+ ,67
+ ,4
+ ,47
+ ,30
+ ,107
+ ,51513
+ ,12198
+ ,69
+ ,69
+ ,1424
+ ,204271
+ ,42
+ ,0
+ ,92
+ ,35
+ ,130
+ ,74163
+ ,26924
+ ,119
+ ,119
+ ,1701
+ ,165543
+ ,65
+ ,1
+ ,70
+ ,32
+ ,115
+ ,51633
+ ,12716
+ ,119
+ ,119
+ ,1249
+ ,141722
+ ,94
+ ,0
+ ,19
+ ,27
+ ,107
+ ,75345
+ ,8172
+ ,75
+ ,65
+ ,946
+ ,116048
+ ,64
+ ,0
+ ,50
+ ,20
+ ,75
+ ,33416
+ ,10855
+ ,63
+ ,61
+ ,1926
+ ,250047
+ ,81
+ ,0
+ ,41
+ ,18
+ ,71
+ ,83305
+ ,11932
+ ,55
+ ,49
+ ,3352
+ ,299775
+ ,95
+ ,9
+ ,91
+ ,31
+ ,120
+ ,98952
+ ,14300
+ ,103
+ ,101
+ ,1641
+ ,195838
+ ,67
+ ,1
+ ,111
+ ,31
+ ,116
+ ,102372
+ ,25515
+ ,197
+ ,196
+ ,2035
+ ,173260
+ ,63
+ ,3
+ ,41
+ ,21
+ ,79
+ ,37238
+ ,2805
+ ,16
+ ,15
+ ,2312
+ ,254488
+ ,83
+ ,10
+ ,120
+ ,39
+ ,150
+ ,103772
+ ,29402
+ ,140
+ ,136
+ ,1369
+ ,104389
+ ,45
+ ,5
+ ,135
+ ,41
+ ,156
+ ,123969
+ ,16440
+ ,89
+ ,89
+ ,1577
+ ,136084
+ ,30
+ ,0
+ ,27
+ ,13
+ ,51
+ ,27142
+ ,11221
+ ,40
+ ,40
+ ,2201
+ ,199476
+ ,70
+ ,2
+ ,87
+ ,32
+ ,118
+ ,135400
+ ,28732
+ ,125
+ ,123
+ ,961
+ ,92499
+ ,32
+ ,0
+ ,25
+ ,18
+ ,71
+ ,21399
+ ,5250
+ ,21
+ ,21
+ ,1900
+ ,224330
+ ,83
+ ,1
+ ,131
+ ,39
+ ,144
+ ,130115
+ ,28608
+ ,167
+ ,163
+ ,1254
+ ,135781
+ ,31
+ ,2
+ ,45
+ ,14
+ ,47
+ ,24874
+ ,8092
+ ,32
+ ,29
+ ,1335
+ ,74408
+ ,67
+ ,4
+ ,29
+ ,7
+ ,28
+ ,34988
+ ,4473
+ ,36
+ ,35
+ ,1597
+ ,81240
+ ,66
+ ,0
+ ,58
+ ,17
+ ,68
+ ,45549
+ ,1572
+ ,13
+ ,13
+ ,207
+ ,14688
+ ,10
+ ,0
+ ,4
+ ,0
+ ,0
+ ,6023
+ ,2065
+ ,5
+ ,5
+ ,1645
+ ,181633
+ ,70
+ ,2
+ ,47
+ ,30
+ ,110
+ ,64466
+ ,14817
+ ,96
+ ,96
+ ,2429
+ ,271856
+ ,103
+ ,1
+ ,109
+ ,37
+ ,147
+ ,54990
+ ,16714
+ ,151
+ ,151
+ ,151
+ ,7199
+ ,5
+ ,0
+ ,7
+ ,0
+ ,0
+ ,1644
+ ,556
+ ,6
+ ,6
+ ,474
+ ,46660
+ ,20
+ ,0
+ ,12
+ ,5
+ ,15
+ ,6179
+ ,2089
+ ,13
+ ,13
+ ,141
+ ,17547
+ ,5
+ ,0
+ ,0
+ ,1
+ ,4
+ ,3926
+ ,2658
+ ,3
+ ,3
+ ,1639
+ ,133368
+ ,36
+ ,1
+ ,37
+ ,16
+ ,64
+ ,32755
+ ,10695
+ ,57
+ ,56
+ ,872
+ ,95227
+ ,34
+ ,0
+ ,37
+ ,32
+ ,111
+ ,34777
+ ,1669
+ ,23
+ ,23
+ ,1318
+ ,152601
+ ,48
+ ,2
+ ,46
+ ,24
+ ,85
+ ,73224
+ ,16267
+ ,61
+ ,57
+ ,1018
+ ,98146
+ ,40
+ ,0
+ ,15
+ ,17
+ ,68
+ ,27114
+ ,7768
+ ,21
+ ,14
+ ,1383
+ ,79619
+ ,43
+ ,3
+ ,42
+ ,11
+ ,40
+ ,20760
+ ,7252
+ ,43
+ ,43
+ ,1314
+ ,59194
+ ,31
+ ,6
+ ,7
+ ,24
+ ,80
+ ,37636
+ ,6387
+ ,20
+ ,20
+ ,1335
+ ,139942
+ ,42
+ ,0
+ ,54
+ ,22
+ ,88
+ ,65461
+ ,18715
+ ,82
+ ,72
+ ,1403
+ ,118612
+ ,46
+ ,2
+ ,54
+ ,12
+ ,48
+ ,30080
+ ,7936
+ ,90
+ ,87
+ ,910
+ ,72880
+ ,33
+ ,0
+ ,14
+ ,19
+ ,76
+ ,24094
+ ,8643
+ ,25
+ ,21
+ ,616
+ ,65475
+ ,18
+ ,2
+ ,16
+ ,13
+ ,51
+ ,69008
+ ,7294
+ ,60
+ ,56
+ ,1407
+ ,99643
+ ,55
+ ,1
+ ,33
+ ,17
+ ,67
+ ,54968
+ ,4570
+ ,61
+ ,59
+ ,771
+ ,71965
+ ,35
+ ,1
+ ,32
+ ,15
+ ,59
+ ,46090
+ ,7185
+ ,85
+ ,82
+ ,766
+ ,77272
+ ,59
+ ,2
+ ,21
+ ,16
+ ,61
+ ,27507
+ ,10058
+ ,43
+ ,43
+ ,473
+ ,49289
+ ,19
+ ,1
+ ,15
+ ,24
+ ,76
+ ,10672
+ ,2342
+ ,25
+ ,25
+ ,1376
+ ,135131
+ ,66
+ ,0
+ ,38
+ ,15
+ ,60
+ ,34029
+ ,8509
+ ,41
+ ,38
+ ,1232
+ ,108446
+ ,60
+ ,1
+ ,22
+ ,17
+ ,68
+ ,46300
+ ,13275
+ ,26
+ ,25
+ ,1521
+ ,89746
+ ,36
+ ,3
+ ,28
+ ,18
+ ,71
+ ,24760
+ ,6816
+ ,38
+ ,38
+ ,572
+ ,44296
+ ,25
+ ,0
+ ,10
+ ,20
+ ,76
+ ,18779
+ ,1930
+ ,12
+ ,12
+ ,1059
+ ,77648
+ ,47
+ ,0
+ ,31
+ ,16
+ ,62
+ ,21280
+ ,8086
+ ,29
+ ,29
+ ,1544
+ ,181528
+ ,54
+ ,0
+ ,32
+ ,16
+ ,61
+ ,40662
+ ,10737
+ ,49
+ ,47
+ ,1230
+ ,134019
+ ,53
+ ,0
+ ,32
+ ,18
+ ,67
+ ,28987
+ ,8033
+ ,46
+ ,45
+ ,1206
+ ,124064
+ ,40
+ ,1
+ ,43
+ ,22
+ ,88
+ ,22827
+ ,7058
+ ,41
+ ,40
+ ,1205
+ ,92630
+ ,40
+ ,4
+ ,27
+ ,8
+ ,30
+ ,18513
+ ,6782
+ ,31
+ ,30
+ ,1255
+ ,121848
+ ,39
+ ,0
+ ,37
+ ,17
+ ,64
+ ,30594
+ ,5401
+ ,41
+ ,41
+ ,613
+ ,52915
+ ,14
+ ,0
+ ,20
+ ,18
+ ,68
+ ,24006
+ ,6521
+ ,26
+ ,25
+ ,721
+ ,81872
+ ,45
+ ,0
+ ,32
+ ,16
+ ,64
+ ,27913
+ ,10856
+ ,23
+ ,23
+ ,1109
+ ,58981
+ ,36
+ ,7
+ ,0
+ ,23
+ ,91
+ ,42744
+ ,2154
+ ,14
+ ,14
+ ,740
+ ,53515
+ ,28
+ ,2
+ ,5
+ ,22
+ ,88
+ ,12934
+ ,6117
+ ,16
+ ,16
+ ,1126
+ ,60812
+ ,44
+ ,0
+ ,26
+ ,13
+ ,52
+ ,22574
+ ,5238
+ ,25
+ ,26
+ ,728
+ ,56375
+ ,30
+ ,7
+ ,10
+ ,13
+ ,49
+ ,41385
+ ,4820
+ ,21
+ ,21
+ ,689
+ ,65490
+ ,22
+ ,3
+ ,27
+ ,16
+ ,62
+ ,18653
+ ,5615
+ ,32
+ ,27
+ ,592
+ ,80949
+ ,17
+ ,0
+ ,11
+ ,16
+ ,61
+ ,18472
+ ,4272
+ ,9
+ ,9
+ ,995
+ ,76302
+ ,31
+ ,0
+ ,29
+ ,20
+ ,76
+ ,30976
+ ,8702
+ ,35
+ ,33
+ ,1613
+ ,104011
+ ,55
+ ,6
+ ,25
+ ,22
+ ,88
+ ,63339
+ ,15340
+ ,42
+ ,42
+ ,2048
+ ,98104
+ ,54
+ ,2
+ ,55
+ ,17
+ ,66
+ ,25568
+ ,8030
+ ,68
+ ,68
+ ,705
+ ,67989
+ ,21
+ ,0
+ ,23
+ ,18
+ ,71
+ ,33747
+ ,9526
+ ,32
+ ,32
+ ,301
+ ,30989
+ ,14
+ ,0
+ ,5
+ ,17
+ ,68
+ ,4154
+ ,1278
+ ,6
+ ,6
+ ,1803
+ ,135458
+ ,81
+ ,3
+ ,43
+ ,12
+ ,48
+ ,19474
+ ,4236
+ ,68
+ ,67
+ ,799
+ ,73504
+ ,35
+ ,0
+ ,23
+ ,7
+ ,25
+ ,35130
+ ,3023
+ ,33
+ ,33
+ ,861
+ ,63123
+ ,43
+ ,1
+ ,34
+ ,17
+ ,68
+ ,39067
+ ,7196
+ ,84
+ ,77
+ ,1186
+ ,61254
+ ,46
+ ,1
+ ,36
+ ,14
+ ,41
+ ,13310
+ ,3394
+ ,46
+ ,46
+ ,1451
+ ,74914
+ ,30
+ ,0
+ ,35
+ ,23
+ ,90
+ ,65892
+ ,6371
+ ,30
+ ,30
+ ,628
+ ,31774
+ ,23
+ ,1
+ ,0
+ ,17
+ ,66
+ ,4143
+ ,1574
+ ,0
+ ,0
+ ,1161
+ ,81437
+ ,38
+ ,0
+ ,37
+ ,14
+ ,54
+ ,28579
+ ,9620
+ ,36
+ ,36
+ ,1463
+ ,87186
+ ,54
+ ,0
+ ,28
+ ,15
+ ,59
+ ,51776
+ ,6978
+ ,47
+ ,46
+ ,742
+ ,50090
+ ,20
+ ,0
+ ,16
+ ,17
+ ,60
+ ,21152
+ ,4911
+ ,20
+ ,18
+ ,979
+ ,65745
+ ,53
+ ,0
+ ,26
+ ,21
+ ,77
+ ,38084
+ ,8645
+ ,50
+ ,48
+ ,675
+ ,56653
+ ,45
+ ,0
+ ,38
+ ,18
+ ,68
+ ,27717
+ ,8987
+ ,30
+ ,29
+ ,1241
+ ,158399
+ ,39
+ ,0
+ ,23
+ ,18
+ ,72
+ ,32928
+ ,5544
+ ,30
+ ,28
+ ,676
+ ,46455
+ ,20
+ ,0
+ ,22
+ ,17
+ ,67
+ ,11342
+ ,3083
+ ,34
+ ,34
+ ,1049
+ ,73624
+ ,24
+ ,0
+ ,30
+ ,17
+ ,64
+ ,19499
+ ,6909
+ ,33
+ ,33
+ ,620
+ ,38395
+ ,31
+ ,0
+ ,16
+ ,16
+ ,63
+ ,16380
+ ,3189
+ ,34
+ ,34
+ ,1081
+ ,91899
+ ,35
+ ,0
+ ,18
+ ,15
+ ,59
+ ,36874
+ ,6745
+ ,37
+ ,33
+ ,1688
+ ,139526
+ ,151
+ ,0
+ ,28
+ ,21
+ ,84
+ ,48259
+ ,16724
+ ,83
+ ,80
+ ,736
+ ,52164
+ ,52
+ ,0
+ ,32
+ ,16
+ ,64
+ ,16734
+ ,4850
+ ,32
+ ,32
+ ,617
+ ,51567
+ ,30
+ ,2
+ ,21
+ ,14
+ ,56
+ ,28207
+ ,7025
+ ,30
+ ,30
+ ,812
+ ,70551
+ ,31
+ ,0
+ ,23
+ ,15
+ ,54
+ ,30143
+ ,6047
+ ,43
+ ,41
+ ,1051
+ ,84856
+ ,29
+ ,1
+ ,29
+ ,17
+ ,67
+ ,41369
+ ,7377
+ ,41
+ ,41
+ ,1656
+ ,102538
+ ,57
+ ,1
+ ,50
+ ,15
+ ,58
+ ,45833
+ ,9078
+ ,51
+ ,51
+ ,705
+ ,86678
+ ,40
+ ,0
+ ,12
+ ,15
+ ,59
+ ,29156
+ ,4605
+ ,19
+ ,18
+ ,945
+ ,85709
+ ,44
+ ,0
+ ,21
+ ,10
+ ,40
+ ,35944
+ ,3238
+ ,37
+ ,34
+ ,554
+ ,34662
+ ,25
+ ,0
+ ,18
+ ,6
+ ,22
+ ,36278
+ ,8100
+ ,33
+ ,31
+ ,1597
+ ,150580
+ ,77
+ ,0
+ ,27
+ ,22
+ ,83
+ ,45588
+ ,9653
+ ,41
+ ,39
+ ,982
+ ,99611
+ ,35
+ ,0
+ ,41
+ ,21
+ ,81
+ ,45097
+ ,8914
+ ,54
+ ,54
+ ,222
+ ,19349
+ ,11
+ ,0
+ ,13
+ ,1
+ ,2
+ ,3895
+ ,786
+ ,14
+ ,14
+ ,1212
+ ,99373
+ ,63
+ ,1
+ ,12
+ ,18
+ ,72
+ ,28394
+ ,6700
+ ,25
+ ,24
+ ,1143
+ ,86230
+ ,44
+ ,0
+ ,21
+ ,17
+ ,61
+ ,18632
+ ,5788
+ ,25
+ ,24
+ ,435
+ ,30837
+ ,19
+ ,0
+ ,8
+ ,4
+ ,15
+ ,2325
+ ,593
+ ,8
+ ,8
+ ,532
+ ,31706
+ ,13
+ ,0
+ ,26
+ ,10
+ ,32
+ ,25139
+ ,4506
+ ,26
+ ,26
+ ,882
+ ,89806
+ ,42
+ ,0
+ ,27
+ ,16
+ ,62
+ ,27975
+ ,6382
+ ,20
+ ,19
+ ,608
+ ,62088
+ ,38
+ ,1
+ ,13
+ ,16
+ ,58
+ ,14483
+ ,5621
+ ,11
+ ,11
+ ,459
+ ,40151
+ ,29
+ ,0
+ ,16
+ ,9
+ ,36
+ ,13127
+ ,3997
+ ,14
+ ,14
+ ,578
+ ,27634
+ ,20
+ ,0
+ ,2
+ ,16
+ ,59
+ ,5839
+ ,520
+ ,3
+ ,1
+ ,826
+ ,76990
+ ,27
+ ,0
+ ,42
+ ,17
+ ,68
+ ,24069
+ ,8891
+ ,40
+ ,39
+ ,509
+ ,37460
+ ,20
+ ,0
+ ,5
+ ,7
+ ,21
+ ,3738
+ ,999
+ ,5
+ ,5
+ ,717
+ ,54157
+ ,19
+ ,0
+ ,37
+ ,15
+ ,55
+ ,18625
+ ,7067
+ ,38
+ ,37
+ ,637
+ ,49862
+ ,37
+ ,0
+ ,17
+ ,14
+ ,54
+ ,36341
+ ,4639
+ ,32
+ ,32
+ ,857
+ ,84337
+ ,26
+ ,0
+ ,38
+ ,14
+ ,55
+ ,24548
+ ,5654
+ ,41
+ ,38
+ ,830
+ ,64175
+ ,42
+ ,0
+ ,37
+ ,18
+ ,72
+ ,21792
+ ,6928
+ ,46
+ ,47
+ ,652
+ ,59382
+ ,49
+ ,0
+ ,29
+ ,12
+ ,41
+ ,26263
+ ,1514
+ ,47
+ ,47
+ ,707
+ ,119308
+ ,30
+ ,0
+ ,32
+ ,16
+ ,61
+ ,23686
+ ,9238
+ ,37
+ ,37
+ ,954
+ ,76702
+ ,49
+ ,0
+ ,35
+ ,21
+ ,67
+ ,49303
+ ,8204
+ ,51
+ ,51
+ ,1461
+ ,103425
+ ,67
+ ,1
+ ,17
+ ,19
+ ,76
+ ,25659
+ ,5926
+ ,49
+ ,45
+ ,672
+ ,70344
+ ,28
+ ,0
+ ,20
+ ,16
+ ,64
+ ,28904
+ ,5785
+ ,21
+ ,21
+ ,778
+ ,43410
+ ,19
+ ,0
+ ,7
+ ,1
+ ,3
+ ,2781
+ ,4
+ ,1
+ ,1
+ ,1141
+ ,104838
+ ,49
+ ,1
+ ,46
+ ,16
+ ,63
+ ,29236
+ ,5930
+ ,44
+ ,42
+ ,680
+ ,62215
+ ,27
+ ,0
+ ,24
+ ,10
+ ,40
+ ,19546
+ ,3710
+ ,26
+ ,26
+ ,1090
+ ,69304
+ ,30
+ ,6
+ ,40
+ ,19
+ ,69
+ ,22818
+ ,705
+ ,21
+ ,21
+ ,616
+ ,53117
+ ,22
+ ,3
+ ,3
+ ,12
+ ,48
+ ,32689
+ ,443
+ ,4
+ ,4
+ ,285
+ ,19764
+ ,12
+ ,1
+ ,10
+ ,2
+ ,8
+ ,5752
+ ,2416
+ ,10
+ ,10
+ ,1145
+ ,86680
+ ,31
+ ,2
+ ,37
+ ,14
+ ,52
+ ,22197
+ ,7747
+ ,43
+ ,43
+ ,733
+ ,84105
+ ,20
+ ,0
+ ,17
+ ,17
+ ,66
+ ,20055
+ ,5432
+ ,34
+ ,34
+ ,888
+ ,77945
+ ,20
+ ,0
+ ,28
+ ,19
+ ,76
+ ,25272
+ ,4913
+ ,32
+ ,31
+ ,849
+ ,89113
+ ,39
+ ,0
+ ,19
+ ,14
+ ,43
+ ,82206
+ ,2650
+ ,20
+ ,19
+ ,1182
+ ,91005
+ ,29
+ ,3
+ ,29
+ ,11
+ ,39
+ ,32073
+ ,2370
+ ,34
+ ,34
+ ,528
+ ,40248
+ ,16
+ ,1
+ ,8
+ ,4
+ ,14
+ ,5444
+ ,775
+ ,6
+ ,6
+ ,642
+ ,64187
+ ,27
+ ,0
+ ,10
+ ,16
+ ,61
+ ,20154
+ ,5576
+ ,12
+ ,11
+ ,947
+ ,50857
+ ,21
+ ,0
+ ,15
+ ,20
+ ,71
+ ,36944
+ ,1352
+ ,24
+ ,24
+ ,819
+ ,56613
+ ,19
+ ,1
+ ,15
+ ,12
+ ,44
+ ,8019
+ ,3080
+ ,16
+ ,16
+ ,757
+ ,62792
+ ,35
+ ,0
+ ,28
+ ,15
+ ,60
+ ,30884
+ ,10205
+ ,72
+ ,72
+ ,894
+ ,72535
+ ,14
+ ,0
+ ,17
+ ,16
+ ,64
+ ,19540
+ ,6095
+ ,27
+ ,21)
+ ,dim=c(11
+ ,289)
+ ,dimnames=list(c('pageviews'
+ ,'time_in_rfc'
+ ,'logins'
+ ,'shared_compendiums'
+ ,'blogged_computations'
+ ,'compendiums_reviewed'
+ ,'feedback_messages_p1'
+ ,'totsize'
+ ,'totrevisions'
+ ,'tothyperlinks'
+ ,'totblogs')
+ ,1:289))
> y <- array(NA,dim=c(11,289),dimnames=list(c('pageviews','time_in_rfc','logins','shared_compendiums','blogged_computations','compendiums_reviewed','feedback_messages_p1','totsize','totrevisions','tothyperlinks','totblogs'),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'
> par3 <- 'No Linear Trend'
> par2 <- 'Do not include Seasonal Dummies'
> par1 <- '1'
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Prof. Dr. P. Wessa
> #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #Technical description: Write here your technical program description (don't use hard returns!)
> library(lattice)
> library(lmtest)
Loading required package: zoo
Attaching package: 'zoo'
The following object(s) are masked from 'package:base':
as.Date, 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
pageviews time_in_rfc logins shared_compendiums blogged_computations
1 1418 210907 56 3 79
2 869 120982 56 4 58
3 1530 176508 54 12 60
4 2172 179321 89 2 108
5 901 123185 40 1 49
6 463 52746 25 3 0
7 3201 385534 92 0 121
8 371 33170 18 0 1
9 1192 101645 63 0 20
10 1583 149061 44 5 43
11 1439 165446 33 0 69
12 1764 237213 84 0 78
13 1495 173326 88 7 86
14 1373 133131 55 7 44
15 2187 258873 60 3 104
16 1491 180083 66 9 63
17 4041 324799 154 0 158
18 1706 230964 53 4 102
19 2152 236785 119 3 77
20 1036 135473 41 0 82
21 1882 202925 61 7 115
22 1929 215147 58 0 101
23 2242 344297 75 1 80
24 1220 153935 33 5 50
25 1289 132943 40 7 83
26 2515 174724 92 0 123
27 2147 174415 100 0 73
28 2352 225548 112 5 81
29 1638 223632 73 0 105
30 1222 124817 40 0 47
31 1812 221698 45 0 105
32 1677 210767 60 3 94
33 1579 170266 62 4 44
34 1731 260561 75 1 114
35 807 84853 31 4 38
36 2452 294424 77 2 107
37 829 101011 34 0 30
38 1940 215641 46 0 71
39 2662 325107 99 0 84
40 186 7176 17 0 0
41 1499 167542 66 2 59
42 865 106408 30 1 33
43 1793 96560 76 0 42
44 2527 265769 146 2 96
45 2747 269651 67 10 106
46 1324 149112 56 6 56
47 2702 175824 107 0 57
48 1383 152871 58 5 59
49 1179 111665 34 4 39
50 2099 116408 61 1 34
51 4308 362301 119 2 76
52 918 78800 42 2 20
53 1831 183167 66 0 91
54 3373 277965 89 8 115
55 1713 150629 44 3 85
56 1438 168809 66 0 76
57 496 24188 24 0 8
58 2253 329267 259 8 79
59 744 65029 17 5 21
60 1161 101097 64 3 30
61 2352 218946 41 1 76
62 2144 244052 68 5 101
63 4691 341570 168 1 94
64 1112 103597 43 1 27
65 2694 233328 132 5 92
66 1973 256462 105 0 123
67 1769 206161 71 12 75
68 3148 311473 112 8 128
69 2474 235800 94 8 105
70 2084 177939 82 8 55
71 1954 207176 70 8 56
72 1226 196553 57 2 41
73 1389 174184 53 0 72
74 1496 143246 103 5 67
75 2269 187559 121 8 75
76 1833 187681 62 2 114
77 1268 119016 52 5 118
78 1943 182192 52 12 77
79 893 73566 32 6 22
80 1762 194979 62 7 66
81 1403 167488 45 2 69
82 1425 143756 46 0 105
83 1857 275541 63 4 116
84 1840 243199 75 3 88
85 1502 182999 88 6 73
86 1441 135649 46 2 99
87 1420 152299 53 0 62
88 1416 120221 37 1 53
89 2970 346485 90 0 118
90 1317 145790 63 5 30
91 1644 193339 78 2 100
92 870 80953 25 0 49
93 1654 122774 45 0 24
94 1054 130585 46 5 67
95 937 112611 41 0 46
96 3004 286468 144 1 57
97 2008 241066 82 0 75
98 2547 148446 91 1 135
99 1885 204713 71 1 68
100 1626 182079 63 2 124
101 1468 140344 53 6 33
102 2445 220516 62 1 98
103 1964 243060 63 4 58
104 1381 162765 32 2 68
105 1369 182613 39 3 81
106 1659 232138 62 0 131
107 2888 265318 117 10 110
108 1290 85574 34 0 37
109 2845 310839 92 9 130
110 1982 225060 93 7 93
111 1904 232317 54 0 118
112 1391 144966 144 0 39
113 602 43287 14 4 13
114 1743 155754 61 4 74
115 1559 164709 109 0 81
116 2014 201940 38 0 109
117 2143 235454 73 0 151
118 2146 220801 75 1 51
119 874 99466 50 0 28
120 1590 92661 61 1 40
121 1590 133328 55 0 56
122 1210 61361 77 0 27
123 2072 125930 75 4 37
124 1281 100750 72 0 83
125 1401 224549 50 4 54
126 834 82316 32 4 27
127 1105 102010 53 3 28
128 1272 101523 42 0 59
129 1944 243511 71 0 133
130 391 22938 10 0 12
131 761 41566 35 5 0
132 1605 152474 65 0 106
133 530 61857 25 4 23
134 1988 99923 66 0 44
135 1386 132487 41 0 71
136 2395 317394 86 1 116
137 387 21054 16 0 4
138 1742 209641 42 5 62
139 620 22648 19 0 12
140 449 31414 19 0 18
141 800 46698 45 0 14
142 1684 131698 65 0 60
143 1050 91735 35 0 7
144 2699 244749 95 2 98
145 1606 184510 49 7 64
146 1502 79863 37 1 29
147 1204 128423 64 8 32
148 1138 97839 38 2 25
149 568 38214 34 0 16
150 1459 151101 32 2 48
151 2158 272458 65 0 100
152 1111 172494 52 0 46
153 1421 108043 62 1 45
154 2833 328107 65 3 129
155 1955 250579 83 0 130
156 2922 351067 95 3 136
157 1002 158015 29 0 59
158 1060 98866 18 0 25
159 956 85439 33 0 32
160 2186 229242 247 4 63
161 3604 351619 139 4 95
162 1035 84207 29 11 14
163 1417 120445 118 0 36
164 3261 324598 110 0 113
165 1587 131069 67 4 47
166 1424 204271 42 0 92
167 1701 165543 65 1 70
168 1249 141722 94 0 19
169 946 116048 64 0 50
170 1926 250047 81 0 41
171 3352 299775 95 9 91
172 1641 195838 67 1 111
173 2035 173260 63 3 41
174 2312 254488 83 10 120
175 1369 104389 45 5 135
176 1577 136084 30 0 27
177 2201 199476 70 2 87
178 961 92499 32 0 25
179 1900 224330 83 1 131
180 1254 135781 31 2 45
181 1335 74408 67 4 29
182 1597 81240 66 0 58
183 207 14688 10 0 4
184 1645 181633 70 2 47
185 2429 271856 103 1 109
186 151 7199 5 0 7
187 474 46660 20 0 12
188 141 17547 5 0 0
189 1639 133368 36 1 37
190 872 95227 34 0 37
191 1318 152601 48 2 46
192 1018 98146 40 0 15
193 1383 79619 43 3 42
194 1314 59194 31 6 7
195 1335 139942 42 0 54
196 1403 118612 46 2 54
197 910 72880 33 0 14
198 616 65475 18 2 16
199 1407 99643 55 1 33
200 771 71965 35 1 32
201 766 77272 59 2 21
202 473 49289 19 1 15
203 1376 135131 66 0 38
204 1232 108446 60 1 22
205 1521 89746 36 3 28
206 572 44296 25 0 10
207 1059 77648 47 0 31
208 1544 181528 54 0 32
209 1230 134019 53 0 32
210 1206 124064 40 1 43
211 1205 92630 40 4 27
212 1255 121848 39 0 37
213 613 52915 14 0 20
214 721 81872 45 0 32
215 1109 58981 36 7 0
216 740 53515 28 2 5
217 1126 60812 44 0 26
218 728 56375 30 7 10
219 689 65490 22 3 27
220 592 80949 17 0 11
221 995 76302 31 0 29
222 1613 104011 55 6 25
223 2048 98104 54 2 55
224 705 67989 21 0 23
225 301 30989 14 0 5
226 1803 135458 81 3 43
227 799 73504 35 0 23
228 861 63123 43 1 34
229 1186 61254 46 1 36
230 1451 74914 30 0 35
231 628 31774 23 1 0
232 1161 81437 38 0 37
233 1463 87186 54 0 28
234 742 50090 20 0 16
235 979 65745 53 0 26
236 675 56653 45 0 38
237 1241 158399 39 0 23
238 676 46455 20 0 22
239 1049 73624 24 0 30
240 620 38395 31 0 16
241 1081 91899 35 0 18
242 1688 139526 151 0 28
243 736 52164 52 0 32
244 617 51567 30 2 21
245 812 70551 31 0 23
246 1051 84856 29 1 29
247 1656 102538 57 1 50
248 705 86678 40 0 12
249 945 85709 44 0 21
250 554 34662 25 0 18
251 1597 150580 77 0 27
252 982 99611 35 0 41
253 222 19349 11 0 13
254 1212 99373 63 1 12
255 1143 86230 44 0 21
256 435 30837 19 0 8
257 532 31706 13 0 26
258 882 89806 42 0 27
259 608 62088 38 1 13
260 459 40151 29 0 16
261 578 27634 20 0 2
262 826 76990 27 0 42
263 509 37460 20 0 5
264 717 54157 19 0 37
265 637 49862 37 0 17
266 857 84337 26 0 38
267 830 64175 42 0 37
268 652 59382 49 0 29
269 707 119308 30 0 32
270 954 76702 49 0 35
271 1461 103425 67 1 17
272 672 70344 28 0 20
273 778 43410 19 0 7
274 1141 104838 49 1 46
275 680 62215 27 0 24
276 1090 69304 30 6 40
277 616 53117 22 3 3
278 285 19764 12 1 10
279 1145 86680 31 2 37
280 733 84105 20 0 17
281 888 77945 20 0 28
282 849 89113 39 0 19
283 1182 91005 29 3 29
284 528 40248 16 1 8
285 642 64187 27 0 10
286 947 50857 21 0 15
287 819 56613 19 1 15
288 757 62792 35 0 28
289 894 72535 14 0 17
compendiums_reviewed feedback_messages_p1 totsize totrevisions
1 30 115 112285 24188
2 28 109 84786 18273
3 38 146 83123 14130
4 30 116 101193 32287
5 22 68 38361 8654
6 26 101 68504 9245
7 25 96 119182 33251
8 18 67 22807 1271
9 11 44 17140 5279
10 26 100 116174 27101
11 25 93 57635 16373
12 38 140 66198 19716
13 44 166 71701 17753
14 30 99 57793 9028
15 40 139 80444 18653
16 34 130 53855 8828
17 47 181 97668 29498
18 30 116 133824 27563
19 31 116 101481 18293
20 23 88 99645 22530
21 36 139 114789 15977
22 36 135 99052 35082
23 30 108 67654 16116
24 25 89 65553 15849
25 39 156 97500 16026
26 34 129 69112 26569
27 31 118 82753 24785
28 31 118 85323 17569
29 33 125 72654 23825
30 25 95 30727 7869
31 33 126 77873 14975
32 35 135 117478 37791
33 42 154 74007 9605
34 43 165 90183 27295
35 30 113 61542 2746
36 33 127 101494 34461
37 13 52 27570 8098
38 32 121 55813 4787
39 36 136 79215 24919
40 0 0 1423 603
41 28 108 55461 16329
42 14 46 31081 12558
43 17 54 22996 7784
44 32 124 83122 28522
45 30 115 70106 22265
46 35 128 60578 14459
47 20 80 39992 14526
48 28 97 79892 22240
49 28 104 49810 11802
50 39 59 71570 7623
51 34 125 100708 11912
52 26 82 33032 7935
53 39 149 82875 18220
54 39 149 139077 19199
55 33 122 71595 19918
56 28 118 72260 21884
57 4 12 5950 2694
58 39 144 115762 15808
59 18 67 32551 3597
60 14 52 31701 5296
61 29 108 80670 25239
62 44 166 143558 29801
63 21 80 117105 18450
64 16 60 23789 7132
65 28 107 120733 34861
66 35 127 105195 35940
67 28 107 73107 16688
68 38 146 132068 24683
69 23 84 149193 46230
70 36 141 46821 10387
71 32 123 87011 21436
72 29 111 95260 30546
73 25 98 55183 19746
74 27 105 106671 15977
75 36 135 73511 22583
76 28 107 92945 17274
77 23 85 78664 16469
78 40 155 70054 14251
79 23 88 22618 3007
80 40 155 74011 16851
81 28 104 83737 21113
82 34 132 69094 17401
83 33 127 93133 23958
84 28 108 95536 23567
85 34 129 225920 13065
86 30 116 62133 15358
87 33 122 61370 14587
88 22 85 43836 12770
89 38 147 106117 24021
90 26 99 38692 9648
91 35 87 84651 20537
92 8 28 56622 7905
93 24 90 15986 4527
94 29 109 95364 30495
95 20 78 26706 7117
96 29 111 89691 17719
97 45 158 67267 27056
98 37 141 126846 33473
99 33 122 41140 9758
100 33 124 102860 21115
101 25 93 51715 7236
102 32 124 55801 13790
103 29 112 111813 32902
104 28 108 120293 25131
105 28 99 138599 30910
106 31 117 161647 35947
107 52 199 115929 29848
108 21 78 24266 6943
109 24 91 162901 42705
110 41 158 109825 31808
111 33 126 129838 26675
112 32 122 37510 8435
113 19 71 43750 7409
114 20 75 40652 14993
115 31 115 87771 36867
116 31 119 85872 33835
117 32 124 89275 24164
118 18 72 44418 12607
119 23 91 192565 22609
120 17 45 35232 5892
121 20 78 40909 17014
122 12 39 13294 5394
123 17 68 32387 9178
124 30 119 140867 6440
125 31 117 120662 21916
126 10 39 21233 4011
127 13 50 44332 5818
128 22 88 61056 18647
129 42 155 101338 20556
130 1 0 1168 238
131 9 36 13497 70
132 32 123 65567 22392
133 11 32 25162 3913
134 25 99 32334 12237
135 36 136 40735 8388
136 31 117 91413 22120
137 0 0 855 338
138 24 88 97068 11727
139 13 39 44339 3704
140 8 25 14116 3988
141 13 52 10288 3030
142 19 75 65622 13520
143 18 71 16563 1421
144 33 124 76643 20923
145 40 151 110681 20237
146 22 71 29011 3219
147 38 145 92696 3769
148 24 87 94785 12252
149 8 27 8773 1888
150 35 131 83209 14497
151 43 162 93815 28864
152 43 165 86687 21721
153 14 54 34553 4821
154 41 159 105547 33644
155 38 147 103487 15923
156 45 170 213688 42935
157 31 119 71220 18864
158 13 49 23517 4977
159 28 104 56926 7785
160 31 120 91721 17939
161 40 150 115168 23436
162 30 112 111194 325
163 16 59 51009 13539
164 37 136 135777 34538
165 30 107 51513 12198
166 35 130 74163 26924
167 32 115 51633 12716
168 27 107 75345 8172
169 20 75 33416 10855
170 18 71 83305 11932
171 31 120 98952 14300
172 31 116 102372 25515
173 21 79 37238 2805
174 39 150 103772 29402
175 41 156 123969 16440
176 13 51 27142 11221
177 32 118 135400 28732
178 18 71 21399 5250
179 39 144 130115 28608
180 14 47 24874 8092
181 7 28 34988 4473
182 17 68 45549 1572
183 0 0 6023 2065
184 30 110 64466 14817
185 37 147 54990 16714
186 0 0 1644 556
187 5 15 6179 2089
188 1 4 3926 2658
189 16 64 32755 10695
190 32 111 34777 1669
191 24 85 73224 16267
192 17 68 27114 7768
193 11 40 20760 7252
194 24 80 37636 6387
195 22 88 65461 18715
196 12 48 30080 7936
197 19 76 24094 8643
198 13 51 69008 7294
199 17 67 54968 4570
200 15 59 46090 7185
201 16 61 27507 10058
202 24 76 10672 2342
203 15 60 34029 8509
204 17 68 46300 13275
205 18 71 24760 6816
206 20 76 18779 1930
207 16 62 21280 8086
208 16 61 40662 10737
209 18 67 28987 8033
210 22 88 22827 7058
211 8 30 18513 6782
212 17 64 30594 5401
213 18 68 24006 6521
214 16 64 27913 10856
215 23 91 42744 2154
216 22 88 12934 6117
217 13 52 22574 5238
218 13 49 41385 4820
219 16 62 18653 5615
220 16 61 18472 4272
221 20 76 30976 8702
222 22 88 63339 15340
223 17 66 25568 8030
224 18 71 33747 9526
225 17 68 4154 1278
226 12 48 19474 4236
227 7 25 35130 3023
228 17 68 39067 7196
229 14 41 13310 3394
230 23 90 65892 6371
231 17 66 4143 1574
232 14 54 28579 9620
233 15 59 51776 6978
234 17 60 21152 4911
235 21 77 38084 8645
236 18 68 27717 8987
237 18 72 32928 5544
238 17 67 11342 3083
239 17 64 19499 6909
240 16 63 16380 3189
241 15 59 36874 6745
242 21 84 48259 16724
243 16 64 16734 4850
244 14 56 28207 7025
245 15 54 30143 6047
246 17 67 41369 7377
247 15 58 45833 9078
248 15 59 29156 4605
249 10 40 35944 3238
250 6 22 36278 8100
251 22 83 45588 9653
252 21 81 45097 8914
253 1 2 3895 786
254 18 72 28394 6700
255 17 61 18632 5788
256 4 15 2325 593
257 10 32 25139 4506
258 16 62 27975 6382
259 16 58 14483 5621
260 9 36 13127 3997
261 16 59 5839 520
262 17 68 24069 8891
263 7 21 3738 999
264 15 55 18625 7067
265 14 54 36341 4639
266 14 55 24548 5654
267 18 72 21792 6928
268 12 41 26263 1514
269 16 61 23686 9238
270 21 67 49303 8204
271 19 76 25659 5926
272 16 64 28904 5785
273 1 3 2781 4
274 16 63 29236 5930
275 10 40 19546 3710
276 19 69 22818 705
277 12 48 32689 443
278 2 8 5752 2416
279 14 52 22197 7747
280 17 66 20055 5432
281 19 76 25272 4913
282 14 43 82206 2650
283 11 39 32073 2370
284 4 14 5444 775
285 16 61 20154 5576
286 20 71 36944 1352
287 12 44 8019 3080
288 15 60 30884 10205
289 16 64 19540 6095
tothyperlinks totblogs
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) time_in_rfc logins
2.810e+02 6.488e-03 5.549e+00
shared_compendiums blogged_computations compendiums_reviewed
1.548e+01 5.570e+00 1.151e+01
feedback_messages_p1 totsize totrevisions
-4.609e+00 -6.998e-04 -2.181e-03
tothyperlinks totblogs
1.142e+00 -3.886e+00
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-1455.81 -171.32 -31.93 135.17 1295.82
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 2.810e+02 4.460e+01 6.299 1.17e-09 ***
time_in_rfc 6.488e-03 4.888e-04 13.275 < 2e-16 ***
logins 5.549e+00 7.489e-01 7.409 1.53e-12 ***
shared_compendiums 1.548e+01 7.286e+00 2.124 0.0345 *
blogged_computations 5.570e+00 1.179e+00 4.724 3.67e-06 ***
compendiums_reviewed 1.151e+01 9.335e+00 1.233 0.2187
feedback_messages_p1 -4.609e+00 2.494e+00 -1.848 0.0657 .
totsize -6.998e-04 8.510e-04 -0.822 0.4116
totrevisions -2.181e-03 4.085e-03 -0.534 0.5939
tothyperlinks 1.142e+00 5.978e+00 0.191 0.8487
totblogs -3.886e+00 6.187e+00 -0.628 0.5304
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 296.3 on 278 degrees of freedom
Multiple R-squared: 0.8436, Adjusted R-squared: 0.838
F-statistic: 149.9 on 10 and 278 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.3722108 7.444216e-01 6.277892e-01
[2,] 0.3710013 7.420027e-01 6.289987e-01
[3,] 0.2386226 4.772452e-01 7.613774e-01
[4,] 0.4349896 8.699791e-01 5.650104e-01
[5,] 0.3426231 6.852462e-01 6.573769e-01
[6,] 0.5355084 9.289833e-01 4.644916e-01
[7,] 0.4749571 9.499143e-01 5.250429e-01
[8,] 0.4368842 8.737685e-01 5.631158e-01
[9,] 0.3759332 7.518664e-01 6.240668e-01
[10,] 0.3333247 6.666493e-01 6.666753e-01
[11,] 0.2630520 5.261041e-01 7.369480e-01
[12,] 0.2006092 4.012185e-01 7.993908e-01
[13,] 0.2013074 4.026149e-01 7.986926e-01
[14,] 0.1634610 3.269220e-01 8.365390e-01
[15,] 0.1398780 2.797561e-01 8.601220e-01
[16,] 0.3184825 6.369650e-01 6.815175e-01
[17,] 0.2740239 5.480477e-01 7.259761e-01
[18,] 0.2207222 4.414444e-01 7.792778e-01
[19,] 0.2373696 4.747391e-01 7.626304e-01
[20,] 0.3300025 6.600051e-01 6.699975e-01
[21,] 0.4694137 9.388273e-01 5.305863e-01
[22,] 0.4105101 8.210202e-01 5.894899e-01
[23,] 0.3526206 7.052412e-01 6.473794e-01
[24,] 0.3072615 6.145230e-01 6.927385e-01
[25,] 0.3234495 6.468989e-01 6.765505e-01
[26,] 0.2770705 5.541410e-01 7.229295e-01
[27,] 0.2342159 4.684317e-01 7.657841e-01
[28,] 0.1933982 3.867964e-01 8.066018e-01
[29,] 0.1600280 3.200561e-01 8.399720e-01
[30,] 0.2011157 4.022314e-01 7.988843e-01
[31,] 0.2333778 4.667555e-01 7.666222e-01
[32,] 0.3024328 6.048657e-01 6.975672e-01
[33,] 0.2612673 5.225345e-01 7.387327e-01
[34,] 0.4533321 9.066643e-01 5.466679e-01
[35,] 0.4134645 8.269290e-01 5.865355e-01
[36,] 0.3910705 7.821410e-01 6.089295e-01
[37,] 0.4218719 8.437437e-01 5.781281e-01
[38,] 0.7963002 4.073995e-01 2.036998e-01
[39,] 0.7630586 4.738829e-01 2.369414e-01
[40,] 0.7530324 4.939351e-01 2.469676e-01
[41,] 0.8749025 2.501950e-01 1.250975e-01
[42,] 0.8987019 2.025963e-01 1.012981e-01
[43,] 0.8802745 2.394509e-01 1.197255e-01
[44,] 0.8566246 2.867508e-01 1.433754e-01
[45,] 0.9998180 3.640261e-04 1.820130e-04
[46,] 0.9997312 5.376233e-04 2.688117e-04
[47,] 0.9996175 7.649161e-04 3.824581e-04
[48,] 0.9998669 2.661961e-04 1.330981e-04
[49,] 0.9998070 3.859710e-04 1.929855e-04
[50,] 0.9999997 5.838564e-07 2.919282e-07
[51,] 0.9999995 9.741227e-07 4.870614e-07
[52,] 1.0000000 7.252717e-08 3.626358e-08
[53,] 1.0000000 3.908778e-08 1.954389e-08
[54,] 1.0000000 3.348605e-08 1.674302e-08
[55,] 1.0000000 5.440655e-08 2.720328e-08
[56,] 1.0000000 6.852266e-08 3.426133e-08
[57,] 1.0000000 4.126860e-08 2.063430e-08
[58,] 1.0000000 5.412400e-08 2.706200e-08
[59,] 1.0000000 6.398804e-08 3.199402e-08
[60,] 1.0000000 8.743933e-08 4.371967e-08
[61,] 0.9999999 1.134763e-07 5.673814e-08
[62,] 0.9999999 1.356674e-07 6.783370e-08
[63,] 0.9999999 2.072890e-07 1.036445e-07
[64,] 0.9999999 1.044833e-07 5.224164e-08
[65,] 0.9999999 1.172723e-07 5.863614e-08
[66,] 0.9999999 1.722232e-07 8.611160e-08
[67,] 0.9999999 2.574124e-07 1.287062e-07
[68,] 0.9999998 4.117052e-07 2.058526e-07
[69,] 0.9999997 6.600898e-07 3.300449e-07
[70,] 0.9999999 1.015619e-07 5.078097e-08
[71,] 1.0000000 9.184309e-08 4.592155e-08
[72,] 1.0000000 3.784003e-08 1.892001e-08
[73,] 1.0000000 6.148729e-08 3.074365e-08
[74,] 0.9999999 1.016939e-07 5.084696e-08
[75,] 0.9999999 1.387735e-07 6.938676e-08
[76,] 0.9999999 2.241423e-07 1.120712e-07
[77,] 0.9999999 2.875193e-07 1.437596e-07
[78,] 0.9999999 1.937463e-07 9.687317e-08
[79,] 0.9999998 3.069292e-07 1.534646e-07
[80,] 0.9999999 1.802157e-07 9.010787e-08
[81,] 0.9999999 2.027043e-07 1.013521e-07
[82,] 0.9999999 2.327675e-07 1.163837e-07
[83,] 0.9999999 2.844615e-07 1.422307e-07
[84,] 0.9999998 3.547124e-07 1.773562e-07
[85,] 1.0000000 1.829892e-08 9.149459e-09
[86,] 1.0000000 2.912662e-08 1.456331e-08
[87,] 1.0000000 4.226312e-08 2.113156e-08
[88,] 1.0000000 6.277149e-08 3.138575e-08
[89,] 1.0000000 5.188252e-08 2.594126e-08
[90,] 1.0000000 8.289389e-08 4.144694e-08
[91,] 0.9999999 1.332129e-07 6.660643e-08
[92,] 0.9999999 1.989038e-07 9.945191e-08
[93,] 0.9999999 1.829704e-07 9.148518e-08
[94,] 0.9999999 2.261578e-07 1.130789e-07
[95,] 0.9999999 2.091529e-07 1.045765e-07
[96,] 0.9999999 2.906650e-07 1.453325e-07
[97,] 0.9999999 2.709493e-07 1.354746e-07
[98,] 0.9999998 4.296310e-07 2.148155e-07
[99,] 0.9999999 2.200903e-07 1.100451e-07
[100,] 0.9999998 3.145114e-07 1.572557e-07
[101,] 0.9999998 4.952340e-07 2.476170e-07
[102,] 0.9999996 7.336029e-07 3.668014e-07
[103,] 0.9999995 9.603610e-07 4.801805e-07
[104,] 0.9999995 1.061592e-06 5.307961e-07
[105,] 0.9999992 1.577095e-06 7.885476e-07
[106,] 0.9999988 2.405884e-06 1.202942e-06
[107,] 0.9999990 2.074450e-06 1.037225e-06
[108,] 0.9999987 2.598944e-06 1.299472e-06
[109,] 0.9999982 3.619487e-06 1.809744e-06
[110,] 0.9999995 9.406366e-07 4.703183e-07
[111,] 0.9999993 1.457601e-06 7.288003e-07
[112,] 0.9999996 8.024513e-07 4.012257e-07
[113,] 0.9999995 9.223510e-07 4.611755e-07
[114,] 0.9999994 1.281024e-06 6.405121e-07
[115,] 0.9999992 1.607577e-06 8.037884e-07
[116,] 0.9999993 1.348966e-06 6.744830e-07
[117,] 0.9999990 1.936611e-06 9.683053e-07
[118,] 0.9999985 2.919267e-06 1.459633e-06
[119,] 0.9999978 4.388049e-06 2.194025e-06
[120,] 0.9999983 3.421467e-06 1.710734e-06
[121,] 0.9999999 1.962622e-07 9.813110e-08
[122,] 0.9999999 1.937442e-07 9.687208e-08
[123,] 0.9999999 1.550257e-07 7.751286e-08
[124,] 0.9999999 2.284357e-07 1.142178e-07
[125,] 0.9999999 2.919143e-07 1.459571e-07
[126,] 0.9999998 3.886354e-07 1.943177e-07
[127,] 0.9999997 5.814200e-07 2.907100e-07
[128,] 0.9999996 8.991257e-07 4.495628e-07
[129,] 0.9999995 9.487743e-07 4.743872e-07
[130,] 0.9999993 1.390621e-06 6.953104e-07
[131,] 0.9999994 1.134400e-06 5.672001e-07
[132,] 0.9999993 1.364461e-06 6.822304e-07
[133,] 0.9999998 4.249875e-07 2.124937e-07
[134,] 0.9999998 3.035437e-07 1.517718e-07
[135,] 0.9999998 3.873339e-07 1.936669e-07
[136,] 0.9999997 5.646145e-07 2.823073e-07
[137,] 0.9999996 8.213975e-07 4.106987e-07
[138,] 0.9999994 1.191955e-06 5.959777e-07
[139,] 0.9999991 1.731612e-06 8.658058e-07
[140,] 0.9999987 2.539669e-06 1.269835e-06
[141,] 0.9999982 3.647111e-06 1.823555e-06
[142,] 0.9999989 2.141036e-06 1.070518e-06
[143,] 0.9999984 3.228071e-06 1.614036e-06
[144,] 0.9999989 2.271317e-06 1.135659e-06
[145,] 0.9999983 3.330492e-06 1.665246e-06
[146,] 0.9999975 5.043385e-06 2.521693e-06
[147,] 0.9999999 1.093613e-07 5.468063e-08
[148,] 1.0000000 7.829399e-08 3.914700e-08
[149,] 1.0000000 4.937416e-08 2.468708e-08
[150,] 1.0000000 6.978213e-08 3.489107e-08
[151,] 1.0000000 2.986366e-08 1.493183e-08
[152,] 1.0000000 4.687373e-08 2.343687e-08
[153,] 1.0000000 5.574620e-08 2.787310e-08
[154,] 1.0000000 8.688180e-08 4.344090e-08
[155,] 1.0000000 8.921572e-08 4.460786e-08
[156,] 1.0000000 5.742224e-08 2.871112e-08
[157,] 1.0000000 4.477995e-08 2.238998e-08
[158,] 1.0000000 3.582252e-08 1.791126e-08
[159,] 1.0000000 5.169617e-08 2.584809e-08
[160,] 1.0000000 6.799225e-08 3.399613e-08
[161,] 1.0000000 3.628376e-08 1.814188e-08
[162,] 1.0000000 1.956073e-08 9.780364e-09
[163,] 1.0000000 5.730488e-09 2.865244e-09
[164,] 1.0000000 4.518820e-09 2.259410e-09
[165,] 1.0000000 7.983050e-09 3.991525e-09
[166,] 1.0000000 2.421045e-09 1.210522e-09
[167,] 1.0000000 3.841263e-09 1.920631e-09
[168,] 1.0000000 6.557068e-09 3.278534e-09
[169,] 1.0000000 9.169112e-09 4.584556e-09
[170,] 1.0000000 1.428382e-08 7.141908e-09
[171,] 1.0000000 2.185132e-08 1.092566e-08
[172,] 1.0000000 8.428313e-09 4.214157e-09
[173,] 1.0000000 1.379933e-08 6.899666e-09
[174,] 1.0000000 2.242498e-08 1.121249e-08
[175,] 1.0000000 3.576083e-08 1.788042e-08
[176,] 1.0000000 1.918004e-08 9.590018e-09
[177,] 1.0000000 1.451005e-08 7.255027e-09
[178,] 1.0000000 1.560167e-08 7.800837e-09
[179,] 1.0000000 2.596734e-08 1.298367e-08
[180,] 1.0000000 3.026583e-08 1.513291e-08
[181,] 1.0000000 1.274833e-08 6.374163e-09
[182,] 1.0000000 2.299459e-08 1.149730e-08
[183,] 1.0000000 4.003817e-08 2.001908e-08
[184,] 1.0000000 5.524481e-08 2.762241e-08
[185,] 1.0000000 8.038796e-08 4.019398e-08
[186,] 0.9999999 1.195336e-07 5.976679e-08
[187,] 0.9999999 1.331632e-07 6.658158e-08
[188,] 0.9999999 1.107031e-07 5.535154e-08
[189,] 0.9999999 1.384225e-07 6.921126e-08
[190,] 0.9999999 2.375324e-07 1.187662e-07
[191,] 0.9999998 3.928072e-07 1.964036e-07
[192,] 0.9999999 1.491997e-07 7.459983e-08
[193,] 0.9999999 2.587372e-07 1.293686e-07
[194,] 0.9999998 4.250984e-07 2.125492e-07
[195,] 0.9999996 7.309792e-07 3.654896e-07
[196,] 0.9999994 1.148635e-06 5.743173e-07
[197,] 0.9999992 1.595948e-06 7.979738e-07
[198,] 0.9999987 2.555783e-06 1.277892e-06
[199,] 0.9999978 4.346053e-06 2.173027e-06
[200,] 0.9999963 7.355729e-06 3.677864e-06
[201,] 0.9999963 7.331383e-06 3.665692e-06
[202,] 0.9999948 1.040146e-05 5.200732e-06
[203,] 0.9999916 1.683387e-05 8.416933e-06
[204,] 0.9999919 1.616026e-05 8.080128e-06
[205,] 0.9999925 1.493640e-05 7.468198e-06
[206,] 0.9999930 1.402761e-05 7.013807e-06
[207,] 0.9999900 2.005733e-05 1.002866e-05
[208,] 0.9999843 3.137272e-05 1.568636e-05
[209,] 0.9999823 3.545944e-05 1.772972e-05
[210,] 0.9999997 5.824852e-07 2.912426e-07
[211,] 0.9999995 1.088693e-06 5.443467e-07
[212,] 0.9999992 1.517302e-06 7.586511e-07
[213,] 0.9999991 1.852452e-06 9.262260e-07
[214,] 0.9999983 3.310685e-06 1.655342e-06
[215,] 0.9999981 3.734311e-06 1.867156e-06
[216,] 0.9999985 3.059882e-06 1.529941e-06
[217,] 0.9999998 3.942929e-07 1.971464e-07
[218,] 0.9999997 6.403503e-07 3.201751e-07
[219,] 0.9999998 4.273186e-07 2.136593e-07
[220,] 1.0000000 4.171591e-08 2.085795e-08
[221,] 1.0000000 8.586827e-08 4.293413e-08
[222,] 0.9999999 1.815475e-07 9.077375e-08
[223,] 0.9999999 2.748571e-07 1.374286e-07
[224,] 0.9999998 4.447852e-07 2.223926e-07
[225,] 0.9999995 9.167322e-07 4.583661e-07
[226,] 0.9999997 6.735302e-07 3.367651e-07
[227,] 0.9999993 1.415988e-06 7.079938e-07
[228,] 0.9999986 2.720844e-06 1.360422e-06
[229,] 0.9999978 4.404225e-06 2.202112e-06
[230,] 0.9999961 7.792359e-06 3.896180e-06
[231,] 0.9999938 1.248943e-05 6.244714e-06
[232,] 0.9999879 2.423052e-05 1.211526e-05
[233,] 0.9999834 3.325146e-05 1.662573e-05
[234,] 0.9999997 6.071247e-07 3.035623e-07
[235,] 0.9999997 6.495275e-07 3.247638e-07
[236,] 0.9999993 1.407431e-06 7.037155e-07
[237,] 0.9999986 2.754152e-06 1.377076e-06
[238,] 0.9999975 5.089866e-06 2.544933e-06
[239,] 0.9999942 1.162438e-05 5.812189e-06
[240,] 0.9999912 1.753321e-05 8.766605e-06
[241,] 0.9999870 2.595655e-05 1.297828e-05
[242,] 0.9999881 2.384029e-05 1.192015e-05
[243,] 0.9999737 5.261056e-05 2.630528e-05
[244,] 0.9999416 1.167675e-04 5.838374e-05
[245,] 0.9998784 2.432008e-04 1.216004e-04
[246,] 0.9997852 4.295842e-04 2.147921e-04
[247,] 0.9995658 8.684883e-04 4.342442e-04
[248,] 0.9991810 1.638010e-03 8.190050e-04
[249,] 0.9984097 3.180517e-03 1.590259e-03
[250,] 0.9971314 5.737256e-03 2.868628e-03
[251,] 0.9943891 1.122177e-02 5.610884e-03
[252,] 0.9892527 2.149454e-02 1.074727e-02
[253,] 0.9815165 3.696701e-02 1.848350e-02
[254,] 0.9700771 5.984580e-02 2.992290e-02
[255,] 0.9954881 9.023823e-03 4.511911e-03
[256,] 0.9962472 7.505512e-03 3.752756e-03
[257,] 0.9935186 1.296289e-02 6.481445e-03
[258,] 0.9853970 2.920599e-02 1.460299e-02
[259,] 0.9681699 6.366028e-02 3.183014e-02
[260,] 0.9696401 6.071975e-02 3.035987e-02
[261,] 0.9234922 1.530156e-01 7.650778e-02
[262,] 0.8233647 3.532706e-01 1.766353e-01
> postscript(file="/var/wessaorg/rcomp/tmp/1fmhp1354796196.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/2l1do1354796196.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/3nwn91354796196.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/4pbno1354796196.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/5k4p01354796196.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
-313.284410 -338.317871 -122.118010 290.751912 -420.049289 -7.764006
7 8 9 10 11 12
-66.714743 -96.494105 -32.827636 306.534413 -50.157078 -314.135596
13 14 15 16 17 18
-428.987641 -94.442156 -95.319353 -319.357162 701.252179 -243.874338
19 20 21 22 23 24
-262.267308 -192.936803 -104.072414 61.841950 -592.254054 -207.043877
25 26 27 28 29 30
16.536069 536.477641 315.160131 45.651393 -493.923872 -25.538194
31 32 33 34 35 36
-149.790987 -175.346402 172.273383 -578.029474 -171.319726 -11.843909
37 38 39 40 41 42
-251.651339 23.831630 -175.214788 -233.523451 -118.137375 -259.513127
43 44 45 46 47 48
445.240918 -207.916516 170.193362 -163.122157 740.016236 -199.667491
49 50 51 52 53 54
72.145294 515.175068 1025.910290 -72.822433 172.789125 822.429217
55 56 57 58 59 60
273.524608 -139.506433 -81.083574 -1455.810618 -57.621074 -136.160986
61 62 63 64 65 66
555.487043 47.018579 1295.824628 -49.579925 539.676790 -466.579379
67 68 69 70 71 72
-405.247864 173.295431 239.351945 218.420294 83.671758 -297.097685
73 74 75 76 77 78
-229.080282 -143.095061 146.102341 -13.340176 -375.717236 151.406608
79 80 81 82 83 84
-27.435709 -34.679711 -57.704687 -3.744423 -644.362473 -329.434627
85 86 87 88 89 90
-427.580385 -39.471310 -24.900853 168.375569 -28.306793 -194.394246
91 92 93 94 95 96
-427.447873 -95.929133 426.895037 -228.361577 -263.000490 178.934346
97 98 99 100 101 102
-241.897594 804.775120 -80.513949 -151.456747 73.166324 366.583653
103 104 105 106 107 108
20.552940 46.879608 -106.771090 -315.652848 269.997873 307.175589
109 110 111 112 113 114
185.956413 -260.127294 -6.474415 -475.840212 118.049772 15.973884
115 116 117 118 119 120
-44.523325 188.457443 -275.590132 57.892846 2.914669 287.428494
121 122 123 124 125 126
195.356898 86.242422 624.838825 -40.396157 -404.910227 -213.977802
127 128 129 130 131 132
-138.344824 213.941281 -351.918708 -141.080222 16.227133 -45.725101
133 134 135 136 137 138
-360.546480 764.161393 90.709939 -261.862098 -129.292983 -107.138747
139 140 141 142 143 144
151.077433 -92.675294 30.697983 270.231161 112.741134 288.680375
145 146 147 148 149 150
-27.622285 487.840241 -165.575180 262.611281 -152.108774 162.475414
151 152 153 154 155 156
-69.532414 -115.485559 57.286287 90.826747 -412.064729 109.903283
157 158 159 160 161 162
-298.071764 139.251821 80.331147 -824.656621 494.322507 81.007337
163 164 165 166 167 168
-198.044004 369.525317 160.162252 -294.152512 132.207285 -158.876296
169 170 171 172 173 174
-393.780837 -323.033277 524.505146 -74.088092 200.213622 -155.805720
175 176 177 178 179 180
-54.205013 334.980077 390.203151 -32.742274 -230.976284 -194.784038
181 182 183 184 185 186
153.753554 288.450527 -224.571829 7.112123 -68.789291 -224.559922
187 188 189 190 191 192
-231.341590 -257.834357 381.080967 -187.092585 -152.381189 -21.102202
193 194 195 196 197 198
272.652158 532.825325 37.503216 122.369347 115.522900 -10.921466
199 200 201 202 203 204
296.347629 4.084141 -235.364180 -176.912403 -112.398290 23.109297
205 206 207 208 209 210
512.214277 -20.354553 54.556863 -116.979023 -130.964990 -64.393227
211 212 213 214 215 216
43.971447 6.060979 4.410485 -301.889568 365.064215 116.395584
217 218 219 220 221 222
251.323447 -75.857973 -140.417827 -225.787323 134.522968 465.405140
223 224 225 226 227 228
824.213720 -9.390097 -146.604489 196.567357 -124.743211 91.207614
229 230 231 232 233 234
207.222006 615.058552 112.691573 162.300023 436.666232 89.501732
235 236 237 238 239 240
120.778893 -211.154108 -177.758943 81.276502 208.726905 46.784861
241 242 243 244 245 246
135.168975 -60.128467 -129.199177 -98.512571 -6.016713 152.313518
247 248 249 250 251 252
386.076928 -249.091301 -61.758315 -32.615772 48.521828 -36.940067
253 254 255 256 257 258
-277.357769 78.345833 117.375302 -148.025948 -40.428945 -178.921042
259 260 261 262 263 264
-238.738812 -213.803789 89.127987 -78.315029 -119.066084 -17.111926
265 266 267 268 269 270
-56.299257 -104.298097 -21.254365 -246.086109 -457.437042 -31.931175
271 272 273 274 275 276
308.550540 -130.783218 78.008463 -111.103767 -125.689811 51.798920
277 278 279 280 281 282
-76.828391 -211.343156 121.616243 -71.539322 78.374695 -180.963695
283 284 285 286 287 288
115.828853 -122.433868 -108.581898 327.762566 87.232026 63.886848
289
158.664345
> postscript(file="/var/wessaorg/rcomp/tmp/6h2kq1354796196.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 -313.284410 NA
1 -338.317871 -313.284410
2 -122.118010 -338.317871
3 290.751912 -122.118010
4 -420.049289 290.751912
5 -7.764006 -420.049289
6 -66.714743 -7.764006
7 -96.494105 -66.714743
8 -32.827636 -96.494105
9 306.534413 -32.827636
10 -50.157078 306.534413
11 -314.135596 -50.157078
12 -428.987641 -314.135596
13 -94.442156 -428.987641
14 -95.319353 -94.442156
15 -319.357162 -95.319353
16 701.252179 -319.357162
17 -243.874338 701.252179
18 -262.267308 -243.874338
19 -192.936803 -262.267308
20 -104.072414 -192.936803
21 61.841950 -104.072414
22 -592.254054 61.841950
23 -207.043877 -592.254054
24 16.536069 -207.043877
25 536.477641 16.536069
26 315.160131 536.477641
27 45.651393 315.160131
28 -493.923872 45.651393
29 -25.538194 -493.923872
30 -149.790987 -25.538194
31 -175.346402 -149.790987
32 172.273383 -175.346402
33 -578.029474 172.273383
34 -171.319726 -578.029474
35 -11.843909 -171.319726
36 -251.651339 -11.843909
37 23.831630 -251.651339
38 -175.214788 23.831630
39 -233.523451 -175.214788
40 -118.137375 -233.523451
41 -259.513127 -118.137375
42 445.240918 -259.513127
43 -207.916516 445.240918
44 170.193362 -207.916516
45 -163.122157 170.193362
46 740.016236 -163.122157
47 -199.667491 740.016236
48 72.145294 -199.667491
49 515.175068 72.145294
50 1025.910290 515.175068
51 -72.822433 1025.910290
52 172.789125 -72.822433
53 822.429217 172.789125
54 273.524608 822.429217
55 -139.506433 273.524608
56 -81.083574 -139.506433
57 -1455.810618 -81.083574
58 -57.621074 -1455.810618
59 -136.160986 -57.621074
60 555.487043 -136.160986
61 47.018579 555.487043
62 1295.824628 47.018579
63 -49.579925 1295.824628
64 539.676790 -49.579925
65 -466.579379 539.676790
66 -405.247864 -466.579379
67 173.295431 -405.247864
68 239.351945 173.295431
69 218.420294 239.351945
70 83.671758 218.420294
71 -297.097685 83.671758
72 -229.080282 -297.097685
73 -143.095061 -229.080282
74 146.102341 -143.095061
75 -13.340176 146.102341
76 -375.717236 -13.340176
77 151.406608 -375.717236
78 -27.435709 151.406608
79 -34.679711 -27.435709
80 -57.704687 -34.679711
81 -3.744423 -57.704687
82 -644.362473 -3.744423
83 -329.434627 -644.362473
84 -427.580385 -329.434627
85 -39.471310 -427.580385
86 -24.900853 -39.471310
87 168.375569 -24.900853
88 -28.306793 168.375569
89 -194.394246 -28.306793
90 -427.447873 -194.394246
91 -95.929133 -427.447873
92 426.895037 -95.929133
93 -228.361577 426.895037
94 -263.000490 -228.361577
95 178.934346 -263.000490
96 -241.897594 178.934346
97 804.775120 -241.897594
98 -80.513949 804.775120
99 -151.456747 -80.513949
100 73.166324 -151.456747
101 366.583653 73.166324
102 20.552940 366.583653
103 46.879608 20.552940
104 -106.771090 46.879608
105 -315.652848 -106.771090
106 269.997873 -315.652848
107 307.175589 269.997873
108 185.956413 307.175589
109 -260.127294 185.956413
110 -6.474415 -260.127294
111 -475.840212 -6.474415
112 118.049772 -475.840212
113 15.973884 118.049772
114 -44.523325 15.973884
115 188.457443 -44.523325
116 -275.590132 188.457443
117 57.892846 -275.590132
118 2.914669 57.892846
119 287.428494 2.914669
120 195.356898 287.428494
121 86.242422 195.356898
122 624.838825 86.242422
123 -40.396157 624.838825
124 -404.910227 -40.396157
125 -213.977802 -404.910227
126 -138.344824 -213.977802
127 213.941281 -138.344824
128 -351.918708 213.941281
129 -141.080222 -351.918708
130 16.227133 -141.080222
131 -45.725101 16.227133
132 -360.546480 -45.725101
133 764.161393 -360.546480
134 90.709939 764.161393
135 -261.862098 90.709939
136 -129.292983 -261.862098
137 -107.138747 -129.292983
138 151.077433 -107.138747
139 -92.675294 151.077433
140 30.697983 -92.675294
141 270.231161 30.697983
142 112.741134 270.231161
143 288.680375 112.741134
144 -27.622285 288.680375
145 487.840241 -27.622285
146 -165.575180 487.840241
147 262.611281 -165.575180
148 -152.108774 262.611281
149 162.475414 -152.108774
150 -69.532414 162.475414
151 -115.485559 -69.532414
152 57.286287 -115.485559
153 90.826747 57.286287
154 -412.064729 90.826747
155 109.903283 -412.064729
156 -298.071764 109.903283
157 139.251821 -298.071764
158 80.331147 139.251821
159 -824.656621 80.331147
160 494.322507 -824.656621
161 81.007337 494.322507
162 -198.044004 81.007337
163 369.525317 -198.044004
164 160.162252 369.525317
165 -294.152512 160.162252
166 132.207285 -294.152512
167 -158.876296 132.207285
168 -393.780837 -158.876296
169 -323.033277 -393.780837
170 524.505146 -323.033277
171 -74.088092 524.505146
172 200.213622 -74.088092
173 -155.805720 200.213622
174 -54.205013 -155.805720
175 334.980077 -54.205013
176 390.203151 334.980077
177 -32.742274 390.203151
178 -230.976284 -32.742274
179 -194.784038 -230.976284
180 153.753554 -194.784038
181 288.450527 153.753554
182 -224.571829 288.450527
183 7.112123 -224.571829
184 -68.789291 7.112123
185 -224.559922 -68.789291
186 -231.341590 -224.559922
187 -257.834357 -231.341590
188 381.080967 -257.834357
189 -187.092585 381.080967
190 -152.381189 -187.092585
191 -21.102202 -152.381189
192 272.652158 -21.102202
193 532.825325 272.652158
194 37.503216 532.825325
195 122.369347 37.503216
196 115.522900 122.369347
197 -10.921466 115.522900
198 296.347629 -10.921466
199 4.084141 296.347629
200 -235.364180 4.084141
201 -176.912403 -235.364180
202 -112.398290 -176.912403
203 23.109297 -112.398290
204 512.214277 23.109297
205 -20.354553 512.214277
206 54.556863 -20.354553
207 -116.979023 54.556863
208 -130.964990 -116.979023
209 -64.393227 -130.964990
210 43.971447 -64.393227
211 6.060979 43.971447
212 4.410485 6.060979
213 -301.889568 4.410485
214 365.064215 -301.889568
215 116.395584 365.064215
216 251.323447 116.395584
217 -75.857973 251.323447
218 -140.417827 -75.857973
219 -225.787323 -140.417827
220 134.522968 -225.787323
221 465.405140 134.522968
222 824.213720 465.405140
223 -9.390097 824.213720
224 -146.604489 -9.390097
225 196.567357 -146.604489
226 -124.743211 196.567357
227 91.207614 -124.743211
228 207.222006 91.207614
229 615.058552 207.222006
230 112.691573 615.058552
231 162.300023 112.691573
232 436.666232 162.300023
233 89.501732 436.666232
234 120.778893 89.501732
235 -211.154108 120.778893
236 -177.758943 -211.154108
237 81.276502 -177.758943
238 208.726905 81.276502
239 46.784861 208.726905
240 135.168975 46.784861
241 -60.128467 135.168975
242 -129.199177 -60.128467
243 -98.512571 -129.199177
244 -6.016713 -98.512571
245 152.313518 -6.016713
246 386.076928 152.313518
247 -249.091301 386.076928
248 -61.758315 -249.091301
249 -32.615772 -61.758315
250 48.521828 -32.615772
251 -36.940067 48.521828
252 -277.357769 -36.940067
253 78.345833 -277.357769
254 117.375302 78.345833
255 -148.025948 117.375302
256 -40.428945 -148.025948
257 -178.921042 -40.428945
258 -238.738812 -178.921042
259 -213.803789 -238.738812
260 89.127987 -213.803789
261 -78.315029 89.127987
262 -119.066084 -78.315029
263 -17.111926 -119.066084
264 -56.299257 -17.111926
265 -104.298097 -56.299257
266 -21.254365 -104.298097
267 -246.086109 -21.254365
268 -457.437042 -246.086109
269 -31.931175 -457.437042
270 308.550540 -31.931175
271 -130.783218 308.550540
272 78.008463 -130.783218
273 -111.103767 78.008463
274 -125.689811 -111.103767
275 51.798920 -125.689811
276 -76.828391 51.798920
277 -211.343156 -76.828391
278 121.616243 -211.343156
279 -71.539322 121.616243
280 78.374695 -71.539322
281 -180.963695 78.374695
282 115.828853 -180.963695
283 -122.433868 115.828853
284 -108.581898 -122.433868
285 327.762566 -108.581898
286 87.232026 327.762566
287 63.886848 87.232026
288 158.664345 63.886848
289 NA 158.664345
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] -338.317871 -313.284410
[2,] -122.118010 -338.317871
[3,] 290.751912 -122.118010
[4,] -420.049289 290.751912
[5,] -7.764006 -420.049289
[6,] -66.714743 -7.764006
[7,] -96.494105 -66.714743
[8,] -32.827636 -96.494105
[9,] 306.534413 -32.827636
[10,] -50.157078 306.534413
[11,] -314.135596 -50.157078
[12,] -428.987641 -314.135596
[13,] -94.442156 -428.987641
[14,] -95.319353 -94.442156
[15,] -319.357162 -95.319353
[16,] 701.252179 -319.357162
[17,] -243.874338 701.252179
[18,] -262.267308 -243.874338
[19,] -192.936803 -262.267308
[20,] -104.072414 -192.936803
[21,] 61.841950 -104.072414
[22,] -592.254054 61.841950
[23,] -207.043877 -592.254054
[24,] 16.536069 -207.043877
[25,] 536.477641 16.536069
[26,] 315.160131 536.477641
[27,] 45.651393 315.160131
[28,] -493.923872 45.651393
[29,] -25.538194 -493.923872
[30,] -149.790987 -25.538194
[31,] -175.346402 -149.790987
[32,] 172.273383 -175.346402
[33,] -578.029474 172.273383
[34,] -171.319726 -578.029474
[35,] -11.843909 -171.319726
[36,] -251.651339 -11.843909
[37,] 23.831630 -251.651339
[38,] -175.214788 23.831630
[39,] -233.523451 -175.214788
[40,] -118.137375 -233.523451
[41,] -259.513127 -118.137375
[42,] 445.240918 -259.513127
[43,] -207.916516 445.240918
[44,] 170.193362 -207.916516
[45,] -163.122157 170.193362
[46,] 740.016236 -163.122157
[47,] -199.667491 740.016236
[48,] 72.145294 -199.667491
[49,] 515.175068 72.145294
[50,] 1025.910290 515.175068
[51,] -72.822433 1025.910290
[52,] 172.789125 -72.822433
[53,] 822.429217 172.789125
[54,] 273.524608 822.429217
[55,] -139.506433 273.524608
[56,] -81.083574 -139.506433
[57,] -1455.810618 -81.083574
[58,] -57.621074 -1455.810618
[59,] -136.160986 -57.621074
[60,] 555.487043 -136.160986
[61,] 47.018579 555.487043
[62,] 1295.824628 47.018579
[63,] -49.579925 1295.824628
[64,] 539.676790 -49.579925
[65,] -466.579379 539.676790
[66,] -405.247864 -466.579379
[67,] 173.295431 -405.247864
[68,] 239.351945 173.295431
[69,] 218.420294 239.351945
[70,] 83.671758 218.420294
[71,] -297.097685 83.671758
[72,] -229.080282 -297.097685
[73,] -143.095061 -229.080282
[74,] 146.102341 -143.095061
[75,] -13.340176 146.102341
[76,] -375.717236 -13.340176
[77,] 151.406608 -375.717236
[78,] -27.435709 151.406608
[79,] -34.679711 -27.435709
[80,] -57.704687 -34.679711
[81,] -3.744423 -57.704687
[82,] -644.362473 -3.744423
[83,] -329.434627 -644.362473
[84,] -427.580385 -329.434627
[85,] -39.471310 -427.580385
[86,] -24.900853 -39.471310
[87,] 168.375569 -24.900853
[88,] -28.306793 168.375569
[89,] -194.394246 -28.306793
[90,] -427.447873 -194.394246
[91,] -95.929133 -427.447873
[92,] 426.895037 -95.929133
[93,] -228.361577 426.895037
[94,] -263.000490 -228.361577
[95,] 178.934346 -263.000490
[96,] -241.897594 178.934346
[97,] 804.775120 -241.897594
[98,] -80.513949 804.775120
[99,] -151.456747 -80.513949
[100,] 73.166324 -151.456747
[101,] 366.583653 73.166324
[102,] 20.552940 366.583653
[103,] 46.879608 20.552940
[104,] -106.771090 46.879608
[105,] -315.652848 -106.771090
[106,] 269.997873 -315.652848
[107,] 307.175589 269.997873
[108,] 185.956413 307.175589
[109,] -260.127294 185.956413
[110,] -6.474415 -260.127294
[111,] -475.840212 -6.474415
[112,] 118.049772 -475.840212
[113,] 15.973884 118.049772
[114,] -44.523325 15.973884
[115,] 188.457443 -44.523325
[116,] -275.590132 188.457443
[117,] 57.892846 -275.590132
[118,] 2.914669 57.892846
[119,] 287.428494 2.914669
[120,] 195.356898 287.428494
[121,] 86.242422 195.356898
[122,] 624.838825 86.242422
[123,] -40.396157 624.838825
[124,] -404.910227 -40.396157
[125,] -213.977802 -404.910227
[126,] -138.344824 -213.977802
[127,] 213.941281 -138.344824
[128,] -351.918708 213.941281
[129,] -141.080222 -351.918708
[130,] 16.227133 -141.080222
[131,] -45.725101 16.227133
[132,] -360.546480 -45.725101
[133,] 764.161393 -360.546480
[134,] 90.709939 764.161393
[135,] -261.862098 90.709939
[136,] -129.292983 -261.862098
[137,] -107.138747 -129.292983
[138,] 151.077433 -107.138747
[139,] -92.675294 151.077433
[140,] 30.697983 -92.675294
[141,] 270.231161 30.697983
[142,] 112.741134 270.231161
[143,] 288.680375 112.741134
[144,] -27.622285 288.680375
[145,] 487.840241 -27.622285
[146,] -165.575180 487.840241
[147,] 262.611281 -165.575180
[148,] -152.108774 262.611281
[149,] 162.475414 -152.108774
[150,] -69.532414 162.475414
[151,] -115.485559 -69.532414
[152,] 57.286287 -115.485559
[153,] 90.826747 57.286287
[154,] -412.064729 90.826747
[155,] 109.903283 -412.064729
[156,] -298.071764 109.903283
[157,] 139.251821 -298.071764
[158,] 80.331147 139.251821
[159,] -824.656621 80.331147
[160,] 494.322507 -824.656621
[161,] 81.007337 494.322507
[162,] -198.044004 81.007337
[163,] 369.525317 -198.044004
[164,] 160.162252 369.525317
[165,] -294.152512 160.162252
[166,] 132.207285 -294.152512
[167,] -158.876296 132.207285
[168,] -393.780837 -158.876296
[169,] -323.033277 -393.780837
[170,] 524.505146 -323.033277
[171,] -74.088092 524.505146
[172,] 200.213622 -74.088092
[173,] -155.805720 200.213622
[174,] -54.205013 -155.805720
[175,] 334.980077 -54.205013
[176,] 390.203151 334.980077
[177,] -32.742274 390.203151
[178,] -230.976284 -32.742274
[179,] -194.784038 -230.976284
[180,] 153.753554 -194.784038
[181,] 288.450527 153.753554
[182,] -224.571829 288.450527
[183,] 7.112123 -224.571829
[184,] -68.789291 7.112123
[185,] -224.559922 -68.789291
[186,] -231.341590 -224.559922
[187,] -257.834357 -231.341590
[188,] 381.080967 -257.834357
[189,] -187.092585 381.080967
[190,] -152.381189 -187.092585
[191,] -21.102202 -152.381189
[192,] 272.652158 -21.102202
[193,] 532.825325 272.652158
[194,] 37.503216 532.825325
[195,] 122.369347 37.503216
[196,] 115.522900 122.369347
[197,] -10.921466 115.522900
[198,] 296.347629 -10.921466
[199,] 4.084141 296.347629
[200,] -235.364180 4.084141
[201,] -176.912403 -235.364180
[202,] -112.398290 -176.912403
[203,] 23.109297 -112.398290
[204,] 512.214277 23.109297
[205,] -20.354553 512.214277
[206,] 54.556863 -20.354553
[207,] -116.979023 54.556863
[208,] -130.964990 -116.979023
[209,] -64.393227 -130.964990
[210,] 43.971447 -64.393227
[211,] 6.060979 43.971447
[212,] 4.410485 6.060979
[213,] -301.889568 4.410485
[214,] 365.064215 -301.889568
[215,] 116.395584 365.064215
[216,] 251.323447 116.395584
[217,] -75.857973 251.323447
[218,] -140.417827 -75.857973
[219,] -225.787323 -140.417827
[220,] 134.522968 -225.787323
[221,] 465.405140 134.522968
[222,] 824.213720 465.405140
[223,] -9.390097 824.213720
[224,] -146.604489 -9.390097
[225,] 196.567357 -146.604489
[226,] -124.743211 196.567357
[227,] 91.207614 -124.743211
[228,] 207.222006 91.207614
[229,] 615.058552 207.222006
[230,] 112.691573 615.058552
[231,] 162.300023 112.691573
[232,] 436.666232 162.300023
[233,] 89.501732 436.666232
[234,] 120.778893 89.501732
[235,] -211.154108 120.778893
[236,] -177.758943 -211.154108
[237,] 81.276502 -177.758943
[238,] 208.726905 81.276502
[239,] 46.784861 208.726905
[240,] 135.168975 46.784861
[241,] -60.128467 135.168975
[242,] -129.199177 -60.128467
[243,] -98.512571 -129.199177
[244,] -6.016713 -98.512571
[245,] 152.313518 -6.016713
[246,] 386.076928 152.313518
[247,] -249.091301 386.076928
[248,] -61.758315 -249.091301
[249,] -32.615772 -61.758315
[250,] 48.521828 -32.615772
[251,] -36.940067 48.521828
[252,] -277.357769 -36.940067
[253,] 78.345833 -277.357769
[254,] 117.375302 78.345833
[255,] -148.025948 117.375302
[256,] -40.428945 -148.025948
[257,] -178.921042 -40.428945
[258,] -238.738812 -178.921042
[259,] -213.803789 -238.738812
[260,] 89.127987 -213.803789
[261,] -78.315029 89.127987
[262,] -119.066084 -78.315029
[263,] -17.111926 -119.066084
[264,] -56.299257 -17.111926
[265,] -104.298097 -56.299257
[266,] -21.254365 -104.298097
[267,] -246.086109 -21.254365
[268,] -457.437042 -246.086109
[269,] -31.931175 -457.437042
[270,] 308.550540 -31.931175
[271,] -130.783218 308.550540
[272,] 78.008463 -130.783218
[273,] -111.103767 78.008463
[274,] -125.689811 -111.103767
[275,] 51.798920 -125.689811
[276,] -76.828391 51.798920
[277,] -211.343156 -76.828391
[278,] 121.616243 -211.343156
[279,] -71.539322 121.616243
[280,] 78.374695 -71.539322
[281,] -180.963695 78.374695
[282,] 115.828853 -180.963695
[283,] -122.433868 115.828853
[284,] -108.581898 -122.433868
[285,] 327.762566 -108.581898
[286,] 87.232026 327.762566
[287,] 63.886848 87.232026
[288,] 158.664345 63.886848
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 -338.317871 -313.284410
2 -122.118010 -338.317871
3 290.751912 -122.118010
4 -420.049289 290.751912
5 -7.764006 -420.049289
6 -66.714743 -7.764006
7 -96.494105 -66.714743
8 -32.827636 -96.494105
9 306.534413 -32.827636
10 -50.157078 306.534413
11 -314.135596 -50.157078
12 -428.987641 -314.135596
13 -94.442156 -428.987641
14 -95.319353 -94.442156
15 -319.357162 -95.319353
16 701.252179 -319.357162
17 -243.874338 701.252179
18 -262.267308 -243.874338
19 -192.936803 -262.267308
20 -104.072414 -192.936803
21 61.841950 -104.072414
22 -592.254054 61.841950
23 -207.043877 -592.254054
24 16.536069 -207.043877
25 536.477641 16.536069
26 315.160131 536.477641
27 45.651393 315.160131
28 -493.923872 45.651393
29 -25.538194 -493.923872
30 -149.790987 -25.538194
31 -175.346402 -149.790987
32 172.273383 -175.346402
33 -578.029474 172.273383
34 -171.319726 -578.029474
35 -11.843909 -171.319726
36 -251.651339 -11.843909
37 23.831630 -251.651339
38 -175.214788 23.831630
39 -233.523451 -175.214788
40 -118.137375 -233.523451
41 -259.513127 -118.137375
42 445.240918 -259.513127
43 -207.916516 445.240918
44 170.193362 -207.916516
45 -163.122157 170.193362
46 740.016236 -163.122157
47 -199.667491 740.016236
48 72.145294 -199.667491
49 515.175068 72.145294
50 1025.910290 515.175068
51 -72.822433 1025.910290
52 172.789125 -72.822433
53 822.429217 172.789125
54 273.524608 822.429217
55 -139.506433 273.524608
56 -81.083574 -139.506433
57 -1455.810618 -81.083574
58 -57.621074 -1455.810618
59 -136.160986 -57.621074
60 555.487043 -136.160986
61 47.018579 555.487043
62 1295.824628 47.018579
63 -49.579925 1295.824628
64 539.676790 -49.579925
65 -466.579379 539.676790
66 -405.247864 -466.579379
67 173.295431 -405.247864
68 239.351945 173.295431
69 218.420294 239.351945
70 83.671758 218.420294
71 -297.097685 83.671758
72 -229.080282 -297.097685
73 -143.095061 -229.080282
74 146.102341 -143.095061
75 -13.340176 146.102341
76 -375.717236 -13.340176
77 151.406608 -375.717236
78 -27.435709 151.406608
79 -34.679711 -27.435709
80 -57.704687 -34.679711
81 -3.744423 -57.704687
82 -644.362473 -3.744423
83 -329.434627 -644.362473
84 -427.580385 -329.434627
85 -39.471310 -427.580385
86 -24.900853 -39.471310
87 168.375569 -24.900853
88 -28.306793 168.375569
89 -194.394246 -28.306793
90 -427.447873 -194.394246
91 -95.929133 -427.447873
92 426.895037 -95.929133
93 -228.361577 426.895037
94 -263.000490 -228.361577
95 178.934346 -263.000490
96 -241.897594 178.934346
97 804.775120 -241.897594
98 -80.513949 804.775120
99 -151.456747 -80.513949
100 73.166324 -151.456747
101 366.583653 73.166324
102 20.552940 366.583653
103 46.879608 20.552940
104 -106.771090 46.879608
105 -315.652848 -106.771090
106 269.997873 -315.652848
107 307.175589 269.997873
108 185.956413 307.175589
109 -260.127294 185.956413
110 -6.474415 -260.127294
111 -475.840212 -6.474415
112 118.049772 -475.840212
113 15.973884 118.049772
114 -44.523325 15.973884
115 188.457443 -44.523325
116 -275.590132 188.457443
117 57.892846 -275.590132
118 2.914669 57.892846
119 287.428494 2.914669
120 195.356898 287.428494
121 86.242422 195.356898
122 624.838825 86.242422
123 -40.396157 624.838825
124 -404.910227 -40.396157
125 -213.977802 -404.910227
126 -138.344824 -213.977802
127 213.941281 -138.344824
128 -351.918708 213.941281
129 -141.080222 -351.918708
130 16.227133 -141.080222
131 -45.725101 16.227133
132 -360.546480 -45.725101
133 764.161393 -360.546480
134 90.709939 764.161393
135 -261.862098 90.709939
136 -129.292983 -261.862098
137 -107.138747 -129.292983
138 151.077433 -107.138747
139 -92.675294 151.077433
140 30.697983 -92.675294
141 270.231161 30.697983
142 112.741134 270.231161
143 288.680375 112.741134
144 -27.622285 288.680375
145 487.840241 -27.622285
146 -165.575180 487.840241
147 262.611281 -165.575180
148 -152.108774 262.611281
149 162.475414 -152.108774
150 -69.532414 162.475414
151 -115.485559 -69.532414
152 57.286287 -115.485559
153 90.826747 57.286287
154 -412.064729 90.826747
155 109.903283 -412.064729
156 -298.071764 109.903283
157 139.251821 -298.071764
158 80.331147 139.251821
159 -824.656621 80.331147
160 494.322507 -824.656621
161 81.007337 494.322507
162 -198.044004 81.007337
163 369.525317 -198.044004
164 160.162252 369.525317
165 -294.152512 160.162252
166 132.207285 -294.152512
167 -158.876296 132.207285
168 -393.780837 -158.876296
169 -323.033277 -393.780837
170 524.505146 -323.033277
171 -74.088092 524.505146
172 200.213622 -74.088092
173 -155.805720 200.213622
174 -54.205013 -155.805720
175 334.980077 -54.205013
176 390.203151 334.980077
177 -32.742274 390.203151
178 -230.976284 -32.742274
179 -194.784038 -230.976284
180 153.753554 -194.784038
181 288.450527 153.753554
182 -224.571829 288.450527
183 7.112123 -224.571829
184 -68.789291 7.112123
185 -224.559922 -68.789291
186 -231.341590 -224.559922
187 -257.834357 -231.341590
188 381.080967 -257.834357
189 -187.092585 381.080967
190 -152.381189 -187.092585
191 -21.102202 -152.381189
192 272.652158 -21.102202
193 532.825325 272.652158
194 37.503216 532.825325
195 122.369347 37.503216
196 115.522900 122.369347
197 -10.921466 115.522900
198 296.347629 -10.921466
199 4.084141 296.347629
200 -235.364180 4.084141
201 -176.912403 -235.364180
202 -112.398290 -176.912403
203 23.109297 -112.398290
204 512.214277 23.109297
205 -20.354553 512.214277
206 54.556863 -20.354553
207 -116.979023 54.556863
208 -130.964990 -116.979023
209 -64.393227 -130.964990
210 43.971447 -64.393227
211 6.060979 43.971447
212 4.410485 6.060979
213 -301.889568 4.410485
214 365.064215 -301.889568
215 116.395584 365.064215
216 251.323447 116.395584
217 -75.857973 251.323447
218 -140.417827 -75.857973
219 -225.787323 -140.417827
220 134.522968 -225.787323
221 465.405140 134.522968
222 824.213720 465.405140
223 -9.390097 824.213720
224 -146.604489 -9.390097
225 196.567357 -146.604489
226 -124.743211 196.567357
227 91.207614 -124.743211
228 207.222006 91.207614
229 615.058552 207.222006
230 112.691573 615.058552
231 162.300023 112.691573
232 436.666232 162.300023
233 89.501732 436.666232
234 120.778893 89.501732
235 -211.154108 120.778893
236 -177.758943 -211.154108
237 81.276502 -177.758943
238 208.726905 81.276502
239 46.784861 208.726905
240 135.168975 46.784861
241 -60.128467 135.168975
242 -129.199177 -60.128467
243 -98.512571 -129.199177
244 -6.016713 -98.512571
245 152.313518 -6.016713
246 386.076928 152.313518
247 -249.091301 386.076928
248 -61.758315 -249.091301
249 -32.615772 -61.758315
250 48.521828 -32.615772
251 -36.940067 48.521828
252 -277.357769 -36.940067
253 78.345833 -277.357769
254 117.375302 78.345833
255 -148.025948 117.375302
256 -40.428945 -148.025948
257 -178.921042 -40.428945
258 -238.738812 -178.921042
259 -213.803789 -238.738812
260 89.127987 -213.803789
261 -78.315029 89.127987
262 -119.066084 -78.315029
263 -17.111926 -119.066084
264 -56.299257 -17.111926
265 -104.298097 -56.299257
266 -21.254365 -104.298097
267 -246.086109 -21.254365
268 -457.437042 -246.086109
269 -31.931175 -457.437042
270 308.550540 -31.931175
271 -130.783218 308.550540
272 78.008463 -130.783218
273 -111.103767 78.008463
274 -125.689811 -111.103767
275 51.798920 -125.689811
276 -76.828391 51.798920
277 -211.343156 -76.828391
278 121.616243 -211.343156
279 -71.539322 121.616243
280 78.374695 -71.539322
281 -180.963695 78.374695
282 115.828853 -180.963695
283 -122.433868 115.828853
284 -108.581898 -122.433868
285 327.762566 -108.581898
286 87.232026 327.762566
287 63.886848 87.232026
288 158.664345 63.886848
> 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/7xz3u1354796196.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/8u7rv1354796196.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/9lmop1354796196.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/1001kk1354796196.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/11a4nd1354796196.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/12yvpa1354796196.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/1328py1354796196.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/14zt3q1354796196.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/15pnxy1354796196.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/164bza1354796196.tab")
+ }
>
> try(system("convert tmp/1fmhp1354796196.ps tmp/1fmhp1354796196.png",intern=TRUE))
character(0)
> try(system("convert tmp/2l1do1354796196.ps tmp/2l1do1354796196.png",intern=TRUE))
character(0)
> try(system("convert tmp/3nwn91354796196.ps tmp/3nwn91354796196.png",intern=TRUE))
character(0)
> try(system("convert tmp/4pbno1354796196.ps tmp/4pbno1354796196.png",intern=TRUE))
character(0)
> try(system("convert tmp/5k4p01354796196.ps tmp/5k4p01354796196.png",intern=TRUE))
character(0)
> try(system("convert tmp/6h2kq1354796196.ps tmp/6h2kq1354796196.png",intern=TRUE))
character(0)
> try(system("convert tmp/7xz3u1354796196.ps tmp/7xz3u1354796196.png",intern=TRUE))
character(0)
> try(system("convert tmp/8u7rv1354796196.ps tmp/8u7rv1354796196.png",intern=TRUE))
character(0)
> try(system("convert tmp/9lmop1354796196.ps tmp/9lmop1354796196.png",intern=TRUE))
character(0)
> try(system("convert tmp/1001kk1354796196.ps tmp/1001kk1354796196.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
14.047 1.131 15.177