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(210907
+ ,79
+ ,81
+ ,94
+ ,56
+ ,3
+ ,1418
+ ,112285
+ ,120982
+ ,58
+ ,55
+ ,103
+ ,56
+ ,4
+ ,869
+ ,84786
+ ,176508
+ ,60
+ ,50
+ ,93
+ ,54
+ ,12
+ ,1530
+ ,83123
+ ,179321
+ ,108
+ ,125
+ ,103
+ ,89
+ ,2
+ ,2172
+ ,101193
+ ,123185
+ ,49
+ ,40
+ ,51
+ ,40
+ ,1
+ ,901
+ ,38361
+ ,52746
+ ,0
+ ,37
+ ,70
+ ,25
+ ,3
+ ,463
+ ,68504
+ ,385534
+ ,121
+ ,63
+ ,91
+ ,92
+ ,0
+ ,3201
+ ,119182
+ ,33170
+ ,1
+ ,44
+ ,22
+ ,18
+ ,0
+ ,371
+ ,22807
+ ,101645
+ ,20
+ ,88
+ ,38
+ ,63
+ ,0
+ ,1192
+ ,17140
+ ,149061
+ ,43
+ ,66
+ ,93
+ ,44
+ ,5
+ ,1583
+ ,116174
+ ,165446
+ ,69
+ ,57
+ ,60
+ ,33
+ ,0
+ ,1439
+ ,57635
+ ,237213
+ ,78
+ ,74
+ ,123
+ ,84
+ ,0
+ ,1764
+ ,66198
+ ,173326
+ ,86
+ ,49
+ ,148
+ ,88
+ ,7
+ ,1495
+ ,71701
+ ,133131
+ ,44
+ ,52
+ ,90
+ ,55
+ ,7
+ ,1373
+ ,57793
+ ,258873
+ ,104
+ ,88
+ ,124
+ ,60
+ ,3
+ ,2187
+ ,80444
+ ,180083
+ ,63
+ ,36
+ ,70
+ ,66
+ ,9
+ ,1491
+ ,53855
+ ,324799
+ ,158
+ ,108
+ ,168
+ ,154
+ ,0
+ ,4041
+ ,97668
+ ,230964
+ ,102
+ ,43
+ ,115
+ ,53
+ ,4
+ ,1706
+ ,133824
+ ,236785
+ ,77
+ ,75
+ ,71
+ ,119
+ ,3
+ ,2152
+ ,101481
+ ,135473
+ ,82
+ ,32
+ ,66
+ ,41
+ ,0
+ ,1036
+ ,99645
+ ,202925
+ ,115
+ ,44
+ ,134
+ ,61
+ ,7
+ ,1882
+ ,114789
+ ,215147
+ ,101
+ ,85
+ ,117
+ ,58
+ ,0
+ ,1929
+ ,99052
+ ,344297
+ ,80
+ ,86
+ ,108
+ ,75
+ ,1
+ ,2242
+ ,67654
+ ,153935
+ ,50
+ ,56
+ ,84
+ ,33
+ ,5
+ ,1220
+ ,65553
+ ,132943
+ ,83
+ ,50
+ ,156
+ ,40
+ ,7
+ ,1289
+ ,97500
+ ,174724
+ ,123
+ ,135
+ ,120
+ ,92
+ ,0
+ ,2515
+ ,69112
+ ,174415
+ ,73
+ ,63
+ ,114
+ ,100
+ ,0
+ ,2147
+ ,82753
+ ,225548
+ ,81
+ ,81
+ ,94
+ ,112
+ ,5
+ ,2352
+ ,85323
+ ,223632
+ ,105
+ ,52
+ ,120
+ ,73
+ ,0
+ ,1638
+ ,72654
+ ,124817
+ ,47
+ ,44
+ ,81
+ ,40
+ ,0
+ ,1222
+ ,30727
+ ,221698
+ ,105
+ ,113
+ ,110
+ ,45
+ ,0
+ ,1812
+ ,77873
+ ,210767
+ ,94
+ ,39
+ ,133
+ ,60
+ ,3
+ ,1677
+ ,117478
+ ,170266
+ ,44
+ ,73
+ ,122
+ ,62
+ ,4
+ ,1579
+ ,74007
+ ,260561
+ ,114
+ ,48
+ ,158
+ ,75
+ ,1
+ ,1731
+ ,90183
+ ,84853
+ ,38
+ ,33
+ ,109
+ ,31
+ ,4
+ ,807
+ ,61542
+ ,294424
+ ,107
+ ,59
+ ,124
+ ,77
+ ,2
+ ,2452
+ ,101494
+ ,101011
+ ,30
+ ,41
+ ,39
+ ,34
+ ,0
+ ,829
+ ,27570
+ ,215641
+ ,71
+ ,69
+ ,92
+ ,46
+ ,0
+ ,1940
+ ,55813
+ ,325107
+ ,84
+ ,64
+ ,126
+ ,99
+ ,0
+ ,2662
+ ,79215
+ ,7176
+ ,0
+ ,1
+ ,0
+ ,17
+ ,0
+ ,186
+ ,1423
+ ,167542
+ ,59
+ ,59
+ ,70
+ ,66
+ ,2
+ ,1499
+ ,55461
+ ,106408
+ ,33
+ ,32
+ ,37
+ ,30
+ ,1
+ ,865
+ ,31081
+ ,96560
+ ,42
+ ,129
+ ,38
+ ,76
+ ,0
+ ,1793
+ ,22996
+ ,265769
+ ,96
+ ,37
+ ,120
+ ,146
+ ,2
+ ,2527
+ ,83122
+ ,269651
+ ,106
+ ,31
+ ,93
+ ,67
+ ,10
+ ,2747
+ ,70106
+ ,149112
+ ,56
+ ,65
+ ,95
+ ,56
+ ,6
+ ,1324
+ ,60578
+ ,175824
+ ,57
+ ,107
+ ,77
+ ,107
+ ,0
+ ,2702
+ ,39992
+ ,152871
+ ,59
+ ,74
+ ,90
+ ,58
+ ,5
+ ,1383
+ ,79892
+ ,111665
+ ,39
+ ,54
+ ,80
+ ,34
+ ,4
+ ,1179
+ ,49810
+ ,116408
+ ,34
+ ,76
+ ,31
+ ,61
+ ,1
+ ,2099
+ ,71570
+ ,362301
+ ,76
+ ,715
+ ,110
+ ,119
+ ,2
+ ,4308
+ ,100708
+ ,78800
+ ,20
+ ,57
+ ,66
+ ,42
+ ,2
+ ,918
+ ,33032
+ ,183167
+ ,91
+ ,66
+ ,138
+ ,66
+ ,0
+ ,1831
+ ,82875
+ ,277965
+ ,115
+ ,106
+ ,133
+ ,89
+ ,8
+ ,3373
+ ,139077
+ ,150629
+ ,85
+ ,54
+ ,113
+ ,44
+ ,3
+ ,1713
+ ,71595
+ ,168809
+ ,76
+ ,32
+ ,100
+ ,66
+ ,0
+ ,1438
+ ,72260
+ ,24188
+ ,8
+ ,20
+ ,7
+ ,24
+ ,0
+ ,496
+ ,5950
+ ,329267
+ ,79
+ ,71
+ ,140
+ ,259
+ ,8
+ ,2253
+ ,115762
+ ,65029
+ ,21
+ ,21
+ ,61
+ ,17
+ ,5
+ ,744
+ ,32551
+ ,101097
+ ,30
+ ,70
+ ,41
+ ,64
+ ,3
+ ,1161
+ ,31701
+ ,218946
+ ,76
+ ,112
+ ,96
+ ,41
+ ,1
+ ,2352
+ ,80670
+ ,244052
+ ,101
+ ,66
+ ,164
+ ,68
+ ,5
+ ,2144
+ ,143558
+ ,341570
+ ,94
+ ,190
+ ,78
+ ,168
+ ,1
+ ,4691
+ ,117105
+ ,103597
+ ,27
+ ,66
+ ,49
+ ,43
+ ,1
+ ,1112
+ ,23789
+ ,233328
+ ,92
+ ,165
+ ,102
+ ,132
+ ,5
+ ,2694
+ ,120733
+ ,256462
+ ,123
+ ,56
+ ,124
+ ,105
+ ,0
+ ,1973
+ ,105195
+ ,206161
+ ,75
+ ,61
+ ,99
+ ,71
+ ,12
+ ,1769
+ ,73107
+ ,311473
+ ,128
+ ,53
+ ,129
+ ,112
+ ,8
+ ,3148
+ ,132068
+ ,235800
+ ,105
+ ,127
+ ,62
+ ,94
+ ,8
+ ,2474
+ ,149193
+ ,177939
+ ,55
+ ,63
+ ,73
+ ,82
+ ,8
+ ,2084
+ ,46821
+ ,207176
+ ,56
+ ,38
+ ,114
+ ,70
+ ,8
+ ,1954
+ ,87011
+ ,196553
+ ,41
+ ,50
+ ,99
+ ,57
+ ,2
+ ,1226
+ ,95260
+ ,174184
+ ,72
+ ,52
+ ,70
+ ,53
+ ,0
+ ,1389
+ ,55183
+ ,143246
+ ,67
+ ,42
+ ,104
+ ,103
+ ,5
+ ,1496
+ ,106671
+ ,187559
+ ,75
+ ,76
+ ,116
+ ,121
+ ,8
+ ,2269
+ ,73511
+ ,187681
+ ,114
+ ,67
+ ,91
+ ,62
+ ,2
+ ,1833
+ ,92945
+ ,119016
+ ,118
+ ,50
+ ,74
+ ,52
+ ,5
+ ,1268
+ ,78664
+ ,182192
+ ,77
+ ,53
+ ,138
+ ,52
+ ,12
+ ,1943
+ ,70054
+ ,73566
+ ,22
+ ,39
+ ,67
+ ,32
+ ,6
+ ,893
+ ,22618
+ ,194979
+ ,66
+ ,50
+ ,151
+ ,62
+ ,7
+ ,1762
+ ,74011
+ ,167488
+ ,69
+ ,77
+ ,72
+ ,45
+ ,2
+ ,1403
+ ,83737
+ ,143756
+ ,105
+ ,57
+ ,120
+ ,46
+ ,0
+ ,1425
+ ,69094
+ ,275541
+ ,116
+ ,73
+ ,115
+ ,63
+ ,4
+ ,1857
+ ,93133
+ ,243199
+ ,88
+ ,34
+ ,105
+ ,75
+ ,3
+ ,1840
+ ,95536
+ ,182999
+ ,73
+ ,39
+ ,104
+ ,88
+ ,6
+ ,1502
+ ,225920
+ ,135649
+ ,99
+ ,46
+ ,108
+ ,46
+ ,2
+ ,1441
+ ,62133
+ ,152299
+ ,62
+ ,63
+ ,98
+ ,53
+ ,0
+ ,1420
+ ,61370
+ ,120221
+ ,53
+ ,35
+ ,69
+ ,37
+ ,1
+ ,1416
+ ,43836
+ ,346485
+ ,118
+ ,106
+ ,111
+ ,90
+ ,0
+ ,2970
+ ,106117
+ ,145790
+ ,30
+ ,43
+ ,99
+ ,63
+ ,5
+ ,1317
+ ,38692
+ ,193339
+ ,100
+ ,47
+ ,71
+ ,78
+ ,2
+ ,1644
+ ,84651
+ ,80953
+ ,49
+ ,31
+ ,27
+ ,25
+ ,0
+ ,870
+ ,56622
+ ,122774
+ ,24
+ ,162
+ ,69
+ ,45
+ ,0
+ ,1654
+ ,15986
+ ,130585
+ ,67
+ ,57
+ ,107
+ ,46
+ ,5
+ ,1054
+ ,95364
+ ,112611
+ ,46
+ ,36
+ ,73
+ ,41
+ ,0
+ ,937
+ ,26706
+ ,286468
+ ,57
+ ,263
+ ,107
+ ,144
+ ,1
+ ,3004
+ ,89691
+ ,241066
+ ,75
+ ,78
+ ,93
+ ,82
+ ,0
+ ,2008
+ ,67267
+ ,148446
+ ,135
+ ,63
+ ,129
+ ,91
+ ,1
+ ,2547
+ ,126846
+ ,204713
+ ,68
+ ,54
+ ,69
+ ,71
+ ,1
+ ,1885
+ ,41140
+ ,182079
+ ,124
+ ,63
+ ,118
+ ,63
+ ,2
+ ,1626
+ ,102860
+ ,140344
+ ,33
+ ,77
+ ,73
+ ,53
+ ,6
+ ,1468
+ ,51715
+ ,220516
+ ,98
+ ,79
+ ,119
+ ,62
+ ,1
+ ,2445
+ ,55801
+ ,243060
+ ,58
+ ,110
+ ,104
+ ,63
+ ,4
+ ,1964
+ ,111813
+ ,162765
+ ,68
+ ,56
+ ,107
+ ,32
+ ,2
+ ,1381
+ ,120293
+ ,182613
+ ,81
+ ,56
+ ,99
+ ,39
+ ,3
+ ,1369
+ ,138599
+ ,232138
+ ,131
+ ,43
+ ,90
+ ,62
+ ,0
+ ,1659
+ ,161647
+ ,265318
+ ,110
+ ,111
+ ,197
+ ,117
+ ,10
+ ,2888
+ ,115929
+ ,85574
+ ,37
+ ,71
+ ,36
+ ,34
+ ,0
+ ,1290
+ ,24266
+ ,310839
+ ,130
+ ,62
+ ,85
+ ,92
+ ,9
+ ,2845
+ ,162901
+ ,225060
+ ,93
+ ,56
+ ,139
+ ,93
+ ,7
+ ,1982
+ ,109825
+ ,232317
+ ,118
+ ,74
+ ,106
+ ,54
+ ,0
+ ,1904
+ ,129838
+ ,144966
+ ,39
+ ,60
+ ,50
+ ,144
+ ,0
+ ,1391
+ ,37510
+ ,43287
+ ,13
+ ,43
+ ,64
+ ,14
+ ,4
+ ,602
+ ,43750
+ ,155754
+ ,74
+ ,68
+ ,31
+ ,61
+ ,4
+ ,1743
+ ,40652
+ ,164709
+ ,81
+ ,53
+ ,63
+ ,109
+ ,0
+ ,1559
+ ,87771
+ ,201940
+ ,109
+ ,87
+ ,92
+ ,38
+ ,0
+ ,2014
+ ,85872
+ ,235454
+ ,151
+ ,46
+ ,106
+ ,73
+ ,0
+ ,2143
+ ,89275
+ ,220801
+ ,51
+ ,105
+ ,63
+ ,75
+ ,1
+ ,2146
+ ,44418
+ ,99466
+ ,28
+ ,32
+ ,69
+ ,50
+ ,0
+ ,874
+ ,192565
+ ,92661
+ ,40
+ ,133
+ ,41
+ ,61
+ ,1
+ ,1590
+ ,35232
+ ,133328
+ ,56
+ ,79
+ ,56
+ ,55
+ ,0
+ ,1590
+ ,40909
+ ,61361
+ ,27
+ ,51
+ ,25
+ ,77
+ ,0
+ ,1210
+ ,13294
+ ,125930
+ ,37
+ ,207
+ ,65
+ ,75
+ ,4
+ ,2072
+ ,32387
+ ,100750
+ ,83
+ ,67
+ ,93
+ ,72
+ ,0
+ ,1281
+ ,140867
+ ,224549
+ ,54
+ ,47
+ ,114
+ ,50
+ ,4
+ ,1401
+ ,120662
+ ,82316
+ ,27
+ ,34
+ ,38
+ ,32
+ ,4
+ ,834
+ ,21233
+ ,102010
+ ,28
+ ,66
+ ,44
+ ,53
+ ,3
+ ,1105
+ ,44332
+ ,101523
+ ,59
+ ,76
+ ,87
+ ,42
+ ,0
+ ,1272
+ ,61056
+ ,243511
+ ,133
+ ,65
+ ,110
+ ,71
+ ,0
+ ,1944
+ ,101338
+ ,22938
+ ,12
+ ,9
+ ,0
+ ,10
+ ,0
+ ,391
+ ,1168
+ ,41566
+ ,0
+ ,42
+ ,27
+ ,35
+ ,5
+ ,761
+ ,13497
+ ,152474
+ ,106
+ ,45
+ ,83
+ ,65
+ ,0
+ ,1605
+ ,65567
+ ,61857
+ ,23
+ ,25
+ ,30
+ ,25
+ ,4
+ ,530
+ ,25162
+ ,99923
+ ,44
+ ,115
+ ,80
+ ,66
+ ,0
+ ,1988
+ ,32334
+ ,132487
+ ,71
+ ,97
+ ,98
+ ,41
+ ,0
+ ,1386
+ ,40735
+ ,317394
+ ,116
+ ,53
+ ,82
+ ,86
+ ,1
+ ,2395
+ ,91413
+ ,21054
+ ,4
+ ,2
+ ,0
+ ,16
+ ,0
+ ,387
+ ,855
+ ,209641
+ ,62
+ ,52
+ ,60
+ ,42
+ ,5
+ ,1742
+ ,97068
+ ,22648
+ ,12
+ ,44
+ ,28
+ ,19
+ ,0
+ ,620
+ ,44339
+ ,31414
+ ,18
+ ,22
+ ,9
+ ,19
+ ,0
+ ,449
+ ,14116
+ ,46698
+ ,14
+ ,35
+ ,33
+ ,45
+ ,0
+ ,800
+ ,10288
+ ,131698
+ ,60
+ ,74
+ ,59
+ ,65
+ ,0
+ ,1684
+ ,65622
+ ,91735
+ ,7
+ ,103
+ ,49
+ ,35
+ ,0
+ ,1050
+ ,16563
+ ,244749
+ ,98
+ ,144
+ ,115
+ ,95
+ ,2
+ ,2699
+ ,76643
+ ,184510
+ ,64
+ ,60
+ ,140
+ ,49
+ ,7
+ ,1606
+ ,110681
+ ,79863
+ ,29
+ ,134
+ ,49
+ ,37
+ ,1
+ ,1502
+ ,29011
+ ,128423
+ ,32
+ ,89
+ ,120
+ ,64
+ ,8
+ ,1204
+ ,92696
+ ,97839
+ ,25
+ ,42
+ ,66
+ ,38
+ ,2
+ ,1138
+ ,94785
+ ,38214
+ ,16
+ ,52
+ ,21
+ ,34
+ ,0
+ ,568
+ ,8773
+ ,151101
+ ,48
+ ,98
+ ,124
+ ,32
+ ,2
+ ,1459
+ ,83209
+ ,272458
+ ,100
+ ,99
+ ,152
+ ,65
+ ,0
+ ,2158
+ ,93815
+ ,172494
+ ,46
+ ,52
+ ,139
+ ,52
+ ,0
+ ,1111
+ ,86687
+ ,108043
+ ,45
+ ,29
+ ,38
+ ,62
+ ,1
+ ,1421
+ ,34553
+ ,328107
+ ,129
+ ,125
+ ,144
+ ,65
+ ,3
+ ,2833
+ ,105547
+ ,250579
+ ,130
+ ,106
+ ,120
+ ,83
+ ,0
+ ,1955
+ ,103487
+ ,351067
+ ,136
+ ,95
+ ,160
+ ,95
+ ,3
+ ,2922
+ ,213688
+ ,158015
+ ,59
+ ,40
+ ,114
+ ,29
+ ,0
+ ,1002
+ ,71220
+ ,98866
+ ,25
+ ,140
+ ,39
+ ,18
+ ,0
+ ,1060
+ ,23517
+ ,85439
+ ,32
+ ,43
+ ,78
+ ,33
+ ,0
+ ,956
+ ,56926
+ ,229242
+ ,63
+ ,128
+ ,119
+ ,247
+ ,4
+ ,2186
+ ,91721
+ ,351619
+ ,95
+ ,142
+ ,141
+ ,139
+ ,4
+ ,3604
+ ,115168
+ ,84207
+ ,14
+ ,73
+ ,101
+ ,29
+ ,11
+ ,1035
+ ,111194
+ ,120445
+ ,36
+ ,72
+ ,56
+ ,118
+ ,0
+ ,1417
+ ,51009
+ ,324598
+ ,113
+ ,128
+ ,133
+ ,110
+ ,0
+ ,3261
+ ,135777
+ ,131069
+ ,47
+ ,61
+ ,83
+ ,67
+ ,4
+ ,1587
+ ,51513
+ ,204271
+ ,92
+ ,73
+ ,116
+ ,42
+ ,0
+ ,1424
+ ,74163
+ ,165543
+ ,70
+ ,148
+ ,90
+ ,65
+ ,1
+ ,1701
+ ,51633
+ ,141722
+ ,19
+ ,64
+ ,36
+ ,94
+ ,0
+ ,1249
+ ,75345
+ ,116048
+ ,50
+ ,45
+ ,50
+ ,64
+ ,0
+ ,946
+ ,33416
+ ,250047
+ ,41
+ ,58
+ ,61
+ ,81
+ ,0
+ ,1926
+ ,83305
+ ,299775
+ ,91
+ ,97
+ ,97
+ ,95
+ ,9
+ ,3352
+ ,98952
+ ,195838
+ ,111
+ ,50
+ ,98
+ ,67
+ ,1
+ ,1641
+ ,102372
+ ,173260
+ ,41
+ ,37
+ ,78
+ ,63
+ ,3
+ ,2035
+ ,37238
+ ,254488
+ ,120
+ ,50
+ ,117
+ ,83
+ ,10
+ ,2312
+ ,103772
+ ,104389
+ ,135
+ ,105
+ ,148
+ ,45
+ ,5
+ ,1369
+ ,123969
+ ,136084
+ ,27
+ ,69
+ ,41
+ ,30
+ ,0
+ ,1577
+ ,27142
+ ,199476
+ ,87
+ ,46
+ ,105
+ ,70
+ ,2
+ ,2201
+ ,135400
+ ,92499
+ ,25
+ ,57
+ ,55
+ ,32
+ ,0
+ ,961
+ ,21399
+ ,224330
+ ,131
+ ,52
+ ,132
+ ,83
+ ,1
+ ,1900
+ ,130115
+ ,135781
+ ,45
+ ,98
+ ,44
+ ,31
+ ,2
+ ,1254
+ ,24874
+ ,74408
+ ,29
+ ,61
+ ,21
+ ,67
+ ,4
+ ,1335
+ ,34988
+ ,81240
+ ,58
+ ,89
+ ,50
+ ,66
+ ,0
+ ,1597
+ ,45549
+ ,14688
+ ,4
+ ,0
+ ,0
+ ,10
+ ,0
+ ,207
+ ,6023
+ ,181633
+ ,47
+ ,48
+ ,73
+ ,70
+ ,2
+ ,1645
+ ,64466
+ ,271856
+ ,109
+ ,91
+ ,86
+ ,103
+ ,1
+ ,2429
+ ,54990
+ ,7199
+ ,7
+ ,0
+ ,0
+ ,5
+ ,0
+ ,151
+ ,1644
+ ,46660
+ ,12
+ ,7
+ ,13
+ ,20
+ ,0
+ ,474
+ ,6179
+ ,17547
+ ,0
+ ,3
+ ,4
+ ,5
+ ,0
+ ,141
+ ,3926
+ ,133368
+ ,37
+ ,54
+ ,57
+ ,36
+ ,1
+ ,1639
+ ,32755
+ ,95227
+ ,37
+ ,70
+ ,48
+ ,34
+ ,0
+ ,872
+ ,34777
+ ,152601
+ ,46
+ ,36
+ ,46
+ ,48
+ ,2
+ ,1318
+ ,73224
+ ,98146
+ ,15
+ ,37
+ ,48
+ ,40
+ ,0
+ ,1018
+ ,27114
+ ,79619
+ ,42
+ ,123
+ ,32
+ ,43
+ ,3
+ ,1383
+ ,20760
+ ,59194
+ ,7
+ ,247
+ ,68
+ ,31
+ ,6
+ ,1314
+ ,37636
+ ,139942
+ ,54
+ ,46
+ ,87
+ ,42
+ ,0
+ ,1335
+ ,65461
+ ,118612
+ ,54
+ ,72
+ ,43
+ ,46
+ ,2
+ ,1403
+ ,30080
+ ,72880
+ ,14
+ ,41
+ ,67
+ ,33
+ ,0
+ ,910
+ ,24094
+ ,65475
+ ,16
+ ,24
+ ,46
+ ,18
+ ,2
+ ,616
+ ,69008
+ ,99643
+ ,33
+ ,45
+ ,46
+ ,55
+ ,1
+ ,1407
+ ,54968
+ ,71965
+ ,32
+ ,33
+ ,56
+ ,35
+ ,1
+ ,771
+ ,46090
+ ,77272
+ ,21
+ ,27
+ ,48
+ ,59
+ ,2
+ ,766
+ ,27507
+ ,49289
+ ,15
+ ,36
+ ,44
+ ,19
+ ,1
+ ,473
+ ,10672
+ ,135131
+ ,38
+ ,87
+ ,60
+ ,66
+ ,0
+ ,1376
+ ,34029
+ ,108446
+ ,22
+ ,90
+ ,65
+ ,60
+ ,1
+ ,1232
+ ,46300
+ ,89746
+ ,28
+ ,114
+ ,55
+ ,36
+ ,3
+ ,1521
+ ,24760
+ ,44296
+ ,10
+ ,31
+ ,38
+ ,25
+ ,0
+ ,572
+ ,18779
+ ,77648
+ ,31
+ ,45
+ ,52
+ ,47
+ ,0
+ ,1059
+ ,21280
+ ,181528
+ ,32
+ ,69
+ ,60
+ ,54
+ ,0
+ ,1544
+ ,40662
+ ,134019
+ ,32
+ ,51
+ ,54
+ ,53
+ ,0
+ ,1230
+ ,28987
+ ,124064
+ ,43
+ ,34
+ ,86
+ ,40
+ ,1
+ ,1206
+ ,22827
+ ,92630
+ ,27
+ ,60
+ ,24
+ ,40
+ ,4
+ ,1205
+ ,18513
+ ,121848
+ ,37
+ ,45
+ ,52
+ ,39
+ ,0
+ ,1255
+ ,30594
+ ,52915
+ ,20
+ ,54
+ ,49
+ ,14
+ ,0
+ ,613
+ ,24006
+ ,81872
+ ,32
+ ,25
+ ,61
+ ,45
+ ,0
+ ,721
+ ,27913
+ ,58981
+ ,0
+ ,38
+ ,61
+ ,36
+ ,7
+ ,1109
+ ,42744
+ ,53515
+ ,5
+ ,52
+ ,81
+ ,28
+ ,2
+ ,740
+ ,12934
+ ,60812
+ ,26
+ ,67
+ ,43
+ ,44
+ ,0
+ ,1126
+ ,22574
+ ,56375
+ ,10
+ ,74
+ ,40
+ ,30
+ ,7
+ ,728
+ ,41385
+ ,65490
+ ,27
+ ,38
+ ,40
+ ,22
+ ,3
+ ,689
+ ,18653
+ ,80949
+ ,11
+ ,30
+ ,56
+ ,17
+ ,0
+ ,592
+ ,18472
+ ,76302
+ ,29
+ ,26
+ ,68
+ ,31
+ ,0
+ ,995
+ ,30976
+ ,104011
+ ,25
+ ,67
+ ,79
+ ,55
+ ,6
+ ,1613
+ ,63339
+ ,98104
+ ,55
+ ,132
+ ,47
+ ,54
+ ,2
+ ,2048
+ ,25568
+ ,67989
+ ,23
+ ,42
+ ,57
+ ,21
+ ,0
+ ,705
+ ,33747
+ ,30989
+ ,5
+ ,35
+ ,41
+ ,14
+ ,0
+ ,301
+ ,4154
+ ,135458
+ ,43
+ ,118
+ ,29
+ ,81
+ ,3
+ ,1803
+ ,19474
+ ,73504
+ ,23
+ ,68
+ ,3
+ ,35
+ ,0
+ ,799
+ ,35130
+ ,63123
+ ,34
+ ,43
+ ,60
+ ,43
+ ,1
+ ,861
+ ,39067
+ ,61254
+ ,36
+ ,76
+ ,30
+ ,46
+ ,1
+ ,1186
+ ,13310
+ ,74914
+ ,35
+ ,64
+ ,79
+ ,30
+ ,0
+ ,1451
+ ,65892
+ ,31774
+ ,0
+ ,48
+ ,47
+ ,23
+ ,1
+ ,628
+ ,4143
+ ,81437
+ ,37
+ ,64
+ ,40
+ ,38
+ ,0
+ ,1161
+ ,28579
+ ,87186
+ ,28
+ ,56
+ ,48
+ ,54
+ ,0
+ ,1463
+ ,51776
+ ,50090
+ ,16
+ ,71
+ ,36
+ ,20
+ ,0
+ ,742
+ ,21152
+ ,65745
+ ,26
+ ,75
+ ,42
+ ,53
+ ,0
+ ,979
+ ,38084
+ ,56653
+ ,38
+ ,39
+ ,49
+ ,45
+ ,0
+ ,675
+ ,27717
+ ,158399
+ ,23
+ ,42
+ ,57
+ ,39
+ ,0
+ ,1241
+ ,32928
+ ,46455
+ ,22
+ ,39
+ ,12
+ ,20
+ ,0
+ ,676
+ ,11342
+ ,73624
+ ,30
+ ,93
+ ,40
+ ,24
+ ,0
+ ,1049
+ ,19499
+ ,38395
+ ,16
+ ,38
+ ,43
+ ,31
+ ,0
+ ,620
+ ,16380
+ ,91899
+ ,18
+ ,60
+ ,33
+ ,35
+ ,0
+ ,1081
+ ,36874
+ ,139526
+ ,28
+ ,71
+ ,77
+ ,151
+ ,0
+ ,1688
+ ,48259
+ ,52164
+ ,32
+ ,52
+ ,43
+ ,52
+ ,0
+ ,736
+ ,16734
+ ,51567
+ ,21
+ ,27
+ ,45
+ ,30
+ ,2
+ ,617
+ ,28207
+ ,70551
+ ,23
+ ,59
+ ,47
+ ,31
+ ,0
+ ,812
+ ,30143
+ ,84856
+ ,29
+ ,40
+ ,43
+ ,29
+ ,1
+ ,1051
+ ,41369
+ ,102538
+ ,50
+ ,79
+ ,45
+ ,57
+ ,1
+ ,1656
+ ,45833
+ ,86678
+ ,12
+ ,44
+ ,50
+ ,40
+ ,0
+ ,705
+ ,29156
+ ,85709
+ ,21
+ ,65
+ ,35
+ ,44
+ ,0
+ ,945
+ ,35944
+ ,34662
+ ,18
+ ,10
+ ,7
+ ,25
+ ,0
+ ,554
+ ,36278
+ ,150580
+ ,27
+ ,124
+ ,71
+ ,77
+ ,0
+ ,1597
+ ,45588
+ ,99611
+ ,41
+ ,81
+ ,67
+ ,35
+ ,0
+ ,982
+ ,45097
+ ,19349
+ ,13
+ ,15
+ ,0
+ ,11
+ ,0
+ ,222
+ ,3895
+ ,99373
+ ,12
+ ,92
+ ,62
+ ,63
+ ,1
+ ,1212
+ ,28394
+ ,86230
+ ,21
+ ,42
+ ,54
+ ,44
+ ,0
+ ,1143
+ ,18632
+ ,30837
+ ,8
+ ,10
+ ,4
+ ,19
+ ,0
+ ,435
+ ,2325
+ ,31706
+ ,26
+ ,24
+ ,25
+ ,13
+ ,0
+ ,532
+ ,25139
+ ,89806
+ ,27
+ ,64
+ ,40
+ ,42
+ ,0
+ ,882
+ ,27975
+ ,62088
+ ,13
+ ,45
+ ,38
+ ,38
+ ,1
+ ,608
+ ,14483
+ ,40151
+ ,16
+ ,22
+ ,19
+ ,29
+ ,0
+ ,459
+ ,13127
+ ,27634
+ ,2
+ ,56
+ ,17
+ ,20
+ ,0
+ ,578
+ ,5839
+ ,76990
+ ,42
+ ,94
+ ,67
+ ,27
+ ,0
+ ,826
+ ,24069
+ ,37460
+ ,5
+ ,19
+ ,14
+ ,20
+ ,0
+ ,509
+ ,3738
+ ,54157
+ ,37
+ ,35
+ ,30
+ ,19
+ ,0
+ ,717
+ ,18625
+ ,49862
+ ,17
+ ,32
+ ,54
+ ,37
+ ,0
+ ,637
+ ,36341
+ ,84337
+ ,38
+ ,35
+ ,35
+ ,26
+ ,0
+ ,857
+ ,24548
+ ,64175
+ ,37
+ ,48
+ ,59
+ ,42
+ ,0
+ ,830
+ ,21792
+ ,59382
+ ,29
+ ,49
+ ,24
+ ,49
+ ,0
+ ,652
+ ,26263
+ ,119308
+ ,32
+ ,48
+ ,58
+ ,30
+ ,0
+ ,707
+ ,23686
+ ,76702
+ ,35
+ ,62
+ ,42
+ ,49
+ ,0
+ ,954
+ ,49303
+ ,103425
+ ,17
+ ,96
+ ,46
+ ,67
+ ,1
+ ,1461
+ ,25659
+ ,70344
+ ,20
+ ,45
+ ,61
+ ,28
+ ,0
+ ,672
+ ,28904
+ ,43410
+ ,7
+ ,63
+ ,3
+ ,19
+ ,0
+ ,778
+ ,2781
+ ,104838
+ ,46
+ ,71
+ ,52
+ ,49
+ ,1
+ ,1141
+ ,29236
+ ,62215
+ ,24
+ ,26
+ ,25
+ ,27
+ ,0
+ ,680
+ ,19546
+ ,69304
+ ,40
+ ,48
+ ,40
+ ,30
+ ,6
+ ,1090
+ ,22818
+ ,53117
+ ,3
+ ,29
+ ,32
+ ,22
+ ,3
+ ,616
+ ,32689
+ ,19764
+ ,10
+ ,19
+ ,4
+ ,12
+ ,1
+ ,285
+ ,5752
+ ,86680
+ ,37
+ ,45
+ ,49
+ ,31
+ ,2
+ ,1145
+ ,22197
+ ,84105
+ ,17
+ ,45
+ ,63
+ ,20
+ ,0
+ ,733
+ ,20055
+ ,77945
+ ,28
+ ,67
+ ,67
+ ,20
+ ,0
+ ,888
+ ,25272
+ ,89113
+ ,19
+ ,30
+ ,32
+ ,39
+ ,0
+ ,849
+ ,82206
+ ,91005
+ ,29
+ ,36
+ ,23
+ ,29
+ ,3
+ ,1182
+ ,32073
+ ,40248
+ ,8
+ ,34
+ ,7
+ ,16
+ ,1
+ ,528
+ ,5444
+ ,64187
+ ,10
+ ,36
+ ,54
+ ,27
+ ,0
+ ,642
+ ,20154
+ ,50857
+ ,15
+ ,34
+ ,37
+ ,21
+ ,0
+ ,947
+ ,36944
+ ,56613
+ ,15
+ ,37
+ ,35
+ ,19
+ ,1
+ ,819
+ ,8019
+ ,62792
+ ,28
+ ,46
+ ,51
+ ,35
+ ,0
+ ,757
+ ,30884
+ ,72535
+ ,17
+ ,44
+ ,39
+ ,14
+ ,0
+ ,894
+ ,19540)
+ ,dim=c(8
+ ,289)
+ ,dimnames=list(c('time_in_rfc'
+ ,'blogged_computations'
+ ,'compendium_views_pr'
+ ,'feedback_messages_p120'
+ ,'logins'
+ ,'shared_compendiums'
+ ,'pageviews'
+ ,'totsize
')
+ ,1:289))
> y <- array(NA,dim=c(8,289),dimnames=list(c('time_in_rfc','blogged_computations','compendium_views_pr','feedback_messages_p120','logins','shared_compendiums','pageviews','totsize
'),1:289))
> for (i in 1:dim(x)[1])
+ {
+ for (j in 1:dim(x)[2])
+ {
+ y[i,j] <- as.numeric(x[i,j])
+ }
+ }
> par3 = 'Linear Trend'
> par2 = 'Do not include Seasonal Dummies'
> par1 = '1'
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Prof. Dr. P. Wessa
> #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #Technical description: Write here your technical program description (don't use hard returns!)
> library(lattice)
> library(lmtest)
Loading required package: zoo
> 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 blogged_computations compendium_views_pr feedback_messages_p120
1 210907 79 81 94
2 120982 58 55 103
3 176508 60 50 93
4 179321 108 125 103
5 123185 49 40 51
6 52746 0 37 70
7 385534 121 63 91
8 33170 1 44 22
9 101645 20 88 38
10 149061 43 66 93
11 165446 69 57 60
12 237213 78 74 123
13 173326 86 49 148
14 133131 44 52 90
15 258873 104 88 124
16 180083 63 36 70
17 324799 158 108 168
18 230964 102 43 115
19 236785 77 75 71
20 135473 82 32 66
21 202925 115 44 134
22 215147 101 85 117
23 344297 80 86 108
24 153935 50 56 84
25 132943 83 50 156
26 174724 123 135 120
27 174415 73 63 114
28 225548 81 81 94
29 223632 105 52 120
30 124817 47 44 81
31 221698 105 113 110
32 210767 94 39 133
33 170266 44 73 122
34 260561 114 48 158
35 84853 38 33 109
36 294424 107 59 124
37 101011 30 41 39
38 215641 71 69 92
39 325107 84 64 126
40 7176 0 1 0
41 167542 59 59 70
42 106408 33 32 37
43 96560 42 129 38
44 265769 96 37 120
45 269651 106 31 93
46 149112 56 65 95
47 175824 57 107 77
48 152871 59 74 90
49 111665 39 54 80
50 116408 34 76 31
51 362301 76 715 110
52 78800 20 57 66
53 183167 91 66 138
54 277965 115 106 133
55 150629 85 54 113
56 168809 76 32 100
57 24188 8 20 7
58 329267 79 71 140
59 65029 21 21 61
60 101097 30 70 41
61 218946 76 112 96
62 244052 101 66 164
63 341570 94 190 78
64 103597 27 66 49
65 233328 92 165 102
66 256462 123 56 124
67 206161 75 61 99
68 311473 128 53 129
69 235800 105 127 62
70 177939 55 63 73
71 207176 56 38 114
72 196553 41 50 99
73 174184 72 52 70
74 143246 67 42 104
75 187559 75 76 116
76 187681 114 67 91
77 119016 118 50 74
78 182192 77 53 138
79 73566 22 39 67
80 194979 66 50 151
81 167488 69 77 72
82 143756 105 57 120
83 275541 116 73 115
84 243199 88 34 105
85 182999 73 39 104
86 135649 99 46 108
87 152299 62 63 98
88 120221 53 35 69
89 346485 118 106 111
90 145790 30 43 99
91 193339 100 47 71
92 80953 49 31 27
93 122774 24 162 69
94 130585 67 57 107
95 112611 46 36 73
96 286468 57 263 107
97 241066 75 78 93
98 148446 135 63 129
99 204713 68 54 69
100 182079 124 63 118
101 140344 33 77 73
102 220516 98 79 119
103 243060 58 110 104
104 162765 68 56 107
105 182613 81 56 99
106 232138 131 43 90
107 265318 110 111 197
108 85574 37 71 36
109 310839 130 62 85
110 225060 93 56 139
111 232317 118 74 106
112 144966 39 60 50
113 43287 13 43 64
114 155754 74 68 31
115 164709 81 53 63
116 201940 109 87 92
117 235454 151 46 106
118 220801 51 105 63
119 99466 28 32 69
120 92661 40 133 41
121 133328 56 79 56
122 61361 27 51 25
123 125930 37 207 65
124 100750 83 67 93
125 224549 54 47 114
126 82316 27 34 38
127 102010 28 66 44
128 101523 59 76 87
129 243511 133 65 110
130 22938 12 9 0
131 41566 0 42 27
132 152474 106 45 83
133 61857 23 25 30
134 99923 44 115 80
135 132487 71 97 98
136 317394 116 53 82
137 21054 4 2 0
138 209641 62 52 60
139 22648 12 44 28
140 31414 18 22 9
141 46698 14 35 33
142 131698 60 74 59
143 91735 7 103 49
144 244749 98 144 115
145 184510 64 60 140
146 79863 29 134 49
147 128423 32 89 120
148 97839 25 42 66
149 38214 16 52 21
150 151101 48 98 124
151 272458 100 99 152
152 172494 46 52 139
153 108043 45 29 38
154 328107 129 125 144
155 250579 130 106 120
156 351067 136 95 160
157 158015 59 40 114
158 98866 25 140 39
159 85439 32 43 78
160 229242 63 128 119
161 351619 95 142 141
162 84207 14 73 101
163 120445 36 72 56
164 324598 113 128 133
165 131069 47 61 83
166 204271 92 73 116
167 165543 70 148 90
168 141722 19 64 36
169 116048 50 45 50
170 250047 41 58 61
171 299775 91 97 97
172 195838 111 50 98
173 173260 41 37 78
174 254488 120 50 117
175 104389 135 105 148
176 136084 27 69 41
177 199476 87 46 105
178 92499 25 57 55
179 224330 131 52 132
180 135781 45 98 44
181 74408 29 61 21
182 81240 58 89 50
183 14688 4 0 0
184 181633 47 48 73
185 271856 109 91 86
186 7199 7 0 0
187 46660 12 7 13
188 17547 0 3 4
189 133368 37 54 57
190 95227 37 70 48
191 152601 46 36 46
192 98146 15 37 48
193 79619 42 123 32
194 59194 7 247 68
195 139942 54 46 87
196 118612 54 72 43
197 72880 14 41 67
198 65475 16 24 46
199 99643 33 45 46
200 71965 32 33 56
201 77272 21 27 48
202 49289 15 36 44
203 135131 38 87 60
204 108446 22 90 65
205 89746 28 114 55
206 44296 10 31 38
207 77648 31 45 52
208 181528 32 69 60
209 134019 32 51 54
210 124064 43 34 86
211 92630 27 60 24
212 121848 37 45 52
213 52915 20 54 49
214 81872 32 25 61
215 58981 0 38 61
216 53515 5 52 81
217 60812 26 67 43
218 56375 10 74 40
219 65490 27 38 40
220 80949 11 30 56
221 76302 29 26 68
222 104011 25 67 79
223 98104 55 132 47
224 67989 23 42 57
225 30989 5 35 41
226 135458 43 118 29
227 73504 23 68 3
228 63123 34 43 60
229 61254 36 76 30
230 74914 35 64 79
231 31774 0 48 47
232 81437 37 64 40
233 87186 28 56 48
234 50090 16 71 36
235 65745 26 75 42
236 56653 38 39 49
237 158399 23 42 57
238 46455 22 39 12
239 73624 30 93 40
240 38395 16 38 43
241 91899 18 60 33
242 139526 28 71 77
243 52164 32 52 43
244 51567 21 27 45
245 70551 23 59 47
246 84856 29 40 43
247 102538 50 79 45
248 86678 12 44 50
249 85709 21 65 35
250 34662 18 10 7
251 150580 27 124 71
252 99611 41 81 67
253 19349 13 15 0
254 99373 12 92 62
255 86230 21 42 54
256 30837 8 10 4
257 31706 26 24 25
258 89806 27 64 40
259 62088 13 45 38
260 40151 16 22 19
261 27634 2 56 17
262 76990 42 94 67
263 37460 5 19 14
264 54157 37 35 30
265 49862 17 32 54
266 84337 38 35 35
267 64175 37 48 59
268 59382 29 49 24
269 119308 32 48 58
270 76702 35 62 42
271 103425 17 96 46
272 70344 20 45 61
273 43410 7 63 3
274 104838 46 71 52
275 62215 24 26 25
276 69304 40 48 40
277 53117 3 29 32
278 19764 10 19 4
279 86680 37 45 49
280 84105 17 45 63
281 77945 28 67 67
282 89113 19 30 32
283 91005 29 36 23
284 40248 8 34 7
285 64187 10 36 54
286 50857 15 34 37
287 56613 15 37 35
288 62792 28 46 51
289 72535 17 44 39
logins shared_compendiums pageviews totsize\r\r\r t
1 56 3 1418 112285 1
2 56 4 869 84786 2
3 54 12 1530 83123 3
4 89 2 2172 101193 4
5 40 1 901 38361 5
6 25 3 463 68504 6
7 92 0 3201 119182 7
8 18 0 371 22807 8
9 63 0 1192 17140 9
10 44 5 1583 116174 10
11 33 0 1439 57635 11
12 84 0 1764 66198 12
13 88 7 1495 71701 13
14 55 7 1373 57793 14
15 60 3 2187 80444 15
16 66 9 1491 53855 16
17 154 0 4041 97668 17
18 53 4 1706 133824 18
19 119 3 2152 101481 19
20 41 0 1036 99645 20
21 61 7 1882 114789 21
22 58 0 1929 99052 22
23 75 1 2242 67654 23
24 33 5 1220 65553 24
25 40 7 1289 97500 25
26 92 0 2515 69112 26
27 100 0 2147 82753 27
28 112 5 2352 85323 28
29 73 0 1638 72654 29
30 40 0 1222 30727 30
31 45 0 1812 77873 31
32 60 3 1677 117478 32
33 62 4 1579 74007 33
34 75 1 1731 90183 34
35 31 4 807 61542 35
36 77 2 2452 101494 36
37 34 0 829 27570 37
38 46 0 1940 55813 38
39 99 0 2662 79215 39
40 17 0 186 1423 40
41 66 2 1499 55461 41
42 30 1 865 31081 42
43 76 0 1793 22996 43
44 146 2 2527 83122 44
45 67 10 2747 70106 45
46 56 6 1324 60578 46
47 107 0 2702 39992 47
48 58 5 1383 79892 48
49 34 4 1179 49810 49
50 61 1 2099 71570 50
51 119 2 4308 100708 51
52 42 2 918 33032 52
53 66 0 1831 82875 53
54 89 8 3373 139077 54
55 44 3 1713 71595 55
56 66 0 1438 72260 56
57 24 0 496 5950 57
58 259 8 2253 115762 58
59 17 5 744 32551 59
60 64 3 1161 31701 60
61 41 1 2352 80670 61
62 68 5 2144 143558 62
63 168 1 4691 117105 63
64 43 1 1112 23789 64
65 132 5 2694 120733 65
66 105 0 1973 105195 66
67 71 12 1769 73107 67
68 112 8 3148 132068 68
69 94 8 2474 149193 69
70 82 8 2084 46821 70
71 70 8 1954 87011 71
72 57 2 1226 95260 72
73 53 0 1389 55183 73
74 103 5 1496 106671 74
75 121 8 2269 73511 75
76 62 2 1833 92945 76
77 52 5 1268 78664 77
78 52 12 1943 70054 78
79 32 6 893 22618 79
80 62 7 1762 74011 80
81 45 2 1403 83737 81
82 46 0 1425 69094 82
83 63 4 1857 93133 83
84 75 3 1840 95536 84
85 88 6 1502 225920 85
86 46 2 1441 62133 86
87 53 0 1420 61370 87
88 37 1 1416 43836 88
89 90 0 2970 106117 89
90 63 5 1317 38692 90
91 78 2 1644 84651 91
92 25 0 870 56622 92
93 45 0 1654 15986 93
94 46 5 1054 95364 94
95 41 0 937 26706 95
96 144 1 3004 89691 96
97 82 0 2008 67267 97
98 91 1 2547 126846 98
99 71 1 1885 41140 99
100 63 2 1626 102860 100
101 53 6 1468 51715 101
102 62 1 2445 55801 102
103 63 4 1964 111813 103
104 32 2 1381 120293 104
105 39 3 1369 138599 105
106 62 0 1659 161647 106
107 117 10 2888 115929 107
108 34 0 1290 24266 108
109 92 9 2845 162901 109
110 93 7 1982 109825 110
111 54 0 1904 129838 111
112 144 0 1391 37510 112
113 14 4 602 43750 113
114 61 4 1743 40652 114
115 109 0 1559 87771 115
116 38 0 2014 85872 116
117 73 0 2143 89275 117
118 75 1 2146 44418 118
119 50 0 874 192565 119
120 61 1 1590 35232 120
121 55 0 1590 40909 121
122 77 0 1210 13294 122
123 75 4 2072 32387 123
124 72 0 1281 140867 124
125 50 4 1401 120662 125
126 32 4 834 21233 126
127 53 3 1105 44332 127
128 42 0 1272 61056 128
129 71 0 1944 101338 129
130 10 0 391 1168 130
131 35 5 761 13497 131
132 65 0 1605 65567 132
133 25 4 530 25162 133
134 66 0 1988 32334 134
135 41 0 1386 40735 135
136 86 1 2395 91413 136
137 16 0 387 855 137
138 42 5 1742 97068 138
139 19 0 620 44339 139
140 19 0 449 14116 140
141 45 0 800 10288 141
142 65 0 1684 65622 142
143 35 0 1050 16563 143
144 95 2 2699 76643 144
145 49 7 1606 110681 145
146 37 1 1502 29011 146
147 64 8 1204 92696 147
148 38 2 1138 94785 148
149 34 0 568 8773 149
150 32 2 1459 83209 150
151 65 0 2158 93815 151
152 52 0 1111 86687 152
153 62 1 1421 34553 153
154 65 3 2833 105547 154
155 83 0 1955 103487 155
156 95 3 2922 213688 156
157 29 0 1002 71220 157
158 18 0 1060 23517 158
159 33 0 956 56926 159
160 247 4 2186 91721 160
161 139 4 3604 115168 161
162 29 11 1035 111194 162
163 118 0 1417 51009 163
164 110 0 3261 135777 164
165 67 4 1587 51513 165
166 42 0 1424 74163 166
167 65 1 1701 51633 167
168 94 0 1249 75345 168
169 64 0 946 33416 169
170 81 0 1926 83305 170
171 95 9 3352 98952 171
172 67 1 1641 102372 172
173 63 3 2035 37238 173
174 83 10 2312 103772 174
175 45 5 1369 123969 175
176 30 0 1577 27142 176
177 70 2 2201 135400 177
178 32 0 961 21399 178
179 83 1 1900 130115 179
180 31 2 1254 24874 180
181 67 4 1335 34988 181
182 66 0 1597 45549 182
183 10 0 207 6023 183
184 70 2 1645 64466 184
185 103 1 2429 54990 185
186 5 0 151 1644 186
187 20 0 474 6179 187
188 5 0 141 3926 188
189 36 1 1639 32755 189
190 34 0 872 34777 190
191 48 2 1318 73224 191
192 40 0 1018 27114 192
193 43 3 1383 20760 193
194 31 6 1314 37636 194
195 42 0 1335 65461 195
196 46 2 1403 30080 196
197 33 0 910 24094 197
198 18 2 616 69008 198
199 55 1 1407 54968 199
200 35 1 771 46090 200
201 59 2 766 27507 201
202 19 1 473 10672 202
203 66 0 1376 34029 203
204 60 1 1232 46300 204
205 36 3 1521 24760 205
206 25 0 572 18779 206
207 47 0 1059 21280 207
208 54 0 1544 40662 208
209 53 0 1230 28987 209
210 40 1 1206 22827 210
211 40 4 1205 18513 211
212 39 0 1255 30594 212
213 14 0 613 24006 213
214 45 0 721 27913 214
215 36 7 1109 42744 215
216 28 2 740 12934 216
217 44 0 1126 22574 217
218 30 7 728 41385 218
219 22 3 689 18653 219
220 17 0 592 18472 220
221 31 0 995 30976 221
222 55 6 1613 63339 222
223 54 2 2048 25568 223
224 21 0 705 33747 224
225 14 0 301 4154 225
226 81 3 1803 19474 226
227 35 0 799 35130 227
228 43 1 861 39067 228
229 46 1 1186 13310 229
230 30 0 1451 65892 230
231 23 1 628 4143 231
232 38 0 1161 28579 232
233 54 0 1463 51776 233
234 20 0 742 21152 234
235 53 0 979 38084 235
236 45 0 675 27717 236
237 39 0 1241 32928 237
238 20 0 676 11342 238
239 24 0 1049 19499 239
240 31 0 620 16380 240
241 35 0 1081 36874 241
242 151 0 1688 48259 242
243 52 0 736 16734 243
244 30 2 617 28207 244
245 31 0 812 30143 245
246 29 1 1051 41369 246
247 57 1 1656 45833 247
248 40 0 705 29156 248
249 44 0 945 35944 249
250 25 0 554 36278 250
251 77 0 1597 45588 251
252 35 0 982 45097 252
253 11 0 222 3895 253
254 63 1 1212 28394 254
255 44 0 1143 18632 255
256 19 0 435 2325 256
257 13 0 532 25139 257
258 42 0 882 27975 258
259 38 1 608 14483 259
260 29 0 459 13127 260
261 20 0 578 5839 261
262 27 0 826 24069 262
263 20 0 509 3738 263
264 19 0 717 18625 264
265 37 0 637 36341 265
266 26 0 857 24548 266
267 42 0 830 21792 267
268 49 0 652 26263 268
269 30 0 707 23686 269
270 49 0 954 49303 270
271 67 1 1461 25659 271
272 28 0 672 28904 272
273 19 0 778 2781 273
274 49 1 1141 29236 274
275 27 0 680 19546 275
276 30 6 1090 22818 276
277 22 3 616 32689 277
278 12 1 285 5752 278
279 31 2 1145 22197 279
280 20 0 733 20055 280
281 20 0 888 25272 281
282 39 0 849 82206 282
283 29 3 1182 32073 283
284 16 1 528 5444 284
285 27 0 642 20154 285
286 21 0 947 36944 286
287 19 1 819 8019 287
288 35 0 757 30884 288
289 14 0 894 19540 289
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) blogged_computations compendium_views_pr
857.0455 303.7885 -105.3452
feedback_messages_p120 logins shared_compendiums
315.8860 145.5679 -1527.7978
pageviews `totsize\r\r\r` t
66.2646 0.2503 -53.5537
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-134526 -13771 341 15298 120427
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 8.570e+02 7.068e+03 0.121 0.903579
blogged_computations 3.038e+02 1.001e+02 3.036 0.002622 **
compendium_views_pr -1.053e+02 4.302e+01 -2.449 0.014949 *
feedback_messages_p120 3.159e+02 8.001e+01 3.948 9.97e-05 ***
logins 1.456e+02 7.874e+01 1.849 0.065560 .
shared_compendiums -1.528e+03 7.263e+02 -2.103 0.036321 *
pageviews 6.626e+01 5.548e+00 11.944 < 2e-16 ***
`totsize\r\r\r` 2.503e-01 7.348e-02 3.406 0.000757 ***
t -5.355e+01 2.521e+01 -2.125 0.034488 *
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 29010 on 280 degrees of freedom
Multiple R-squared: 0.8793, Adjusted R-squared: 0.8759
F-statistic: 255 on 8 and 280 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.8026480 3.947039e-01 1.973520e-01
[2,] 0.6813588 6.372824e-01 3.186412e-01
[3,] 0.5664261 8.671477e-01 4.335739e-01
[4,] 0.4426078 8.852155e-01 5.573922e-01
[5,] 0.5426066 9.147868e-01 4.573934e-01
[6,] 0.9097224 1.805553e-01 9.027764e-02
[7,] 0.8708865 2.582269e-01 1.291135e-01
[8,] 0.8579714 2.840571e-01 1.420286e-01
[9,] 0.8895884 2.208232e-01 1.104116e-01
[10,] 0.8722575 2.554851e-01 1.277425e-01
[11,] 0.8457848 3.084305e-01 1.542152e-01
[12,] 0.9992363 1.527413e-03 7.637066e-04
[13,] 0.9987298 2.540438e-03 1.270219e-03
[14,] 0.9987319 2.536156e-03 1.268078e-03
[15,] 0.9995995 8.010495e-04 4.005247e-04
[16,] 0.9997623 4.753800e-04 2.376900e-04
[17,] 0.9996037 7.925020e-04 3.962510e-04
[18,] 0.9995992 8.016203e-04 4.008102e-04
[19,] 0.9994042 1.191692e-03 5.958458e-04
[20,] 0.9994487 1.102540e-03 5.512700e-04
[21,] 0.9991086 1.782775e-03 8.913877e-04
[22,] 0.9987348 2.530480e-03 1.265240e-03
[23,] 0.9991637 1.672566e-03 8.362831e-04
[24,] 0.9990330 1.933975e-03 9.669873e-04
[25,] 0.9988027 2.394662e-03 1.197331e-03
[26,] 0.9982586 3.482897e-03 1.741449e-03
[27,] 0.9975478 4.904351e-03 2.452176e-03
[28,] 0.9985692 2.861645e-03 1.430822e-03
[29,] 0.9987660 2.467985e-03 1.233993e-03
[30,] 0.9983005 3.398958e-03 1.699479e-03
[31,] 0.9976672 4.665542e-03 2.332771e-03
[32,] 0.9975619 4.876225e-03 2.438112e-03
[33,] 0.9964846 7.030890e-03 3.515445e-03
[34,] 0.9958735 8.253051e-03 4.126525e-03
[35,] 0.9950342 9.931613e-03 4.965807e-03
[36,] 0.9967258 6.548342e-03 3.274171e-03
[37,] 0.9958564 8.287150e-03 4.143575e-03
[38,] 0.9943510 1.129800e-02 5.649001e-03
[39,] 0.9974769 5.046204e-03 2.523102e-03
[40,] 0.9998609 2.782776e-04 1.391388e-04
[41,] 0.9997920 4.159889e-04 2.079944e-04
[42,] 0.9997657 4.686139e-04 2.343070e-04
[43,] 0.9998531 2.937282e-04 1.468641e-04
[44,] 0.9998607 2.785973e-04 1.392986e-04
[45,] 0.9998016 3.967474e-04 1.983737e-04
[46,] 0.9997141 5.718778e-04 2.859389e-04
[47,] 0.9999529 9.422995e-05 4.711498e-05
[48,] 0.9999289 1.421343e-04 7.106713e-05
[49,] 0.9998931 2.138768e-04 1.069384e-04
[50,] 0.9998481 3.038168e-04 1.519084e-04
[51,] 0.9997793 4.414870e-04 2.207435e-04
[52,] 0.9998595 2.810431e-04 1.405216e-04
[53,] 0.9998038 3.923056e-04 1.961528e-04
[54,] 0.9997866 4.267453e-04 2.133727e-04
[55,] 0.9997214 5.571724e-04 2.785862e-04
[56,] 0.9997459 5.081415e-04 2.540707e-04
[57,] 0.9996424 7.151755e-04 3.575877e-04
[58,] 0.9995081 9.838168e-04 4.919084e-04
[59,] 0.9993089 1.382137e-03 6.910685e-04
[60,] 0.9991501 1.699775e-03 8.498877e-04
[61,] 0.9995896 8.208572e-04 4.104286e-04
[62,] 0.9995477 9.046319e-04 4.523159e-04
[63,] 0.9996263 7.473914e-04 3.736957e-04
[64,] 0.9996474 7.051547e-04 3.525773e-04
[65,] 0.9995450 9.100231e-04 4.550115e-04
[66,] 0.9995827 8.346358e-04 4.173179e-04
[67,] 0.9994288 1.142332e-03 5.711659e-04
[68,] 0.9992091 1.581849e-03 7.909244e-04
[69,] 0.9989326 2.134753e-03 1.067377e-03
[70,] 0.9987713 2.457411e-03 1.228706e-03
[71,] 0.9987509 2.498203e-03 1.249102e-03
[72,] 0.9996941 6.118130e-04 3.059065e-04
[73,] 0.9997663 4.674123e-04 2.337062e-04
[74,] 0.9997411 5.177935e-04 2.588967e-04
[75,] 0.9997490 5.019089e-04 2.509545e-04
[76,] 0.9996488 7.024251e-04 3.512126e-04
[77,] 0.9995535 8.930740e-04 4.465370e-04
[78,] 0.9998254 3.492482e-04 1.746241e-04
[79,] 0.9997824 4.352234e-04 2.176117e-04
[80,] 0.9997125 5.750584e-04 2.875292e-04
[81,] 0.9996078 7.844997e-04 3.922498e-04
[82,] 0.9994588 1.082489e-03 5.412447e-04
[83,] 0.9992682 1.463690e-03 7.318452e-04
[84,] 0.9990499 1.900273e-03 9.501365e-04
[85,] 0.9991119 1.776229e-03 8.881146e-04
[86,] 0.9993359 1.328119e-03 6.640594e-04
[87,] 0.9999996 8.626920e-07 4.313460e-07
[88,] 0.9999996 8.429039e-07 4.214520e-07
[89,] 0.9999994 1.109019e-06 5.545095e-07
[90,] 0.9999992 1.545155e-06 7.725777e-07
[91,] 0.9999990 1.964914e-06 9.824572e-07
[92,] 0.9999996 7.124137e-07 3.562068e-07
[93,] 0.9999994 1.114142e-06 5.570709e-07
[94,] 0.9999992 1.539787e-06 7.698937e-07
[95,] 0.9999990 2.042600e-06 1.021300e-06
[96,] 0.9999992 1.574929e-06 7.874643e-07
[97,] 0.9999990 2.051676e-06 1.025838e-06
[98,] 0.9999990 2.082021e-06 1.041010e-06
[99,] 0.9999984 3.172007e-06 1.586003e-06
[100,] 0.9999978 4.473893e-06 2.236946e-06
[101,] 0.9999967 6.506268e-06 3.253134e-06
[102,] 0.9999956 8.718464e-06 4.359232e-06
[103,] 0.9999938 1.246910e-05 6.234552e-06
[104,] 0.9999910 1.792110e-05 8.960551e-06
[105,] 0.9999870 2.596574e-05 1.298287e-05
[106,] 0.9999830 3.408323e-05 1.704162e-05
[107,] 0.9999888 2.234907e-05 1.117453e-05
[108,] 0.9999887 2.259257e-05 1.129628e-05
[109,] 0.9999885 2.294458e-05 1.147229e-05
[110,] 0.9999838 3.249072e-05 1.624536e-05
[111,] 0.9999871 2.573186e-05 1.286593e-05
[112,] 0.9999843 3.141846e-05 1.570923e-05
[113,] 0.9999976 4.778673e-06 2.389337e-06
[114,] 0.9999995 1.043523e-06 5.217616e-07
[115,] 0.9999993 1.416762e-06 7.083811e-07
[116,] 0.9999990 2.073775e-06 1.036887e-06
[117,] 0.9999991 1.759447e-06 8.797236e-07
[118,] 0.9999988 2.310344e-06 1.155172e-06
[119,] 0.9999983 3.494644e-06 1.747322e-06
[120,] 0.9999974 5.266982e-06 2.633491e-06
[121,] 0.9999977 4.626703e-06 2.313352e-06
[122,] 0.9999971 5.853882e-06 2.926941e-06
[123,] 0.9999997 6.952090e-07 3.476045e-07
[124,] 0.9999995 9.791189e-07 4.895595e-07
[125,] 1.0000000 7.289551e-08 3.644776e-08
[126,] 0.9999999 1.160952e-07 5.804759e-08
[127,] 1.0000000 3.933056e-08 1.966528e-08
[128,] 1.0000000 3.485228e-08 1.742614e-08
[129,] 1.0000000 5.633664e-08 2.816832e-08
[130,] 1.0000000 6.698680e-08 3.349340e-08
[131,] 1.0000000 6.009510e-08 3.004755e-08
[132,] 1.0000000 8.399503e-08 4.199752e-08
[133,] 0.9999999 1.268164e-07 6.340821e-08
[134,] 0.9999999 1.954728e-07 9.773639e-08
[135,] 0.9999999 1.602550e-07 8.012750e-08
[136,] 0.9999999 2.421748e-07 1.210874e-07
[137,] 0.9999998 3.018233e-07 1.509116e-07
[138,] 0.9999998 4.562874e-07 2.281437e-07
[139,] 0.9999997 6.957098e-07 3.478549e-07
[140,] 0.9999997 5.602047e-07 2.801024e-07
[141,] 0.9999997 6.202345e-07 3.101172e-07
[142,] 0.9999997 6.627639e-07 3.313820e-07
[143,] 0.9999999 2.894049e-07 1.447025e-07
[144,] 0.9999998 3.013428e-07 1.506714e-07
[145,] 0.9999999 2.879696e-07 1.439848e-07
[146,] 0.9999999 2.675885e-07 1.337942e-07
[147,] 0.9999999 2.702299e-07 1.351150e-07
[148,] 0.9999998 3.339899e-07 1.669950e-07
[149,] 0.9999997 5.338862e-07 2.669431e-07
[150,] 0.9999997 6.050862e-07 3.025431e-07
[151,] 0.9999996 8.648947e-07 4.324474e-07
[152,] 0.9999995 1.007450e-06 5.037248e-07
[153,] 0.9999992 1.520466e-06 7.602328e-07
[154,] 0.9999990 2.010977e-06 1.005488e-06
[155,] 0.9999995 9.632174e-07 4.816087e-07
[156,] 0.9999994 1.251632e-06 6.258162e-07
[157,] 0.9999993 1.452232e-06 7.261162e-07
[158,] 0.9999990 1.962459e-06 9.812295e-07
[159,] 1.0000000 6.857280e-08 3.428640e-08
[160,] 1.0000000 7.153306e-08 3.576653e-08
[161,] 0.9999999 1.004740e-07 5.023702e-08
[162,] 0.9999999 1.639423e-07 8.197115e-08
[163,] 1.0000000 7.391583e-08 3.695792e-08
[164,] 1.0000000 5.291284e-09 2.645642e-09
[165,] 1.0000000 6.609671e-09 3.304836e-09
[166,] 1.0000000 8.109397e-09 4.054698e-09
[167,] 1.0000000 1.346114e-08 6.730571e-09
[168,] 1.0000000 2.120307e-08 1.060154e-08
[169,] 1.0000000 8.352994e-09 4.176497e-09
[170,] 1.0000000 9.696424e-09 4.848212e-09
[171,] 1.0000000 6.381423e-10 3.190712e-10
[172,] 1.0000000 1.134803e-09 5.674017e-10
[173,] 1.0000000 8.634431e-10 4.317215e-10
[174,] 1.0000000 1.033885e-10 5.169427e-11
[175,] 1.0000000 1.752413e-10 8.762067e-11
[176,] 1.0000000 3.188824e-10 1.594412e-10
[177,] 1.0000000 5.757686e-10 2.878843e-10
[178,] 1.0000000 1.003423e-09 5.017114e-10
[179,] 1.0000000 1.482370e-09 7.411850e-10
[180,] 1.0000000 5.329343e-10 2.664671e-10
[181,] 1.0000000 8.454807e-10 4.227404e-10
[182,] 1.0000000 1.425859e-09 7.129293e-10
[183,] 1.0000000 2.407034e-09 1.203517e-09
[184,] 1.0000000 3.723410e-09 1.861705e-09
[185,] 1.0000000 5.503356e-09 2.751678e-09
[186,] 1.0000000 9.123764e-09 4.561882e-09
[187,] 1.0000000 1.556371e-08 7.781855e-09
[188,] 1.0000000 2.387216e-08 1.193608e-08
[189,] 1.0000000 4.237646e-08 2.118823e-08
[190,] 1.0000000 7.402676e-08 3.701338e-08
[191,] 0.9999999 1.276420e-07 6.382100e-08
[192,] 0.9999999 1.426127e-07 7.130637e-08
[193,] 0.9999999 2.358677e-07 1.179338e-07
[194,] 0.9999998 3.666916e-07 1.833458e-07
[195,] 0.9999997 6.026002e-07 3.013001e-07
[196,] 0.9999996 8.779492e-07 4.389746e-07
[197,] 1.0000000 2.889177e-08 1.444589e-08
[198,] 1.0000000 1.109052e-08 5.545262e-09
[199,] 1.0000000 1.132810e-08 5.664050e-09
[200,] 1.0000000 1.138417e-08 5.692084e-09
[201,] 1.0000000 5.360419e-09 2.680209e-09
[202,] 1.0000000 1.043743e-08 5.218714e-09
[203,] 1.0000000 1.551091e-08 7.755454e-09
[204,] 1.0000000 2.396202e-08 1.198101e-08
[205,] 1.0000000 3.809559e-08 1.904780e-08
[206,] 1.0000000 4.551524e-08 2.275762e-08
[207,] 1.0000000 8.436733e-08 4.218367e-08
[208,] 0.9999999 1.084483e-07 5.422415e-08
[209,] 1.0000000 6.103934e-08 3.051967e-08
[210,] 0.9999999 1.135204e-07 5.676019e-08
[211,] 0.9999999 1.918536e-07 9.592681e-08
[212,] 1.0000000 9.082595e-08 4.541297e-08
[213,] 0.9999999 1.575523e-07 7.877615e-08
[214,] 0.9999999 2.927092e-07 1.463546e-07
[215,] 0.9999998 3.425875e-07 1.712938e-07
[216,] 0.9999999 2.978407e-07 1.489204e-07
[217,] 0.9999997 5.062589e-07 2.531294e-07
[218,] 0.9999996 7.100118e-07 3.550059e-07
[219,] 1.0000000 7.400682e-08 3.700341e-08
[220,] 1.0000000 7.932638e-08 3.966319e-08
[221,] 0.9999999 1.550363e-07 7.751813e-08
[222,] 1.0000000 7.501654e-08 3.750827e-08
[223,] 0.9999999 1.274761e-07 6.373807e-08
[224,] 0.9999999 1.674849e-07 8.374247e-08
[225,] 0.9999999 2.909156e-07 1.454578e-07
[226,] 1.0000000 1.191653e-09 5.958263e-10
[227,] 1.0000000 2.828302e-09 1.414151e-09
[228,] 1.0000000 6.159047e-09 3.079524e-09
[229,] 1.0000000 6.812864e-09 3.406432e-09
[230,] 1.0000000 1.232355e-08 6.161773e-09
[231,] 1.0000000 2.698261e-08 1.349131e-08
[232,] 1.0000000 2.874906e-08 1.437453e-08
[233,] 1.0000000 5.809263e-08 2.904631e-08
[234,] 0.9999999 1.340533e-07 6.702665e-08
[235,] 0.9999998 3.034887e-07 1.517444e-07
[236,] 0.9999998 3.094724e-07 1.547362e-07
[237,] 0.9999998 3.819331e-07 1.909666e-07
[238,] 0.9999996 8.078442e-07 4.039221e-07
[239,] 0.9999992 1.649223e-06 8.246115e-07
[240,] 0.9999991 1.894189e-06 9.470945e-07
[241,] 0.9999979 4.112367e-06 2.056183e-06
[242,] 0.9999954 9.179102e-06 4.589551e-06
[243,] 0.9999899 2.027985e-05 1.013992e-05
[244,] 0.9999779 4.417161e-05 2.208581e-05
[245,] 0.9999531 9.379897e-05 4.689949e-05
[246,] 0.9999403 1.193483e-04 5.967417e-05
[247,] 0.9999150 1.700955e-04 8.504774e-05
[248,] 0.9998470 3.060888e-04 1.530444e-04
[249,] 0.9996777 6.445862e-04 3.222931e-04
[250,] 0.9994319 1.136196e-03 5.680981e-04
[251,] 0.9990279 1.944137e-03 9.720683e-04
[252,] 0.9980765 3.846987e-03 1.923494e-03
[253,] 0.9971994 5.601195e-03 2.800598e-03
[254,] 0.9974298 5.140481e-03 2.570241e-03
[255,] 0.9949007 1.019857e-02 5.099283e-03
[256,] 0.9970569 5.886180e-03 2.943090e-03
[257,] 0.9945902 1.081958e-02 5.409790e-03
[258,] 0.9994968 1.006438e-03 5.032188e-04
[259,] 0.9991689 1.662265e-03 8.311323e-04
[260,] 0.9978952 4.209536e-03 2.104768e-03
[261,] 0.9948092 1.038153e-02 5.190763e-03
[262,] 0.9911931 1.761377e-02 8.806885e-03
[263,] 0.9790116 4.197670e-02 2.098835e-02
[264,] 0.9497541 1.004917e-01 5.024587e-02
[265,] 0.8925233 2.149535e-01 1.074767e-01
[266,] 0.8759831 2.480338e-01 1.240169e-01
> postscript(file="/var/wessaorg/rcomp/tmp/1de2p1324466480.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/264xe1324466480.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/3q5ap1324466480.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/4795c1324466480.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/5juzk1324466480.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> qqnorm(mysum$resid, main='Residual Normal Q-Q Plot')
> qqline(mysum$resid)
> grid()
> dev.off()
null device
1
> (myerror <- as.ts(mysum$resid))
Time Series:
Start = 1
End = 289
Frequency = 1
1 2 3 4 5
39311.89863 -4972.98917 21759.92248 -52650.25936 22214.18449
6 7 8 9 10
-12884.05949 70852.89528 -2788.76485 13.24984 -19484.53103
11 12 13 14 15
16685.95385 36559.53426 -13674.68602 -6050.96424 28122.54285
16 17 18 19 20
34489.00676 -79472.81624 20145.41417 18289.56003 -6256.57013
21 22 23 24 25
-21058.69191 -4274.86990 120426.53329 24125.90655 -40744.09666
26 27 28 29 30
-83134.99463 -54083.51388 -5448.57143 22651.63412 -4150.84535
31 32 33 34 35
21649.41276 486.78876 894.87108 35376.51515 -23907.22038
36 37 38 39 40
34001.26317 18239.07291 24240.14611 57128.89000 -6589.53964
41 42 43 44 45
15298.32456 21522.02447 -48798.04363 -2353.94603 18841.98645
46 47 48 49 50
8664.97771 -57514.33608 3585.19881 -7427.40490 -58238.70703
51 52 53 54 55
56723.59676 -12347.54869 -30814.14673 -44831.05397 -36361.26423
56 57 58 59 60
-3334.23601 -14001.00292 67027.86879 -8387.45094 -837.26567
61 62 63 64 65
-741.83980 -9277.54637 -52174.27831 5068.19553 -27142.39072
66 67 68 69 70
16152.25866 33741.72169 -5526.75993 -2202.30423 -1827.58717
71 72 73 74 75
11877.50789 50769.15196 25162.84293 -35611.31504 -34844.72236
76 77 78 79 80
-16118.33614 -35313.21680 -11408.28481 -7125.58562 2313.89018
81 82 83 84 85
17955.01964 -34923.67530 65832.39974 38350.94122 -23937.45809
86 87 88 89 90
-34624.47622 -4223.08798 -18792.41659 54186.68826 15410.96423
91 92 93 94 95
11077.39931 -10585.54936 -5276.53501 -6152.71367 8858.39833
96 97 98 99 100
26403.29203 39628.97342 -134525.99604 28380.69994 -21333.35190
101 102 103 104 105
11155.85303 -17397.07988 47649.56301 -4298.99524 10903.14281
106 107 108 109 110
13849.58930 -35899.70759 -21135.24029 27075.53204 2168.30495
111 112 113 114 115
9347.35468 6262.56626 -17921.48319 7451.54542 -10052.69102
116 117 118 119 120
-6192.99932 -8620.75945 39220.86161 -35334.19029 -34391.25173
121 122 123 124 125
-11033.35493 -38405.16575 -28518.05125 -71619.23916 58721.09947
126 127 128 129 130
12457.08959 5053.40988 -35560.79787 16744.00235 -1311.85449
131 132 133 134 135
-7640.92877 -27218.89093 15347.24417 -69714.11314 -11452.95326
136 137 138 139 140
75689.71901 -1658.10440 45664.29384 -33566.25443 -3990.35253
141 142 143 144 145
-19735.06652 -28097.47854 12964.00702 -8126.90361 3514.30084
146 147 148 149 150
-33995.14486 -2888.06724 -20716.72490 -5462.93781 -4257.33307
151 152 153 154 155
35784.21577 24488.90835 -17546.08380 44968.81845 24262.52367
156 157 158 159 160
20366.97621 27403.10064 22558.45456 -19132.46386 -3943.91461
161 162 163 164 165
19181.04977 -20268.57120 -16563.36847 3586.21046 -16717.60936
166 167 168 169 170
36368.12609 5953.73709 24156.52882 17633.01975 72415.25487
171 172 173 174 175
13047.05478 2197.08783 -283.91364 18825.45159 -84449.68761
176 177 178 179 180
15108.76463 -33521.45373 8517.00026 -11975.80454 36540.21971
181 182 183 184 185
-26634.41095 -60739.68987 5736.37107 26076.02674 42029.49517
186 187 188 189 190
3031.21024 12935.80988 14756.88250 -1441.99545 14081.31360
191 192 193 194 195
27666.71469 11684.07754 -19327.53844 -20696.39693 -472.17236
196 197 198 199 200
1711.10049 -9659.48215 704.98200 -23841.08339 -8307.29664
201 202 203 204 205
5306.01056 9334.50984 14509.62400 347.77189 -21644.91521
206 207 208 209 210
-3546.99712 -15567.82447 50055.30204 26473.38992 7883.18184
211 212 213 214 215
9416.53186 12922.97848 -1066.52060 4805.98484 -24358.30973
216 217 218 219 220
-10695.15857 -29515.63779 7044.81888 10583.11348 27676.75548
221 222 223 224 225
-18468.72380 -32024.12292 -55373.60902 340.87940 8375.18731
226 227 228 229 230
5355.00092 17200.42427 -21838.24561 -26835.06367 -59478.25015
231 232 233 234 235
-10973.46033 -13746.04441 -36725.50802 -4361.65297 -17910.76376
236 237 238 239 240
-12694.87786 53513.99412 1433.55821 -4270.69738 -13745.65009
241 242 243 244 245
8421.71025 -19633.49856 -14033.98363 -3228.80810 1333.17937
246 247 248 249 250
-3696.73029 -34147.18670 24461.88553 9578.34396 -8861.25266
251 252 253 254 255
17155.42935 5710.10336 12385.22920 3518.42688 -6792.00279
256 257 258 259 260
8876.56821 -12091.52252 17109.94531 15971.82214 6750.99521
261 262 263 264 265
-1997.38968 1454.61676 9172.29191 -4530.12013 -12345.74724
266 267 268 269 270
12095.46304 -13771.19778 4738.22989 52727.10689 -9751.74114
271 272 273 274 275
-3960.24266 7610.07914 5720.07428 7204.32321 9754.30241
276 277 278 279 280
-9641.70626 11510.80817 10951.10329 -4098.58943 21416.05188
281 282 283 284 285
1445.36551 8129.50701 7031.44229 16389.03335 10773.53438
286 287 288 289
-22401.63529 1895.28429 -5398.14417 8137.64420
> postscript(file="/var/wessaorg/rcomp/tmp/643dw1324466480.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> dum <- cbind(lag(myerror,k=1),myerror)
> dum
Time Series:
Start = 0
End = 289
Frequency = 1
lag(myerror, k = 1) myerror
0 39311.89863 NA
1 -4972.98917 39311.89863
2 21759.92248 -4972.98917
3 -52650.25936 21759.92248
4 22214.18449 -52650.25936
5 -12884.05949 22214.18449
6 70852.89528 -12884.05949
7 -2788.76485 70852.89528
8 13.24984 -2788.76485
9 -19484.53103 13.24984
10 16685.95385 -19484.53103
11 36559.53426 16685.95385
12 -13674.68602 36559.53426
13 -6050.96424 -13674.68602
14 28122.54285 -6050.96424
15 34489.00676 28122.54285
16 -79472.81624 34489.00676
17 20145.41417 -79472.81624
18 18289.56003 20145.41417
19 -6256.57013 18289.56003
20 -21058.69191 -6256.57013
21 -4274.86990 -21058.69191
22 120426.53329 -4274.86990
23 24125.90655 120426.53329
24 -40744.09666 24125.90655
25 -83134.99463 -40744.09666
26 -54083.51388 -83134.99463
27 -5448.57143 -54083.51388
28 22651.63412 -5448.57143
29 -4150.84535 22651.63412
30 21649.41276 -4150.84535
31 486.78876 21649.41276
32 894.87108 486.78876
33 35376.51515 894.87108
34 -23907.22038 35376.51515
35 34001.26317 -23907.22038
36 18239.07291 34001.26317
37 24240.14611 18239.07291
38 57128.89000 24240.14611
39 -6589.53964 57128.89000
40 15298.32456 -6589.53964
41 21522.02447 15298.32456
42 -48798.04363 21522.02447
43 -2353.94603 -48798.04363
44 18841.98645 -2353.94603
45 8664.97771 18841.98645
46 -57514.33608 8664.97771
47 3585.19881 -57514.33608
48 -7427.40490 3585.19881
49 -58238.70703 -7427.40490
50 56723.59676 -58238.70703
51 -12347.54869 56723.59676
52 -30814.14673 -12347.54869
53 -44831.05397 -30814.14673
54 -36361.26423 -44831.05397
55 -3334.23601 -36361.26423
56 -14001.00292 -3334.23601
57 67027.86879 -14001.00292
58 -8387.45094 67027.86879
59 -837.26567 -8387.45094
60 -741.83980 -837.26567
61 -9277.54637 -741.83980
62 -52174.27831 -9277.54637
63 5068.19553 -52174.27831
64 -27142.39072 5068.19553
65 16152.25866 -27142.39072
66 33741.72169 16152.25866
67 -5526.75993 33741.72169
68 -2202.30423 -5526.75993
69 -1827.58717 -2202.30423
70 11877.50789 -1827.58717
71 50769.15196 11877.50789
72 25162.84293 50769.15196
73 -35611.31504 25162.84293
74 -34844.72236 -35611.31504
75 -16118.33614 -34844.72236
76 -35313.21680 -16118.33614
77 -11408.28481 -35313.21680
78 -7125.58562 -11408.28481
79 2313.89018 -7125.58562
80 17955.01964 2313.89018
81 -34923.67530 17955.01964
82 65832.39974 -34923.67530
83 38350.94122 65832.39974
84 -23937.45809 38350.94122
85 -34624.47622 -23937.45809
86 -4223.08798 -34624.47622
87 -18792.41659 -4223.08798
88 54186.68826 -18792.41659
89 15410.96423 54186.68826
90 11077.39931 15410.96423
91 -10585.54936 11077.39931
92 -5276.53501 -10585.54936
93 -6152.71367 -5276.53501
94 8858.39833 -6152.71367
95 26403.29203 8858.39833
96 39628.97342 26403.29203
97 -134525.99604 39628.97342
98 28380.69994 -134525.99604
99 -21333.35190 28380.69994
100 11155.85303 -21333.35190
101 -17397.07988 11155.85303
102 47649.56301 -17397.07988
103 -4298.99524 47649.56301
104 10903.14281 -4298.99524
105 13849.58930 10903.14281
106 -35899.70759 13849.58930
107 -21135.24029 -35899.70759
108 27075.53204 -21135.24029
109 2168.30495 27075.53204
110 9347.35468 2168.30495
111 6262.56626 9347.35468
112 -17921.48319 6262.56626
113 7451.54542 -17921.48319
114 -10052.69102 7451.54542
115 -6192.99932 -10052.69102
116 -8620.75945 -6192.99932
117 39220.86161 -8620.75945
118 -35334.19029 39220.86161
119 -34391.25173 -35334.19029
120 -11033.35493 -34391.25173
121 -38405.16575 -11033.35493
122 -28518.05125 -38405.16575
123 -71619.23916 -28518.05125
124 58721.09947 -71619.23916
125 12457.08959 58721.09947
126 5053.40988 12457.08959
127 -35560.79787 5053.40988
128 16744.00235 -35560.79787
129 -1311.85449 16744.00235
130 -7640.92877 -1311.85449
131 -27218.89093 -7640.92877
132 15347.24417 -27218.89093
133 -69714.11314 15347.24417
134 -11452.95326 -69714.11314
135 75689.71901 -11452.95326
136 -1658.10440 75689.71901
137 45664.29384 -1658.10440
138 -33566.25443 45664.29384
139 -3990.35253 -33566.25443
140 -19735.06652 -3990.35253
141 -28097.47854 -19735.06652
142 12964.00702 -28097.47854
143 -8126.90361 12964.00702
144 3514.30084 -8126.90361
145 -33995.14486 3514.30084
146 -2888.06724 -33995.14486
147 -20716.72490 -2888.06724
148 -5462.93781 -20716.72490
149 -4257.33307 -5462.93781
150 35784.21577 -4257.33307
151 24488.90835 35784.21577
152 -17546.08380 24488.90835
153 44968.81845 -17546.08380
154 24262.52367 44968.81845
155 20366.97621 24262.52367
156 27403.10064 20366.97621
157 22558.45456 27403.10064
158 -19132.46386 22558.45456
159 -3943.91461 -19132.46386
160 19181.04977 -3943.91461
161 -20268.57120 19181.04977
162 -16563.36847 -20268.57120
163 3586.21046 -16563.36847
164 -16717.60936 3586.21046
165 36368.12609 -16717.60936
166 5953.73709 36368.12609
167 24156.52882 5953.73709
168 17633.01975 24156.52882
169 72415.25487 17633.01975
170 13047.05478 72415.25487
171 2197.08783 13047.05478
172 -283.91364 2197.08783
173 18825.45159 -283.91364
174 -84449.68761 18825.45159
175 15108.76463 -84449.68761
176 -33521.45373 15108.76463
177 8517.00026 -33521.45373
178 -11975.80454 8517.00026
179 36540.21971 -11975.80454
180 -26634.41095 36540.21971
181 -60739.68987 -26634.41095
182 5736.37107 -60739.68987
183 26076.02674 5736.37107
184 42029.49517 26076.02674
185 3031.21024 42029.49517
186 12935.80988 3031.21024
187 14756.88250 12935.80988
188 -1441.99545 14756.88250
189 14081.31360 -1441.99545
190 27666.71469 14081.31360
191 11684.07754 27666.71469
192 -19327.53844 11684.07754
193 -20696.39693 -19327.53844
194 -472.17236 -20696.39693
195 1711.10049 -472.17236
196 -9659.48215 1711.10049
197 704.98200 -9659.48215
198 -23841.08339 704.98200
199 -8307.29664 -23841.08339
200 5306.01056 -8307.29664
201 9334.50984 5306.01056
202 14509.62400 9334.50984
203 347.77189 14509.62400
204 -21644.91521 347.77189
205 -3546.99712 -21644.91521
206 -15567.82447 -3546.99712
207 50055.30204 -15567.82447
208 26473.38992 50055.30204
209 7883.18184 26473.38992
210 9416.53186 7883.18184
211 12922.97848 9416.53186
212 -1066.52060 12922.97848
213 4805.98484 -1066.52060
214 -24358.30973 4805.98484
215 -10695.15857 -24358.30973
216 -29515.63779 -10695.15857
217 7044.81888 -29515.63779
218 10583.11348 7044.81888
219 27676.75548 10583.11348
220 -18468.72380 27676.75548
221 -32024.12292 -18468.72380
222 -55373.60902 -32024.12292
223 340.87940 -55373.60902
224 8375.18731 340.87940
225 5355.00092 8375.18731
226 17200.42427 5355.00092
227 -21838.24561 17200.42427
228 -26835.06367 -21838.24561
229 -59478.25015 -26835.06367
230 -10973.46033 -59478.25015
231 -13746.04441 -10973.46033
232 -36725.50802 -13746.04441
233 -4361.65297 -36725.50802
234 -17910.76376 -4361.65297
235 -12694.87786 -17910.76376
236 53513.99412 -12694.87786
237 1433.55821 53513.99412
238 -4270.69738 1433.55821
239 -13745.65009 -4270.69738
240 8421.71025 -13745.65009
241 -19633.49856 8421.71025
242 -14033.98363 -19633.49856
243 -3228.80810 -14033.98363
244 1333.17937 -3228.80810
245 -3696.73029 1333.17937
246 -34147.18670 -3696.73029
247 24461.88553 -34147.18670
248 9578.34396 24461.88553
249 -8861.25266 9578.34396
250 17155.42935 -8861.25266
251 5710.10336 17155.42935
252 12385.22920 5710.10336
253 3518.42688 12385.22920
254 -6792.00279 3518.42688
255 8876.56821 -6792.00279
256 -12091.52252 8876.56821
257 17109.94531 -12091.52252
258 15971.82214 17109.94531
259 6750.99521 15971.82214
260 -1997.38968 6750.99521
261 1454.61676 -1997.38968
262 9172.29191 1454.61676
263 -4530.12013 9172.29191
264 -12345.74724 -4530.12013
265 12095.46304 -12345.74724
266 -13771.19778 12095.46304
267 4738.22989 -13771.19778
268 52727.10689 4738.22989
269 -9751.74114 52727.10689
270 -3960.24266 -9751.74114
271 7610.07914 -3960.24266
272 5720.07428 7610.07914
273 7204.32321 5720.07428
274 9754.30241 7204.32321
275 -9641.70626 9754.30241
276 11510.80817 -9641.70626
277 10951.10329 11510.80817
278 -4098.58943 10951.10329
279 21416.05188 -4098.58943
280 1445.36551 21416.05188
281 8129.50701 1445.36551
282 7031.44229 8129.50701
283 16389.03335 7031.44229
284 10773.53438 16389.03335
285 -22401.63529 10773.53438
286 1895.28429 -22401.63529
287 -5398.14417 1895.28429
288 8137.64420 -5398.14417
289 NA 8137.64420
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] -4972.98917 39311.89863
[2,] 21759.92248 -4972.98917
[3,] -52650.25936 21759.92248
[4,] 22214.18449 -52650.25936
[5,] -12884.05949 22214.18449
[6,] 70852.89528 -12884.05949
[7,] -2788.76485 70852.89528
[8,] 13.24984 -2788.76485
[9,] -19484.53103 13.24984
[10,] 16685.95385 -19484.53103
[11,] 36559.53426 16685.95385
[12,] -13674.68602 36559.53426
[13,] -6050.96424 -13674.68602
[14,] 28122.54285 -6050.96424
[15,] 34489.00676 28122.54285
[16,] -79472.81624 34489.00676
[17,] 20145.41417 -79472.81624
[18,] 18289.56003 20145.41417
[19,] -6256.57013 18289.56003
[20,] -21058.69191 -6256.57013
[21,] -4274.86990 -21058.69191
[22,] 120426.53329 -4274.86990
[23,] 24125.90655 120426.53329
[24,] -40744.09666 24125.90655
[25,] -83134.99463 -40744.09666
[26,] -54083.51388 -83134.99463
[27,] -5448.57143 -54083.51388
[28,] 22651.63412 -5448.57143
[29,] -4150.84535 22651.63412
[30,] 21649.41276 -4150.84535
[31,] 486.78876 21649.41276
[32,] 894.87108 486.78876
[33,] 35376.51515 894.87108
[34,] -23907.22038 35376.51515
[35,] 34001.26317 -23907.22038
[36,] 18239.07291 34001.26317
[37,] 24240.14611 18239.07291
[38,] 57128.89000 24240.14611
[39,] -6589.53964 57128.89000
[40,] 15298.32456 -6589.53964
[41,] 21522.02447 15298.32456
[42,] -48798.04363 21522.02447
[43,] -2353.94603 -48798.04363
[44,] 18841.98645 -2353.94603
[45,] 8664.97771 18841.98645
[46,] -57514.33608 8664.97771
[47,] 3585.19881 -57514.33608
[48,] -7427.40490 3585.19881
[49,] -58238.70703 -7427.40490
[50,] 56723.59676 -58238.70703
[51,] -12347.54869 56723.59676
[52,] -30814.14673 -12347.54869
[53,] -44831.05397 -30814.14673
[54,] -36361.26423 -44831.05397
[55,] -3334.23601 -36361.26423
[56,] -14001.00292 -3334.23601
[57,] 67027.86879 -14001.00292
[58,] -8387.45094 67027.86879
[59,] -837.26567 -8387.45094
[60,] -741.83980 -837.26567
[61,] -9277.54637 -741.83980
[62,] -52174.27831 -9277.54637
[63,] 5068.19553 -52174.27831
[64,] -27142.39072 5068.19553
[65,] 16152.25866 -27142.39072
[66,] 33741.72169 16152.25866
[67,] -5526.75993 33741.72169
[68,] -2202.30423 -5526.75993
[69,] -1827.58717 -2202.30423
[70,] 11877.50789 -1827.58717
[71,] 50769.15196 11877.50789
[72,] 25162.84293 50769.15196
[73,] -35611.31504 25162.84293
[74,] -34844.72236 -35611.31504
[75,] -16118.33614 -34844.72236
[76,] -35313.21680 -16118.33614
[77,] -11408.28481 -35313.21680
[78,] -7125.58562 -11408.28481
[79,] 2313.89018 -7125.58562
[80,] 17955.01964 2313.89018
[81,] -34923.67530 17955.01964
[82,] 65832.39974 -34923.67530
[83,] 38350.94122 65832.39974
[84,] -23937.45809 38350.94122
[85,] -34624.47622 -23937.45809
[86,] -4223.08798 -34624.47622
[87,] -18792.41659 -4223.08798
[88,] 54186.68826 -18792.41659
[89,] 15410.96423 54186.68826
[90,] 11077.39931 15410.96423
[91,] -10585.54936 11077.39931
[92,] -5276.53501 -10585.54936
[93,] -6152.71367 -5276.53501
[94,] 8858.39833 -6152.71367
[95,] 26403.29203 8858.39833
[96,] 39628.97342 26403.29203
[97,] -134525.99604 39628.97342
[98,] 28380.69994 -134525.99604
[99,] -21333.35190 28380.69994
[100,] 11155.85303 -21333.35190
[101,] -17397.07988 11155.85303
[102,] 47649.56301 -17397.07988
[103,] -4298.99524 47649.56301
[104,] 10903.14281 -4298.99524
[105,] 13849.58930 10903.14281
[106,] -35899.70759 13849.58930
[107,] -21135.24029 -35899.70759
[108,] 27075.53204 -21135.24029
[109,] 2168.30495 27075.53204
[110,] 9347.35468 2168.30495
[111,] 6262.56626 9347.35468
[112,] -17921.48319 6262.56626
[113,] 7451.54542 -17921.48319
[114,] -10052.69102 7451.54542
[115,] -6192.99932 -10052.69102
[116,] -8620.75945 -6192.99932
[117,] 39220.86161 -8620.75945
[118,] -35334.19029 39220.86161
[119,] -34391.25173 -35334.19029
[120,] -11033.35493 -34391.25173
[121,] -38405.16575 -11033.35493
[122,] -28518.05125 -38405.16575
[123,] -71619.23916 -28518.05125
[124,] 58721.09947 -71619.23916
[125,] 12457.08959 58721.09947
[126,] 5053.40988 12457.08959
[127,] -35560.79787 5053.40988
[128,] 16744.00235 -35560.79787
[129,] -1311.85449 16744.00235
[130,] -7640.92877 -1311.85449
[131,] -27218.89093 -7640.92877
[132,] 15347.24417 -27218.89093
[133,] -69714.11314 15347.24417
[134,] -11452.95326 -69714.11314
[135,] 75689.71901 -11452.95326
[136,] -1658.10440 75689.71901
[137,] 45664.29384 -1658.10440
[138,] -33566.25443 45664.29384
[139,] -3990.35253 -33566.25443
[140,] -19735.06652 -3990.35253
[141,] -28097.47854 -19735.06652
[142,] 12964.00702 -28097.47854
[143,] -8126.90361 12964.00702
[144,] 3514.30084 -8126.90361
[145,] -33995.14486 3514.30084
[146,] -2888.06724 -33995.14486
[147,] -20716.72490 -2888.06724
[148,] -5462.93781 -20716.72490
[149,] -4257.33307 -5462.93781
[150,] 35784.21577 -4257.33307
[151,] 24488.90835 35784.21577
[152,] -17546.08380 24488.90835
[153,] 44968.81845 -17546.08380
[154,] 24262.52367 44968.81845
[155,] 20366.97621 24262.52367
[156,] 27403.10064 20366.97621
[157,] 22558.45456 27403.10064
[158,] -19132.46386 22558.45456
[159,] -3943.91461 -19132.46386
[160,] 19181.04977 -3943.91461
[161,] -20268.57120 19181.04977
[162,] -16563.36847 -20268.57120
[163,] 3586.21046 -16563.36847
[164,] -16717.60936 3586.21046
[165,] 36368.12609 -16717.60936
[166,] 5953.73709 36368.12609
[167,] 24156.52882 5953.73709
[168,] 17633.01975 24156.52882
[169,] 72415.25487 17633.01975
[170,] 13047.05478 72415.25487
[171,] 2197.08783 13047.05478
[172,] -283.91364 2197.08783
[173,] 18825.45159 -283.91364
[174,] -84449.68761 18825.45159
[175,] 15108.76463 -84449.68761
[176,] -33521.45373 15108.76463
[177,] 8517.00026 -33521.45373
[178,] -11975.80454 8517.00026
[179,] 36540.21971 -11975.80454
[180,] -26634.41095 36540.21971
[181,] -60739.68987 -26634.41095
[182,] 5736.37107 -60739.68987
[183,] 26076.02674 5736.37107
[184,] 42029.49517 26076.02674
[185,] 3031.21024 42029.49517
[186,] 12935.80988 3031.21024
[187,] 14756.88250 12935.80988
[188,] -1441.99545 14756.88250
[189,] 14081.31360 -1441.99545
[190,] 27666.71469 14081.31360
[191,] 11684.07754 27666.71469
[192,] -19327.53844 11684.07754
[193,] -20696.39693 -19327.53844
[194,] -472.17236 -20696.39693
[195,] 1711.10049 -472.17236
[196,] -9659.48215 1711.10049
[197,] 704.98200 -9659.48215
[198,] -23841.08339 704.98200
[199,] -8307.29664 -23841.08339
[200,] 5306.01056 -8307.29664
[201,] 9334.50984 5306.01056
[202,] 14509.62400 9334.50984
[203,] 347.77189 14509.62400
[204,] -21644.91521 347.77189
[205,] -3546.99712 -21644.91521
[206,] -15567.82447 -3546.99712
[207,] 50055.30204 -15567.82447
[208,] 26473.38992 50055.30204
[209,] 7883.18184 26473.38992
[210,] 9416.53186 7883.18184
[211,] 12922.97848 9416.53186
[212,] -1066.52060 12922.97848
[213,] 4805.98484 -1066.52060
[214,] -24358.30973 4805.98484
[215,] -10695.15857 -24358.30973
[216,] -29515.63779 -10695.15857
[217,] 7044.81888 -29515.63779
[218,] 10583.11348 7044.81888
[219,] 27676.75548 10583.11348
[220,] -18468.72380 27676.75548
[221,] -32024.12292 -18468.72380
[222,] -55373.60902 -32024.12292
[223,] 340.87940 -55373.60902
[224,] 8375.18731 340.87940
[225,] 5355.00092 8375.18731
[226,] 17200.42427 5355.00092
[227,] -21838.24561 17200.42427
[228,] -26835.06367 -21838.24561
[229,] -59478.25015 -26835.06367
[230,] -10973.46033 -59478.25015
[231,] -13746.04441 -10973.46033
[232,] -36725.50802 -13746.04441
[233,] -4361.65297 -36725.50802
[234,] -17910.76376 -4361.65297
[235,] -12694.87786 -17910.76376
[236,] 53513.99412 -12694.87786
[237,] 1433.55821 53513.99412
[238,] -4270.69738 1433.55821
[239,] -13745.65009 -4270.69738
[240,] 8421.71025 -13745.65009
[241,] -19633.49856 8421.71025
[242,] -14033.98363 -19633.49856
[243,] -3228.80810 -14033.98363
[244,] 1333.17937 -3228.80810
[245,] -3696.73029 1333.17937
[246,] -34147.18670 -3696.73029
[247,] 24461.88553 -34147.18670
[248,] 9578.34396 24461.88553
[249,] -8861.25266 9578.34396
[250,] 17155.42935 -8861.25266
[251,] 5710.10336 17155.42935
[252,] 12385.22920 5710.10336
[253,] 3518.42688 12385.22920
[254,] -6792.00279 3518.42688
[255,] 8876.56821 -6792.00279
[256,] -12091.52252 8876.56821
[257,] 17109.94531 -12091.52252
[258,] 15971.82214 17109.94531
[259,] 6750.99521 15971.82214
[260,] -1997.38968 6750.99521
[261,] 1454.61676 -1997.38968
[262,] 9172.29191 1454.61676
[263,] -4530.12013 9172.29191
[264,] -12345.74724 -4530.12013
[265,] 12095.46304 -12345.74724
[266,] -13771.19778 12095.46304
[267,] 4738.22989 -13771.19778
[268,] 52727.10689 4738.22989
[269,] -9751.74114 52727.10689
[270,] -3960.24266 -9751.74114
[271,] 7610.07914 -3960.24266
[272,] 5720.07428 7610.07914
[273,] 7204.32321 5720.07428
[274,] 9754.30241 7204.32321
[275,] -9641.70626 9754.30241
[276,] 11510.80817 -9641.70626
[277,] 10951.10329 11510.80817
[278,] -4098.58943 10951.10329
[279,] 21416.05188 -4098.58943
[280,] 1445.36551 21416.05188
[281,] 8129.50701 1445.36551
[282,] 7031.44229 8129.50701
[283,] 16389.03335 7031.44229
[284,] 10773.53438 16389.03335
[285,] -22401.63529 10773.53438
[286,] 1895.28429 -22401.63529
[287,] -5398.14417 1895.28429
[288,] 8137.64420 -5398.14417
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 -4972.98917 39311.89863
2 21759.92248 -4972.98917
3 -52650.25936 21759.92248
4 22214.18449 -52650.25936
5 -12884.05949 22214.18449
6 70852.89528 -12884.05949
7 -2788.76485 70852.89528
8 13.24984 -2788.76485
9 -19484.53103 13.24984
10 16685.95385 -19484.53103
11 36559.53426 16685.95385
12 -13674.68602 36559.53426
13 -6050.96424 -13674.68602
14 28122.54285 -6050.96424
15 34489.00676 28122.54285
16 -79472.81624 34489.00676
17 20145.41417 -79472.81624
18 18289.56003 20145.41417
19 -6256.57013 18289.56003
20 -21058.69191 -6256.57013
21 -4274.86990 -21058.69191
22 120426.53329 -4274.86990
23 24125.90655 120426.53329
24 -40744.09666 24125.90655
25 -83134.99463 -40744.09666
26 -54083.51388 -83134.99463
27 -5448.57143 -54083.51388
28 22651.63412 -5448.57143
29 -4150.84535 22651.63412
30 21649.41276 -4150.84535
31 486.78876 21649.41276
32 894.87108 486.78876
33 35376.51515 894.87108
34 -23907.22038 35376.51515
35 34001.26317 -23907.22038
36 18239.07291 34001.26317
37 24240.14611 18239.07291
38 57128.89000 24240.14611
39 -6589.53964 57128.89000
40 15298.32456 -6589.53964
41 21522.02447 15298.32456
42 -48798.04363 21522.02447
43 -2353.94603 -48798.04363
44 18841.98645 -2353.94603
45 8664.97771 18841.98645
46 -57514.33608 8664.97771
47 3585.19881 -57514.33608
48 -7427.40490 3585.19881
49 -58238.70703 -7427.40490
50 56723.59676 -58238.70703
51 -12347.54869 56723.59676
52 -30814.14673 -12347.54869
53 -44831.05397 -30814.14673
54 -36361.26423 -44831.05397
55 -3334.23601 -36361.26423
56 -14001.00292 -3334.23601
57 67027.86879 -14001.00292
58 -8387.45094 67027.86879
59 -837.26567 -8387.45094
60 -741.83980 -837.26567
61 -9277.54637 -741.83980
62 -52174.27831 -9277.54637
63 5068.19553 -52174.27831
64 -27142.39072 5068.19553
65 16152.25866 -27142.39072
66 33741.72169 16152.25866
67 -5526.75993 33741.72169
68 -2202.30423 -5526.75993
69 -1827.58717 -2202.30423
70 11877.50789 -1827.58717
71 50769.15196 11877.50789
72 25162.84293 50769.15196
73 -35611.31504 25162.84293
74 -34844.72236 -35611.31504
75 -16118.33614 -34844.72236
76 -35313.21680 -16118.33614
77 -11408.28481 -35313.21680
78 -7125.58562 -11408.28481
79 2313.89018 -7125.58562
80 17955.01964 2313.89018
81 -34923.67530 17955.01964
82 65832.39974 -34923.67530
83 38350.94122 65832.39974
84 -23937.45809 38350.94122
85 -34624.47622 -23937.45809
86 -4223.08798 -34624.47622
87 -18792.41659 -4223.08798
88 54186.68826 -18792.41659
89 15410.96423 54186.68826
90 11077.39931 15410.96423
91 -10585.54936 11077.39931
92 -5276.53501 -10585.54936
93 -6152.71367 -5276.53501
94 8858.39833 -6152.71367
95 26403.29203 8858.39833
96 39628.97342 26403.29203
97 -134525.99604 39628.97342
98 28380.69994 -134525.99604
99 -21333.35190 28380.69994
100 11155.85303 -21333.35190
101 -17397.07988 11155.85303
102 47649.56301 -17397.07988
103 -4298.99524 47649.56301
104 10903.14281 -4298.99524
105 13849.58930 10903.14281
106 -35899.70759 13849.58930
107 -21135.24029 -35899.70759
108 27075.53204 -21135.24029
109 2168.30495 27075.53204
110 9347.35468 2168.30495
111 6262.56626 9347.35468
112 -17921.48319 6262.56626
113 7451.54542 -17921.48319
114 -10052.69102 7451.54542
115 -6192.99932 -10052.69102
116 -8620.75945 -6192.99932
117 39220.86161 -8620.75945
118 -35334.19029 39220.86161
119 -34391.25173 -35334.19029
120 -11033.35493 -34391.25173
121 -38405.16575 -11033.35493
122 -28518.05125 -38405.16575
123 -71619.23916 -28518.05125
124 58721.09947 -71619.23916
125 12457.08959 58721.09947
126 5053.40988 12457.08959
127 -35560.79787 5053.40988
128 16744.00235 -35560.79787
129 -1311.85449 16744.00235
130 -7640.92877 -1311.85449
131 -27218.89093 -7640.92877
132 15347.24417 -27218.89093
133 -69714.11314 15347.24417
134 -11452.95326 -69714.11314
135 75689.71901 -11452.95326
136 -1658.10440 75689.71901
137 45664.29384 -1658.10440
138 -33566.25443 45664.29384
139 -3990.35253 -33566.25443
140 -19735.06652 -3990.35253
141 -28097.47854 -19735.06652
142 12964.00702 -28097.47854
143 -8126.90361 12964.00702
144 3514.30084 -8126.90361
145 -33995.14486 3514.30084
146 -2888.06724 -33995.14486
147 -20716.72490 -2888.06724
148 -5462.93781 -20716.72490
149 -4257.33307 -5462.93781
150 35784.21577 -4257.33307
151 24488.90835 35784.21577
152 -17546.08380 24488.90835
153 44968.81845 -17546.08380
154 24262.52367 44968.81845
155 20366.97621 24262.52367
156 27403.10064 20366.97621
157 22558.45456 27403.10064
158 -19132.46386 22558.45456
159 -3943.91461 -19132.46386
160 19181.04977 -3943.91461
161 -20268.57120 19181.04977
162 -16563.36847 -20268.57120
163 3586.21046 -16563.36847
164 -16717.60936 3586.21046
165 36368.12609 -16717.60936
166 5953.73709 36368.12609
167 24156.52882 5953.73709
168 17633.01975 24156.52882
169 72415.25487 17633.01975
170 13047.05478 72415.25487
171 2197.08783 13047.05478
172 -283.91364 2197.08783
173 18825.45159 -283.91364
174 -84449.68761 18825.45159
175 15108.76463 -84449.68761
176 -33521.45373 15108.76463
177 8517.00026 -33521.45373
178 -11975.80454 8517.00026
179 36540.21971 -11975.80454
180 -26634.41095 36540.21971
181 -60739.68987 -26634.41095
182 5736.37107 -60739.68987
183 26076.02674 5736.37107
184 42029.49517 26076.02674
185 3031.21024 42029.49517
186 12935.80988 3031.21024
187 14756.88250 12935.80988
188 -1441.99545 14756.88250
189 14081.31360 -1441.99545
190 27666.71469 14081.31360
191 11684.07754 27666.71469
192 -19327.53844 11684.07754
193 -20696.39693 -19327.53844
194 -472.17236 -20696.39693
195 1711.10049 -472.17236
196 -9659.48215 1711.10049
197 704.98200 -9659.48215
198 -23841.08339 704.98200
199 -8307.29664 -23841.08339
200 5306.01056 -8307.29664
201 9334.50984 5306.01056
202 14509.62400 9334.50984
203 347.77189 14509.62400
204 -21644.91521 347.77189
205 -3546.99712 -21644.91521
206 -15567.82447 -3546.99712
207 50055.30204 -15567.82447
208 26473.38992 50055.30204
209 7883.18184 26473.38992
210 9416.53186 7883.18184
211 12922.97848 9416.53186
212 -1066.52060 12922.97848
213 4805.98484 -1066.52060
214 -24358.30973 4805.98484
215 -10695.15857 -24358.30973
216 -29515.63779 -10695.15857
217 7044.81888 -29515.63779
218 10583.11348 7044.81888
219 27676.75548 10583.11348
220 -18468.72380 27676.75548
221 -32024.12292 -18468.72380
222 -55373.60902 -32024.12292
223 340.87940 -55373.60902
224 8375.18731 340.87940
225 5355.00092 8375.18731
226 17200.42427 5355.00092
227 -21838.24561 17200.42427
228 -26835.06367 -21838.24561
229 -59478.25015 -26835.06367
230 -10973.46033 -59478.25015
231 -13746.04441 -10973.46033
232 -36725.50802 -13746.04441
233 -4361.65297 -36725.50802
234 -17910.76376 -4361.65297
235 -12694.87786 -17910.76376
236 53513.99412 -12694.87786
237 1433.55821 53513.99412
238 -4270.69738 1433.55821
239 -13745.65009 -4270.69738
240 8421.71025 -13745.65009
241 -19633.49856 8421.71025
242 -14033.98363 -19633.49856
243 -3228.80810 -14033.98363
244 1333.17937 -3228.80810
245 -3696.73029 1333.17937
246 -34147.18670 -3696.73029
247 24461.88553 -34147.18670
248 9578.34396 24461.88553
249 -8861.25266 9578.34396
250 17155.42935 -8861.25266
251 5710.10336 17155.42935
252 12385.22920 5710.10336
253 3518.42688 12385.22920
254 -6792.00279 3518.42688
255 8876.56821 -6792.00279
256 -12091.52252 8876.56821
257 17109.94531 -12091.52252
258 15971.82214 17109.94531
259 6750.99521 15971.82214
260 -1997.38968 6750.99521
261 1454.61676 -1997.38968
262 9172.29191 1454.61676
263 -4530.12013 9172.29191
264 -12345.74724 -4530.12013
265 12095.46304 -12345.74724
266 -13771.19778 12095.46304
267 4738.22989 -13771.19778
268 52727.10689 4738.22989
269 -9751.74114 52727.10689
270 -3960.24266 -9751.74114
271 7610.07914 -3960.24266
272 5720.07428 7610.07914
273 7204.32321 5720.07428
274 9754.30241 7204.32321
275 -9641.70626 9754.30241
276 11510.80817 -9641.70626
277 10951.10329 11510.80817
278 -4098.58943 10951.10329
279 21416.05188 -4098.58943
280 1445.36551 21416.05188
281 8129.50701 1445.36551
282 7031.44229 8129.50701
283 16389.03335 7031.44229
284 10773.53438 16389.03335
285 -22401.63529 10773.53438
286 1895.28429 -22401.63529
287 -5398.14417 1895.28429
288 8137.64420 -5398.14417
> 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/7l93k1324466480.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/880xl1324466480.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/9e0ak1324466480.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/1048zv1324466480.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/11d9kl1324466480.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/12iqx11324466480.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/13kwmp1324466480.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/14zzml1324466480.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/159p6h1324466480.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/16q4241324466480.tab")
+ }
>
> try(system("convert tmp/1de2p1324466480.ps tmp/1de2p1324466480.png",intern=TRUE))
character(0)
> try(system("convert tmp/264xe1324466480.ps tmp/264xe1324466480.png",intern=TRUE))
character(0)
> try(system("convert tmp/3q5ap1324466480.ps tmp/3q5ap1324466480.png",intern=TRUE))
character(0)
> try(system("convert tmp/4795c1324466480.ps tmp/4795c1324466480.png",intern=TRUE))
character(0)
> try(system("convert tmp/5juzk1324466480.ps tmp/5juzk1324466480.png",intern=TRUE))
character(0)
> try(system("convert tmp/643dw1324466480.ps tmp/643dw1324466480.png",intern=TRUE))
character(0)
> try(system("convert tmp/7l93k1324466480.ps tmp/7l93k1324466480.png",intern=TRUE))
character(0)
> try(system("convert tmp/880xl1324466480.ps tmp/880xl1324466480.png",intern=TRUE))
character(0)
> try(system("convert tmp/9e0ak1324466480.ps tmp/9e0ak1324466480.png",intern=TRUE))
character(0)
> try(system("convert tmp/1048zv1324466480.ps tmp/1048zv1324466480.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
10.162 0.704 10.878