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(1818
+ ,279055
+ ,73
+ ,504
+ ,95
+ ,3
+ ,96
+ ,42
+ ,159
+ ,130
+ ,140824
+ ,32033
+ ,186099
+ ,165
+ ,165
+ ,1433
+ ,212408
+ ,75
+ ,510
+ ,68
+ ,4
+ ,75
+ ,38
+ ,149
+ ,143
+ ,110459
+ ,20654
+ ,113854
+ ,135
+ ,132
+ ,2059
+ ,233939
+ ,83
+ ,710
+ ,64
+ ,16
+ ,70
+ ,46
+ ,178
+ ,118
+ ,105079
+ ,16346
+ ,99776
+ ,121
+ ,121
+ ,2733
+ ,222117
+ ,106
+ ,1154
+ ,139
+ ,2
+ ,134
+ ,42
+ ,164
+ ,146
+ ,112098
+ ,35926
+ ,106194
+ ,148
+ ,145
+ ,1399
+ ,189911
+ ,56
+ ,415
+ ,51
+ ,1
+ ,83
+ ,30
+ ,100
+ ,73
+ ,43929
+ ,10621
+ ,100792
+ ,73
+ ,71
+ ,631
+ ,70849
+ ,28
+ ,179
+ ,46
+ ,3
+ ,8
+ ,35
+ ,129
+ ,89
+ ,76173
+ ,10024
+ ,47552
+ ,49
+ ,47
+ ,5460
+ ,605767
+ ,135
+ ,2563
+ ,118
+ ,0
+ ,173
+ ,40
+ ,156
+ ,146
+ ,187326
+ ,43068
+ ,250931
+ ,185
+ ,177
+ ,381
+ ,33186
+ ,19
+ ,111
+ ,46
+ ,0
+ ,1
+ ,18
+ ,67
+ ,22
+ ,22807
+ ,1271
+ ,6853
+ ,5
+ ,5
+ ,2150
+ ,227332
+ ,62
+ ,763
+ ,79
+ ,7
+ ,88
+ ,38
+ ,148
+ ,132
+ ,144408
+ ,34416
+ ,115466
+ ,125
+ ,124
+ ,2042
+ ,267925
+ ,49
+ ,661
+ ,76
+ ,0
+ ,104
+ ,37
+ ,132
+ ,92
+ ,66485
+ ,20318
+ ,110896
+ ,93
+ ,92
+ ,2536
+ ,371987
+ ,122
+ ,981
+ ,82
+ ,0
+ ,114
+ ,46
+ ,169
+ ,147
+ ,79089
+ ,24409
+ ,169351
+ ,154
+ ,149
+ ,2377
+ ,264989
+ ,131
+ ,733
+ ,66
+ ,7
+ ,125
+ ,60
+ ,230
+ ,203
+ ,81625
+ ,20648
+ ,94853
+ ,98
+ ,93
+ ,2100
+ ,212638
+ ,87
+ ,785
+ ,60
+ ,10
+ ,57
+ ,37
+ ,122
+ ,113
+ ,68788
+ ,12347
+ ,72591
+ ,70
+ ,70
+ ,3020
+ ,368577
+ ,85
+ ,1186
+ ,117
+ ,4
+ ,139
+ ,55
+ ,191
+ ,171
+ ,103297
+ ,21857
+ ,101345
+ ,148
+ ,148
+ ,2265
+ ,269455
+ ,88
+ ,724
+ ,50
+ ,10
+ ,87
+ ,44
+ ,162
+ ,87
+ ,69446
+ ,11034
+ ,113713
+ ,100
+ ,100
+ ,5139
+ ,398124
+ ,191
+ ,1774
+ ,133
+ ,0
+ ,176
+ ,63
+ ,237
+ ,208
+ ,114948
+ ,33433
+ ,165354
+ ,150
+ ,142
+ ,2363
+ ,335567
+ ,77
+ ,845
+ ,63
+ ,8
+ ,114
+ ,40
+ ,156
+ ,153
+ ,167949
+ ,35902
+ ,164263
+ ,197
+ ,194
+ ,3564
+ ,432711
+ ,173
+ ,1390
+ ,100
+ ,4
+ ,121
+ ,43
+ ,157
+ ,97
+ ,125081
+ ,22355
+ ,135213
+ ,114
+ ,113
+ ,1477
+ ,182016
+ ,58
+ ,514
+ ,44
+ ,3
+ ,103
+ ,32
+ ,123
+ ,95
+ ,125818
+ ,31219
+ ,111669
+ ,169
+ ,162
+ ,2398
+ ,267365
+ ,89
+ ,692
+ ,65
+ ,8
+ ,135
+ ,52
+ ,203
+ ,197
+ ,136588
+ ,21983
+ ,134163
+ ,200
+ ,186
+ ,2546
+ ,279428
+ ,73
+ ,847
+ ,103
+ ,0
+ ,123
+ ,49
+ ,187
+ ,160
+ ,112431
+ ,40085
+ ,140303
+ ,148
+ ,147
+ ,3150
+ ,508849
+ ,111
+ ,1397
+ ,103
+ ,1
+ ,99
+ ,41
+ ,152
+ ,148
+ ,103037
+ ,18507
+ ,150773
+ ,140
+ ,137
+ ,1705
+ ,220142
+ ,49
+ ,569
+ ,62
+ ,5
+ ,77
+ ,25
+ ,89
+ ,84
+ ,82317
+ ,16278
+ ,111848
+ ,74
+ ,71
+ ,1787
+ ,200004
+ ,58
+ ,636
+ ,70
+ ,9
+ ,103
+ ,57
+ ,227
+ ,227
+ ,118906
+ ,24662
+ ,102509
+ ,128
+ ,123
+ ,3792
+ ,257139
+ ,133
+ ,1370
+ ,159
+ ,1
+ ,158
+ ,45
+ ,165
+ ,154
+ ,83515
+ ,31452
+ ,96785
+ ,140
+ ,134
+ ,3108
+ ,270941
+ ,138
+ ,1092
+ ,78
+ ,0
+ ,116
+ ,42
+ ,162
+ ,151
+ ,104581
+ ,32580
+ ,116136
+ ,116
+ ,115
+ ,3230
+ ,324969
+ ,134
+ ,1201
+ ,101
+ ,5
+ ,114
+ ,45
+ ,174
+ ,142
+ ,103129
+ ,22883
+ ,158376
+ ,147
+ ,138
+ ,2348
+ ,329962
+ ,92
+ ,763
+ ,73
+ ,0
+ ,150
+ ,43
+ ,154
+ ,148
+ ,83243
+ ,27652
+ ,153990
+ ,132
+ ,125
+ ,1780
+ ,190867
+ ,60
+ ,652
+ ,58
+ ,0
+ ,64
+ ,36
+ ,129
+ ,110
+ ,37110
+ ,9845
+ ,64057
+ ,70
+ ,66
+ ,3218
+ ,393860
+ ,79
+ ,1213
+ ,147
+ ,0
+ ,150
+ ,45
+ ,174
+ ,149
+ ,113344
+ ,20190
+ ,230054
+ ,144
+ ,137
+ ,2692
+ ,327660
+ ,89
+ ,1111
+ ,54
+ ,3
+ ,143
+ ,50
+ ,195
+ ,179
+ ,139165
+ ,46201
+ ,184531
+ ,155
+ ,152
+ ,2187
+ ,269239
+ ,83
+ ,758
+ ,84
+ ,6
+ ,50
+ ,50
+ ,186
+ ,149
+ ,86652
+ ,10971
+ ,114198
+ ,165
+ ,159
+ ,2577
+ ,396136
+ ,106
+ ,906
+ ,56
+ ,1
+ ,145
+ ,51
+ ,197
+ ,187
+ ,112302
+ ,34811
+ ,198299
+ ,161
+ ,159
+ ,1293
+ ,130446
+ ,49
+ ,456
+ ,45
+ ,4
+ ,56
+ ,42
+ ,157
+ ,153
+ ,69652
+ ,3029
+ ,33750
+ ,31
+ ,31
+ ,3567
+ ,430118
+ ,104
+ ,1293
+ ,87
+ ,4
+ ,141
+ ,44
+ ,168
+ ,163
+ ,119442
+ ,38941
+ ,189723
+ ,199
+ ,185
+ ,2764
+ ,273950
+ ,56
+ ,1186
+ ,87
+ ,0
+ ,83
+ ,42
+ ,159
+ ,127
+ ,69867
+ ,4958
+ ,100826
+ ,78
+ ,78
+ ,3755
+ ,428077
+ ,128
+ ,1348
+ ,77
+ ,0
+ ,112
+ ,44
+ ,161
+ ,151
+ ,101629
+ ,32344
+ ,188355
+ ,121
+ ,117
+ ,2075
+ ,254312
+ ,93
+ ,695
+ ,72
+ ,2
+ ,79
+ ,40
+ ,153
+ ,100
+ ,70168
+ ,19433
+ ,104470
+ ,112
+ ,109
+ ,995
+ ,120351
+ ,35
+ ,306
+ ,36
+ ,1
+ ,33
+ ,17
+ ,55
+ ,46
+ ,31081
+ ,12558
+ ,58391
+ ,41
+ ,41
+ ,3750
+ ,395658
+ ,212
+ ,1319
+ ,51
+ ,2
+ ,152
+ ,43
+ ,166
+ ,156
+ ,103925
+ ,36524
+ ,164808
+ ,158
+ ,149
+ ,3413
+ ,345875
+ ,86
+ ,1566
+ ,44
+ ,10
+ ,126
+ ,41
+ ,151
+ ,128
+ ,92622
+ ,26041
+ ,134097
+ ,123
+ ,123
+ ,2053
+ ,216827
+ ,82
+ ,784
+ ,75
+ ,10
+ ,97
+ ,41
+ ,148
+ ,111
+ ,79011
+ ,16637
+ ,80238
+ ,104
+ ,103
+ ,1984
+ ,224524
+ ,83
+ ,730
+ ,87
+ ,5
+ ,84
+ ,40
+ ,129
+ ,119
+ ,93487
+ ,28395
+ ,133252
+ ,94
+ ,87
+ ,1825
+ ,182485
+ ,69
+ ,488
+ ,97
+ ,6
+ ,68
+ ,49
+ ,181
+ ,148
+ ,64520
+ ,16747
+ ,54518
+ ,73
+ ,71
+ ,2599
+ ,157164
+ ,85
+ ,1051
+ ,90
+ ,1
+ ,50
+ ,52
+ ,93
+ ,65
+ ,93473
+ ,9105
+ ,121850
+ ,52
+ ,51
+ ,5572
+ ,459455
+ ,157
+ ,2089
+ ,860
+ ,2
+ ,101
+ ,42
+ ,150
+ ,134
+ ,114360
+ ,11941
+ ,79367
+ ,71
+ ,70
+ ,918
+ ,78800
+ ,42
+ ,330
+ ,57
+ ,2
+ ,20
+ ,26
+ ,82
+ ,66
+ ,33032
+ ,7935
+ ,56968
+ ,21
+ ,21
+ ,2685
+ ,255072
+ ,85
+ ,764
+ ,99
+ ,1
+ ,107
+ ,59
+ ,229
+ ,201
+ ,96125
+ ,19499
+ ,106314
+ ,155
+ ,155
+ ,4145
+ ,368086
+ ,123
+ ,1410
+ ,120
+ ,10
+ ,150
+ ,50
+ ,193
+ ,177
+ ,151911
+ ,22938
+ ,191889
+ ,174
+ ,172
+ ,2841
+ ,230299
+ ,70
+ ,1187
+ ,76
+ ,3
+ ,129
+ ,50
+ ,176
+ ,156
+ ,89256
+ ,25314
+ ,104864
+ ,136
+ ,133
+ ,2175
+ ,244782
+ ,81
+ ,691
+ ,56
+ ,0
+ ,99
+ ,47
+ ,179
+ ,158
+ ,95676
+ ,28527
+ ,160792
+ ,128
+ ,125
+ ,496
+ ,24188
+ ,24
+ ,218
+ ,20
+ ,0
+ ,8
+ ,4
+ ,12
+ ,7
+ ,5950
+ ,2694
+ ,15049
+ ,7
+ ,7
+ ,2699
+ ,400109
+ ,334
+ ,865
+ ,94
+ ,8
+ ,88
+ ,51
+ ,181
+ ,175
+ ,149695
+ ,20867
+ ,191179
+ ,165
+ ,158
+ ,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
+ ,3333
+ ,309810
+ ,67
+ ,1229
+ ,133
+ ,1
+ ,102
+ ,41
+ ,148
+ ,133
+ ,100087
+ ,32982
+ ,129711
+ ,137
+ ,133
+ ,2970
+ ,375638
+ ,91
+ ,790
+ ,86
+ ,5
+ ,166
+ ,61
+ ,230
+ ,228
+ ,169707
+ ,38975
+ ,210012
+ ,174
+ ,169
+ ,3968
+ ,367127
+ ,204
+ ,1208
+ ,224
+ ,6
+ ,132
+ ,40
+ ,148
+ ,140
+ ,150491
+ ,42721
+ ,194679
+ ,257
+ ,256
+ ,2878
+ ,381998
+ ,155
+ ,1102
+ ,65
+ ,0
+ ,161
+ ,44
+ ,160
+ ,155
+ ,120192
+ ,41455
+ ,197680
+ ,207
+ ,190
+ ,2399
+ ,280106
+ ,90
+ ,919
+ ,86
+ ,12
+ ,90
+ ,40
+ ,155
+ ,141
+ ,95893
+ ,23923
+ ,81180
+ ,103
+ ,100
+ ,4121
+ ,400971
+ ,153
+ ,1352
+ ,70
+ ,10
+ ,160
+ ,51
+ ,198
+ ,181
+ ,151715
+ ,26719
+ ,197765
+ ,171
+ ,171
+ ,3294
+ ,315924
+ ,122
+ ,1190
+ ,148
+ ,12
+ ,139
+ ,29
+ ,104
+ ,75
+ ,176225
+ ,53405
+ ,214738
+ ,279
+ ,267
+ ,3132
+ ,291391
+ ,124
+ ,1257
+ ,72
+ ,11
+ ,104
+ ,43
+ ,169
+ ,97
+ ,59900
+ ,12526
+ ,96252
+ ,83
+ ,80
+ ,2868
+ ,295075
+ ,93
+ ,1030
+ ,59
+ ,8
+ ,103
+ ,42
+ ,163
+ ,142
+ ,104767
+ ,26584
+ ,124527
+ ,130
+ ,126
+ ,1778
+ ,280018
+ ,81
+ ,669
+ ,67
+ ,3
+ ,66
+ ,41
+ ,151
+ ,136
+ ,114799
+ ,37062
+ ,153242
+ ,131
+ ,132
+ ,2109
+ ,267432
+ ,71
+ ,542
+ ,58
+ ,0
+ ,163
+ ,30
+ ,116
+ ,87
+ ,72128
+ ,25696
+ ,145707
+ ,126
+ ,121
+ ,2148
+ ,217181
+ ,141
+ ,652
+ ,60
+ ,6
+ ,93
+ ,39
+ ,153
+ ,140
+ ,143592
+ ,24634
+ ,113963
+ ,158
+ ,156
+ ,3009
+ ,258166
+ ,159
+ ,894
+ ,105
+ ,10
+ ,85
+ ,51
+ ,195
+ ,169
+ ,89626
+ ,27269
+ ,134904
+ ,138
+ ,133
+ ,2562
+ ,264771
+ ,88
+ ,917
+ ,84
+ ,2
+ ,154
+ ,40
+ ,149
+ ,129
+ ,131072
+ ,25270
+ ,114268
+ ,200
+ ,199
+ ,1737
+ ,182961
+ ,73
+ ,637
+ ,63
+ ,5
+ ,143
+ ,29
+ ,106
+ ,92
+ ,126817
+ ,24634
+ ,94333
+ ,104
+ ,98
+ ,2680
+ ,256967
+ ,74
+ ,900
+ ,67
+ ,13
+ ,107
+ ,47
+ ,179
+ ,160
+ ,81351
+ ,17828
+ ,102204
+ ,111
+ ,109
+ ,893
+ ,73566
+ ,32
+ ,385
+ ,39
+ ,6
+ ,22
+ ,23
+ ,88
+ ,67
+ ,22618
+ ,3007
+ ,23824
+ ,26
+ ,25
+ ,2389
+ ,272362
+ ,93
+ ,784
+ ,60
+ ,7
+ ,85
+ ,48
+ ,185
+ ,179
+ ,88977
+ ,20065
+ ,111563
+ ,115
+ ,113
+ ,2197
+ ,229056
+ ,62
+ ,910
+ ,94
+ ,2
+ ,101
+ ,38
+ ,133
+ ,90
+ ,92059
+ ,24648
+ ,91313
+ ,127
+ ,126
+ ,2227
+ ,229851
+ ,70
+ ,781
+ ,67
+ ,5
+ ,131
+ ,42
+ ,164
+ ,144
+ ,81897
+ ,21588
+ ,89770
+ ,140
+ ,137
+ ,2370
+ ,371391
+ ,91
+ ,1001
+ ,96
+ ,4
+ ,140
+ ,46
+ ,169
+ ,144
+ ,108146
+ ,25217
+ ,100125
+ ,121
+ ,121
+ ,3226
+ ,398210
+ ,104
+ ,1265
+ ,54
+ ,3
+ ,156
+ ,40
+ ,153
+ ,144
+ ,126372
+ ,30927
+ ,165278
+ ,183
+ ,178
+ ,1978
+ ,220419
+ ,111
+ ,587
+ ,54
+ ,6
+ ,81
+ ,45
+ ,166
+ ,134
+ ,249771
+ ,18487
+ ,181712
+ ,68
+ ,63
+ ,2516
+ ,231884
+ ,72
+ ,767
+ ,62
+ ,2
+ ,137
+ ,42
+ ,164
+ ,146
+ ,71154
+ ,18050
+ ,80906
+ ,112
+ ,109
+ ,2147
+ ,219381
+ ,73
+ ,746
+ ,71
+ ,0
+ ,102
+ ,41
+ ,146
+ ,121
+ ,71571
+ ,17696
+ ,75881
+ ,103
+ ,101
+ ,2150
+ ,206169
+ ,54
+ ,795
+ ,50
+ ,1
+ ,74
+ ,37
+ ,141
+ ,112
+ ,55918
+ ,17326
+ ,83963
+ ,63
+ ,61
+ ,4229
+ ,483074
+ ,132
+ ,1272
+ ,117
+ ,1
+ ,161
+ ,47
+ ,183
+ ,145
+ ,160141
+ ,39361
+ ,175721
+ ,166
+ ,157
+ ,1380
+ ,146100
+ ,72
+ ,657
+ ,45
+ ,5
+ ,30
+ ,26
+ ,99
+ ,99
+ ,38692
+ ,9648
+ ,68580
+ ,38
+ ,38
+ ,2449
+ ,295224
+ ,109
+ ,703
+ ,61
+ ,2
+ ,120
+ ,48
+ ,134
+ ,96
+ ,102812
+ ,26759
+ ,136323
+ ,163
+ ,159
+ ,870
+ ,80953
+ ,25
+ ,437
+ ,31
+ ,0
+ ,49
+ ,8
+ ,28
+ ,27
+ ,56622
+ ,7905
+ ,55792
+ ,59
+ ,58
+ ,2700
+ ,217384
+ ,63
+ ,1060
+ ,175
+ ,0
+ ,121
+ ,27
+ ,101
+ ,77
+ ,15986
+ ,4527
+ ,25157
+ ,27
+ ,27
+ ,1574
+ ,179344
+ ,62
+ ,459
+ ,70
+ ,6
+ ,76
+ ,38
+ ,139
+ ,137
+ ,123534
+ ,41517
+ ,100922
+ ,108
+ ,108
+ ,4046
+ ,415550
+ ,222
+ ,1586
+ ,284
+ ,1
+ ,85
+ ,41
+ ,159
+ ,151
+ ,108535
+ ,21261
+ ,118845
+ ,88
+ ,83
+ ,3259
+ ,389059
+ ,129
+ ,1084
+ ,95
+ ,4
+ ,151
+ ,61
+ ,222
+ ,126
+ ,93879
+ ,36099
+ ,170492
+ ,92
+ ,88
+ ,3098
+ ,180679
+ ,106
+ ,1051
+ ,72
+ ,1
+ ,165
+ ,45
+ ,171
+ ,159
+ ,144551
+ ,39039
+ ,81716
+ ,170
+ ,164
+ ,2615
+ ,299505
+ ,104
+ ,846
+ ,63
+ ,1
+ ,89
+ ,41
+ ,154
+ ,101
+ ,56750
+ ,13841
+ ,115750
+ ,98
+ ,96
+ ,2404
+ ,292260
+ ,84
+ ,732
+ ,75
+ ,3
+ ,168
+ ,42
+ ,154
+ ,144
+ ,127654
+ ,23841
+ ,105590
+ ,205
+ ,192
+ ,1932
+ ,199481
+ ,68
+ ,632
+ ,90
+ ,10
+ ,48
+ ,35
+ ,129
+ ,102
+ ,65594
+ ,8589
+ ,92795
+ ,96
+ ,94
+ ,3147
+ ,282361
+ ,78
+ ,1128
+ ,89
+ ,1
+ ,149
+ ,36
+ ,140
+ ,135
+ ,59938
+ ,15049
+ ,82390
+ ,107
+ ,107
+ ,2598
+ ,329281
+ ,89
+ ,971
+ ,138
+ ,4
+ ,75
+ ,40
+ ,156
+ ,147
+ ,146975
+ ,39038
+ ,135599
+ ,150
+ ,144
+ ,2108
+ ,234577
+ ,48
+ ,711
+ ,68
+ ,5
+ ,107
+ ,40
+ ,156
+ ,155
+ ,165904
+ ,36774
+ ,127667
+ ,138
+ ,136
+ ,2193
+ ,297995
+ ,67
+ ,738
+ ,80
+ ,7
+ ,116
+ ,38
+ ,138
+ ,138
+ ,169265
+ ,40076
+ ,163073
+ ,177
+ ,171
+ ,2478
+ ,342490
+ ,90
+ ,898
+ ,65
+ ,0
+ ,181
+ ,43
+ ,153
+ ,113
+ ,183500
+ ,43840
+ ,211381
+ ,213
+ ,210
+ ,4198
+ ,416463
+ ,163
+ ,1369
+ ,130
+ ,12
+ ,155
+ ,65
+ ,251
+ ,248
+ ,165986
+ ,43146
+ ,189944
+ ,208
+ ,193
+ ,4165
+ ,429565
+ ,120
+ ,1538
+ ,85
+ ,13
+ ,165
+ ,33
+ ,126
+ ,116
+ ,184923
+ ,50099
+ ,226168
+ ,307
+ ,297
+ ,2842
+ ,297080
+ ,142
+ ,893
+ ,83
+ ,9
+ ,121
+ ,51
+ ,198
+ ,176
+ ,140358
+ ,40312
+ ,117495
+ ,125
+ ,125
+ ,2562
+ ,331792
+ ,71
+ ,926
+ ,89
+ ,0
+ ,176
+ ,45
+ ,168
+ ,140
+ ,149959
+ ,32616
+ ,195894
+ ,208
+ ,204
+ ,2449
+ ,229772
+ ,202
+ ,800
+ ,116
+ ,0
+ ,86
+ ,36
+ ,138
+ ,59
+ ,57224
+ ,11338
+ ,80684
+ ,73
+ ,70
+ ,602
+ ,43287
+ ,14
+ ,214
+ ,43
+ ,4
+ ,13
+ ,19
+ ,71
+ ,64
+ ,43750
+ ,7409
+ ,19630
+ ,49
+ ,49
+ ,2579
+ ,238089
+ ,87
+ ,833
+ ,87
+ ,4
+ ,120
+ ,25
+ ,90
+ ,40
+ ,48029
+ ,18213
+ ,88634
+ ,82
+ ,82
+ ,2591
+ ,263322
+ ,160
+ ,906
+ ,80
+ ,0
+ ,117
+ ,44
+ ,167
+ ,98
+ ,104978
+ ,45873
+ ,139292
+ ,206
+ ,205
+ ,2957
+ ,302082
+ ,61
+ ,1288
+ ,132
+ ,0
+ ,133
+ ,45
+ ,172
+ ,139
+ ,100046
+ ,39844
+ ,128602
+ ,112
+ ,111
+ ,2786
+ ,321797
+ ,95
+ ,1079
+ ,59
+ ,0
+ ,169
+ ,44
+ ,162
+ ,135
+ ,101047
+ ,28317
+ ,135848
+ ,139
+ ,135
+ ,1477
+ ,193926
+ ,96
+ ,490
+ ,50
+ ,0
+ ,39
+ ,35
+ ,129
+ ,97
+ ,197426
+ ,24797
+ ,178377
+ ,60
+ ,59
+ ,3350
+ ,175138
+ ,105
+ ,990
+ ,87
+ ,0
+ ,125
+ ,46
+ ,179
+ ,142
+ ,160902
+ ,7471
+ ,106330
+ ,70
+ ,70
+ ,2107
+ ,354041
+ ,78
+ ,677
+ ,62
+ ,5
+ ,82
+ ,44
+ ,163
+ ,155
+ ,147172
+ ,27259
+ ,178303
+ ,112
+ ,108
+ ,2332
+ ,303273
+ ,91
+ ,696
+ ,70
+ ,1
+ ,148
+ ,45
+ ,164
+ ,115
+ ,109432
+ ,23201
+ ,116938
+ ,142
+ ,141
+ ,400
+ ,23668
+ ,13
+ ,156
+ ,9
+ ,0
+ ,12
+ ,1
+ ,0
+ ,0
+ ,1168
+ ,238
+ ,5841
+ ,11
+ ,11
+ ,2233
+ ,196743
+ ,79
+ ,785
+ ,54
+ ,0
+ ,146
+ ,40
+ ,155
+ ,103
+ ,83248
+ ,28830
+ ,106020
+ ,130
+ ,130
+ ,530
+ ,61857
+ ,25
+ ,192
+ ,25
+ ,4
+ ,23
+ ,11
+ ,32
+ ,30
+ ,25162
+ ,3913
+ ,24610
+ ,31
+ ,28
+ ,2033
+ ,217543
+ ,54
+ ,641
+ ,113
+ ,0
+ ,87
+ ,51
+ ,189
+ ,130
+ ,45724
+ ,9935
+ ,74151
+ ,132
+ ,101
+ ,3246
+ ,440711
+ ,128
+ ,1251
+ ,63
+ ,1
+ ,164
+ ,38
+ ,140
+ ,102
+ ,110529
+ ,27738
+ ,232241
+ ,219
+ ,216
+ ,387
+ ,21054
+ ,16
+ ,146
+ ,2
+ ,0
+ ,4
+ ,0
+ ,0
+ ,0
+ ,855
+ ,338
+ ,6622
+ ,4
+ ,4
+ ,2137
+ ,252805
+ ,52
+ ,866
+ ,67
+ ,5
+ ,81
+ ,30
+ ,111
+ ,77
+ ,101382
+ ,13326
+ ,127097
+ ,102
+ ,97
+ ,492
+ ,31961
+ ,22
+ ,200
+ ,22
+ ,0
+ ,18
+ ,8
+ ,25
+ ,9
+ ,14116
+ ,3988
+ ,13155
+ ,39
+ ,39
+ ,3838
+ ,360436
+ ,125
+ ,1351
+ ,157
+ ,3
+ ,118
+ ,43
+ ,159
+ ,150
+ ,89506
+ ,24347
+ ,160501
+ ,125
+ ,119
+ ,2193
+ ,251948
+ ,77
+ ,740
+ ,79
+ ,7
+ ,76
+ ,48
+ ,183
+ ,163
+ ,135356
+ ,27111
+ ,91502
+ ,121
+ ,118
+ ,1796
+ ,187320
+ ,97
+ ,524
+ ,113
+ ,14
+ ,55
+ ,49
+ ,184
+ ,148
+ ,116066
+ ,3938
+ ,24469
+ ,42
+ ,41
+ ,1907
+ ,180842
+ ,58
+ ,724
+ ,50
+ ,3
+ ,62
+ ,32
+ ,119
+ ,94
+ ,144244
+ ,17416
+ ,88229
+ ,111
+ ,107
+ ,568
+ ,38214
+ ,34
+ ,276
+ ,52
+ ,0
+ ,16
+ ,8
+ ,27
+ ,21
+ ,8773
+ ,1888
+ ,13983
+ ,16
+ ,16
+ ,2602
+ ,280392
+ ,56
+ ,862
+ ,113
+ ,3
+ ,98
+ ,43
+ ,163
+ ,151
+ ,102153
+ ,18700
+ ,80716
+ ,70
+ ,69
+ ,2819
+ ,358276
+ ,84
+ ,1031
+ ,115
+ ,0
+ ,137
+ ,52
+ ,198
+ ,187
+ ,117440
+ ,36809
+ ,157384
+ ,162
+ ,160
+ ,1464
+ ,211775
+ ,67
+ ,511
+ ,78
+ ,0
+ ,50
+ ,53
+ ,205
+ ,171
+ ,104128
+ ,24959
+ ,122975
+ ,173
+ ,158
+ ,3946
+ ,447335
+ ,90
+ ,1716
+ ,135
+ ,4
+ ,152
+ ,49
+ ,191
+ ,170
+ ,134238
+ ,37343
+ ,191469
+ ,171
+ ,161
+ ,2554
+ ,348017
+ ,99
+ ,884
+ ,120
+ ,0
+ ,163
+ ,48
+ ,187
+ ,145
+ ,134047
+ ,21849
+ ,231257
+ ,172
+ ,165
+ ,3506
+ ,441946
+ ,133
+ ,1201
+ ,122
+ ,3
+ ,142
+ ,56
+ ,210
+ ,198
+ ,279488
+ ,49809
+ ,258287
+ ,254
+ ,246
+ ,1552
+ ,215177
+ ,43
+ ,575
+ ,54
+ ,0
+ ,80
+ ,45
+ ,166
+ ,152
+ ,79756
+ ,21654
+ ,122531
+ ,90
+ ,89
+ ,1389
+ ,130177
+ ,47
+ ,481
+ ,63
+ ,0
+ ,59
+ ,40
+ ,145
+ ,112
+ ,66089
+ ,8728
+ ,61394
+ ,50
+ ,49
+ ,3101
+ ,318037
+ ,365
+ ,1031
+ ,162
+ ,4
+ ,94
+ ,48
+ ,187
+ ,173
+ ,102070
+ ,20920
+ ,86480
+ ,113
+ ,107
+ ,4541
+ ,466139
+ ,198
+ ,1574
+ ,162
+ ,5
+ ,128
+ ,50
+ ,186
+ ,177
+ ,146760
+ ,27195
+ ,195791
+ ,187
+ ,182
+ ,1872
+ ,162279
+ ,62
+ ,575
+ ,107
+ ,16
+ ,63
+ ,43
+ ,164
+ ,153
+ ,154771
+ ,1037
+ ,18284
+ ,16
+ ,16
+ ,4403
+ ,416643
+ ,140
+ ,1827
+ ,146
+ ,6
+ ,127
+ ,46
+ ,172
+ ,161
+ ,165933
+ ,42570
+ ,147581
+ ,175
+ ,173
+ ,2113
+ ,178322
+ ,86
+ ,790
+ ,77
+ ,5
+ ,60
+ ,40
+ ,147
+ ,115
+ ,64593
+ ,17672
+ ,72558
+ ,90
+ ,90
+ ,2046
+ ,292443
+ ,54
+ ,668
+ ,87
+ ,2
+ ,118
+ ,45
+ ,167
+ ,147
+ ,92280
+ ,34245
+ ,147341
+ ,140
+ ,140
+ ,2564
+ ,283913
+ ,100
+ ,905
+ ,192
+ ,1
+ ,110
+ ,46
+ ,158
+ ,124
+ ,67150
+ ,16786
+ ,114651
+ ,145
+ ,142
+ ,2073
+ ,244931
+ ,127
+ ,689
+ ,75
+ ,2
+ ,46
+ ,37
+ ,144
+ ,57
+ ,128692
+ ,20954
+ ,100187
+ ,141
+ ,126
+ ,4112
+ ,387072
+ ,125
+ ,1613
+ ,131
+ ,9
+ ,96
+ ,45
+ ,169
+ ,144
+ ,124089
+ ,16378
+ ,130332
+ ,125
+ ,123
+ ,2340
+ ,246963
+ ,93
+ ,811
+ ,67
+ ,1
+ ,128
+ ,39
+ ,145
+ ,126
+ ,125386
+ ,31852
+ ,134218
+ ,241
+ ,239
+ ,2035
+ ,173260
+ ,63
+ ,716
+ ,37
+ ,3
+ ,41
+ ,21
+ ,79
+ ,78
+ ,37238
+ ,2805
+ ,10901
+ ,16
+ ,15
+ ,3241
+ ,346748
+ ,108
+ ,1034
+ ,61
+ ,11
+ ,146
+ ,50
+ ,194
+ ,153
+ ,140015
+ ,38086
+ ,145758
+ ,175
+ ,170
+ ,1991
+ ,178402
+ ,60
+ ,739
+ ,127
+ ,5
+ ,147
+ ,55
+ ,212
+ ,196
+ ,150047
+ ,21166
+ ,75767
+ ,132
+ ,123
+ ,2828
+ ,268750
+ ,96
+ ,1086
+ ,58
+ ,2
+ ,121
+ ,40
+ ,148
+ ,130
+ ,154451
+ ,34672
+ ,134969
+ ,154
+ ,151
+ ,2748
+ ,314070
+ ,112
+ ,852
+ ,71
+ ,1
+ ,185
+ ,48
+ ,171
+ ,159
+ ,156349
+ ,36171
+ ,169216
+ ,198
+ ,194
+ ,2
+ ,1
+ ,0
+ ,0
+ ,0
+ ,9
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,207
+ ,14688
+ ,10
+ ,85
+ ,0
+ ,0
+ ,4
+ ,0
+ ,0
+ ,0
+ ,6023
+ ,2065
+ ,7953
+ ,5
+ ,5
+ ,5
+ ,98
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,8
+ ,455
+ ,2
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,2449
+ ,291847
+ ,95
+ ,816
+ ,72
+ ,2
+ ,85
+ ,46
+ ,141
+ ,94
+ ,84601
+ ,19354
+ ,105406
+ ,125
+ ,122
+ ,3490
+ ,415421
+ ,168
+ ,1142
+ ,123
+ ,3
+ ,164
+ ,52
+ ,204
+ ,129
+ ,68946
+ ,22124
+ ,174586
+ ,174
+ ,173
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,4
+ ,203
+ ,4
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,151
+ ,7199
+ ,5
+ ,74
+ ,0
+ ,0
+ ,7
+ ,0
+ ,0
+ ,0
+ ,1644
+ ,556
+ ,4245
+ ,6
+ ,6
+ ,475
+ ,46660
+ ,21
+ ,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
+ ,1145
+ ,121550
+ ,46
+ ,309
+ ,106
+ ,0
+ ,37
+ ,48
+ ,172
+ ,89
+ ,52789
+ ,1813
+ ,15673
+ ,35
+ ,35
+ ,29
+ ,969
+ ,2
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,2080
+ ,242774
+ ,75
+ ,695
+ ,53
+ ,2
+ ,62
+ ,34
+ ,125
+ ,71
+ ,100350
+ ,17372
+ ,75882
+ ,80
+ ,72)
+ ,dim=c(15
+ ,164)
+ ,dimnames=list(c('#Pageviews'
+ ,'TotalTimeRFC'
+ ,'#Logins'
+ ,'
#CompendiumViews'
+ ,'#CompendiumViews(PR)'
+ ,'#SharedCOmp'
+ ,'#BloggedComputations'
+ ,'#ReviewedCompendiums'
+ ,'#submittedFeedbackMsgs'
+ ,'#LongPeerRev'
+ ,'#characters'
+ ,'#revisions'
+ ,'
#seconds(comp)'
+ ,'#hyperlinks'
+ ,'#blogs')
+ ,1:164))
> y <- array(NA,dim=c(15,164),dimnames=list(c('#Pageviews','TotalTimeRFC','#Logins','
#CompendiumViews','#CompendiumViews(PR)','#SharedCOmp','#BloggedComputations','#ReviewedCompendiums','#submittedFeedbackMsgs','#LongPeerRev','#characters','#revisions','
#seconds(comp)','#hyperlinks','#blogs'),1:164))
> 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 = '11'
> 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
#characters #Pageviews TotalTimeRFC #Logins \r#CompendiumViews
1 140824 1818 279055 73 504
2 110459 1433 212408 75 510
3 105079 2059 233939 83 710
4 112098 2733 222117 106 1154
5 43929 1399 189911 56 415
6 76173 631 70849 28 179
7 187326 5460 605767 135 2563
8 22807 381 33186 19 111
9 144408 2150 227332 62 763
10 66485 2042 267925 49 661
11 79089 2536 371987 122 981
12 81625 2377 264989 131 733
13 68788 2100 212638 87 785
14 103297 3020 368577 85 1186
15 69446 2265 269455 88 724
16 114948 5139 398124 191 1774
17 167949 2363 335567 77 845
18 125081 3564 432711 173 1390
19 125818 1477 182016 58 514
20 136588 2398 267365 89 692
21 112431 2546 279428 73 847
22 103037 3150 508849 111 1397
23 82317 1705 220142 49 569
24 118906 1787 200004 58 636
25 83515 3792 257139 133 1370
26 104581 3108 270941 138 1092
27 103129 3230 324969 134 1201
28 83243 2348 329962 92 763
29 37110 1780 190867 60 652
30 113344 3218 393860 79 1213
31 139165 2692 327660 89 1111
32 86652 2187 269239 83 758
33 112302 2577 396136 106 906
34 69652 1293 130446 49 456
35 119442 3567 430118 104 1293
36 69867 2764 273950 56 1186
37 101629 3755 428077 128 1348
38 70168 2075 254312 93 695
39 31081 995 120351 35 306
40 103925 3750 395658 212 1319
41 92622 3413 345875 86 1566
42 79011 2053 216827 82 784
43 93487 1984 224524 83 730
44 64520 1825 182485 69 488
45 93473 2599 157164 85 1051
46 114360 5572 459455 157 2089
47 33032 918 78800 42 330
48 96125 2685 255072 85 764
49 151911 4145 368086 123 1410
50 89256 2841 230299 70 1187
51 95676 2175 244782 81 691
52 5950 496 24188 24 218
53 149695 2699 400109 334 865
54 32551 744 65029 17 255
55 31701 1161 101097 64 454
56 100087 3333 309810 67 1229
57 169707 2970 375638 91 790
58 150491 3968 367127 204 1208
59 120192 2878 381998 155 1102
60 95893 2399 280106 90 919
61 151715 4121 400971 153 1352
62 176225 3294 315924 122 1190
63 59900 3132 291391 124 1257
64 104767 2868 295075 93 1030
65 114799 1778 280018 81 669
66 72128 2109 267432 71 542
67 143592 2148 217181 141 652
68 89626 3009 258166 159 894
69 131072 2562 264771 88 917
70 126817 1737 182961 73 637
71 81351 2680 256967 74 900
72 22618 893 73566 32 385
73 88977 2389 272362 93 784
74 92059 2197 229056 62 910
75 81897 2227 229851 70 781
76 108146 2370 371391 91 1001
77 126372 3226 398210 104 1265
78 249771 1978 220419 111 587
79 71154 2516 231884 72 767
80 71571 2147 219381 73 746
81 55918 2150 206169 54 795
82 160141 4229 483074 132 1272
83 38692 1380 146100 72 657
84 102812 2449 295224 109 703
85 56622 870 80953 25 437
86 15986 2700 217384 63 1060
87 123534 1574 179344 62 459
88 108535 4046 415550 222 1586
89 93879 3259 389059 129 1084
90 144551 3098 180679 106 1051
91 56750 2615 299505 104 846
92 127654 2404 292260 84 732
93 65594 1932 199481 68 632
94 59938 3147 282361 78 1128
95 146975 2598 329281 89 971
96 165904 2108 234577 48 711
97 169265 2193 297995 67 738
98 183500 2478 342490 90 898
99 165986 4198 416463 163 1369
100 184923 4165 429565 120 1538
101 140358 2842 297080 142 893
102 149959 2562 331792 71 926
103 57224 2449 229772 202 800
104 43750 602 43287 14 214
105 48029 2579 238089 87 833
106 104978 2591 263322 160 906
107 100046 2957 302082 61 1288
108 101047 2786 321797 95 1079
109 197426 1477 193926 96 490
110 160902 3350 175138 105 990
111 147172 2107 354041 78 677
112 109432 2332 303273 91 696
113 1168 400 23668 13 156
114 83248 2233 196743 79 785
115 25162 530 61857 25 192
116 45724 2033 217543 54 641
117 110529 3246 440711 128 1251
118 855 387 21054 16 146
119 101382 2137 252805 52 866
120 14116 492 31961 22 200
121 89506 3838 360436 125 1351
122 135356 2193 251948 77 740
123 116066 1796 187320 97 524
124 144244 1907 180842 58 724
125 8773 568 38214 34 276
126 102153 2602 280392 56 862
127 117440 2819 358276 84 1031
128 104128 1464 211775 67 511
129 134238 3946 447335 90 1716
130 134047 2554 348017 99 884
131 279488 3506 441946 133 1201
132 79756 1552 215177 43 575
133 66089 1389 130177 47 481
134 102070 3101 318037 365 1031
135 146760 4541 466139 198 1574
136 154771 1872 162279 62 575
137 165933 4403 416643 140 1827
138 64593 2113 178322 86 790
139 92280 2046 292443 54 668
140 67150 2564 283913 100 905
141 128692 2073 244931 127 689
142 124089 4112 387072 125 1613
143 125386 2340 246963 93 811
144 37238 2035 173260 63 716
145 140015 3241 346748 108 1034
146 150047 1991 178402 60 739
147 154451 2828 268750 96 1086
148 156349 2748 314070 112 852
149 0 2 1 0 0
150 6023 207 14688 10 85
151 0 5 98 1 0
152 0 8 455 2 0
153 0 0 0 0 0
154 0 0 0 0 0
155 84601 2449 291847 95 816
156 68946 3490 415421 168 1142
157 0 0 0 0 0
158 0 4 203 4 0
159 1644 151 7199 5 74
160 6179 475 46660 21 259
161 3926 141 17547 5 69
162 52789 1145 121550 46 309
163 0 29 969 2 0
164 100350 2080 242774 75 695
#CompendiumViews(PR) #SharedCOmp #BloggedComputations #ReviewedCompendiums
1 95 3 96 42
2 68 4 75 38
3 64 16 70 46
4 139 2 134 42
5 51 1 83 30
6 46 3 8 35
7 118 0 173 40
8 46 0 1 18
9 79 7 88 38
10 76 0 104 37
11 82 0 114 46
12 66 7 125 60
13 60 10 57 37
14 117 4 139 55
15 50 10 87 44
16 133 0 176 63
17 63 8 114 40
18 100 4 121 43
19 44 3 103 32
20 65 8 135 52
21 103 0 123 49
22 103 1 99 41
23 62 5 77 25
24 70 9 103 57
25 159 1 158 45
26 78 0 116 42
27 101 5 114 45
28 73 0 150 43
29 58 0 64 36
30 147 0 150 45
31 54 3 143 50
32 84 6 50 50
33 56 1 145 51
34 45 4 56 42
35 87 4 141 44
36 87 0 83 42
37 77 0 112 44
38 72 2 79 40
39 36 1 33 17
40 51 2 152 43
41 44 10 126 41
42 75 10 97 41
43 87 5 84 40
44 97 6 68 49
45 90 1 50 52
46 860 2 101 42
47 57 2 20 26
48 99 1 107 59
49 120 10 150 50
50 76 3 129 50
51 56 0 99 47
52 20 0 8 4
53 94 8 88 51
54 21 5 21 18
55 70 3 30 14
56 133 1 102 41
57 86 5 166 61
58 224 6 132 40
59 65 0 161 44
60 86 12 90 40
61 70 10 160 51
62 148 12 139 29
63 72 11 104 43
64 59 8 103 42
65 67 3 66 41
66 58 0 163 30
67 60 6 93 39
68 105 10 85 51
69 84 2 154 40
70 63 5 143 29
71 67 13 107 47
72 39 6 22 23
73 60 7 85 48
74 94 2 101 38
75 67 5 131 42
76 96 4 140 46
77 54 3 156 40
78 54 6 81 45
79 62 2 137 42
80 71 0 102 41
81 50 1 74 37
82 117 1 161 47
83 45 5 30 26
84 61 2 120 48
85 31 0 49 8
86 175 0 121 27
87 70 6 76 38
88 284 1 85 41
89 95 4 151 61
90 72 1 165 45
91 63 1 89 41
92 75 3 168 42
93 90 10 48 35
94 89 1 149 36
95 138 4 75 40
96 68 5 107 40
97 80 7 116 38
98 65 0 181 43
99 130 12 155 65
100 85 13 165 33
101 83 9 121 51
102 89 0 176 45
103 116 0 86 36
104 43 4 13 19
105 87 4 120 25
106 80 0 117 44
107 132 0 133 45
108 59 0 169 44
109 50 0 39 35
110 87 0 125 46
111 62 5 82 44
112 70 1 148 45
113 9 0 12 1
114 54 0 146 40
115 25 4 23 11
116 113 0 87 51
117 63 1 164 38
118 2 0 4 0
119 67 5 81 30
120 22 0 18 8
121 157 3 118 43
122 79 7 76 48
123 113 14 55 49
124 50 3 62 32
125 52 0 16 8
126 113 3 98 43
127 115 0 137 52
128 78 0 50 53
129 135 4 152 49
130 120 0 163 48
131 122 3 142 56
132 54 0 80 45
133 63 0 59 40
134 162 4 94 48
135 162 5 128 50
136 107 16 63 43
137 146 6 127 46
138 77 5 60 40
139 87 2 118 45
140 192 1 110 46
141 75 2 46 37
142 131 9 96 45
143 67 1 128 39
144 37 3 41 21
145 61 11 146 50
146 127 5 147 55
147 58 2 121 40
148 71 1 185 48
149 0 9 0 0
150 0 0 4 0
151 0 0 0 0
152 0 0 0 0
153 0 1 0 0
154 0 0 0 0
155 72 2 85 46
156 123 3 164 52
157 0 0 0 0
158 0 0 0 0
159 0 0 7 0
160 7 0 12 5
161 3 0 0 1
162 106 0 37 48
163 0 0 0 0
164 53 2 62 34
#submittedFeedbackMsgs #LongPeerRev #revisions \r#seconds(comp) #hyperlinks
1 159 130 32033 186099 165
2 149 143 20654 113854 135
3 178 118 16346 99776 121
4 164 146 35926 106194 148
5 100 73 10621 100792 73
6 129 89 10024 47552 49
7 156 146 43068 250931 185
8 67 22 1271 6853 5
9 148 132 34416 115466 125
10 132 92 20318 110896 93
11 169 147 24409 169351 154
12 230 203 20648 94853 98
13 122 113 12347 72591 70
14 191 171 21857 101345 148
15 162 87 11034 113713 100
16 237 208 33433 165354 150
17 156 153 35902 164263 197
18 157 97 22355 135213 114
19 123 95 31219 111669 169
20 203 197 21983 134163 200
21 187 160 40085 140303 148
22 152 148 18507 150773 140
23 89 84 16278 111848 74
24 227 227 24662 102509 128
25 165 154 31452 96785 140
26 162 151 32580 116136 116
27 174 142 22883 158376 147
28 154 148 27652 153990 132
29 129 110 9845 64057 70
30 174 149 20190 230054 144
31 195 179 46201 184531 155
32 186 149 10971 114198 165
33 197 187 34811 198299 161
34 157 153 3029 33750 31
35 168 163 38941 189723 199
36 159 127 4958 100826 78
37 161 151 32344 188355 121
38 153 100 19433 104470 112
39 55 46 12558 58391 41
40 166 156 36524 164808 158
41 151 128 26041 134097 123
42 148 111 16637 80238 104
43 129 119 28395 133252 94
44 181 148 16747 54518 73
45 93 65 9105 121850 52
46 150 134 11941 79367 71
47 82 66 7935 56968 21
48 229 201 19499 106314 155
49 193 177 22938 191889 174
50 176 156 25314 104864 136
51 179 158 28527 160792 128
52 12 7 2694 15049 7
53 181 175 20867 191179 165
54 67 61 3597 25109 21
55 52 41 5296 45824 35
56 148 133 32982 129711 137
57 230 228 38975 210012 174
58 148 140 42721 194679 257
59 160 155 41455 197680 207
60 155 141 23923 81180 103
61 198 181 26719 197765 171
62 104 75 53405 214738 279
63 169 97 12526 96252 83
64 163 142 26584 124527 130
65 151 136 37062 153242 131
66 116 87 25696 145707 126
67 153 140 24634 113963 158
68 195 169 27269 134904 138
69 149 129 25270 114268 200
70 106 92 24634 94333 104
71 179 160 17828 102204 111
72 88 67 3007 23824 26
73 185 179 20065 111563 115
74 133 90 24648 91313 127
75 164 144 21588 89770 140
76 169 144 25217 100125 121
77 153 144 30927 165278 183
78 166 134 18487 181712 68
79 164 146 18050 80906 112
80 146 121 17696 75881 103
81 141 112 17326 83963 63
82 183 145 39361 175721 166
83 99 99 9648 68580 38
84 134 96 26759 136323 163
85 28 27 7905 55792 59
86 101 77 4527 25157 27
87 139 137 41517 100922 108
88 159 151 21261 118845 88
89 222 126 36099 170492 92
90 171 159 39039 81716 170
91 154 101 13841 115750 98
92 154 144 23841 105590 205
93 129 102 8589 92795 96
94 140 135 15049 82390 107
95 156 147 39038 135599 150
96 156 155 36774 127667 138
97 138 138 40076 163073 177
98 153 113 43840 211381 213
99 251 248 43146 189944 208
100 126 116 50099 226168 307
101 198 176 40312 117495 125
102 168 140 32616 195894 208
103 138 59 11338 80684 73
104 71 64 7409 19630 49
105 90 40 18213 88634 82
106 167 98 45873 139292 206
107 172 139 39844 128602 112
108 162 135 28317 135848 139
109 129 97 24797 178377 60
110 179 142 7471 106330 70
111 163 155 27259 178303 112
112 164 115 23201 116938 142
113 0 0 238 5841 11
114 155 103 28830 106020 130
115 32 30 3913 24610 31
116 189 130 9935 74151 132
117 140 102 27738 232241 219
118 0 0 338 6622 4
119 111 77 13326 127097 102
120 25 9 3988 13155 39
121 159 150 24347 160501 125
122 183 163 27111 91502 121
123 184 148 3938 24469 42
124 119 94 17416 88229 111
125 27 21 1888 13983 16
126 163 151 18700 80716 70
127 198 187 36809 157384 162
128 205 171 24959 122975 173
129 191 170 37343 191469 171
130 187 145 21849 231257 172
131 210 198 49809 258287 254
132 166 152 21654 122531 90
133 145 112 8728 61394 50
134 187 173 20920 86480 113
135 186 177 27195 195791 187
136 164 153 1037 18284 16
137 172 161 42570 147581 175
138 147 115 17672 72558 90
139 167 147 34245 147341 140
140 158 124 16786 114651 145
141 144 57 20954 100187 141
142 169 144 16378 130332 125
143 145 126 31852 134218 241
144 79 78 2805 10901 16
145 194 153 38086 145758 175
146 212 196 21166 75767 132
147 148 130 34672 134969 154
148 171 159 36171 169216 198
149 0 0 0 0 0
150 0 0 2065 7953 5
151 0 0 0 0 0
152 0 0 0 0 0
153 0 0 0 0 0
154 0 0 0 0 0
155 141 94 19354 105406 125
156 204 129 22124 174586 174
157 0 0 0 0 0
158 0 0 0 0 0
159 0 0 556 4245 6
160 15 13 2089 21509 13
161 4 4 2658 7670 3
162 172 89 1813 15673 35
163 0 0 0 0 0
164 125 71 17372 75882 80
#blogs
1 165
2 132
3 121
4 145
5 71
6 47
7 177
8 5
9 124
10 92
11 149
12 93
13 70
14 148
15 100
16 142
17 194
18 113
19 162
20 186
21 147
22 137
23 71
24 123
25 134
26 115
27 138
28 125
29 66
30 137
31 152
32 159
33 159
34 31
35 185
36 78
37 117
38 109
39 41
40 149
41 123
42 103
43 87
44 71
45 51
46 70
47 21
48 155
49 172
50 133
51 125
52 7
53 158
54 21
55 35
56 133
57 169
58 256
59 190
60 100
61 171
62 267
63 80
64 126
65 132
66 121
67 156
68 133
69 199
70 98
71 109
72 25
73 113
74 126
75 137
76 121
77 178
78 63
79 109
80 101
81 61
82 157
83 38
84 159
85 58
86 27
87 108
88 83
89 88
90 164
91 96
92 192
93 94
94 107
95 144
96 136
97 171
98 210
99 193
100 297
101 125
102 204
103 70
104 49
105 82
106 205
107 111
108 135
109 59
110 70
111 108
112 141
113 11
114 130
115 28
116 101
117 216
118 4
119 97
120 39
121 119
122 118
123 41
124 107
125 16
126 69
127 160
128 158
129 161
130 165
131 246
132 89
133 49
134 107
135 182
136 16
137 173
138 90
139 140
140 142
141 126
142 123
143 239
144 15
145 170
146 123
147 151
148 194
149 0
150 5
151 0
152 0
153 0
154 0
155 122
156 173
157 0
158 0
159 6
160 13
161 3
162 35
163 0
164 72
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) `#Pageviews` TotalTimeRFC
2007.5482 1.3033 -0.1346
`#Logins` `\r#CompendiumViews` `#CompendiumViews(PR)`
11.4375 0.7662 83.7354
`#SharedCOmp` `#BloggedComputations` `#ReviewedCompendiums`
2215.9504 -57.5784 180.9099
`#submittedFeedbackMsgs` `#LongPeerRev` `#revisions`
1.7201 195.6183 1.0818
`\r#seconds(comp)` `#hyperlinks` `#blogs`
0.4894 46.2764 1.6596
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-54870 -17966 -3622 12340 113158
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 2007.54820 5720.49379 0.351 0.726129
`#Pageviews` 1.30332 11.09686 0.117 0.906662
TotalTimeRFC -0.13457 0.06023 -2.234 0.026951 *
`#Logins` 11.43751 63.04208 0.181 0.856279
`\r#CompendiumViews` 0.76622 22.75897 0.034 0.973188
`#CompendiumViews(PR)` 83.73545 39.07438 2.143 0.033739 *
`#SharedCOmp` 2215.95039 632.10599 3.506 0.000602 ***
`#BloggedComputations` -57.57844 101.18089 -0.569 0.570169
`#ReviewedCompendiums` 180.90987 789.68163 0.229 0.819111
`#submittedFeedbackMsgs` 1.72008 245.24958 0.007 0.994413
`#LongPeerRev` 195.61833 119.96372 1.631 0.105077
`#revisions` 1.08179 0.34483 3.137 0.002056 **
`\r#seconds(comp)` 0.48939 0.09437 5.186 6.9e-07 ***
`#hyperlinks` 46.27643 565.52359 0.082 0.934892
`#blogs` 1.65961 588.08108 0.003 0.997752
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 26870 on 149 degrees of freedom
Multiple R-squared: 0.7575, Adjusted R-squared: 0.7347
F-statistic: 33.25 on 14 and 149 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,] 4.035083e-01 8.070167e-01 5.964917e-01
[2,] 2.918978e-01 5.837956e-01 7.081022e-01
[3,] 2.053465e-01 4.106931e-01 7.946535e-01
[4,] 1.776715e-01 3.553429e-01 8.223285e-01
[5,] 1.819106e-01 3.638212e-01 8.180894e-01
[6,] 1.127139e-01 2.254277e-01 8.872861e-01
[7,] 7.569600e-02 1.513920e-01 9.243040e-01
[8,] 6.076837e-02 1.215367e-01 9.392316e-01
[9,] 3.490342e-02 6.980684e-02 9.650966e-01
[10,] 2.451188e-02 4.902376e-02 9.754881e-01
[11,] 1.440132e-02 2.880264e-02 9.855987e-01
[12,] 1.042873e-02 2.085746e-02 9.895713e-01
[13,] 8.459113e-03 1.691823e-02 9.915409e-01
[14,] 5.248684e-03 1.049737e-02 9.947513e-01
[15,] 3.955222e-03 7.910444e-03 9.960448e-01
[16,] 2.618767e-03 5.237534e-03 9.973812e-01
[17,] 5.493682e-03 1.098736e-02 9.945063e-01
[18,] 6.695414e-03 1.339083e-02 9.933046e-01
[19,] 4.136331e-03 8.272662e-03 9.958637e-01
[20,] 2.513750e-03 5.027500e-03 9.974862e-01
[21,] 1.704871e-03 3.409742e-03 9.982951e-01
[22,] 1.645629e-03 3.291259e-03 9.983544e-01
[23,] 1.034165e-03 2.068329e-03 9.989658e-01
[24,] 1.893750e-03 3.787499e-03 9.981063e-01
[25,] 1.130375e-03 2.260751e-03 9.988696e-01
[26,] 9.447444e-04 1.889489e-03 9.990553e-01
[27,] 6.338986e-04 1.267797e-03 9.993661e-01
[28,] 1.916879e-03 3.833758e-03 9.980831e-01
[29,] 1.165253e-03 2.330505e-03 9.988347e-01
[30,] 9.818672e-04 1.963734e-03 9.990181e-01
[31,] 5.927755e-04 1.185551e-03 9.994072e-01
[32,] 4.080640e-04 8.161280e-04 9.995919e-01
[33,] 3.401573e-04 6.803147e-04 9.996598e-01
[34,] 2.932781e-04 5.865562e-04 9.997067e-01
[35,] 1.931911e-04 3.863823e-04 9.998068e-01
[36,] 1.267269e-04 2.534537e-04 9.998733e-01
[37,] 7.138579e-05 1.427716e-04 9.999286e-01
[38,] 4.512553e-05 9.025106e-05 9.999549e-01
[39,] 2.880172e-05 5.760345e-05 9.999712e-01
[40,] 4.208158e-05 8.416317e-05 9.999579e-01
[41,] 5.167582e-05 1.033516e-04 9.999483e-01
[42,] 3.869586e-05 7.739173e-05 9.999613e-01
[43,] 2.241254e-05 4.482508e-05 9.999776e-01
[44,] 1.480506e-05 2.961013e-05 9.999852e-01
[45,] 8.948988e-06 1.789798e-05 9.999911e-01
[46,] 7.189854e-06 1.437971e-05 9.999928e-01
[47,] 4.094589e-06 8.189179e-06 9.999959e-01
[48,] 3.152853e-06 6.305706e-06 9.999968e-01
[49,] 2.091895e-06 4.183789e-06 9.999979e-01
[50,] 4.728956e-06 9.457911e-06 9.999953e-01
[51,] 2.762782e-05 5.525564e-05 9.999724e-01
[52,] 2.801046e-05 5.602092e-05 9.999720e-01
[53,] 1.005778e-04 2.011555e-04 9.998994e-01
[54,] 1.482880e-04 2.965759e-04 9.998517e-01
[55,] 1.226380e-04 2.452761e-04 9.998774e-01
[56,] 1.172913e-04 2.345827e-04 9.998827e-01
[57,] 7.406862e-05 1.481372e-04 9.999259e-01
[58,] 5.582646e-05 1.116529e-04 9.999442e-01
[59,] 5.832367e-05 1.166473e-04 9.999417e-01
[60,] 3.735419e-05 7.470839e-05 9.999626e-01
[61,] 1.432667e-01 2.865334e-01 8.567333e-01
[62,] 1.265579e-01 2.531157e-01 8.734421e-01
[63,] 1.053001e-01 2.106002e-01 8.946999e-01
[64,] 9.401342e-02 1.880268e-01 9.059866e-01
[65,] 1.636647e-01 3.273295e-01 8.363353e-01
[66,] 1.835500e-01 3.671001e-01 8.164500e-01
[67,] 1.547664e-01 3.095329e-01 8.452336e-01
[68,] 1.390123e-01 2.780245e-01 8.609877e-01
[69,] 1.167495e-01 2.334989e-01 8.832505e-01
[70,] 1.047155e-01 2.094310e-01 8.952845e-01
[71,] 1.017600e-01 2.035200e-01 8.982400e-01
[72,] 1.263713e-01 2.527426e-01 8.736287e-01
[73,] 1.424096e-01 2.848191e-01 8.575904e-01
[74,] 1.311511e-01 2.623021e-01 8.688489e-01
[75,] 1.613315e-01 3.226631e-01 8.386685e-01
[76,] 1.908061e-01 3.816122e-01 8.091939e-01
[77,] 1.605230e-01 3.210460e-01 8.394770e-01
[78,] 1.625296e-01 3.250592e-01 8.374704e-01
[79,] 1.804789e-01 3.609579e-01 8.195211e-01
[80,] 1.684351e-01 3.368703e-01 8.315649e-01
[81,] 1.961057e-01 3.922114e-01 8.038943e-01
[82,] 3.119578e-01 6.239156e-01 6.880422e-01
[83,] 2.766338e-01 5.532677e-01 7.233662e-01
[84,] 2.803195e-01 5.606391e-01 7.196805e-01
[85,] 2.679005e-01 5.358009e-01 7.320995e-01
[86,] 2.293976e-01 4.587952e-01 7.706024e-01
[87,] 2.002153e-01 4.004307e-01 7.997847e-01
[88,] 1.762169e-01 3.524338e-01 8.237831e-01
[89,] 1.934570e-01 3.869140e-01 8.065430e-01
[90,] 1.679263e-01 3.358526e-01 8.320737e-01
[91,] 1.415901e-01 2.831801e-01 8.584099e-01
[92,] 2.659712e-01 5.319425e-01 7.340288e-01
[93,] 4.236885e-01 8.473771e-01 5.763115e-01
[94,] 3.809812e-01 7.619623e-01 6.190188e-01
[95,] 3.885803e-01 7.771605e-01 6.114197e-01
[96,] 3.384486e-01 6.768972e-01 6.615514e-01
[97,] 3.061096e-01 6.122191e-01 6.938904e-01
[98,] 2.654736e-01 5.309471e-01 7.345264e-01
[99,] 3.639549e-01 7.279098e-01 6.360451e-01
[100,] 3.390958e-01 6.781917e-01 6.609042e-01
[101,] 2.888327e-01 5.776654e-01 7.111673e-01
[102,] 2.542341e-01 5.084682e-01 7.457659e-01
[103,] 2.097799e-01 4.195597e-01 7.902201e-01
[104,] 4.664905e-01 9.329810e-01 5.335095e-01
[105,] 4.345698e-01 8.691397e-01 5.654302e-01
[106,] 4.664315e-01 9.328631e-01 5.335685e-01
[107,] 6.401324e-01 7.197352e-01 3.598676e-01
[108,] 5.817913e-01 8.364173e-01 4.182087e-01
[109,] 5.277537e-01 9.444927e-01 4.722463e-01
[110,] 4.662552e-01 9.325104e-01 5.337448e-01
[111,] 6.940184e-01 6.119632e-01 3.059816e-01
[112,] 6.689446e-01 6.621107e-01 3.310554e-01
[113,] 6.170634e-01 7.658733e-01 3.829366e-01
[114,] 9.282299e-01 1.435402e-01 7.177008e-02
[115,] 9.052814e-01 1.894373e-01 9.471865e-02
[116,] 8.718180e-01 2.563640e-01 1.281820e-01
[117,] 8.415257e-01 3.169485e-01 1.584743e-01
[118,] 8.002699e-01 3.994602e-01 1.997301e-01
[119,] 9.998536e-01 2.928360e-04 1.464180e-04
[120,] 9.999660e-01 6.794797e-05 3.397398e-05
[121,] 9.999952e-01 9.571989e-06 4.785995e-06
[122,] 9.999818e-01 3.648765e-05 1.824383e-05
[123,] 9.999546e-01 9.082343e-05 4.541172e-05
[124,] 9.999755e-01 4.892810e-05 2.446405e-05
[125,] 9.998991e-01 2.017248e-04 1.008624e-04
[126,] 9.999532e-01 9.367411e-05 4.683706e-05
[127,] 1.000000e+00 5.215413e-12 2.607706e-12
[128,] 1.000000e+00 1.333313e-09 6.666567e-10
[129,] 1.000000e+00 1.501842e-46 7.509209e-47
> postscript(file="/var/wessaorg/rcomp/tmp/1pbws1324657832.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/2s8cx1324657832.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/3f4u01324657832.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/4etha1324657832.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/5ei351324657832.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 = 164
Frequency = 1
1 2 3 4 5 6
-3235.1158 4047.4910 -10430.5494 -8397.3459 -21202.3649 11956.4062
7 8 9 10 11 12
41981.9440 8028.6073 19157.8235 -9195.9164 -32639.5942 -28127.7267
13 14 15 16 17 18
-14085.1628 10350.3198 -19787.5285 -21487.6604 23953.4317 40503.6678
19 20 21 22 23 24
19825.1591 3193.8197 -17922.0156 20777.8686 -2118.1646 -17621.5119
25 26 27 28 29 30
-24475.0051 -2337.6128 -20293.6050 -25119.3849 -17489.7168 -24265.3880
31 32 33 34 35 36
-19170.6450 -15121.5707 -28781.2640 14119.0885 -24167.4302 5907.5960
37 38 39 40 41 42
-20897.4226 -12402.6701 -16287.2906 -19417.2188 -20469.2994 -12520.6634
43 44 45 46 47 48
-27088.8250 -20394.8776 6322.3814 3666.0250 -23579.7426 -11898.0234
49 50 51 52 53 54
-3084.0860 -17839.0395 -32069.0919 -7832.8588 5195.9761 -6184.1358
55 56 57 58 59 60
-10551.7231 -13168.3895 -5294.4597 -23677.7237 -22988.0830 -7516.2415
61 62 63 64 65 66
-2156.8653 -17162.9960 -25447.1584 -11441.2579 -17390.4330 -21119.9499
67 68 69 70 71 72
28237.6477 -54870.0228 31477.5602 36507.6020 -34091.9041 -20184.5272
73 74 75 76 77 78
-22628.0621 7972.2649 -12668.9292 23518.5414 10628.2012 113157.6419
79 80 81 82 83 84
-6998.8404 2258.4201 -16211.1541 38702.7557 -29961.3973 1763.3304
85 86 87 88 89 90
18520.2476 -8129.6766 -5625.3689 9718.8950 -33406.4071 32157.7438
91 92 93 94 95 96
-16318.5161 33764.4116 -26090.1981 -5910.0604 15998.5331 34312.2844
97 98 99 100 101 102
21896.4492 36388.6623 -27194.9116 -2139.0280 1292.9298 13490.0784
103 104 105 106 107 108
1723.9163 -1306.4031 -15684.1017 -22540.6041 -17277.0476 3171.7521
109 110 111 112 113 114
73652.0190 76168.3423 16430.8060 24016.6658 -2330.8594 -9770.1490
115 116 117 118 119 120
-4949.4417 -23513.3420 -18096.1716 -2853.7538 15147.9727 -1314.7478
121 122 123 124 125 126
-32726.4485 24409.4087 41717.4212 63889.8100 -8119.7858 22036.7329
127 128 129 130 131 132
-14451.5092 -14941.1620 -11201.9780 -9341.8261 79583.7451 -22028.3782
133 134 135 136 137 138
5727.6532 3631.6381 2635.9444 81397.4709 26272.3958 -20758.6662
139 140 141 142 143 144
-32139.1625 -27852.5453 50376.1336 16671.1648 7778.0820 18912.2340
145 146 147 148 149 150
-2880.4969 40004.9537 36813.9923 22074.8555 -21953.5738 -592.5941
151 152 153 154 155 156
-2012.3140 -1979.6184 -4223.4986 -2007.5482 5937.1394 -44712.3679
157 158 159 160 161 162
-2007.5482 -2031.1929 -2268.9377 -7384.9250 -4008.2983 20439.4194
163 164
-1937.8169 41607.1494
> postscript(file="/var/wessaorg/rcomp/tmp/6jht11324657832.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 = 164
Frequency = 1
lag(myerror, k = 1) myerror
0 -3235.1158 NA
1 4047.4910 -3235.1158
2 -10430.5494 4047.4910
3 -8397.3459 -10430.5494
4 -21202.3649 -8397.3459
5 11956.4062 -21202.3649
6 41981.9440 11956.4062
7 8028.6073 41981.9440
8 19157.8235 8028.6073
9 -9195.9164 19157.8235
10 -32639.5942 -9195.9164
11 -28127.7267 -32639.5942
12 -14085.1628 -28127.7267
13 10350.3198 -14085.1628
14 -19787.5285 10350.3198
15 -21487.6604 -19787.5285
16 23953.4317 -21487.6604
17 40503.6678 23953.4317
18 19825.1591 40503.6678
19 3193.8197 19825.1591
20 -17922.0156 3193.8197
21 20777.8686 -17922.0156
22 -2118.1646 20777.8686
23 -17621.5119 -2118.1646
24 -24475.0051 -17621.5119
25 -2337.6128 -24475.0051
26 -20293.6050 -2337.6128
27 -25119.3849 -20293.6050
28 -17489.7168 -25119.3849
29 -24265.3880 -17489.7168
30 -19170.6450 -24265.3880
31 -15121.5707 -19170.6450
32 -28781.2640 -15121.5707
33 14119.0885 -28781.2640
34 -24167.4302 14119.0885
35 5907.5960 -24167.4302
36 -20897.4226 5907.5960
37 -12402.6701 -20897.4226
38 -16287.2906 -12402.6701
39 -19417.2188 -16287.2906
40 -20469.2994 -19417.2188
41 -12520.6634 -20469.2994
42 -27088.8250 -12520.6634
43 -20394.8776 -27088.8250
44 6322.3814 -20394.8776
45 3666.0250 6322.3814
46 -23579.7426 3666.0250
47 -11898.0234 -23579.7426
48 -3084.0860 -11898.0234
49 -17839.0395 -3084.0860
50 -32069.0919 -17839.0395
51 -7832.8588 -32069.0919
52 5195.9761 -7832.8588
53 -6184.1358 5195.9761
54 -10551.7231 -6184.1358
55 -13168.3895 -10551.7231
56 -5294.4597 -13168.3895
57 -23677.7237 -5294.4597
58 -22988.0830 -23677.7237
59 -7516.2415 -22988.0830
60 -2156.8653 -7516.2415
61 -17162.9960 -2156.8653
62 -25447.1584 -17162.9960
63 -11441.2579 -25447.1584
64 -17390.4330 -11441.2579
65 -21119.9499 -17390.4330
66 28237.6477 -21119.9499
67 -54870.0228 28237.6477
68 31477.5602 -54870.0228
69 36507.6020 31477.5602
70 -34091.9041 36507.6020
71 -20184.5272 -34091.9041
72 -22628.0621 -20184.5272
73 7972.2649 -22628.0621
74 -12668.9292 7972.2649
75 23518.5414 -12668.9292
76 10628.2012 23518.5414
77 113157.6419 10628.2012
78 -6998.8404 113157.6419
79 2258.4201 -6998.8404
80 -16211.1541 2258.4201
81 38702.7557 -16211.1541
82 -29961.3973 38702.7557
83 1763.3304 -29961.3973
84 18520.2476 1763.3304
85 -8129.6766 18520.2476
86 -5625.3689 -8129.6766
87 9718.8950 -5625.3689
88 -33406.4071 9718.8950
89 32157.7438 -33406.4071
90 -16318.5161 32157.7438
91 33764.4116 -16318.5161
92 -26090.1981 33764.4116
93 -5910.0604 -26090.1981
94 15998.5331 -5910.0604
95 34312.2844 15998.5331
96 21896.4492 34312.2844
97 36388.6623 21896.4492
98 -27194.9116 36388.6623
99 -2139.0280 -27194.9116
100 1292.9298 -2139.0280
101 13490.0784 1292.9298
102 1723.9163 13490.0784
103 -1306.4031 1723.9163
104 -15684.1017 -1306.4031
105 -22540.6041 -15684.1017
106 -17277.0476 -22540.6041
107 3171.7521 -17277.0476
108 73652.0190 3171.7521
109 76168.3423 73652.0190
110 16430.8060 76168.3423
111 24016.6658 16430.8060
112 -2330.8594 24016.6658
113 -9770.1490 -2330.8594
114 -4949.4417 -9770.1490
115 -23513.3420 -4949.4417
116 -18096.1716 -23513.3420
117 -2853.7538 -18096.1716
118 15147.9727 -2853.7538
119 -1314.7478 15147.9727
120 -32726.4485 -1314.7478
121 24409.4087 -32726.4485
122 41717.4212 24409.4087
123 63889.8100 41717.4212
124 -8119.7858 63889.8100
125 22036.7329 -8119.7858
126 -14451.5092 22036.7329
127 -14941.1620 -14451.5092
128 -11201.9780 -14941.1620
129 -9341.8261 -11201.9780
130 79583.7451 -9341.8261
131 -22028.3782 79583.7451
132 5727.6532 -22028.3782
133 3631.6381 5727.6532
134 2635.9444 3631.6381
135 81397.4709 2635.9444
136 26272.3958 81397.4709
137 -20758.6662 26272.3958
138 -32139.1625 -20758.6662
139 -27852.5453 -32139.1625
140 50376.1336 -27852.5453
141 16671.1648 50376.1336
142 7778.0820 16671.1648
143 18912.2340 7778.0820
144 -2880.4969 18912.2340
145 40004.9537 -2880.4969
146 36813.9923 40004.9537
147 22074.8555 36813.9923
148 -21953.5738 22074.8555
149 -592.5941 -21953.5738
150 -2012.3140 -592.5941
151 -1979.6184 -2012.3140
152 -4223.4986 -1979.6184
153 -2007.5482 -4223.4986
154 5937.1394 -2007.5482
155 -44712.3679 5937.1394
156 -2007.5482 -44712.3679
157 -2031.1929 -2007.5482
158 -2268.9377 -2031.1929
159 -7384.9250 -2268.9377
160 -4008.2983 -7384.9250
161 20439.4194 -4008.2983
162 -1937.8169 20439.4194
163 41607.1494 -1937.8169
164 NA 41607.1494
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] 4047.4910 -3235.1158
[2,] -10430.5494 4047.4910
[3,] -8397.3459 -10430.5494
[4,] -21202.3649 -8397.3459
[5,] 11956.4062 -21202.3649
[6,] 41981.9440 11956.4062
[7,] 8028.6073 41981.9440
[8,] 19157.8235 8028.6073
[9,] -9195.9164 19157.8235
[10,] -32639.5942 -9195.9164
[11,] -28127.7267 -32639.5942
[12,] -14085.1628 -28127.7267
[13,] 10350.3198 -14085.1628
[14,] -19787.5285 10350.3198
[15,] -21487.6604 -19787.5285
[16,] 23953.4317 -21487.6604
[17,] 40503.6678 23953.4317
[18,] 19825.1591 40503.6678
[19,] 3193.8197 19825.1591
[20,] -17922.0156 3193.8197
[21,] 20777.8686 -17922.0156
[22,] -2118.1646 20777.8686
[23,] -17621.5119 -2118.1646
[24,] -24475.0051 -17621.5119
[25,] -2337.6128 -24475.0051
[26,] -20293.6050 -2337.6128
[27,] -25119.3849 -20293.6050
[28,] -17489.7168 -25119.3849
[29,] -24265.3880 -17489.7168
[30,] -19170.6450 -24265.3880
[31,] -15121.5707 -19170.6450
[32,] -28781.2640 -15121.5707
[33,] 14119.0885 -28781.2640
[34,] -24167.4302 14119.0885
[35,] 5907.5960 -24167.4302
[36,] -20897.4226 5907.5960
[37,] -12402.6701 -20897.4226
[38,] -16287.2906 -12402.6701
[39,] -19417.2188 -16287.2906
[40,] -20469.2994 -19417.2188
[41,] -12520.6634 -20469.2994
[42,] -27088.8250 -12520.6634
[43,] -20394.8776 -27088.8250
[44,] 6322.3814 -20394.8776
[45,] 3666.0250 6322.3814
[46,] -23579.7426 3666.0250
[47,] -11898.0234 -23579.7426
[48,] -3084.0860 -11898.0234
[49,] -17839.0395 -3084.0860
[50,] -32069.0919 -17839.0395
[51,] -7832.8588 -32069.0919
[52,] 5195.9761 -7832.8588
[53,] -6184.1358 5195.9761
[54,] -10551.7231 -6184.1358
[55,] -13168.3895 -10551.7231
[56,] -5294.4597 -13168.3895
[57,] -23677.7237 -5294.4597
[58,] -22988.0830 -23677.7237
[59,] -7516.2415 -22988.0830
[60,] -2156.8653 -7516.2415
[61,] -17162.9960 -2156.8653
[62,] -25447.1584 -17162.9960
[63,] -11441.2579 -25447.1584
[64,] -17390.4330 -11441.2579
[65,] -21119.9499 -17390.4330
[66,] 28237.6477 -21119.9499
[67,] -54870.0228 28237.6477
[68,] 31477.5602 -54870.0228
[69,] 36507.6020 31477.5602
[70,] -34091.9041 36507.6020
[71,] -20184.5272 -34091.9041
[72,] -22628.0621 -20184.5272
[73,] 7972.2649 -22628.0621
[74,] -12668.9292 7972.2649
[75,] 23518.5414 -12668.9292
[76,] 10628.2012 23518.5414
[77,] 113157.6419 10628.2012
[78,] -6998.8404 113157.6419
[79,] 2258.4201 -6998.8404
[80,] -16211.1541 2258.4201
[81,] 38702.7557 -16211.1541
[82,] -29961.3973 38702.7557
[83,] 1763.3304 -29961.3973
[84,] 18520.2476 1763.3304
[85,] -8129.6766 18520.2476
[86,] -5625.3689 -8129.6766
[87,] 9718.8950 -5625.3689
[88,] -33406.4071 9718.8950
[89,] 32157.7438 -33406.4071
[90,] -16318.5161 32157.7438
[91,] 33764.4116 -16318.5161
[92,] -26090.1981 33764.4116
[93,] -5910.0604 -26090.1981
[94,] 15998.5331 -5910.0604
[95,] 34312.2844 15998.5331
[96,] 21896.4492 34312.2844
[97,] 36388.6623 21896.4492
[98,] -27194.9116 36388.6623
[99,] -2139.0280 -27194.9116
[100,] 1292.9298 -2139.0280
[101,] 13490.0784 1292.9298
[102,] 1723.9163 13490.0784
[103,] -1306.4031 1723.9163
[104,] -15684.1017 -1306.4031
[105,] -22540.6041 -15684.1017
[106,] -17277.0476 -22540.6041
[107,] 3171.7521 -17277.0476
[108,] 73652.0190 3171.7521
[109,] 76168.3423 73652.0190
[110,] 16430.8060 76168.3423
[111,] 24016.6658 16430.8060
[112,] -2330.8594 24016.6658
[113,] -9770.1490 -2330.8594
[114,] -4949.4417 -9770.1490
[115,] -23513.3420 -4949.4417
[116,] -18096.1716 -23513.3420
[117,] -2853.7538 -18096.1716
[118,] 15147.9727 -2853.7538
[119,] -1314.7478 15147.9727
[120,] -32726.4485 -1314.7478
[121,] 24409.4087 -32726.4485
[122,] 41717.4212 24409.4087
[123,] 63889.8100 41717.4212
[124,] -8119.7858 63889.8100
[125,] 22036.7329 -8119.7858
[126,] -14451.5092 22036.7329
[127,] -14941.1620 -14451.5092
[128,] -11201.9780 -14941.1620
[129,] -9341.8261 -11201.9780
[130,] 79583.7451 -9341.8261
[131,] -22028.3782 79583.7451
[132,] 5727.6532 -22028.3782
[133,] 3631.6381 5727.6532
[134,] 2635.9444 3631.6381
[135,] 81397.4709 2635.9444
[136,] 26272.3958 81397.4709
[137,] -20758.6662 26272.3958
[138,] -32139.1625 -20758.6662
[139,] -27852.5453 -32139.1625
[140,] 50376.1336 -27852.5453
[141,] 16671.1648 50376.1336
[142,] 7778.0820 16671.1648
[143,] 18912.2340 7778.0820
[144,] -2880.4969 18912.2340
[145,] 40004.9537 -2880.4969
[146,] 36813.9923 40004.9537
[147,] 22074.8555 36813.9923
[148,] -21953.5738 22074.8555
[149,] -592.5941 -21953.5738
[150,] -2012.3140 -592.5941
[151,] -1979.6184 -2012.3140
[152,] -4223.4986 -1979.6184
[153,] -2007.5482 -4223.4986
[154,] 5937.1394 -2007.5482
[155,] -44712.3679 5937.1394
[156,] -2007.5482 -44712.3679
[157,] -2031.1929 -2007.5482
[158,] -2268.9377 -2031.1929
[159,] -7384.9250 -2268.9377
[160,] -4008.2983 -7384.9250
[161,] 20439.4194 -4008.2983
[162,] -1937.8169 20439.4194
[163,] 41607.1494 -1937.8169
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 4047.4910 -3235.1158
2 -10430.5494 4047.4910
3 -8397.3459 -10430.5494
4 -21202.3649 -8397.3459
5 11956.4062 -21202.3649
6 41981.9440 11956.4062
7 8028.6073 41981.9440
8 19157.8235 8028.6073
9 -9195.9164 19157.8235
10 -32639.5942 -9195.9164
11 -28127.7267 -32639.5942
12 -14085.1628 -28127.7267
13 10350.3198 -14085.1628
14 -19787.5285 10350.3198
15 -21487.6604 -19787.5285
16 23953.4317 -21487.6604
17 40503.6678 23953.4317
18 19825.1591 40503.6678
19 3193.8197 19825.1591
20 -17922.0156 3193.8197
21 20777.8686 -17922.0156
22 -2118.1646 20777.8686
23 -17621.5119 -2118.1646
24 -24475.0051 -17621.5119
25 -2337.6128 -24475.0051
26 -20293.6050 -2337.6128
27 -25119.3849 -20293.6050
28 -17489.7168 -25119.3849
29 -24265.3880 -17489.7168
30 -19170.6450 -24265.3880
31 -15121.5707 -19170.6450
32 -28781.2640 -15121.5707
33 14119.0885 -28781.2640
34 -24167.4302 14119.0885
35 5907.5960 -24167.4302
36 -20897.4226 5907.5960
37 -12402.6701 -20897.4226
38 -16287.2906 -12402.6701
39 -19417.2188 -16287.2906
40 -20469.2994 -19417.2188
41 -12520.6634 -20469.2994
42 -27088.8250 -12520.6634
43 -20394.8776 -27088.8250
44 6322.3814 -20394.8776
45 3666.0250 6322.3814
46 -23579.7426 3666.0250
47 -11898.0234 -23579.7426
48 -3084.0860 -11898.0234
49 -17839.0395 -3084.0860
50 -32069.0919 -17839.0395
51 -7832.8588 -32069.0919
52 5195.9761 -7832.8588
53 -6184.1358 5195.9761
54 -10551.7231 -6184.1358
55 -13168.3895 -10551.7231
56 -5294.4597 -13168.3895
57 -23677.7237 -5294.4597
58 -22988.0830 -23677.7237
59 -7516.2415 -22988.0830
60 -2156.8653 -7516.2415
61 -17162.9960 -2156.8653
62 -25447.1584 -17162.9960
63 -11441.2579 -25447.1584
64 -17390.4330 -11441.2579
65 -21119.9499 -17390.4330
66 28237.6477 -21119.9499
67 -54870.0228 28237.6477
68 31477.5602 -54870.0228
69 36507.6020 31477.5602
70 -34091.9041 36507.6020
71 -20184.5272 -34091.9041
72 -22628.0621 -20184.5272
73 7972.2649 -22628.0621
74 -12668.9292 7972.2649
75 23518.5414 -12668.9292
76 10628.2012 23518.5414
77 113157.6419 10628.2012
78 -6998.8404 113157.6419
79 2258.4201 -6998.8404
80 -16211.1541 2258.4201
81 38702.7557 -16211.1541
82 -29961.3973 38702.7557
83 1763.3304 -29961.3973
84 18520.2476 1763.3304
85 -8129.6766 18520.2476
86 -5625.3689 -8129.6766
87 9718.8950 -5625.3689
88 -33406.4071 9718.8950
89 32157.7438 -33406.4071
90 -16318.5161 32157.7438
91 33764.4116 -16318.5161
92 -26090.1981 33764.4116
93 -5910.0604 -26090.1981
94 15998.5331 -5910.0604
95 34312.2844 15998.5331
96 21896.4492 34312.2844
97 36388.6623 21896.4492
98 -27194.9116 36388.6623
99 -2139.0280 -27194.9116
100 1292.9298 -2139.0280
101 13490.0784 1292.9298
102 1723.9163 13490.0784
103 -1306.4031 1723.9163
104 -15684.1017 -1306.4031
105 -22540.6041 -15684.1017
106 -17277.0476 -22540.6041
107 3171.7521 -17277.0476
108 73652.0190 3171.7521
109 76168.3423 73652.0190
110 16430.8060 76168.3423
111 24016.6658 16430.8060
112 -2330.8594 24016.6658
113 -9770.1490 -2330.8594
114 -4949.4417 -9770.1490
115 -23513.3420 -4949.4417
116 -18096.1716 -23513.3420
117 -2853.7538 -18096.1716
118 15147.9727 -2853.7538
119 -1314.7478 15147.9727
120 -32726.4485 -1314.7478
121 24409.4087 -32726.4485
122 41717.4212 24409.4087
123 63889.8100 41717.4212
124 -8119.7858 63889.8100
125 22036.7329 -8119.7858
126 -14451.5092 22036.7329
127 -14941.1620 -14451.5092
128 -11201.9780 -14941.1620
129 -9341.8261 -11201.9780
130 79583.7451 -9341.8261
131 -22028.3782 79583.7451
132 5727.6532 -22028.3782
133 3631.6381 5727.6532
134 2635.9444 3631.6381
135 81397.4709 2635.9444
136 26272.3958 81397.4709
137 -20758.6662 26272.3958
138 -32139.1625 -20758.6662
139 -27852.5453 -32139.1625
140 50376.1336 -27852.5453
141 16671.1648 50376.1336
142 7778.0820 16671.1648
143 18912.2340 7778.0820
144 -2880.4969 18912.2340
145 40004.9537 -2880.4969
146 36813.9923 40004.9537
147 22074.8555 36813.9923
148 -21953.5738 22074.8555
149 -592.5941 -21953.5738
150 -2012.3140 -592.5941
151 -1979.6184 -2012.3140
152 -4223.4986 -1979.6184
153 -2007.5482 -4223.4986
154 5937.1394 -2007.5482
155 -44712.3679 5937.1394
156 -2007.5482 -44712.3679
157 -2031.1929 -2007.5482
158 -2268.9377 -2031.1929
159 -7384.9250 -2268.9377
160 -4008.2983 -7384.9250
161 20439.4194 -4008.2983
162 -1937.8169 20439.4194
163 41607.1494 -1937.8169
> 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/7cvkv1324657832.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/8fwl51324657832.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/98iy11324657832.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/10qay61324657832.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/11et7t1324657832.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/12mo411324657832.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/13677q1324657832.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/14p0r11324657832.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/15wcwc1324657832.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/169hhi1324657832.tab")
+ }
>
> try(system("convert tmp/1pbws1324657832.ps tmp/1pbws1324657832.png",intern=TRUE))
character(0)
> try(system("convert tmp/2s8cx1324657832.ps tmp/2s8cx1324657832.png",intern=TRUE))
character(0)
> try(system("convert tmp/3f4u01324657832.ps tmp/3f4u01324657832.png",intern=TRUE))
character(0)
> try(system("convert tmp/4etha1324657832.ps tmp/4etha1324657832.png",intern=TRUE))
character(0)
> try(system("convert tmp/5ei351324657832.ps tmp/5ei351324657832.png",intern=TRUE))
character(0)
> try(system("convert tmp/6jht11324657832.ps tmp/6jht11324657832.png",intern=TRUE))
character(0)
> try(system("convert tmp/7cvkv1324657832.ps tmp/7cvkv1324657832.png",intern=TRUE))
character(0)
> try(system("convert tmp/8fwl51324657832.ps tmp/8fwl51324657832.png",intern=TRUE))
character(0)
> try(system("convert tmp/98iy11324657832.ps tmp/98iy11324657832.png",intern=TRUE))
character(0)
> try(system("convert tmp/10qay61324657832.ps tmp/10qay61324657832.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
6.545 0.721 7.296