R version 2.13.0 (2011-04-13)
Copyright (C) 2011 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i486-pc-linux-gnu (32-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> x <- array(list(1418
+ ,210907
+ ,56
+ ,396
+ ,81
+ ,3
+ ,79
+ ,30
+ ,115
+ ,94
+ ,112285
+ ,24188
+ ,146283
+ ,144
+ ,145
+ ,869
+ ,120982
+ ,56
+ ,297
+ ,55
+ ,4
+ ,58
+ ,28
+ ,109
+ ,103
+ ,84786
+ ,18273
+ ,98364
+ ,103
+ ,101
+ ,1530
+ ,176508
+ ,54
+ ,559
+ ,50
+ ,12
+ ,60
+ ,38
+ ,146
+ ,93
+ ,83123
+ ,14130
+ ,86146
+ ,98
+ ,98
+ ,2172
+ ,179321
+ ,89
+ ,967
+ ,125
+ ,2
+ ,108
+ ,30
+ ,116
+ ,103
+ ,101193
+ ,32287
+ ,96933
+ ,135
+ ,132
+ ,901
+ ,123185
+ ,40
+ ,270
+ ,40
+ ,1
+ ,49
+ ,22
+ ,68
+ ,51
+ ,38361
+ ,8654
+ ,79234
+ ,61
+ ,60
+ ,463
+ ,52746
+ ,25
+ ,143
+ ,37
+ ,3
+ ,0
+ ,26
+ ,101
+ ,70
+ ,68504
+ ,9245
+ ,42551
+ ,39
+ ,38
+ ,3201
+ ,385534
+ ,92
+ ,1562
+ ,63
+ ,0
+ ,121
+ ,25
+ ,96
+ ,91
+ ,119182
+ ,33251
+ ,195663
+ ,150
+ ,144
+ ,371
+ ,33170
+ ,18
+ ,109
+ ,44
+ ,0
+ ,1
+ ,18
+ ,67
+ ,22
+ ,22807
+ ,1271
+ ,6853
+ ,5
+ ,5
+ ,1192
+ ,101645
+ ,63
+ ,371
+ ,88
+ ,0
+ ,20
+ ,11
+ ,44
+ ,38
+ ,17140
+ ,5279
+ ,21529
+ ,28
+ ,28
+ ,1583
+ ,149061
+ ,44
+ ,656
+ ,66
+ ,5
+ ,43
+ ,26
+ ,100
+ ,93
+ ,116174
+ ,27101
+ ,95757
+ ,84
+ ,84
+ ,1439
+ ,165446
+ ,33
+ ,511
+ ,57
+ ,0
+ ,69
+ ,25
+ ,93
+ ,60
+ ,57635
+ ,16373
+ ,85584
+ ,80
+ ,79
+ ,1764
+ ,237213
+ ,84
+ ,655
+ ,74
+ ,0
+ ,78
+ ,38
+ ,140
+ ,123
+ ,66198
+ ,19716
+ ,143983
+ ,130
+ ,127
+ ,1495
+ ,173326
+ ,88
+ ,465
+ ,49
+ ,7
+ ,86
+ ,44
+ ,166
+ ,148
+ ,71701
+ ,17753
+ ,75851
+ ,82
+ ,78
+ ,1373
+ ,133131
+ ,55
+ ,525
+ ,52
+ ,7
+ ,44
+ ,30
+ ,99
+ ,90
+ ,57793
+ ,9028
+ ,59238
+ ,60
+ ,60
+ ,2187
+ ,258873
+ ,60
+ ,885
+ ,88
+ ,3
+ ,104
+ ,40
+ ,139
+ ,124
+ ,80444
+ ,18653
+ ,93163
+ ,131
+ ,131
+ ,1491
+ ,180083
+ ,66
+ ,497
+ ,36
+ ,9
+ ,63
+ ,34
+ ,130
+ ,70
+ ,53855
+ ,8828
+ ,96037
+ ,84
+ ,84
+ ,4041
+ ,324799
+ ,154
+ ,1436
+ ,108
+ ,0
+ ,158
+ ,47
+ ,181
+ ,168
+ ,97668
+ ,29498
+ ,151511
+ ,140
+ ,133
+ ,1706
+ ,230964
+ ,53
+ ,612
+ ,43
+ ,4
+ ,102
+ ,30
+ ,116
+ ,115
+ ,133824
+ ,27563
+ ,136368
+ ,151
+ ,150
+ ,2152
+ ,236785
+ ,119
+ ,865
+ ,75
+ ,3
+ ,77
+ ,31
+ ,116
+ ,71
+ ,101481
+ ,18293
+ ,112642
+ ,91
+ ,91
+ ,1036
+ ,135473
+ ,41
+ ,385
+ ,32
+ ,0
+ ,82
+ ,23
+ ,88
+ ,66
+ ,99645
+ ,22530
+ ,94728
+ ,138
+ ,132
+ ,1882
+ ,202925
+ ,61
+ ,567
+ ,44
+ ,7
+ ,115
+ ,36
+ ,139
+ ,134
+ ,114789
+ ,15977
+ ,105499
+ ,150
+ ,136
+ ,1929
+ ,215147
+ ,58
+ ,639
+ ,85
+ ,0
+ ,101
+ ,36
+ ,135
+ ,117
+ ,99052
+ ,35082
+ ,121527
+ ,124
+ ,124
+ ,2242
+ ,344297
+ ,75
+ ,963
+ ,86
+ ,1
+ ,80
+ ,30
+ ,108
+ ,108
+ ,67654
+ ,16116
+ ,127766
+ ,119
+ ,118
+ ,1220
+ ,153935
+ ,33
+ ,398
+ ,56
+ ,5
+ ,50
+ ,25
+ ,89
+ ,84
+ ,65553
+ ,15849
+ ,98958
+ ,73
+ ,70
+ ,1289
+ ,132943
+ ,40
+ ,410
+ ,50
+ ,7
+ ,83
+ ,39
+ ,156
+ ,156
+ ,97500
+ ,16026
+ ,77900
+ ,110
+ ,107
+ ,2515
+ ,174724
+ ,92
+ ,966
+ ,135
+ ,0
+ ,123
+ ,34
+ ,129
+ ,120
+ ,69112
+ ,26569
+ ,85646
+ ,123
+ ,119
+ ,2147
+ ,174415
+ ,100
+ ,801
+ ,63
+ ,0
+ ,73
+ ,31
+ ,118
+ ,114
+ ,82753
+ ,24785
+ ,98579
+ ,90
+ ,89
+ ,2352
+ ,225548
+ ,112
+ ,892
+ ,81
+ ,5
+ ,81
+ ,31
+ ,118
+ ,94
+ ,85323
+ ,17569
+ ,130767
+ ,116
+ ,112
+ ,1638
+ ,223632
+ ,73
+ ,513
+ ,52
+ ,0
+ ,105
+ ,33
+ ,125
+ ,120
+ ,72654
+ ,23825
+ ,131741
+ ,113
+ ,108
+ ,1222
+ ,124817
+ ,40
+ ,469
+ ,44
+ ,0
+ ,47
+ ,25
+ ,95
+ ,81
+ ,30727
+ ,7869
+ ,53907
+ ,56
+ ,52
+ ,1812
+ ,221698
+ ,45
+ ,683
+ ,113
+ ,0
+ ,105
+ ,33
+ ,126
+ ,110
+ ,77873
+ ,14975
+ ,178812
+ ,115
+ ,112
+ ,1677
+ ,210767
+ ,60
+ ,643
+ ,39
+ ,3
+ ,94
+ ,35
+ ,135
+ ,133
+ ,117478
+ ,37791
+ ,146761
+ ,119
+ ,116
+ ,1579
+ ,170266
+ ,62
+ ,535
+ ,73
+ ,4
+ ,44
+ ,42
+ ,154
+ ,122
+ ,74007
+ ,9605
+ ,82036
+ ,129
+ ,123
+ ,1731
+ ,260561
+ ,75
+ ,625
+ ,48
+ ,1
+ ,114
+ ,43
+ ,165
+ ,158
+ ,90183
+ ,27295
+ ,163253
+ ,127
+ ,125
+ ,807
+ ,84853
+ ,31
+ ,264
+ ,33
+ ,4
+ ,38
+ ,30
+ ,113
+ ,109
+ ,61542
+ ,2746
+ ,27032
+ ,27
+ ,27
+ ,2452
+ ,294424
+ ,77
+ ,992
+ ,59
+ ,2
+ ,107
+ ,33
+ ,127
+ ,124
+ ,101494
+ ,34461
+ ,171975
+ ,175
+ ,162
+ ,829
+ ,101011
+ ,34
+ ,238
+ ,41
+ ,0
+ ,30
+ ,13
+ ,52
+ ,39
+ ,27570
+ ,8098
+ ,65990
+ ,35
+ ,32
+ ,1940
+ ,215641
+ ,46
+ ,818
+ ,69
+ ,0
+ ,71
+ ,32
+ ,121
+ ,92
+ ,55813
+ ,4787
+ ,86572
+ ,64
+ ,64
+ ,2662
+ ,325107
+ ,99
+ ,937
+ ,64
+ ,0
+ ,84
+ ,36
+ ,136
+ ,126
+ ,79215
+ ,24919
+ ,159676
+ ,96
+ ,92
+ ,186
+ ,7176
+ ,17
+ ,70
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,1423
+ ,603
+ ,1929
+ ,0
+ ,0
+ ,1499
+ ,167542
+ ,66
+ ,507
+ ,59
+ ,2
+ ,59
+ ,28
+ ,108
+ ,70
+ ,55461
+ ,16329
+ ,85371
+ ,84
+ ,83
+ ,865
+ ,106408
+ ,30
+ ,260
+ ,32
+ ,1
+ ,33
+ ,14
+ ,46
+ ,37
+ ,31081
+ ,12558
+ ,58391
+ ,41
+ ,41
+ ,1793
+ ,96560
+ ,76
+ ,503
+ ,129
+ ,0
+ ,42
+ ,17
+ ,54
+ ,38
+ ,22996
+ ,7784
+ ,31580
+ ,47
+ ,47
+ ,2527
+ ,265769
+ ,146
+ ,927
+ ,37
+ ,2
+ ,96
+ ,32
+ ,124
+ ,120
+ ,83122
+ ,28522
+ ,136815
+ ,126
+ ,120
+ ,2747
+ ,269651
+ ,67
+ ,1269
+ ,31
+ ,10
+ ,106
+ ,30
+ ,115
+ ,93
+ ,70106
+ ,22265
+ ,120642
+ ,105
+ ,105
+ ,1324
+ ,149112
+ ,56
+ ,537
+ ,65
+ ,6
+ ,56
+ ,35
+ ,128
+ ,95
+ ,60578
+ ,14459
+ ,69107
+ ,80
+ ,79
+ ,2702
+ ,175824
+ ,107
+ ,910
+ ,107
+ ,0
+ ,57
+ ,20
+ ,80
+ ,77
+ ,39992
+ ,14526
+ ,50495
+ ,70
+ ,65
+ ,1383
+ ,152871
+ ,58
+ ,532
+ ,74
+ ,5
+ ,59
+ ,28
+ ,97
+ ,90
+ ,79892
+ ,22240
+ ,108016
+ ,73
+ ,70
+ ,1179
+ ,111665
+ ,34
+ ,345
+ ,54
+ ,4
+ ,39
+ ,28
+ ,104
+ ,80
+ ,49810
+ ,11802
+ ,46341
+ ,57
+ ,55
+ ,2099
+ ,116408
+ ,61
+ ,918
+ ,76
+ ,1
+ ,34
+ ,39
+ ,59
+ ,31
+ ,71570
+ ,7623
+ ,78348
+ ,40
+ ,39
+ ,4308
+ ,362301
+ ,119
+ ,1635
+ ,715
+ ,2
+ ,76
+ ,34
+ ,125
+ ,110
+ ,100708
+ ,11912
+ ,79336
+ ,68
+ ,67
+ ,918
+ ,78800
+ ,42
+ ,330
+ ,57
+ ,2
+ ,20
+ ,26
+ ,82
+ ,66
+ ,33032
+ ,7935
+ ,56968
+ ,21
+ ,21
+ ,1831
+ ,183167
+ ,66
+ ,557
+ ,66
+ ,0
+ ,91
+ ,39
+ ,149
+ ,138
+ ,82875
+ ,18220
+ ,93176
+ ,127
+ ,127
+ ,3373
+ ,277965
+ ,89
+ ,1178
+ ,106
+ ,8
+ ,115
+ ,39
+ ,149
+ ,133
+ ,139077
+ ,19199
+ ,161632
+ ,154
+ ,152
+ ,1713
+ ,150629
+ ,44
+ ,740
+ ,54
+ ,3
+ ,85
+ ,33
+ ,122
+ ,113
+ ,71595
+ ,19918
+ ,87850
+ ,116
+ ,113
+ ,1438
+ ,168809
+ ,66
+ ,452
+ ,32
+ ,0
+ ,76
+ ,28
+ ,118
+ ,100
+ ,72260
+ ,21884
+ ,127969
+ ,102
+ ,99
+ ,496
+ ,24188
+ ,24
+ ,218
+ ,20
+ ,0
+ ,8
+ ,4
+ ,12
+ ,7
+ ,5950
+ ,2694
+ ,15049
+ ,7
+ ,7
+ ,2253
+ ,329267
+ ,259
+ ,764
+ ,71
+ ,8
+ ,79
+ ,39
+ ,144
+ ,140
+ ,115762
+ ,15808
+ ,155135
+ ,148
+ ,141
+ ,744
+ ,65029
+ ,17
+ ,255
+ ,21
+ ,5
+ ,21
+ ,18
+ ,67
+ ,61
+ ,32551
+ ,3597
+ ,25109
+ ,21
+ ,21
+ ,1161
+ ,101097
+ ,64
+ ,454
+ ,70
+ ,3
+ ,30
+ ,14
+ ,52
+ ,41
+ ,31701
+ ,5296
+ ,45824
+ ,35
+ ,35
+ ,2352
+ ,218946
+ ,41
+ ,866
+ ,112
+ ,1
+ ,76
+ ,29
+ ,108
+ ,96
+ ,80670
+ ,25239
+ ,102996
+ ,112
+ ,109
+ ,2144
+ ,244052
+ ,68
+ ,574
+ ,66
+ ,5
+ ,101
+ ,44
+ ,166
+ ,164
+ ,143558
+ ,29801
+ ,160604
+ ,137
+ ,133
+ ,4691
+ ,341570
+ ,168
+ ,1276
+ ,190
+ ,1
+ ,94
+ ,21
+ ,80
+ ,78
+ ,117105
+ ,18450
+ ,158051
+ ,135
+ ,123
+ ,1112
+ ,103597
+ ,43
+ ,379
+ ,66
+ ,1
+ ,27
+ ,16
+ ,60
+ ,49
+ ,23789
+ ,7132
+ ,44547
+ ,26
+ ,26
+ ,2694
+ ,233328
+ ,132
+ ,825
+ ,165
+ ,5
+ ,92
+ ,28
+ ,107
+ ,102
+ ,120733
+ ,34861
+ ,162647
+ ,230
+ ,230
+ ,1973
+ ,256462
+ ,105
+ ,798
+ ,56
+ ,0
+ ,123
+ ,35
+ ,127
+ ,124
+ ,105195
+ ,35940
+ ,174141
+ ,181
+ ,166
+ ,1769
+ ,206161
+ ,71
+ ,663
+ ,61
+ ,12
+ ,75
+ ,28
+ ,107
+ ,99
+ ,73107
+ ,16688
+ ,60622
+ ,71
+ ,68
+ ,3148
+ ,311473
+ ,112
+ ,1069
+ ,53
+ ,8
+ ,128
+ ,38
+ ,146
+ ,129
+ ,132068
+ ,24683
+ ,179566
+ ,147
+ ,147
+ ,2474
+ ,235800
+ ,94
+ ,921
+ ,127
+ ,8
+ ,105
+ ,23
+ ,84
+ ,62
+ ,149193
+ ,46230
+ ,184301
+ ,190
+ ,179
+ ,2084
+ ,177939
+ ,82
+ ,858
+ ,63
+ ,8
+ ,55
+ ,36
+ ,141
+ ,73
+ ,46821
+ ,10387
+ ,75661
+ ,64
+ ,61
+ ,1954
+ ,207176
+ ,70
+ ,711
+ ,38
+ ,8
+ ,56
+ ,32
+ ,123
+ ,114
+ ,87011
+ ,21436
+ ,96144
+ ,105
+ ,101
+ ,1226
+ ,196553
+ ,57
+ ,503
+ ,50
+ ,2
+ ,41
+ ,29
+ ,111
+ ,99
+ ,95260
+ ,30546
+ ,129847
+ ,107
+ ,108
+ ,1389
+ ,174184
+ ,53
+ ,382
+ ,52
+ ,0
+ ,72
+ ,25
+ ,98
+ ,70
+ ,55183
+ ,19746
+ ,117286
+ ,94
+ ,90
+ ,1496
+ ,143246
+ ,103
+ ,464
+ ,42
+ ,5
+ ,67
+ ,27
+ ,105
+ ,104
+ ,106671
+ ,15977
+ ,71180
+ ,116
+ ,114
+ ,2269
+ ,187559
+ ,121
+ ,717
+ ,76
+ ,8
+ ,75
+ ,36
+ ,135
+ ,116
+ ,73511
+ ,22583
+ ,109377
+ ,106
+ ,103
+ ,1833
+ ,187681
+ ,62
+ ,690
+ ,67
+ ,2
+ ,114
+ ,28
+ ,107
+ ,91
+ ,92945
+ ,17274
+ ,85298
+ ,143
+ ,142
+ ,1268
+ ,119016
+ ,52
+ ,462
+ ,50
+ ,5
+ ,118
+ ,23
+ ,85
+ ,74
+ ,78664
+ ,16469
+ ,73631
+ ,81
+ ,79
+ ,1943
+ ,182192
+ ,52
+ ,657
+ ,53
+ ,12
+ ,77
+ ,40
+ ,155
+ ,138
+ ,70054
+ ,14251
+ ,86767
+ ,89
+ ,88
+ ,893
+ ,73566
+ ,32
+ ,385
+ ,39
+ ,6
+ ,22
+ ,23
+ ,88
+ ,67
+ ,22618
+ ,3007
+ ,23824
+ ,26
+ ,25
+ ,1762
+ ,194979
+ ,62
+ ,577
+ ,50
+ ,7
+ ,66
+ ,40
+ ,155
+ ,151
+ ,74011
+ ,16851
+ ,93487
+ ,84
+ ,83
+ ,1403
+ ,167488
+ ,45
+ ,619
+ ,77
+ ,2
+ ,69
+ ,28
+ ,104
+ ,72
+ ,83737
+ ,21113
+ ,82981
+ ,113
+ ,113
+ ,1425
+ ,143756
+ ,46
+ ,479
+ ,57
+ ,0
+ ,105
+ ,34
+ ,132
+ ,120
+ ,69094
+ ,17401
+ ,73815
+ ,120
+ ,118
+ ,1857
+ ,275541
+ ,63
+ ,817
+ ,73
+ ,4
+ ,116
+ ,33
+ ,127
+ ,115
+ ,93133
+ ,23958
+ ,94552
+ ,110
+ ,110
+ ,1840
+ ,243199
+ ,75
+ ,752
+ ,34
+ ,3
+ ,88
+ ,28
+ ,108
+ ,105
+ ,95536
+ ,23567
+ ,132190
+ ,134
+ ,129
+ ,1502
+ ,182999
+ ,88
+ ,430
+ ,39
+ ,6
+ ,73
+ ,34
+ ,129
+ ,104
+ ,225920
+ ,13065
+ ,128754
+ ,54
+ ,51
+ ,1441
+ ,135649
+ ,46
+ ,451
+ ,46
+ ,2
+ ,99
+ ,30
+ ,116
+ ,108
+ ,62133
+ ,15358
+ ,66363
+ ,96
+ ,93
+ ,1420
+ ,152299
+ ,53
+ ,537
+ ,63
+ ,0
+ ,62
+ ,33
+ ,122
+ ,98
+ ,61370
+ ,14587
+ ,67808
+ ,78
+ ,76
+ ,1416
+ ,120221
+ ,37
+ ,519
+ ,35
+ ,1
+ ,53
+ ,22
+ ,85
+ ,69
+ ,43836
+ ,12770
+ ,61724
+ ,51
+ ,49
+ ,2970
+ ,346485
+ ,90
+ ,1000
+ ,106
+ ,0
+ ,118
+ ,38
+ ,147
+ ,111
+ ,106117
+ ,24021
+ ,131722
+ ,121
+ ,118
+ ,1317
+ ,145790
+ ,63
+ ,637
+ ,43
+ ,5
+ ,30
+ ,26
+ ,99
+ ,99
+ ,38692
+ ,9648
+ ,68580
+ ,38
+ ,38
+ ,1644
+ ,193339
+ ,78
+ ,465
+ ,47
+ ,2
+ ,100
+ ,35
+ ,87
+ ,71
+ ,84651
+ ,20537
+ ,106175
+ ,145
+ ,141
+ ,870
+ ,80953
+ ,25
+ ,437
+ ,31
+ ,0
+ ,49
+ ,8
+ ,28
+ ,27
+ ,56622
+ ,7905
+ ,55792
+ ,59
+ ,58
+ ,1654
+ ,122774
+ ,45
+ ,711
+ ,162
+ ,0
+ ,24
+ ,24
+ ,90
+ ,69
+ ,15986
+ ,4527
+ ,25157
+ ,27
+ ,27
+ ,1054
+ ,130585
+ ,46
+ ,299
+ ,57
+ ,5
+ ,67
+ ,29
+ ,109
+ ,107
+ ,95364
+ ,30495
+ ,76669
+ ,91
+ ,91
+ ,937
+ ,112611
+ ,41
+ ,248
+ ,36
+ ,0
+ ,46
+ ,20
+ ,78
+ ,73
+ ,26706
+ ,7117
+ ,57283
+ ,48
+ ,48
+ ,3004
+ ,286468
+ ,144
+ ,1162
+ ,263
+ ,1
+ ,57
+ ,29
+ ,111
+ ,107
+ ,89691
+ ,17719
+ ,105805
+ ,68
+ ,63
+ ,2008
+ ,241066
+ ,82
+ ,714
+ ,78
+ ,0
+ ,75
+ ,45
+ ,158
+ ,93
+ ,67267
+ ,27056
+ ,129484
+ ,58
+ ,56
+ ,2547
+ ,148446
+ ,91
+ ,905
+ ,63
+ ,1
+ ,135
+ ,37
+ ,141
+ ,129
+ ,126846
+ ,33473
+ ,72413
+ ,150
+ ,144
+ ,1885
+ ,204713
+ ,71
+ ,649
+ ,54
+ ,1
+ ,68
+ ,33
+ ,122
+ ,69
+ ,41140
+ ,9758
+ ,87831
+ ,74
+ ,73
+ ,1626
+ ,182079
+ ,63
+ ,512
+ ,63
+ ,2
+ ,124
+ ,33
+ ,124
+ ,118
+ ,102860
+ ,21115
+ ,96971
+ ,181
+ ,168
+ ,1468
+ ,140344
+ ,53
+ ,472
+ ,77
+ ,6
+ ,33
+ ,25
+ ,93
+ ,73
+ ,51715
+ ,7236
+ ,71299
+ ,65
+ ,64
+ ,2445
+ ,220516
+ ,62
+ ,905
+ ,79
+ ,1
+ ,98
+ ,32
+ ,124
+ ,119
+ ,55801
+ ,13790
+ ,77494
+ ,97
+ ,97
+ ,1964
+ ,243060
+ ,63
+ ,786
+ ,110
+ ,4
+ ,58
+ ,29
+ ,112
+ ,104
+ ,111813
+ ,32902
+ ,120336
+ ,121
+ ,117
+ ,1381
+ ,162765
+ ,32
+ ,489
+ ,56
+ ,2
+ ,68
+ ,28
+ ,108
+ ,107
+ ,120293
+ ,25131
+ ,93913
+ ,99
+ ,100
+ ,1369
+ ,182613
+ ,39
+ ,479
+ ,56
+ ,3
+ ,81
+ ,28
+ ,99
+ ,99
+ ,138599
+ ,30910
+ ,136048
+ ,152
+ ,149
+ ,1659
+ ,232138
+ ,62
+ ,617
+ ,43
+ ,0
+ ,131
+ ,31
+ ,117
+ ,90
+ ,161647
+ ,35947
+ ,181248
+ ,188
+ ,187
+ ,2888
+ ,265318
+ ,117
+ ,925
+ ,111
+ ,10
+ ,110
+ ,52
+ ,199
+ ,197
+ ,115929
+ ,29848
+ ,146123
+ ,138
+ ,127
+ ,1290
+ ,85574
+ ,34
+ ,351
+ ,71
+ ,0
+ ,37
+ ,21
+ ,78
+ ,36
+ ,24266
+ ,6943
+ ,32036
+ ,40
+ ,37
+ ,2845
+ ,310839
+ ,92
+ ,1144
+ ,62
+ ,9
+ ,130
+ ,24
+ ,91
+ ,85
+ ,162901
+ ,42705
+ ,186646
+ ,254
+ ,245
+ ,1982
+ ,225060
+ ,93
+ ,669
+ ,56
+ ,7
+ ,93
+ ,41
+ ,158
+ ,139
+ ,109825
+ ,31808
+ ,102255
+ ,87
+ ,87
+ ,1904
+ ,232317
+ ,54
+ ,707
+ ,74
+ ,0
+ ,118
+ ,33
+ ,126
+ ,106
+ ,129838
+ ,26675
+ ,168237
+ ,178
+ ,177
+ ,1391
+ ,144966
+ ,144
+ ,458
+ ,60
+ ,0
+ ,39
+ ,32
+ ,122
+ ,50
+ ,37510
+ ,8435
+ ,64219
+ ,51
+ ,49
+ ,602
+ ,43287
+ ,14
+ ,214
+ ,43
+ ,4
+ ,13
+ ,19
+ ,71
+ ,64
+ ,43750
+ ,7409
+ ,19630
+ ,49
+ ,49
+ ,1743
+ ,155754
+ ,61
+ ,599
+ ,68
+ ,4
+ ,74
+ ,20
+ ,75
+ ,31
+ ,40652
+ ,14993
+ ,76825
+ ,73
+ ,73
+ ,1559
+ ,164709
+ ,109
+ ,572
+ ,53
+ ,0
+ ,81
+ ,31
+ ,115
+ ,63
+ ,87771
+ ,36867
+ ,115338
+ ,176
+ ,177
+ ,2014
+ ,201940
+ ,38
+ ,897
+ ,87
+ ,0
+ ,109
+ ,31
+ ,119
+ ,92
+ ,85872
+ ,33835
+ ,109427
+ ,94
+ ,94
+ ,2143
+ ,235454
+ ,73
+ ,819
+ ,46
+ ,0
+ ,151
+ ,32
+ ,124
+ ,106
+ ,89275
+ ,24164
+ ,118168
+ ,120
+ ,117
+ ,2146
+ ,220801
+ ,75
+ ,720
+ ,105
+ ,1
+ ,51
+ ,18
+ ,72
+ ,63
+ ,44418
+ ,12607
+ ,84845
+ ,66
+ ,60
+ ,874
+ ,99466
+ ,50
+ ,273
+ ,32
+ ,0
+ ,28
+ ,23
+ ,91
+ ,69
+ ,192565
+ ,22609
+ ,153197
+ ,56
+ ,55
+ ,1590
+ ,92661
+ ,61
+ ,508
+ ,133
+ ,1
+ ,40
+ ,17
+ ,45
+ ,41
+ ,35232
+ ,5892
+ ,29877
+ ,39
+ ,39
+ ,1590
+ ,133328
+ ,55
+ ,506
+ ,79
+ ,0
+ ,56
+ ,20
+ ,78
+ ,56
+ ,40909
+ ,17014
+ ,63506
+ ,66
+ ,64
+ ,1210
+ ,61361
+ ,77
+ ,451
+ ,51
+ ,0
+ ,27
+ ,12
+ ,39
+ ,25
+ ,13294
+ ,5394
+ ,22445
+ ,27
+ ,26
+ ,2072
+ ,125930
+ ,75
+ ,699
+ ,207
+ ,4
+ ,37
+ ,17
+ ,68
+ ,65
+ ,32387
+ ,9178
+ ,47695
+ ,65
+ ,64
+ ,1281
+ ,100750
+ ,72
+ ,407
+ ,67
+ ,0
+ ,83
+ ,30
+ ,119
+ ,93
+ ,140867
+ ,6440
+ ,68370
+ ,58
+ ,58
+ ,1401
+ ,224549
+ ,50
+ ,465
+ ,47
+ ,4
+ ,54
+ ,31
+ ,117
+ ,114
+ ,120662
+ ,21916
+ ,146304
+ ,98
+ ,95
+ ,834
+ ,82316
+ ,32
+ ,245
+ ,34
+ ,4
+ ,27
+ ,10
+ ,39
+ ,38
+ ,21233
+ ,4011
+ ,38233
+ ,25
+ ,25
+ ,1105
+ ,102010
+ ,53
+ ,370
+ ,66
+ ,3
+ ,28
+ ,13
+ ,50
+ ,44
+ ,44332
+ ,5818
+ ,42071
+ ,26
+ ,26
+ ,1272
+ ,101523
+ ,42
+ ,316
+ ,76
+ ,0
+ ,59
+ ,22
+ ,88
+ ,87
+ ,61056
+ ,18647
+ ,50517
+ ,77
+ ,76
+ ,1944
+ ,243511
+ ,71
+ ,603
+ ,65
+ ,0
+ ,133
+ ,42
+ ,155
+ ,110
+ ,101338
+ ,20556
+ ,103950
+ ,130
+ ,129
+ ,391
+ ,22938
+ ,10
+ ,154
+ ,9
+ ,0
+ ,12
+ ,1
+ ,0
+ ,0
+ ,1168
+ ,238
+ ,5841
+ ,11
+ ,11
+ ,761
+ ,41566
+ ,35
+ ,229
+ ,42
+ ,5
+ ,0
+ ,9
+ ,36
+ ,27
+ ,13497
+ ,70
+ ,2341
+ ,2
+ ,2
+ ,1605
+ ,152474
+ ,65
+ ,577
+ ,45
+ ,0
+ ,106
+ ,32
+ ,123
+ ,83
+ ,65567
+ ,22392
+ ,84396
+ ,101
+ ,101
+ ,530
+ ,61857
+ ,25
+ ,192
+ ,25
+ ,4
+ ,23
+ ,11
+ ,32
+ ,30
+ ,25162
+ ,3913
+ ,24610
+ ,31
+ ,28
+ ,1988
+ ,99923
+ ,66
+ ,617
+ ,115
+ ,0
+ ,44
+ ,25
+ ,99
+ ,80
+ ,32334
+ ,12237
+ ,35753
+ ,36
+ ,36
+ ,1386
+ ,132487
+ ,41
+ ,411
+ ,97
+ ,0
+ ,71
+ ,36
+ ,136
+ ,98
+ ,40735
+ ,8388
+ ,55515
+ ,120
+ ,89
+ ,2395
+ ,317394
+ ,86
+ ,975
+ ,53
+ ,1
+ ,116
+ ,31
+ ,117
+ ,82
+ ,91413
+ ,22120
+ ,209056
+ ,195
+ ,193
+ ,387
+ ,21054
+ ,16
+ ,146
+ ,2
+ ,0
+ ,4
+ ,0
+ ,0
+ ,0
+ ,855
+ ,338
+ ,6622
+ ,4
+ ,4
+ ,1742
+ ,209641
+ ,42
+ ,705
+ ,52
+ ,5
+ ,62
+ ,24
+ ,88
+ ,60
+ ,97068
+ ,11727
+ ,115814
+ ,89
+ ,84
+ ,620
+ ,22648
+ ,19
+ ,184
+ ,44
+ ,0
+ ,12
+ ,13
+ ,39
+ ,28
+ ,44339
+ ,3704
+ ,11609
+ ,24
+ ,23
+ ,449
+ ,31414
+ ,19
+ ,200
+ ,22
+ ,0
+ ,18
+ ,8
+ ,25
+ ,9
+ ,14116
+ ,3988
+ ,13155
+ ,39
+ ,39
+ ,800
+ ,46698
+ ,45
+ ,274
+ ,35
+ ,0
+ ,14
+ ,13
+ ,52
+ ,33
+ ,10288
+ ,3030
+ ,18274
+ ,14
+ ,14
+ ,1684
+ ,131698
+ ,65
+ ,502
+ ,74
+ ,0
+ ,60
+ ,19
+ ,75
+ ,59
+ ,65622
+ ,13520
+ ,72875
+ ,78
+ ,78
+ ,1050
+ ,91735
+ ,35
+ ,382
+ ,103
+ ,0
+ ,7
+ ,18
+ ,71
+ ,49
+ ,16563
+ ,1421
+ ,10112
+ ,15
+ ,14
+ ,2699
+ ,244749
+ ,95
+ ,964
+ ,144
+ ,2
+ ,98
+ ,33
+ ,124
+ ,115
+ ,76643
+ ,20923
+ ,142775
+ ,106
+ ,101
+ ,1606
+ ,184510
+ ,49
+ ,537
+ ,60
+ ,7
+ ,64
+ ,40
+ ,151
+ ,140
+ ,110681
+ ,20237
+ ,68847
+ ,83
+ ,82
+ ,1502
+ ,79863
+ ,37
+ ,438
+ ,134
+ ,1
+ ,29
+ ,22
+ ,71
+ ,49
+ ,29011
+ ,3219
+ ,17659
+ ,24
+ ,24
+ ,1204
+ ,128423
+ ,64
+ ,369
+ ,89
+ ,8
+ ,32
+ ,38
+ ,145
+ ,120
+ ,92696
+ ,3769
+ ,20112
+ ,37
+ ,36
+ ,1138
+ ,97839
+ ,38
+ ,417
+ ,42
+ ,2
+ ,25
+ ,24
+ ,87
+ ,66
+ ,94785
+ ,12252
+ ,61023
+ ,77
+ ,75
+ ,568
+ ,38214
+ ,34
+ ,276
+ ,52
+ ,0
+ ,16
+ ,8
+ ,27
+ ,21
+ ,8773
+ ,1888
+ ,13983
+ ,16
+ ,16
+ ,1459
+ ,151101
+ ,32
+ ,514
+ ,98
+ ,2
+ ,48
+ ,35
+ ,131
+ ,124
+ ,83209
+ ,14497
+ ,65176
+ ,56
+ ,55
+ ,2158
+ ,272458
+ ,65
+ ,822
+ ,99
+ ,0
+ ,100
+ ,43
+ ,162
+ ,152
+ ,93815
+ ,28864
+ ,132432
+ ,132
+ ,131
+ ,1111
+ ,172494
+ ,52
+ ,389
+ ,52
+ ,0
+ ,46
+ ,43
+ ,165
+ ,139
+ ,86687
+ ,21721
+ ,112494
+ ,144
+ ,131
+ ,1421
+ ,108043
+ ,62
+ ,466
+ ,29
+ ,1
+ ,45
+ ,14
+ ,54
+ ,38
+ ,34553
+ ,4821
+ ,45109
+ ,40
+ ,39
+ ,2833
+ ,328107
+ ,65
+ ,1255
+ ,125
+ ,3
+ ,129
+ ,41
+ ,159
+ ,144
+ ,105547
+ ,33644
+ ,170875
+ ,153
+ ,144
+ ,1955
+ ,250579
+ ,83
+ ,694
+ ,106
+ ,0
+ ,130
+ ,38
+ ,147
+ ,120
+ ,103487
+ ,15923
+ ,180759
+ ,143
+ ,139
+ ,2922
+ ,351067
+ ,95
+ ,1024
+ ,95
+ ,3
+ ,136
+ ,45
+ ,170
+ ,160
+ ,213688
+ ,42935
+ ,214921
+ ,220
+ ,211
+ ,1002
+ ,158015
+ ,29
+ ,400
+ ,40
+ ,0
+ ,59
+ ,31
+ ,119
+ ,114
+ ,71220
+ ,18864
+ ,100226
+ ,79
+ ,78
+ ,1060
+ ,98866
+ ,18
+ ,397
+ ,140
+ ,0
+ ,25
+ ,13
+ ,49
+ ,39
+ ,23517
+ ,4977
+ ,32043
+ ,50
+ ,50
+ ,956
+ ,85439
+ ,33
+ ,350
+ ,43
+ ,0
+ ,32
+ ,28
+ ,104
+ ,78
+ ,56926
+ ,7785
+ ,54454
+ ,39
+ ,39
+ ,2186
+ ,229242
+ ,247
+ ,719
+ ,128
+ ,4
+ ,63
+ ,31
+ ,120
+ ,119
+ ,91721
+ ,17939
+ ,78876
+ ,95
+ ,90
+ ,3604
+ ,351619
+ ,139
+ ,1277
+ ,142
+ ,4
+ ,95
+ ,40
+ ,150
+ ,141
+ ,115168
+ ,23436
+ ,170745
+ ,169
+ ,166
+ ,1035
+ ,84207
+ ,29
+ ,356
+ ,73
+ ,11
+ ,14
+ ,30
+ ,112
+ ,101
+ ,111194
+ ,325
+ ,6940
+ ,12
+ ,12
+ ,1417
+ ,120445
+ ,118
+ ,457
+ ,72
+ ,0
+ ,36
+ ,16
+ ,59
+ ,56
+ ,51009
+ ,13539
+ ,49025
+ ,63
+ ,57
+ ,3261
+ ,324598
+ ,110
+ ,1402
+ ,128
+ ,0
+ ,113
+ ,37
+ ,136
+ ,133
+ ,135777
+ ,34538
+ ,122037
+ ,134
+ ,133
+ ,1587
+ ,131069
+ ,67
+ ,600
+ ,61
+ ,4
+ ,47
+ ,30
+ ,107
+ ,83
+ ,51513
+ ,12198
+ ,53782
+ ,69
+ ,69
+ ,1424
+ ,204271
+ ,42
+ ,480
+ ,73
+ ,0
+ ,92
+ ,35
+ ,130
+ ,116
+ ,74163
+ ,26924
+ ,127748
+ ,119
+ ,119
+ ,1701
+ ,165543
+ ,65
+ ,595
+ ,148
+ ,1
+ ,70
+ ,32
+ ,115
+ ,90
+ ,51633
+ ,12716
+ ,86839
+ ,119
+ ,119
+ ,1249
+ ,141722
+ ,94
+ ,436
+ ,64
+ ,0
+ ,19
+ ,27
+ ,107
+ ,36
+ ,75345
+ ,8172
+ ,44830
+ ,75
+ ,65
+ ,946
+ ,116048
+ ,64
+ ,230
+ ,45
+ ,0
+ ,50
+ ,20
+ ,75
+ ,50
+ ,33416
+ ,10855
+ ,77395
+ ,63
+ ,61
+ ,1926
+ ,250047
+ ,81
+ ,651
+ ,58
+ ,0
+ ,41
+ ,18
+ ,71
+ ,61
+ ,83305
+ ,11932
+ ,89324
+ ,55
+ ,49
+ ,3352
+ ,299775
+ ,95
+ ,1367
+ ,97
+ ,9
+ ,91
+ ,31
+ ,120
+ ,97
+ ,98952
+ ,14300
+ ,103300
+ ,103
+ ,101
+ ,1641
+ ,195838
+ ,67
+ ,564
+ ,50
+ ,1
+ ,111
+ ,31
+ ,116
+ ,98
+ ,102372
+ ,25515
+ ,112283
+ ,197
+ ,196
+ ,2035
+ ,173260
+ ,63
+ ,716
+ ,37
+ ,3
+ ,41
+ ,21
+ ,79
+ ,78
+ ,37238
+ ,2805
+ ,10901
+ ,16
+ ,15
+ ,2312
+ ,254488
+ ,83
+ ,747
+ ,50
+ ,10
+ ,120
+ ,39
+ ,150
+ ,117
+ ,103772
+ ,29402
+ ,120691
+ ,140
+ ,136
+ ,1369
+ ,104389
+ ,45
+ ,467
+ ,105
+ ,5
+ ,135
+ ,41
+ ,156
+ ,148
+ ,123969
+ ,16440
+ ,58106
+ ,89
+ ,89
+ ,1577
+ ,136084
+ ,30
+ ,671
+ ,69
+ ,0
+ ,27
+ ,13
+ ,51
+ ,41
+ ,27142
+ ,11221
+ ,57140
+ ,40
+ ,40
+ ,2201
+ ,199476
+ ,70
+ ,861
+ ,46
+ ,2
+ ,87
+ ,32
+ ,118
+ ,105
+ ,135400
+ ,28732
+ ,122422
+ ,125
+ ,123
+ ,961
+ ,92499
+ ,32
+ ,319
+ ,57
+ ,0
+ ,25
+ ,18
+ ,71
+ ,55
+ ,21399
+ ,5250
+ ,25899
+ ,21
+ ,21
+ ,1900
+ ,224330
+ ,83
+ ,612
+ ,52
+ ,1
+ ,131
+ ,39
+ ,144
+ ,132
+ ,130115
+ ,28608
+ ,139296
+ ,167
+ ,163
+ ,1254
+ ,135781
+ ,31
+ ,433
+ ,98
+ ,2
+ ,45
+ ,14
+ ,47
+ ,44
+ ,24874
+ ,8092
+ ,52678
+ ,32
+ ,29
+ ,1335
+ ,74408
+ ,67
+ ,434
+ ,61
+ ,4
+ ,29
+ ,7
+ ,28
+ ,21
+ ,34988
+ ,4473
+ ,23853
+ ,36
+ ,35
+ ,1597
+ ,81240
+ ,66
+ ,503
+ ,89
+ ,0
+ ,58
+ ,17
+ ,68
+ ,50
+ ,45549
+ ,1572
+ ,17306
+ ,13
+ ,13
+ ,207
+ ,14688
+ ,10
+ ,85
+ ,0
+ ,0
+ ,4
+ ,0
+ ,0
+ ,0
+ ,6023
+ ,2065
+ ,7953
+ ,5
+ ,5
+ ,1645
+ ,181633
+ ,70
+ ,564
+ ,48
+ ,2
+ ,47
+ ,30
+ ,110
+ ,73
+ ,64466
+ ,14817
+ ,89455
+ ,96
+ ,96
+ ,2429
+ ,271856
+ ,103
+ ,824
+ ,91
+ ,1
+ ,109
+ ,37
+ ,147
+ ,86
+ ,54990
+ ,16714
+ ,147866
+ ,151
+ ,151
+ ,151
+ ,7199
+ ,5
+ ,74
+ ,0
+ ,0
+ ,7
+ ,0
+ ,0
+ ,0
+ ,1644
+ ,556
+ ,4245
+ ,6
+ ,6
+ ,474
+ ,46660
+ ,20
+ ,259
+ ,7
+ ,0
+ ,12
+ ,5
+ ,15
+ ,13
+ ,6179
+ ,2089
+ ,21509
+ ,13
+ ,13
+ ,141
+ ,17547
+ ,5
+ ,69
+ ,3
+ ,0
+ ,0
+ ,1
+ ,4
+ ,4
+ ,3926
+ ,2658
+ ,7670
+ ,3
+ ,3
+ ,1639
+ ,133368
+ ,36
+ ,535
+ ,54
+ ,1
+ ,37
+ ,16
+ ,64
+ ,57
+ ,32755
+ ,10695
+ ,66675
+ ,57
+ ,56
+ ,872
+ ,95227
+ ,34
+ ,239
+ ,70
+ ,0
+ ,37
+ ,32
+ ,111
+ ,48
+ ,34777
+ ,1669
+ ,14336
+ ,23
+ ,23
+ ,1318
+ ,152601
+ ,48
+ ,438
+ ,36
+ ,2
+ ,46
+ ,24
+ ,85
+ ,46
+ ,73224
+ ,16267
+ ,53608
+ ,61
+ ,57
+ ,1018
+ ,98146
+ ,40
+ ,459
+ ,37
+ ,0
+ ,15
+ ,17
+ ,68
+ ,48
+ ,27114
+ ,7768
+ ,30059
+ ,21
+ ,14
+ ,1383
+ ,79619
+ ,43
+ ,426
+ ,123
+ ,3
+ ,42
+ ,11
+ ,40
+ ,32
+ ,20760
+ ,7252
+ ,29668
+ ,43
+ ,43
+ ,1314
+ ,59194
+ ,31
+ ,288
+ ,247
+ ,6
+ ,7
+ ,24
+ ,80
+ ,68
+ ,37636
+ ,6387
+ ,22097
+ ,20
+ ,20
+ ,1335
+ ,139942
+ ,42
+ ,498
+ ,46
+ ,0
+ ,54
+ ,22
+ ,88
+ ,87
+ ,65461
+ ,18715
+ ,96841
+ ,82
+ ,72
+ ,1403
+ ,118612
+ ,46
+ ,454
+ ,72
+ ,2
+ ,54
+ ,12
+ ,48
+ ,43
+ ,30080
+ ,7936
+ ,41907
+ ,90
+ ,87
+ ,910
+ ,72880
+ ,33
+ ,376
+ ,41
+ ,0
+ ,14
+ ,19
+ ,76
+ ,67
+ ,24094
+ ,8643
+ ,27080
+ ,25
+ ,21)
+ ,dim=c(15
+ ,197)
+ ,dimnames=list(c('pageviews'
+ ,'time_in_rfc'
+ ,'logins'
+ ,'compendium_views_info'
+ ,'compendium_views_pr'
+ ,'shared_compendiums'
+ ,'blogged_computations'
+ ,'compendiums_reviewed'
+ ,'feedback_messages_p1'
+ ,'feedback_messages_p120'
+ ,'totsize'
+ ,'totrevisions'
+ ,'totseconds'
+ ,'tothyperlinks'
+ ,'totblogs
')
+ ,1:197))
> y <- array(NA,dim=c(15,197),dimnames=list(c('pageviews','time_in_rfc','logins','compendium_views_info','compendium_views_pr','shared_compendiums','blogged_computations','compendiums_reviewed','feedback_messages_p1','feedback_messages_p120','totsize','totrevisions','totseconds','tothyperlinks','totblogs
'),1:197))
> 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 = '2'
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Prof. Dr. P. Wessa
> #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #Technical description: Write here your technical program description (don't use hard returns!)
> library(lattice)
> library(lmtest)
Loading required package: zoo
> n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test
> par1 <- as.numeric(par1)
> x <- t(y)
> k <- length(x[1,])
> n <- length(x[,1])
> x1 <- cbind(x[,par1], x[,1:k!=par1])
> mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1])
> colnames(x1) <- mycolnames #colnames(x)[par1]
> x <- x1
> if (par3 == 'First Differences'){
+ x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep='')))
+ for (i in 1:n-1) {
+ for (j in 1:k) {
+ x2[i,j] <- x[i+1,j] - x[i,j]
+ }
+ }
+ x <- x2
+ }
> if (par2 == 'Include Monthly Dummies'){
+ x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep ='')))
+ for (i in 1:11){
+ x2[seq(i,n,12),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> if (par2 == 'Include Quarterly Dummies'){
+ x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep ='')))
+ for (i in 1:3){
+ x2[seq(i,n,4),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> k <- length(x[1,])
> if (par3 == 'Linear Trend'){
+ x <- cbind(x, c(1:n))
+ colnames(x)[k+1] <- 't'
+ }
> x
time_in_rfc pageviews logins compendium_views_info compendium_views_pr
1 210907 1418 56 396 81
2 120982 869 56 297 55
3 176508 1530 54 559 50
4 179321 2172 89 967 125
5 123185 901 40 270 40
6 52746 463 25 143 37
7 385534 3201 92 1562 63
8 33170 371 18 109 44
9 101645 1192 63 371 88
10 149061 1583 44 656 66
11 165446 1439 33 511 57
12 237213 1764 84 655 74
13 173326 1495 88 465 49
14 133131 1373 55 525 52
15 258873 2187 60 885 88
16 180083 1491 66 497 36
17 324799 4041 154 1436 108
18 230964 1706 53 612 43
19 236785 2152 119 865 75
20 135473 1036 41 385 32
21 202925 1882 61 567 44
22 215147 1929 58 639 85
23 344297 2242 75 963 86
24 153935 1220 33 398 56
25 132943 1289 40 410 50
26 174724 2515 92 966 135
27 174415 2147 100 801 63
28 225548 2352 112 892 81
29 223632 1638 73 513 52
30 124817 1222 40 469 44
31 221698 1812 45 683 113
32 210767 1677 60 643 39
33 170266 1579 62 535 73
34 260561 1731 75 625 48
35 84853 807 31 264 33
36 294424 2452 77 992 59
37 101011 829 34 238 41
38 215641 1940 46 818 69
39 325107 2662 99 937 64
40 7176 186 17 70 1
41 167542 1499 66 507 59
42 106408 865 30 260 32
43 96560 1793 76 503 129
44 265769 2527 146 927 37
45 269651 2747 67 1269 31
46 149112 1324 56 537 65
47 175824 2702 107 910 107
48 152871 1383 58 532 74
49 111665 1179 34 345 54
50 116408 2099 61 918 76
51 362301 4308 119 1635 715
52 78800 918 42 330 57
53 183167 1831 66 557 66
54 277965 3373 89 1178 106
55 150629 1713 44 740 54
56 168809 1438 66 452 32
57 24188 496 24 218 20
58 329267 2253 259 764 71
59 65029 744 17 255 21
60 101097 1161 64 454 70
61 218946 2352 41 866 112
62 244052 2144 68 574 66
63 341570 4691 168 1276 190
64 103597 1112 43 379 66
65 233328 2694 132 825 165
66 256462 1973 105 798 56
67 206161 1769 71 663 61
68 311473 3148 112 1069 53
69 235800 2474 94 921 127
70 177939 2084 82 858 63
71 207176 1954 70 711 38
72 196553 1226 57 503 50
73 174184 1389 53 382 52
74 143246 1496 103 464 42
75 187559 2269 121 717 76
76 187681 1833 62 690 67
77 119016 1268 52 462 50
78 182192 1943 52 657 53
79 73566 893 32 385 39
80 194979 1762 62 577 50
81 167488 1403 45 619 77
82 143756 1425 46 479 57
83 275541 1857 63 817 73
84 243199 1840 75 752 34
85 182999 1502 88 430 39
86 135649 1441 46 451 46
87 152299 1420 53 537 63
88 120221 1416 37 519 35
89 346485 2970 90 1000 106
90 145790 1317 63 637 43
91 193339 1644 78 465 47
92 80953 870 25 437 31
93 122774 1654 45 711 162
94 130585 1054 46 299 57
95 112611 937 41 248 36
96 286468 3004 144 1162 263
97 241066 2008 82 714 78
98 148446 2547 91 905 63
99 204713 1885 71 649 54
100 182079 1626 63 512 63
101 140344 1468 53 472 77
102 220516 2445 62 905 79
103 243060 1964 63 786 110
104 162765 1381 32 489 56
105 182613 1369 39 479 56
106 232138 1659 62 617 43
107 265318 2888 117 925 111
108 85574 1290 34 351 71
109 310839 2845 92 1144 62
110 225060 1982 93 669 56
111 232317 1904 54 707 74
112 144966 1391 144 458 60
113 43287 602 14 214 43
114 155754 1743 61 599 68
115 164709 1559 109 572 53
116 201940 2014 38 897 87
117 235454 2143 73 819 46
118 220801 2146 75 720 105
119 99466 874 50 273 32
120 92661 1590 61 508 133
121 133328 1590 55 506 79
122 61361 1210 77 451 51
123 125930 2072 75 699 207
124 100750 1281 72 407 67
125 224549 1401 50 465 47
126 82316 834 32 245 34
127 102010 1105 53 370 66
128 101523 1272 42 316 76
129 243511 1944 71 603 65
130 22938 391 10 154 9
131 41566 761 35 229 42
132 152474 1605 65 577 45
133 61857 530 25 192 25
134 99923 1988 66 617 115
135 132487 1386 41 411 97
136 317394 2395 86 975 53
137 21054 387 16 146 2
138 209641 1742 42 705 52
139 22648 620 19 184 44
140 31414 449 19 200 22
141 46698 800 45 274 35
142 131698 1684 65 502 74
143 91735 1050 35 382 103
144 244749 2699 95 964 144
145 184510 1606 49 537 60
146 79863 1502 37 438 134
147 128423 1204 64 369 89
148 97839 1138 38 417 42
149 38214 568 34 276 52
150 151101 1459 32 514 98
151 272458 2158 65 822 99
152 172494 1111 52 389 52
153 108043 1421 62 466 29
154 328107 2833 65 1255 125
155 250579 1955 83 694 106
156 351067 2922 95 1024 95
157 158015 1002 29 400 40
158 98866 1060 18 397 140
159 85439 956 33 350 43
160 229242 2186 247 719 128
161 351619 3604 139 1277 142
162 84207 1035 29 356 73
163 120445 1417 118 457 72
164 324598 3261 110 1402 128
165 131069 1587 67 600 61
166 204271 1424 42 480 73
167 165543 1701 65 595 148
168 141722 1249 94 436 64
169 116048 946 64 230 45
170 250047 1926 81 651 58
171 299775 3352 95 1367 97
172 195838 1641 67 564 50
173 173260 2035 63 716 37
174 254488 2312 83 747 50
175 104389 1369 45 467 105
176 136084 1577 30 671 69
177 199476 2201 70 861 46
178 92499 961 32 319 57
179 224330 1900 83 612 52
180 135781 1254 31 433 98
181 74408 1335 67 434 61
182 81240 1597 66 503 89
183 14688 207 10 85 0
184 181633 1645 70 564 48
185 271856 2429 103 824 91
186 7199 151 5 74 0
187 46660 474 20 259 7
188 17547 141 5 69 3
189 133368 1639 36 535 54
190 95227 872 34 239 70
191 152601 1318 48 438 36
192 98146 1018 40 459 37
193 79619 1383 43 426 123
194 59194 1314 31 288 247
195 139942 1335 42 498 46
196 118612 1403 46 454 72
197 72880 910 33 376 41
shared_compendiums blogged_computations compendiums_reviewed
1 3 79 30
2 4 58 28
3 12 60 38
4 2 108 30
5 1 49 22
6 3 0 26
7 0 121 25
8 0 1 18
9 0 20 11
10 5 43 26
11 0 69 25
12 0 78 38
13 7 86 44
14 7 44 30
15 3 104 40
16 9 63 34
17 0 158 47
18 4 102 30
19 3 77 31
20 0 82 23
21 7 115 36
22 0 101 36
23 1 80 30
24 5 50 25
25 7 83 39
26 0 123 34
27 0 73 31
28 5 81 31
29 0 105 33
30 0 47 25
31 0 105 33
32 3 94 35
33 4 44 42
34 1 114 43
35 4 38 30
36 2 107 33
37 0 30 13
38 0 71 32
39 0 84 36
40 0 0 0
41 2 59 28
42 1 33 14
43 0 42 17
44 2 96 32
45 10 106 30
46 6 56 35
47 0 57 20
48 5 59 28
49 4 39 28
50 1 34 39
51 2 76 34
52 2 20 26
53 0 91 39
54 8 115 39
55 3 85 33
56 0 76 28
57 0 8 4
58 8 79 39
59 5 21 18
60 3 30 14
61 1 76 29
62 5 101 44
63 1 94 21
64 1 27 16
65 5 92 28
66 0 123 35
67 12 75 28
68 8 128 38
69 8 105 23
70 8 55 36
71 8 56 32
72 2 41 29
73 0 72 25
74 5 67 27
75 8 75 36
76 2 114 28
77 5 118 23
78 12 77 40
79 6 22 23
80 7 66 40
81 2 69 28
82 0 105 34
83 4 116 33
84 3 88 28
85 6 73 34
86 2 99 30
87 0 62 33
88 1 53 22
89 0 118 38
90 5 30 26
91 2 100 35
92 0 49 8
93 0 24 24
94 5 67 29
95 0 46 20
96 1 57 29
97 0 75 45
98 1 135 37
99 1 68 33
100 2 124 33
101 6 33 25
102 1 98 32
103 4 58 29
104 2 68 28
105 3 81 28
106 0 131 31
107 10 110 52
108 0 37 21
109 9 130 24
110 7 93 41
111 0 118 33
112 0 39 32
113 4 13 19
114 4 74 20
115 0 81 31
116 0 109 31
117 0 151 32
118 1 51 18
119 0 28 23
120 1 40 17
121 0 56 20
122 0 27 12
123 4 37 17
124 0 83 30
125 4 54 31
126 4 27 10
127 3 28 13
128 0 59 22
129 0 133 42
130 0 12 1
131 5 0 9
132 0 106 32
133 4 23 11
134 0 44 25
135 0 71 36
136 1 116 31
137 0 4 0
138 5 62 24
139 0 12 13
140 0 18 8
141 0 14 13
142 0 60 19
143 0 7 18
144 2 98 33
145 7 64 40
146 1 29 22
147 8 32 38
148 2 25 24
149 0 16 8
150 2 48 35
151 0 100 43
152 0 46 43
153 1 45 14
154 3 129 41
155 0 130 38
156 3 136 45
157 0 59 31
158 0 25 13
159 0 32 28
160 4 63 31
161 4 95 40
162 11 14 30
163 0 36 16
164 0 113 37
165 4 47 30
166 0 92 35
167 1 70 32
168 0 19 27
169 0 50 20
170 0 41 18
171 9 91 31
172 1 111 31
173 3 41 21
174 10 120 39
175 5 135 41
176 0 27 13
177 2 87 32
178 0 25 18
179 1 131 39
180 2 45 14
181 4 29 7
182 0 58 17
183 0 4 0
184 2 47 30
185 1 109 37
186 0 7 0
187 0 12 5
188 0 0 1
189 1 37 16
190 0 37 32
191 2 46 24
192 0 15 17
193 3 42 11
194 6 7 24
195 0 54 22
196 2 54 12
197 0 14 19
feedback_messages_p1 feedback_messages_p120 totsize totrevisions totseconds
1 115 94 112285 24188 146283
2 109 103 84786 18273 98364
3 146 93 83123 14130 86146
4 116 103 101193 32287 96933
5 68 51 38361 8654 79234
6 101 70 68504 9245 42551
7 96 91 119182 33251 195663
8 67 22 22807 1271 6853
9 44 38 17140 5279 21529
10 100 93 116174 27101 95757
11 93 60 57635 16373 85584
12 140 123 66198 19716 143983
13 166 148 71701 17753 75851
14 99 90 57793 9028 59238
15 139 124 80444 18653 93163
16 130 70 53855 8828 96037
17 181 168 97668 29498 151511
18 116 115 133824 27563 136368
19 116 71 101481 18293 112642
20 88 66 99645 22530 94728
21 139 134 114789 15977 105499
22 135 117 99052 35082 121527
23 108 108 67654 16116 127766
24 89 84 65553 15849 98958
25 156 156 97500 16026 77900
26 129 120 69112 26569 85646
27 118 114 82753 24785 98579
28 118 94 85323 17569 130767
29 125 120 72654 23825 131741
30 95 81 30727 7869 53907
31 126 110 77873 14975 178812
32 135 133 117478 37791 146761
33 154 122 74007 9605 82036
34 165 158 90183 27295 163253
35 113 109 61542 2746 27032
36 127 124 101494 34461 171975
37 52 39 27570 8098 65990
38 121 92 55813 4787 86572
39 136 126 79215 24919 159676
40 0 0 1423 603 1929
41 108 70 55461 16329 85371
42 46 37 31081 12558 58391
43 54 38 22996 7784 31580
44 124 120 83122 28522 136815
45 115 93 70106 22265 120642
46 128 95 60578 14459 69107
47 80 77 39992 14526 50495
48 97 90 79892 22240 108016
49 104 80 49810 11802 46341
50 59 31 71570 7623 78348
51 125 110 100708 11912 79336
52 82 66 33032 7935 56968
53 149 138 82875 18220 93176
54 149 133 139077 19199 161632
55 122 113 71595 19918 87850
56 118 100 72260 21884 127969
57 12 7 5950 2694 15049
58 144 140 115762 15808 155135
59 67 61 32551 3597 25109
60 52 41 31701 5296 45824
61 108 96 80670 25239 102996
62 166 164 143558 29801 160604
63 80 78 117105 18450 158051
64 60 49 23789 7132 44547
65 107 102 120733 34861 162647
66 127 124 105195 35940 174141
67 107 99 73107 16688 60622
68 146 129 132068 24683 179566
69 84 62 149193 46230 184301
70 141 73 46821 10387 75661
71 123 114 87011 21436 96144
72 111 99 95260 30546 129847
73 98 70 55183 19746 117286
74 105 104 106671 15977 71180
75 135 116 73511 22583 109377
76 107 91 92945 17274 85298
77 85 74 78664 16469 73631
78 155 138 70054 14251 86767
79 88 67 22618 3007 23824
80 155 151 74011 16851 93487
81 104 72 83737 21113 82981
82 132 120 69094 17401 73815
83 127 115 93133 23958 94552
84 108 105 95536 23567 132190
85 129 104 225920 13065 128754
86 116 108 62133 15358 66363
87 122 98 61370 14587 67808
88 85 69 43836 12770 61724
89 147 111 106117 24021 131722
90 99 99 38692 9648 68580
91 87 71 84651 20537 106175
92 28 27 56622 7905 55792
93 90 69 15986 4527 25157
94 109 107 95364 30495 76669
95 78 73 26706 7117 57283
96 111 107 89691 17719 105805
97 158 93 67267 27056 129484
98 141 129 126846 33473 72413
99 122 69 41140 9758 87831
100 124 118 102860 21115 96971
101 93 73 51715 7236 71299
102 124 119 55801 13790 77494
103 112 104 111813 32902 120336
104 108 107 120293 25131 93913
105 99 99 138599 30910 136048
106 117 90 161647 35947 181248
107 199 197 115929 29848 146123
108 78 36 24266 6943 32036
109 91 85 162901 42705 186646
110 158 139 109825 31808 102255
111 126 106 129838 26675 168237
112 122 50 37510 8435 64219
113 71 64 43750 7409 19630
114 75 31 40652 14993 76825
115 115 63 87771 36867 115338
116 119 92 85872 33835 109427
117 124 106 89275 24164 118168
118 72 63 44418 12607 84845
119 91 69 192565 22609 153197
120 45 41 35232 5892 29877
121 78 56 40909 17014 63506
122 39 25 13294 5394 22445
123 68 65 32387 9178 47695
124 119 93 140867 6440 68370
125 117 114 120662 21916 146304
126 39 38 21233 4011 38233
127 50 44 44332 5818 42071
128 88 87 61056 18647 50517
129 155 110 101338 20556 103950
130 0 0 1168 238 5841
131 36 27 13497 70 2341
132 123 83 65567 22392 84396
133 32 30 25162 3913 24610
134 99 80 32334 12237 35753
135 136 98 40735 8388 55515
136 117 82 91413 22120 209056
137 0 0 855 338 6622
138 88 60 97068 11727 115814
139 39 28 44339 3704 11609
140 25 9 14116 3988 13155
141 52 33 10288 3030 18274
142 75 59 65622 13520 72875
143 71 49 16563 1421 10112
144 124 115 76643 20923 142775
145 151 140 110681 20237 68847
146 71 49 29011 3219 17659
147 145 120 92696 3769 20112
148 87 66 94785 12252 61023
149 27 21 8773 1888 13983
150 131 124 83209 14497 65176
151 162 152 93815 28864 132432
152 165 139 86687 21721 112494
153 54 38 34553 4821 45109
154 159 144 105547 33644 170875
155 147 120 103487 15923 180759
156 170 160 213688 42935 214921
157 119 114 71220 18864 100226
158 49 39 23517 4977 32043
159 104 78 56926 7785 54454
160 120 119 91721 17939 78876
161 150 141 115168 23436 170745
162 112 101 111194 325 6940
163 59 56 51009 13539 49025
164 136 133 135777 34538 122037
165 107 83 51513 12198 53782
166 130 116 74163 26924 127748
167 115 90 51633 12716 86839
168 107 36 75345 8172 44830
169 75 50 33416 10855 77395
170 71 61 83305 11932 89324
171 120 97 98952 14300 103300
172 116 98 102372 25515 112283
173 79 78 37238 2805 10901
174 150 117 103772 29402 120691
175 156 148 123969 16440 58106
176 51 41 27142 11221 57140
177 118 105 135400 28732 122422
178 71 55 21399 5250 25899
179 144 132 130115 28608 139296
180 47 44 24874 8092 52678
181 28 21 34988 4473 23853
182 68 50 45549 1572 17306
183 0 0 6023 2065 7953
184 110 73 64466 14817 89455
185 147 86 54990 16714 147866
186 0 0 1644 556 4245
187 15 13 6179 2089 21509
188 4 4 3926 2658 7670
189 64 57 32755 10695 66675
190 111 48 34777 1669 14336
191 85 46 73224 16267 53608
192 68 48 27114 7768 30059
193 40 32 20760 7252 29668
194 80 68 37636 6387 22097
195 88 87 65461 18715 96841
196 48 43 30080 7936 41907
197 76 67 24094 8643 27080
tothyperlinks totblogs\r
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
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) pageviews logins
-7939.4657 12.1105 137.6262
compendium_views_info compendium_views_pr shared_compendiums
102.3639 16.3926 36.1221
blogged_computations compendiums_reviewed feedback_messages_p1
-6.1631 -1191.3912 514.8623
feedback_messages_p120 totsize totrevisions
126.7349 -0.1172 -0.7415
totseconds tothyperlinks `totblogs\r`
0.8515 -3.8620 38.2951
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-59568 -13428 737 10629 88702
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -7.939e+03 5.137e+03 -1.545 0.1240
pageviews 1.211e+01 9.818e+00 1.234 0.2190
logins 1.376e+02 7.258e+01 1.896 0.0595 .
compendium_views_info 1.024e+02 2.104e+01 4.865 2.47e-06 ***
compendium_views_pr 1.639e+01 3.849e+01 0.426 0.6707
shared_compendiums 3.612e+01 6.679e+02 0.054 0.9569
blogged_computations -6.163e+00 1.103e+02 -0.056 0.9555
compendiums_reviewed -1.191e+03 8.075e+02 -1.475 0.1418
feedback_messages_p1 5.149e+02 2.456e+02 2.096 0.0375 *
feedback_messages_p120 1.267e+02 1.278e+02 0.992 0.3225
totsize -1.172e-01 7.762e-02 -1.510 0.1329
totrevisions -7.415e-01 3.720e-01 -1.993 0.0477 *
totseconds 8.515e-01 8.384e-02 10.157 < 2e-16 ***
tothyperlinks -3.862e+00 5.063e+02 -0.008 0.9939
`totblogs\r` 3.830e+01 5.244e+02 0.073 0.9419
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 23950 on 182 degrees of freedom
Multiple R-squared: 0.9215, Adjusted R-squared: 0.9155
F-statistic: 152.6 on 14 and 182 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.6396043 0.72079133 0.360395666
[2,] 0.4962176 0.99243526 0.503782368
[3,] 0.4634835 0.92696694 0.536516530
[4,] 0.4279532 0.85590641 0.572046794
[5,] 0.4985760 0.99715204 0.501423982
[6,] 0.8309480 0.33810401 0.169052004
[7,] 0.7650753 0.46984941 0.234924705
[8,] 0.8364381 0.32712380 0.163561899
[9,] 0.8671995 0.26560104 0.132800518
[10,] 0.8703069 0.25938623 0.129693116
[11,] 0.8688771 0.26224573 0.131122864
[12,] 0.8513319 0.29733624 0.148668118
[13,] 0.8153572 0.36928565 0.184642823
[14,] 0.8984873 0.20302536 0.101512679
[15,] 0.8673095 0.26538093 0.132690463
[16,] 0.8304574 0.33908525 0.169542627
[17,] 0.7869548 0.42609030 0.213045151
[18,] 0.7352704 0.52945913 0.264729565
[19,] 0.6894657 0.62106854 0.310534272
[20,] 0.6302082 0.73958366 0.369791832
[21,] 0.5698416 0.86031671 0.430158353
[22,] 0.6495865 0.70082696 0.350413480
[23,] 0.6808025 0.63839510 0.319197548
[24,] 0.6288684 0.74226320 0.371131600
[25,] 0.5835048 0.83299042 0.416495212
[26,] 0.5300534 0.93989316 0.469946580
[27,] 0.4767479 0.95349588 0.523252060
[28,] 0.5125580 0.97488407 0.487442036
[29,] 0.4598070 0.91961407 0.540192966
[30,] 0.4122923 0.82458457 0.587707715
[31,] 0.3626077 0.72521538 0.637392308
[32,] 0.3175698 0.63513969 0.682430155
[33,] 0.5319544 0.93609122 0.468045612
[34,] 0.8392447 0.32151064 0.160755321
[35,] 0.8366084 0.32678324 0.163391619
[36,] 0.8156620 0.36867591 0.184337954
[37,] 0.9107623 0.17847548 0.089237740
[38,] 0.9417017 0.11659658 0.058298290
[39,] 0.9425742 0.11485154 0.057425770
[40,] 0.9355042 0.12899154 0.064495768
[41,] 0.9435060 0.11298802 0.056494008
[42,] 0.9285740 0.14285195 0.071425977
[43,] 0.9166462 0.16670753 0.083353766
[44,] 0.9003719 0.19925628 0.099628141
[45,] 0.8813819 0.23723618 0.118618089
[46,] 0.8567825 0.28643497 0.143217485
[47,] 0.8284745 0.34305097 0.171525484
[48,] 0.8921308 0.21573849 0.107869247
[49,] 0.8732195 0.25356095 0.126780474
[50,] 0.9354290 0.12914194 0.064570971
[51,] 0.9263860 0.14722797 0.073613983
[52,] 0.9253814 0.14923719 0.074618597
[53,] 0.9304088 0.13918237 0.069591186
[54,] 0.9175315 0.16493705 0.082468527
[55,] 0.9049489 0.19010222 0.095051110
[56,] 0.8878940 0.22421203 0.112106015
[57,] 0.8648295 0.27034103 0.135170513
[58,] 0.8942003 0.21159946 0.105799729
[59,] 0.8773731 0.24525371 0.122626854
[60,] 0.8579668 0.28406635 0.142033173
[61,] 0.8489347 0.30213052 0.151065262
[62,] 0.8390905 0.32181899 0.160909493
[63,] 0.8115056 0.37698876 0.188494381
[64,] 0.7910664 0.41786729 0.208933647
[65,] 0.7593955 0.48120896 0.240604479
[66,] 0.9514610 0.09707802 0.048539012
[67,] 0.9442973 0.11140538 0.055702692
[68,] 0.9309788 0.13804246 0.069021228
[69,] 0.9155268 0.16894633 0.084473167
[70,] 0.8977377 0.20452464 0.102262321
[71,] 0.8859264 0.22814724 0.114073622
[72,] 0.9841733 0.03165339 0.015826695
[73,] 0.9837781 0.03244389 0.016221946
[74,] 0.9872543 0.02549137 0.012745687
[75,] 0.9853202 0.02935959 0.014679793
[76,] 0.9815136 0.03697280 0.018486398
[77,] 0.9791882 0.04162355 0.020811775
[78,] 0.9741886 0.05162286 0.025811430
[79,] 0.9739777 0.05204469 0.026022346
[80,] 0.9684984 0.06300313 0.031501565
[81,] 0.9907666 0.01846690 0.009233448
[82,] 0.9884246 0.02315076 0.011575381
[83,] 0.9853246 0.02935078 0.014675392
[84,] 0.9810873 0.03782535 0.018912677
[85,] 0.9755217 0.04895654 0.024478268
[86,] 0.9827386 0.03452276 0.017261378
[87,] 0.9792782 0.04144352 0.020721761
[88,] 0.9734532 0.05309352 0.026546760
[89,] 0.9674742 0.06505165 0.032525823
[90,] 0.9859834 0.02803311 0.014016553
[91,] 0.9820715 0.03585692 0.017928462
[92,] 0.9777421 0.04451590 0.022257948
[93,] 0.9747007 0.05059859 0.025299295
[94,] 0.9707533 0.05849340 0.029246701
[95,] 0.9635682 0.07286352 0.036431760
[96,] 0.9564376 0.08712487 0.043562433
[97,] 0.9452206 0.10955870 0.054779351
[98,] 0.9417137 0.11657265 0.058286323
[99,] 0.9288176 0.14236479 0.071182396
[100,] 0.9133421 0.17331586 0.086657932
[101,] 0.9354184 0.12916318 0.064581590
[102,] 0.9779452 0.04410965 0.022054825
[103,] 0.9724840 0.05503208 0.027516038
[104,] 0.9644005 0.07119899 0.035599495
[105,] 0.9724882 0.05502369 0.027511843
[106,] 0.9706236 0.05875274 0.029376371
[107,] 0.9736877 0.05262451 0.026312255
[108,] 0.9712189 0.05756221 0.028781107
[109,] 0.9636783 0.07264339 0.036321693
[110,] 0.9543643 0.09127131 0.045635657
[111,] 0.9418316 0.11633687 0.058168437
[112,] 0.9699936 0.06001285 0.030006423
[113,] 0.9610653 0.07786942 0.038934710
[114,] 0.9505456 0.09890890 0.049454449
[115,] 0.9436856 0.11262888 0.056314442
[116,] 0.9319795 0.13604091 0.068020456
[117,] 0.9663623 0.06727541 0.033637704
[118,] 0.9567190 0.08656202 0.043281010
[119,] 0.9464212 0.10715768 0.053578839
[120,] 0.9311360 0.13772795 0.068863973
[121,] 0.9192856 0.16142882 0.080714408
[122,] 0.9061296 0.18774090 0.093870450
[123,] 0.8825585 0.23488295 0.117441474
[124,] 0.8797396 0.24052079 0.120260394
[125,] 0.8713149 0.25737025 0.128685126
[126,] 0.8549396 0.29012072 0.145060361
[127,] 0.8831307 0.23373859 0.116869294
[128,] 0.8777742 0.24445170 0.122225848
[129,] 0.8678096 0.26438084 0.132190419
[130,] 0.8828236 0.23435289 0.117176446
[131,] 0.8795583 0.24088345 0.120441726
[132,] 0.8517718 0.29645649 0.148228246
[133,] 0.8184119 0.36317622 0.181588109
[134,] 0.8351282 0.32974356 0.164871782
[135,] 0.8004462 0.39910761 0.199553806
[136,] 0.7817648 0.43647041 0.218235207
[137,] 0.7416725 0.51665500 0.258327500
[138,] 0.7163129 0.56737425 0.283687125
[139,] 0.6682145 0.66357109 0.331785546
[140,] 0.6462343 0.70753142 0.353765708
[141,] 0.6891198 0.62176041 0.310880204
[142,] 0.6989147 0.60217070 0.301085349
[143,] 0.7135172 0.57296560 0.286482798
[144,] 0.6742470 0.65150608 0.325753041
[145,] 0.6840275 0.63194494 0.315972472
[146,] 0.6327578 0.73448438 0.367242188
[147,] 0.7054043 0.58919148 0.294595742
[148,] 0.6837321 0.63253587 0.316267934
[149,] 0.7546363 0.49072732 0.245363661
[150,] 0.6940491 0.61190189 0.305950947
[151,] 0.6319082 0.73618353 0.368091766
[152,] 0.5657087 0.86858269 0.434291346
[153,] 0.9928024 0.01439519 0.007197593
[154,] 0.9896305 0.02073892 0.010369458
[155,] 0.9813828 0.03723430 0.018617152
[156,] 0.9824644 0.03507111 0.017535555
[157,] 0.9694695 0.06106103 0.030530514
[158,] 0.9418167 0.11636653 0.058183267
[159,] 0.9056212 0.18875763 0.094378816
[160,] 0.8412857 0.31742862 0.158714309
[161,] 0.8568520 0.28629600 0.143147999
[162,] 0.7223970 0.55520592 0.277602962
> postscript(file="/var/wessaorg/rcomp/tmp/1z52h1324643131.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/218hu1324643131.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/3zyms1324643131.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/40gol1324643131.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/5nqwy1324643131.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 = 197
Frequency = 1
1 2 3 4 5 6
18629.8426 -19962.3868 2184.1463 -40048.9234 12794.7695 -16268.2434
7 8 9 10 11 12
17761.2362 4034.6184 19443.0666 -19396.8272 16196.9733 -2541.3402
13 14 15 16 17 18
5425.4600 -3381.2409 40355.6716 2368.7121 -43926.2356 24358.4350
19 20 21 22 23 24
7490.1493 2718.0440 6544.7869 13980.7583 88702.3619 7110.1992
25 26 27 28 29 30
-23088.4456 -51127.6947 -37185.0578 -30631.9639 19016.0202 -3624.6022
31 32 33 34 35 36
-44999.6443 -7681.7119 -8271.9471 6863.3007 324.6237 3832.0597
37 38 39 40 41 42
5105.7930 7660.7811 36272.6738 2312.6218 7607.6430 22896.9178
43 44 45 46 47 48
-13427.8281 -2617.7557 -11891.3992 -2584.6184 -15119.5528 -15924.3746
49 50 51 52 53 54
7456.4408 -46368.5505 35982.0311 -23908.3645 -5924.3297 -48531.8440
55 56 57 58 59 60
-37799.2675 -23817.5431 -12437.3571 38039.4125 -1284.3325 -9250.3179
61 62 63 64 65 66
6206.1226 1033.5098 -1038.7070 3611.3823 -37376.8315 -10165.2951
67 68 69 70 71 72
46776.9707 -14359.9140 -30164.9411 -31742.6557 9575.6570 14044.1224
73 74 75 76 77 78
7102.5924 -4318.3766 -35733.9604 6704.4889 -12845.5964 -16410.6471
79 80 81 82 83 84
-16505.1135 736.9085 9572.9482 -9788.5944 71791.5362 17559.6607
85 86 87 88 89 90
1900.6016 -5279.1195 1866.3606 -13422.6638 72243.0732 -17358.5891
91 92 93 94 95 96
40610.3780 -15336.7579 -14409.9581 18355.0574 10258.9850 10506.8156
97 98 99 100 101 102
15737.5982 -51638.6012 15628.7027 8822.9597 -5218.1691 3580.1255
103 104 105 106 107 108
31214.1987 12453.1885 7076.8849 -1495.2913 -32526.3699 -3907.4730
109 110 111 112 113 114
8400.5732 23384.0109 -16440.4778 -8124.7723 -10478.7656 1301.8540
115 116 117 118 119 120
-17108.6625 -9579.6335 8317.7902 33935.5337 -59568.4654 -7419.6903
121 122 123 124 125 126
-597.4792 -26184.5607 -31616.2747 -35230.7374 25250.7260 9199.2714
127 128 129 130 131 132
7337.5456 275.1696 46619.1231 5080.4902 -501.2017 -13899.2061
133 134 135 136 137 138
16639.2735 -40217.1281 -4002.1755 -6744.8806 1725.5863 12143.4851
139 140 141 142 143 144
-9890.0683 -1836.4771 -17837.2662 -12415.7544 13978.4754 -33224.2253
145 146 147 148 149 150
29141.9994 -9001.2417 24354.5818 -15436.9291 -11472.5724 4615.0270
151 152 153 154 155 156
23842.3100 -5895.7631 -5761.9076 -1420.6442 -28939.4022 17334.6605
157 158 159 160 161 162
4111.4414 11467.1983 -24417.9029 15131.1734 -3495.3261 9896.1910
163 164 165 166 167 168
1162.9119 23217.8904 -14665.4880 15258.8436 -14057.2167 23064.5841
169 170 171 172 173 174
2592.4714 74217.4063 6379.7798 10122.7420 45237.1484 27024.1064
175 176 177 178 179 180
-35374.2115 -3354.5249 -24540.9160 12504.2328 2948.5116 28517.2492
181 182 183 184 185 186
-11267.2279 -21610.0609 5360.4479 10629.4041 -3779.8793 1873.3622
187 188 189 190 191 192
-350.4889 10399.6776 -6999.8265 29746.8729 44803.6374 2797.8458
193 194 195 196 197
-11060.3048 -18322.4501 -18945.7572 10709.7038 -13337.4288
> postscript(file="/var/wessaorg/rcomp/tmp/6r1sv1324643131.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 = 197
Frequency = 1
lag(myerror, k = 1) myerror
0 18629.8426 NA
1 -19962.3868 18629.8426
2 2184.1463 -19962.3868
3 -40048.9234 2184.1463
4 12794.7695 -40048.9234
5 -16268.2434 12794.7695
6 17761.2362 -16268.2434
7 4034.6184 17761.2362
8 19443.0666 4034.6184
9 -19396.8272 19443.0666
10 16196.9733 -19396.8272
11 -2541.3402 16196.9733
12 5425.4600 -2541.3402
13 -3381.2409 5425.4600
14 40355.6716 -3381.2409
15 2368.7121 40355.6716
16 -43926.2356 2368.7121
17 24358.4350 -43926.2356
18 7490.1493 24358.4350
19 2718.0440 7490.1493
20 6544.7869 2718.0440
21 13980.7583 6544.7869
22 88702.3619 13980.7583
23 7110.1992 88702.3619
24 -23088.4456 7110.1992
25 -51127.6947 -23088.4456
26 -37185.0578 -51127.6947
27 -30631.9639 -37185.0578
28 19016.0202 -30631.9639
29 -3624.6022 19016.0202
30 -44999.6443 -3624.6022
31 -7681.7119 -44999.6443
32 -8271.9471 -7681.7119
33 6863.3007 -8271.9471
34 324.6237 6863.3007
35 3832.0597 324.6237
36 5105.7930 3832.0597
37 7660.7811 5105.7930
38 36272.6738 7660.7811
39 2312.6218 36272.6738
40 7607.6430 2312.6218
41 22896.9178 7607.6430
42 -13427.8281 22896.9178
43 -2617.7557 -13427.8281
44 -11891.3992 -2617.7557
45 -2584.6184 -11891.3992
46 -15119.5528 -2584.6184
47 -15924.3746 -15119.5528
48 7456.4408 -15924.3746
49 -46368.5505 7456.4408
50 35982.0311 -46368.5505
51 -23908.3645 35982.0311
52 -5924.3297 -23908.3645
53 -48531.8440 -5924.3297
54 -37799.2675 -48531.8440
55 -23817.5431 -37799.2675
56 -12437.3571 -23817.5431
57 38039.4125 -12437.3571
58 -1284.3325 38039.4125
59 -9250.3179 -1284.3325
60 6206.1226 -9250.3179
61 1033.5098 6206.1226
62 -1038.7070 1033.5098
63 3611.3823 -1038.7070
64 -37376.8315 3611.3823
65 -10165.2951 -37376.8315
66 46776.9707 -10165.2951
67 -14359.9140 46776.9707
68 -30164.9411 -14359.9140
69 -31742.6557 -30164.9411
70 9575.6570 -31742.6557
71 14044.1224 9575.6570
72 7102.5924 14044.1224
73 -4318.3766 7102.5924
74 -35733.9604 -4318.3766
75 6704.4889 -35733.9604
76 -12845.5964 6704.4889
77 -16410.6471 -12845.5964
78 -16505.1135 -16410.6471
79 736.9085 -16505.1135
80 9572.9482 736.9085
81 -9788.5944 9572.9482
82 71791.5362 -9788.5944
83 17559.6607 71791.5362
84 1900.6016 17559.6607
85 -5279.1195 1900.6016
86 1866.3606 -5279.1195
87 -13422.6638 1866.3606
88 72243.0732 -13422.6638
89 -17358.5891 72243.0732
90 40610.3780 -17358.5891
91 -15336.7579 40610.3780
92 -14409.9581 -15336.7579
93 18355.0574 -14409.9581
94 10258.9850 18355.0574
95 10506.8156 10258.9850
96 15737.5982 10506.8156
97 -51638.6012 15737.5982
98 15628.7027 -51638.6012
99 8822.9597 15628.7027
100 -5218.1691 8822.9597
101 3580.1255 -5218.1691
102 31214.1987 3580.1255
103 12453.1885 31214.1987
104 7076.8849 12453.1885
105 -1495.2913 7076.8849
106 -32526.3699 -1495.2913
107 -3907.4730 -32526.3699
108 8400.5732 -3907.4730
109 23384.0109 8400.5732
110 -16440.4778 23384.0109
111 -8124.7723 -16440.4778
112 -10478.7656 -8124.7723
113 1301.8540 -10478.7656
114 -17108.6625 1301.8540
115 -9579.6335 -17108.6625
116 8317.7902 -9579.6335
117 33935.5337 8317.7902
118 -59568.4654 33935.5337
119 -7419.6903 -59568.4654
120 -597.4792 -7419.6903
121 -26184.5607 -597.4792
122 -31616.2747 -26184.5607
123 -35230.7374 -31616.2747
124 25250.7260 -35230.7374
125 9199.2714 25250.7260
126 7337.5456 9199.2714
127 275.1696 7337.5456
128 46619.1231 275.1696
129 5080.4902 46619.1231
130 -501.2017 5080.4902
131 -13899.2061 -501.2017
132 16639.2735 -13899.2061
133 -40217.1281 16639.2735
134 -4002.1755 -40217.1281
135 -6744.8806 -4002.1755
136 1725.5863 -6744.8806
137 12143.4851 1725.5863
138 -9890.0683 12143.4851
139 -1836.4771 -9890.0683
140 -17837.2662 -1836.4771
141 -12415.7544 -17837.2662
142 13978.4754 -12415.7544
143 -33224.2253 13978.4754
144 29141.9994 -33224.2253
145 -9001.2417 29141.9994
146 24354.5818 -9001.2417
147 -15436.9291 24354.5818
148 -11472.5724 -15436.9291
149 4615.0270 -11472.5724
150 23842.3100 4615.0270
151 -5895.7631 23842.3100
152 -5761.9076 -5895.7631
153 -1420.6442 -5761.9076
154 -28939.4022 -1420.6442
155 17334.6605 -28939.4022
156 4111.4414 17334.6605
157 11467.1983 4111.4414
158 -24417.9029 11467.1983
159 15131.1734 -24417.9029
160 -3495.3261 15131.1734
161 9896.1910 -3495.3261
162 1162.9119 9896.1910
163 23217.8904 1162.9119
164 -14665.4880 23217.8904
165 15258.8436 -14665.4880
166 -14057.2167 15258.8436
167 23064.5841 -14057.2167
168 2592.4714 23064.5841
169 74217.4063 2592.4714
170 6379.7798 74217.4063
171 10122.7420 6379.7798
172 45237.1484 10122.7420
173 27024.1064 45237.1484
174 -35374.2115 27024.1064
175 -3354.5249 -35374.2115
176 -24540.9160 -3354.5249
177 12504.2328 -24540.9160
178 2948.5116 12504.2328
179 28517.2492 2948.5116
180 -11267.2279 28517.2492
181 -21610.0609 -11267.2279
182 5360.4479 -21610.0609
183 10629.4041 5360.4479
184 -3779.8793 10629.4041
185 1873.3622 -3779.8793
186 -350.4889 1873.3622
187 10399.6776 -350.4889
188 -6999.8265 10399.6776
189 29746.8729 -6999.8265
190 44803.6374 29746.8729
191 2797.8458 44803.6374
192 -11060.3048 2797.8458
193 -18322.4501 -11060.3048
194 -18945.7572 -18322.4501
195 10709.7038 -18945.7572
196 -13337.4288 10709.7038
197 NA -13337.4288
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] -19962.3868 18629.8426
[2,] 2184.1463 -19962.3868
[3,] -40048.9234 2184.1463
[4,] 12794.7695 -40048.9234
[5,] -16268.2434 12794.7695
[6,] 17761.2362 -16268.2434
[7,] 4034.6184 17761.2362
[8,] 19443.0666 4034.6184
[9,] -19396.8272 19443.0666
[10,] 16196.9733 -19396.8272
[11,] -2541.3402 16196.9733
[12,] 5425.4600 -2541.3402
[13,] -3381.2409 5425.4600
[14,] 40355.6716 -3381.2409
[15,] 2368.7121 40355.6716
[16,] -43926.2356 2368.7121
[17,] 24358.4350 -43926.2356
[18,] 7490.1493 24358.4350
[19,] 2718.0440 7490.1493
[20,] 6544.7869 2718.0440
[21,] 13980.7583 6544.7869
[22,] 88702.3619 13980.7583
[23,] 7110.1992 88702.3619
[24,] -23088.4456 7110.1992
[25,] -51127.6947 -23088.4456
[26,] -37185.0578 -51127.6947
[27,] -30631.9639 -37185.0578
[28,] 19016.0202 -30631.9639
[29,] -3624.6022 19016.0202
[30,] -44999.6443 -3624.6022
[31,] -7681.7119 -44999.6443
[32,] -8271.9471 -7681.7119
[33,] 6863.3007 -8271.9471
[34,] 324.6237 6863.3007
[35,] 3832.0597 324.6237
[36,] 5105.7930 3832.0597
[37,] 7660.7811 5105.7930
[38,] 36272.6738 7660.7811
[39,] 2312.6218 36272.6738
[40,] 7607.6430 2312.6218
[41,] 22896.9178 7607.6430
[42,] -13427.8281 22896.9178
[43,] -2617.7557 -13427.8281
[44,] -11891.3992 -2617.7557
[45,] -2584.6184 -11891.3992
[46,] -15119.5528 -2584.6184
[47,] -15924.3746 -15119.5528
[48,] 7456.4408 -15924.3746
[49,] -46368.5505 7456.4408
[50,] 35982.0311 -46368.5505
[51,] -23908.3645 35982.0311
[52,] -5924.3297 -23908.3645
[53,] -48531.8440 -5924.3297
[54,] -37799.2675 -48531.8440
[55,] -23817.5431 -37799.2675
[56,] -12437.3571 -23817.5431
[57,] 38039.4125 -12437.3571
[58,] -1284.3325 38039.4125
[59,] -9250.3179 -1284.3325
[60,] 6206.1226 -9250.3179
[61,] 1033.5098 6206.1226
[62,] -1038.7070 1033.5098
[63,] 3611.3823 -1038.7070
[64,] -37376.8315 3611.3823
[65,] -10165.2951 -37376.8315
[66,] 46776.9707 -10165.2951
[67,] -14359.9140 46776.9707
[68,] -30164.9411 -14359.9140
[69,] -31742.6557 -30164.9411
[70,] 9575.6570 -31742.6557
[71,] 14044.1224 9575.6570
[72,] 7102.5924 14044.1224
[73,] -4318.3766 7102.5924
[74,] -35733.9604 -4318.3766
[75,] 6704.4889 -35733.9604
[76,] -12845.5964 6704.4889
[77,] -16410.6471 -12845.5964
[78,] -16505.1135 -16410.6471
[79,] 736.9085 -16505.1135
[80,] 9572.9482 736.9085
[81,] -9788.5944 9572.9482
[82,] 71791.5362 -9788.5944
[83,] 17559.6607 71791.5362
[84,] 1900.6016 17559.6607
[85,] -5279.1195 1900.6016
[86,] 1866.3606 -5279.1195
[87,] -13422.6638 1866.3606
[88,] 72243.0732 -13422.6638
[89,] -17358.5891 72243.0732
[90,] 40610.3780 -17358.5891
[91,] -15336.7579 40610.3780
[92,] -14409.9581 -15336.7579
[93,] 18355.0574 -14409.9581
[94,] 10258.9850 18355.0574
[95,] 10506.8156 10258.9850
[96,] 15737.5982 10506.8156
[97,] -51638.6012 15737.5982
[98,] 15628.7027 -51638.6012
[99,] 8822.9597 15628.7027
[100,] -5218.1691 8822.9597
[101,] 3580.1255 -5218.1691
[102,] 31214.1987 3580.1255
[103,] 12453.1885 31214.1987
[104,] 7076.8849 12453.1885
[105,] -1495.2913 7076.8849
[106,] -32526.3699 -1495.2913
[107,] -3907.4730 -32526.3699
[108,] 8400.5732 -3907.4730
[109,] 23384.0109 8400.5732
[110,] -16440.4778 23384.0109
[111,] -8124.7723 -16440.4778
[112,] -10478.7656 -8124.7723
[113,] 1301.8540 -10478.7656
[114,] -17108.6625 1301.8540
[115,] -9579.6335 -17108.6625
[116,] 8317.7902 -9579.6335
[117,] 33935.5337 8317.7902
[118,] -59568.4654 33935.5337
[119,] -7419.6903 -59568.4654
[120,] -597.4792 -7419.6903
[121,] -26184.5607 -597.4792
[122,] -31616.2747 -26184.5607
[123,] -35230.7374 -31616.2747
[124,] 25250.7260 -35230.7374
[125,] 9199.2714 25250.7260
[126,] 7337.5456 9199.2714
[127,] 275.1696 7337.5456
[128,] 46619.1231 275.1696
[129,] 5080.4902 46619.1231
[130,] -501.2017 5080.4902
[131,] -13899.2061 -501.2017
[132,] 16639.2735 -13899.2061
[133,] -40217.1281 16639.2735
[134,] -4002.1755 -40217.1281
[135,] -6744.8806 -4002.1755
[136,] 1725.5863 -6744.8806
[137,] 12143.4851 1725.5863
[138,] -9890.0683 12143.4851
[139,] -1836.4771 -9890.0683
[140,] -17837.2662 -1836.4771
[141,] -12415.7544 -17837.2662
[142,] 13978.4754 -12415.7544
[143,] -33224.2253 13978.4754
[144,] 29141.9994 -33224.2253
[145,] -9001.2417 29141.9994
[146,] 24354.5818 -9001.2417
[147,] -15436.9291 24354.5818
[148,] -11472.5724 -15436.9291
[149,] 4615.0270 -11472.5724
[150,] 23842.3100 4615.0270
[151,] -5895.7631 23842.3100
[152,] -5761.9076 -5895.7631
[153,] -1420.6442 -5761.9076
[154,] -28939.4022 -1420.6442
[155,] 17334.6605 -28939.4022
[156,] 4111.4414 17334.6605
[157,] 11467.1983 4111.4414
[158,] -24417.9029 11467.1983
[159,] 15131.1734 -24417.9029
[160,] -3495.3261 15131.1734
[161,] 9896.1910 -3495.3261
[162,] 1162.9119 9896.1910
[163,] 23217.8904 1162.9119
[164,] -14665.4880 23217.8904
[165,] 15258.8436 -14665.4880
[166,] -14057.2167 15258.8436
[167,] 23064.5841 -14057.2167
[168,] 2592.4714 23064.5841
[169,] 74217.4063 2592.4714
[170,] 6379.7798 74217.4063
[171,] 10122.7420 6379.7798
[172,] 45237.1484 10122.7420
[173,] 27024.1064 45237.1484
[174,] -35374.2115 27024.1064
[175,] -3354.5249 -35374.2115
[176,] -24540.9160 -3354.5249
[177,] 12504.2328 -24540.9160
[178,] 2948.5116 12504.2328
[179,] 28517.2492 2948.5116
[180,] -11267.2279 28517.2492
[181,] -21610.0609 -11267.2279
[182,] 5360.4479 -21610.0609
[183,] 10629.4041 5360.4479
[184,] -3779.8793 10629.4041
[185,] 1873.3622 -3779.8793
[186,] -350.4889 1873.3622
[187,] 10399.6776 -350.4889
[188,] -6999.8265 10399.6776
[189,] 29746.8729 -6999.8265
[190,] 44803.6374 29746.8729
[191,] 2797.8458 44803.6374
[192,] -11060.3048 2797.8458
[193,] -18322.4501 -11060.3048
[194,] -18945.7572 -18322.4501
[195,] 10709.7038 -18945.7572
[196,] -13337.4288 10709.7038
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 -19962.3868 18629.8426
2 2184.1463 -19962.3868
3 -40048.9234 2184.1463
4 12794.7695 -40048.9234
5 -16268.2434 12794.7695
6 17761.2362 -16268.2434
7 4034.6184 17761.2362
8 19443.0666 4034.6184
9 -19396.8272 19443.0666
10 16196.9733 -19396.8272
11 -2541.3402 16196.9733
12 5425.4600 -2541.3402
13 -3381.2409 5425.4600
14 40355.6716 -3381.2409
15 2368.7121 40355.6716
16 -43926.2356 2368.7121
17 24358.4350 -43926.2356
18 7490.1493 24358.4350
19 2718.0440 7490.1493
20 6544.7869 2718.0440
21 13980.7583 6544.7869
22 88702.3619 13980.7583
23 7110.1992 88702.3619
24 -23088.4456 7110.1992
25 -51127.6947 -23088.4456
26 -37185.0578 -51127.6947
27 -30631.9639 -37185.0578
28 19016.0202 -30631.9639
29 -3624.6022 19016.0202
30 -44999.6443 -3624.6022
31 -7681.7119 -44999.6443
32 -8271.9471 -7681.7119
33 6863.3007 -8271.9471
34 324.6237 6863.3007
35 3832.0597 324.6237
36 5105.7930 3832.0597
37 7660.7811 5105.7930
38 36272.6738 7660.7811
39 2312.6218 36272.6738
40 7607.6430 2312.6218
41 22896.9178 7607.6430
42 -13427.8281 22896.9178
43 -2617.7557 -13427.8281
44 -11891.3992 -2617.7557
45 -2584.6184 -11891.3992
46 -15119.5528 -2584.6184
47 -15924.3746 -15119.5528
48 7456.4408 -15924.3746
49 -46368.5505 7456.4408
50 35982.0311 -46368.5505
51 -23908.3645 35982.0311
52 -5924.3297 -23908.3645
53 -48531.8440 -5924.3297
54 -37799.2675 -48531.8440
55 -23817.5431 -37799.2675
56 -12437.3571 -23817.5431
57 38039.4125 -12437.3571
58 -1284.3325 38039.4125
59 -9250.3179 -1284.3325
60 6206.1226 -9250.3179
61 1033.5098 6206.1226
62 -1038.7070 1033.5098
63 3611.3823 -1038.7070
64 -37376.8315 3611.3823
65 -10165.2951 -37376.8315
66 46776.9707 -10165.2951
67 -14359.9140 46776.9707
68 -30164.9411 -14359.9140
69 -31742.6557 -30164.9411
70 9575.6570 -31742.6557
71 14044.1224 9575.6570
72 7102.5924 14044.1224
73 -4318.3766 7102.5924
74 -35733.9604 -4318.3766
75 6704.4889 -35733.9604
76 -12845.5964 6704.4889
77 -16410.6471 -12845.5964
78 -16505.1135 -16410.6471
79 736.9085 -16505.1135
80 9572.9482 736.9085
81 -9788.5944 9572.9482
82 71791.5362 -9788.5944
83 17559.6607 71791.5362
84 1900.6016 17559.6607
85 -5279.1195 1900.6016
86 1866.3606 -5279.1195
87 -13422.6638 1866.3606
88 72243.0732 -13422.6638
89 -17358.5891 72243.0732
90 40610.3780 -17358.5891
91 -15336.7579 40610.3780
92 -14409.9581 -15336.7579
93 18355.0574 -14409.9581
94 10258.9850 18355.0574
95 10506.8156 10258.9850
96 15737.5982 10506.8156
97 -51638.6012 15737.5982
98 15628.7027 -51638.6012
99 8822.9597 15628.7027
100 -5218.1691 8822.9597
101 3580.1255 -5218.1691
102 31214.1987 3580.1255
103 12453.1885 31214.1987
104 7076.8849 12453.1885
105 -1495.2913 7076.8849
106 -32526.3699 -1495.2913
107 -3907.4730 -32526.3699
108 8400.5732 -3907.4730
109 23384.0109 8400.5732
110 -16440.4778 23384.0109
111 -8124.7723 -16440.4778
112 -10478.7656 -8124.7723
113 1301.8540 -10478.7656
114 -17108.6625 1301.8540
115 -9579.6335 -17108.6625
116 8317.7902 -9579.6335
117 33935.5337 8317.7902
118 -59568.4654 33935.5337
119 -7419.6903 -59568.4654
120 -597.4792 -7419.6903
121 -26184.5607 -597.4792
122 -31616.2747 -26184.5607
123 -35230.7374 -31616.2747
124 25250.7260 -35230.7374
125 9199.2714 25250.7260
126 7337.5456 9199.2714
127 275.1696 7337.5456
128 46619.1231 275.1696
129 5080.4902 46619.1231
130 -501.2017 5080.4902
131 -13899.2061 -501.2017
132 16639.2735 -13899.2061
133 -40217.1281 16639.2735
134 -4002.1755 -40217.1281
135 -6744.8806 -4002.1755
136 1725.5863 -6744.8806
137 12143.4851 1725.5863
138 -9890.0683 12143.4851
139 -1836.4771 -9890.0683
140 -17837.2662 -1836.4771
141 -12415.7544 -17837.2662
142 13978.4754 -12415.7544
143 -33224.2253 13978.4754
144 29141.9994 -33224.2253
145 -9001.2417 29141.9994
146 24354.5818 -9001.2417
147 -15436.9291 24354.5818
148 -11472.5724 -15436.9291
149 4615.0270 -11472.5724
150 23842.3100 4615.0270
151 -5895.7631 23842.3100
152 -5761.9076 -5895.7631
153 -1420.6442 -5761.9076
154 -28939.4022 -1420.6442
155 17334.6605 -28939.4022
156 4111.4414 17334.6605
157 11467.1983 4111.4414
158 -24417.9029 11467.1983
159 15131.1734 -24417.9029
160 -3495.3261 15131.1734
161 9896.1910 -3495.3261
162 1162.9119 9896.1910
163 23217.8904 1162.9119
164 -14665.4880 23217.8904
165 15258.8436 -14665.4880
166 -14057.2167 15258.8436
167 23064.5841 -14057.2167
168 2592.4714 23064.5841
169 74217.4063 2592.4714
170 6379.7798 74217.4063
171 10122.7420 6379.7798
172 45237.1484 10122.7420
173 27024.1064 45237.1484
174 -35374.2115 27024.1064
175 -3354.5249 -35374.2115
176 -24540.9160 -3354.5249
177 12504.2328 -24540.9160
178 2948.5116 12504.2328
179 28517.2492 2948.5116
180 -11267.2279 28517.2492
181 -21610.0609 -11267.2279
182 5360.4479 -21610.0609
183 10629.4041 5360.4479
184 -3779.8793 10629.4041
185 1873.3622 -3779.8793
186 -350.4889 1873.3622
187 10399.6776 -350.4889
188 -6999.8265 10399.6776
189 29746.8729 -6999.8265
190 44803.6374 29746.8729
191 2797.8458 44803.6374
192 -11060.3048 2797.8458
193 -18322.4501 -11060.3048
194 -18945.7572 -18322.4501
195 10709.7038 -18945.7572
196 -13337.4288 10709.7038
> 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/7ldo51324643131.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/88o371324643131.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/9t6ov1324643131.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/10f5qe1324643131.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/114c741324643131.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/12tvqh1324643131.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/13mmfd1324643131.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/14n8991324643131.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/15lxft1324643131.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/16mj0g1324643131.tab")
+ }
>
> try(system("convert tmp/1z52h1324643131.ps tmp/1z52h1324643131.png",intern=TRUE))
character(0)
> try(system("convert tmp/218hu1324643131.ps tmp/218hu1324643131.png",intern=TRUE))
character(0)
> try(system("convert tmp/3zyms1324643131.ps tmp/3zyms1324643131.png",intern=TRUE))
character(0)
> try(system("convert tmp/40gol1324643131.ps tmp/40gol1324643131.png",intern=TRUE))
character(0)
> try(system("convert tmp/5nqwy1324643131.ps tmp/5nqwy1324643131.png",intern=TRUE))
character(0)
> try(system("convert tmp/6r1sv1324643131.ps tmp/6r1sv1324643131.png",intern=TRUE))
character(0)
> try(system("convert tmp/7ldo51324643131.ps tmp/7ldo51324643131.png",intern=TRUE))
character(0)
> try(system("convert tmp/88o371324643131.ps tmp/88o371324643131.png",intern=TRUE))
character(0)
> try(system("convert tmp/9t6ov1324643131.ps tmp/9t6ov1324643131.png",intern=TRUE))
character(0)
> try(system("convert tmp/10f5qe1324643131.ps tmp/10f5qe1324643131.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
7.630 0.703 8.342