R version 2.15.2 (2012-10-26) -- "Trick or Treat"
Copyright (C) 2012 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i686-pc-linux-gnu (32-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> x <- array(list(1418
+ ,210907
+ ,56
+ ,396
+ ,81
+ ,3
+ ,79
+ ,30
+ ,115
+ ,94
+ ,112285
+ ,24188
+ ,869
+ ,120982
+ ,56
+ ,297
+ ,55
+ ,4
+ ,58
+ ,28
+ ,109
+ ,103
+ ,84786
+ ,18273
+ ,1530
+ ,176508
+ ,54
+ ,559
+ ,50
+ ,12
+ ,60
+ ,38
+ ,146
+ ,93
+ ,83123
+ ,14130
+ ,2172
+ ,179321
+ ,89
+ ,967
+ ,125
+ ,2
+ ,108
+ ,30
+ ,116
+ ,103
+ ,101193
+ ,32287
+ ,901
+ ,123185
+ ,40
+ ,270
+ ,40
+ ,1
+ ,49
+ ,22
+ ,68
+ ,51
+ ,38361
+ ,8654
+ ,463
+ ,52746
+ ,25
+ ,143
+ ,37
+ ,3
+ ,0
+ ,26
+ ,101
+ ,70
+ ,68504
+ ,9245
+ ,3201
+ ,385534
+ ,92
+ ,1562
+ ,63
+ ,0
+ ,121
+ ,25
+ ,96
+ ,91
+ ,119182
+ ,33251
+ ,371
+ ,33170
+ ,18
+ ,109
+ ,44
+ ,0
+ ,1
+ ,18
+ ,67
+ ,22
+ ,22807
+ ,1271
+ ,1192
+ ,101645
+ ,63
+ ,371
+ ,88
+ ,0
+ ,20
+ ,11
+ ,44
+ ,38
+ ,17140
+ ,5279
+ ,1583
+ ,149061
+ ,44
+ ,656
+ ,66
+ ,5
+ ,43
+ ,26
+ ,100
+ ,93
+ ,116174
+ ,27101
+ ,1439
+ ,165446
+ ,33
+ ,511
+ ,57
+ ,0
+ ,69
+ ,25
+ ,93
+ ,60
+ ,57635
+ ,16373
+ ,1764
+ ,237213
+ ,84
+ ,655
+ ,74
+ ,0
+ ,78
+ ,38
+ ,140
+ ,123
+ ,66198
+ ,19716
+ ,1495
+ ,173326
+ ,88
+ ,465
+ ,49
+ ,7
+ ,86
+ ,44
+ ,166
+ ,148
+ ,71701
+ ,17753
+ ,1373
+ ,133131
+ ,55
+ ,525
+ ,52
+ ,7
+ ,44
+ ,30
+ ,99
+ ,90
+ ,57793
+ ,9028
+ ,2187
+ ,258873
+ ,60
+ ,885
+ ,88
+ ,3
+ ,104
+ ,40
+ ,139
+ ,124
+ ,80444
+ ,18653
+ ,1491
+ ,180083
+ ,66
+ ,497
+ ,36
+ ,9
+ ,63
+ ,34
+ ,130
+ ,70
+ ,53855
+ ,8828
+ ,4041
+ ,324799
+ ,154
+ ,1436
+ ,108
+ ,0
+ ,158
+ ,47
+ ,181
+ ,168
+ ,97668
+ ,29498
+ ,1706
+ ,230964
+ ,53
+ ,612
+ ,43
+ ,4
+ ,102
+ ,30
+ ,116
+ ,115
+ ,133824
+ ,27563
+ ,2152
+ ,236785
+ ,119
+ ,865
+ ,75
+ ,3
+ ,77
+ ,31
+ ,116
+ ,71
+ ,101481
+ ,18293
+ ,1036
+ ,135473
+ ,41
+ ,385
+ ,32
+ ,0
+ ,82
+ ,23
+ ,88
+ ,66
+ ,99645
+ ,22530
+ ,1882
+ ,202925
+ ,61
+ ,567
+ ,44
+ ,7
+ ,115
+ ,36
+ ,139
+ ,134
+ ,114789
+ ,15977
+ ,1929
+ ,215147
+ ,58
+ ,639
+ ,85
+ ,0
+ ,101
+ ,36
+ ,135
+ ,117
+ ,99052
+ ,35082
+ ,2242
+ ,344297
+ ,75
+ ,963
+ ,86
+ ,1
+ ,80
+ ,30
+ ,108
+ ,108
+ ,67654
+ ,16116
+ ,1220
+ ,153935
+ ,33
+ ,398
+ ,56
+ ,5
+ ,50
+ ,25
+ ,89
+ ,84
+ ,65553
+ ,15849
+ ,1289
+ ,132943
+ ,40
+ ,410
+ ,50
+ ,7
+ ,83
+ ,39
+ ,156
+ ,156
+ ,97500
+ ,16026
+ ,2515
+ ,174724
+ ,92
+ ,966
+ ,135
+ ,0
+ ,123
+ ,34
+ ,129
+ ,120
+ ,69112
+ ,26569
+ ,2147
+ ,174415
+ ,100
+ ,801
+ ,63
+ ,0
+ ,73
+ ,31
+ ,118
+ ,114
+ ,82753
+ ,24785
+ ,2352
+ ,225548
+ ,112
+ ,892
+ ,81
+ ,5
+ ,81
+ ,31
+ ,118
+ ,94
+ ,85323
+ ,17569
+ ,1638
+ ,223632
+ ,73
+ ,513
+ ,52
+ ,0
+ ,105
+ ,33
+ ,125
+ ,120
+ ,72654
+ ,23825
+ ,1222
+ ,124817
+ ,40
+ ,469
+ ,44
+ ,0
+ ,47
+ ,25
+ ,95
+ ,81
+ ,30727
+ ,7869
+ ,1812
+ ,221698
+ ,45
+ ,683
+ ,113
+ ,0
+ ,105
+ ,33
+ ,126
+ ,110
+ ,77873
+ ,14975
+ ,1677
+ ,210767
+ ,60
+ ,643
+ ,39
+ ,3
+ ,94
+ ,35
+ ,135
+ ,133
+ ,117478
+ ,37791
+ ,1579
+ ,170266
+ ,62
+ ,535
+ ,73
+ ,4
+ ,44
+ ,42
+ ,154
+ ,122
+ ,74007
+ ,9605
+ ,1731
+ ,260561
+ ,75
+ ,625
+ ,48
+ ,1
+ ,114
+ ,43
+ ,165
+ ,158
+ ,90183
+ ,27295
+ ,807
+ ,84853
+ ,31
+ ,264
+ ,33
+ ,4
+ ,38
+ ,30
+ ,113
+ ,109
+ ,61542
+ ,2746
+ ,2452
+ ,294424
+ ,77
+ ,992
+ ,59
+ ,2
+ ,107
+ ,33
+ ,127
+ ,124
+ ,101494
+ ,34461
+ ,829
+ ,101011
+ ,34
+ ,238
+ ,41
+ ,0
+ ,30
+ ,13
+ ,52
+ ,39
+ ,27570
+ ,8098
+ ,1940
+ ,215641
+ ,46
+ ,818
+ ,69
+ ,0
+ ,71
+ ,32
+ ,121
+ ,92
+ ,55813
+ ,4787
+ ,2662
+ ,325107
+ ,99
+ ,937
+ ,64
+ ,0
+ ,84
+ ,36
+ ,136
+ ,126
+ ,79215
+ ,24919
+ ,186
+ ,7176
+ ,17
+ ,70
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,1423
+ ,603
+ ,1499
+ ,167542
+ ,66
+ ,507
+ ,59
+ ,2
+ ,59
+ ,28
+ ,108
+ ,70
+ ,55461
+ ,16329
+ ,865
+ ,106408
+ ,30
+ ,260
+ ,32
+ ,1
+ ,33
+ ,14
+ ,46
+ ,37
+ ,31081
+ ,12558
+ ,1793
+ ,96560
+ ,76
+ ,503
+ ,129
+ ,0
+ ,42
+ ,17
+ ,54
+ ,38
+ ,22996
+ ,7784
+ ,2527
+ ,265769
+ ,146
+ ,927
+ ,37
+ ,2
+ ,96
+ ,32
+ ,124
+ ,120
+ ,83122
+ ,28522
+ ,2747
+ ,269651
+ ,67
+ ,1269
+ ,31
+ ,10
+ ,106
+ ,30
+ ,115
+ ,93
+ ,70106
+ ,22265
+ ,1324
+ ,149112
+ ,56
+ ,537
+ ,65
+ ,6
+ ,56
+ ,35
+ ,128
+ ,95
+ ,60578
+ ,14459
+ ,2702
+ ,175824
+ ,107
+ ,910
+ ,107
+ ,0
+ ,57
+ ,20
+ ,80
+ ,77
+ ,39992
+ ,14526
+ ,1383
+ ,152871
+ ,58
+ ,532
+ ,74
+ ,5
+ ,59
+ ,28
+ ,97
+ ,90
+ ,79892
+ ,22240
+ ,1179
+ ,111665
+ ,34
+ ,345
+ ,54
+ ,4
+ ,39
+ ,28
+ ,104
+ ,80
+ ,49810
+ ,11802
+ ,2099
+ ,116408
+ ,61
+ ,918
+ ,76
+ ,1
+ ,34
+ ,39
+ ,59
+ ,31
+ ,71570
+ ,7623
+ ,4308
+ ,362301
+ ,119
+ ,1635
+ ,715
+ ,2
+ ,76
+ ,34
+ ,125
+ ,110
+ ,100708
+ ,11912
+ ,918
+ ,78800
+ ,42
+ ,330
+ ,57
+ ,2
+ ,20
+ ,26
+ ,82
+ ,66
+ ,33032
+ ,7935
+ ,1831
+ ,183167
+ ,66
+ ,557
+ ,66
+ ,0
+ ,91
+ ,39
+ ,149
+ ,138
+ ,82875
+ ,18220
+ ,3373
+ ,277965
+ ,89
+ ,1178
+ ,106
+ ,8
+ ,115
+ ,39
+ ,149
+ ,133
+ ,139077
+ ,19199
+ ,1713
+ ,150629
+ ,44
+ ,740
+ ,54
+ ,3
+ ,85
+ ,33
+ ,122
+ ,113
+ ,71595
+ ,19918
+ ,1438
+ ,168809
+ ,66
+ ,452
+ ,32
+ ,0
+ ,76
+ ,28
+ ,118
+ ,100
+ ,72260
+ ,21884
+ ,496
+ ,24188
+ ,24
+ ,218
+ ,20
+ ,0
+ ,8
+ ,4
+ ,12
+ ,7
+ ,5950
+ ,2694
+ ,2253
+ ,329267
+ ,259
+ ,764
+ ,71
+ ,8
+ ,79
+ ,39
+ ,144
+ ,140
+ ,115762
+ ,15808
+ ,744
+ ,65029
+ ,17
+ ,255
+ ,21
+ ,5
+ ,21
+ ,18
+ ,67
+ ,61
+ ,32551
+ ,3597
+ ,1161
+ ,101097
+ ,64
+ ,454
+ ,70
+ ,3
+ ,30
+ ,14
+ ,52
+ ,41
+ ,31701
+ ,5296
+ ,2352
+ ,218946
+ ,41
+ ,866
+ ,112
+ ,1
+ ,76
+ ,29
+ ,108
+ ,96
+ ,80670
+ ,25239
+ ,2144
+ ,244052
+ ,68
+ ,574
+ ,66
+ ,5
+ ,101
+ ,44
+ ,166
+ ,164
+ ,143558
+ ,29801
+ ,4691
+ ,341570
+ ,168
+ ,1276
+ ,190
+ ,1
+ ,94
+ ,21
+ ,80
+ ,78
+ ,117105
+ ,18450
+ ,1112
+ ,103597
+ ,43
+ ,379
+ ,66
+ ,1
+ ,27
+ ,16
+ ,60
+ ,49
+ ,23789
+ ,7132
+ ,2694
+ ,233328
+ ,132
+ ,825
+ ,165
+ ,5
+ ,92
+ ,28
+ ,107
+ ,102
+ ,120733
+ ,34861
+ ,1973
+ ,256462
+ ,105
+ ,798
+ ,56
+ ,0
+ ,123
+ ,35
+ ,127
+ ,124
+ ,105195
+ ,35940
+ ,1769
+ ,206161
+ ,71
+ ,663
+ ,61
+ ,12
+ ,75
+ ,28
+ ,107
+ ,99
+ ,73107
+ ,16688
+ ,3148
+ ,311473
+ ,112
+ ,1069
+ ,53
+ ,8
+ ,128
+ ,38
+ ,146
+ ,129
+ ,132068
+ ,24683
+ ,2474
+ ,235800
+ ,94
+ ,921
+ ,127
+ ,8
+ ,105
+ ,23
+ ,84
+ ,62
+ ,149193
+ ,46230
+ ,2084
+ ,177939
+ ,82
+ ,858
+ ,63
+ ,8
+ ,55
+ ,36
+ ,141
+ ,73
+ ,46821
+ ,10387
+ ,1954
+ ,207176
+ ,70
+ ,711
+ ,38
+ ,8
+ ,56
+ ,32
+ ,123
+ ,114
+ ,87011
+ ,21436
+ ,1226
+ ,196553
+ ,57
+ ,503
+ ,50
+ ,2
+ ,41
+ ,29
+ ,111
+ ,99
+ ,95260
+ ,30546
+ ,1389
+ ,174184
+ ,53
+ ,382
+ ,52
+ ,0
+ ,72
+ ,25
+ ,98
+ ,70
+ ,55183
+ ,19746
+ ,1496
+ ,143246
+ ,103
+ ,464
+ ,42
+ ,5
+ ,67
+ ,27
+ ,105
+ ,104
+ ,106671
+ ,15977
+ ,2269
+ ,187559
+ ,121
+ ,717
+ ,76
+ ,8
+ ,75
+ ,36
+ ,135
+ ,116
+ ,73511
+ ,22583
+ ,1833
+ ,187681
+ ,62
+ ,690
+ ,67
+ ,2
+ ,114
+ ,28
+ ,107
+ ,91
+ ,92945
+ ,17274
+ ,1268
+ ,119016
+ ,52
+ ,462
+ ,50
+ ,5
+ ,118
+ ,23
+ ,85
+ ,74
+ ,78664
+ ,16469
+ ,1943
+ ,182192
+ ,52
+ ,657
+ ,53
+ ,12
+ ,77
+ ,40
+ ,155
+ ,138
+ ,70054
+ ,14251
+ ,893
+ ,73566
+ ,32
+ ,385
+ ,39
+ ,6
+ ,22
+ ,23
+ ,88
+ ,67
+ ,22618
+ ,3007
+ ,1762
+ ,194979
+ ,62
+ ,577
+ ,50
+ ,7
+ ,66
+ ,40
+ ,155
+ ,151
+ ,74011
+ ,16851
+ ,1403
+ ,167488
+ ,45
+ ,619
+ ,77
+ ,2
+ ,69
+ ,28
+ ,104
+ ,72
+ ,83737
+ ,21113
+ ,1425
+ ,143756
+ ,46
+ ,479
+ ,57
+ ,0
+ ,105
+ ,34
+ ,132
+ ,120
+ ,69094
+ ,17401
+ ,1857
+ ,275541
+ ,63
+ ,817
+ ,73
+ ,4
+ ,116
+ ,33
+ ,127
+ ,115
+ ,93133
+ ,23958
+ ,1840
+ ,243199
+ ,75
+ ,752
+ ,34
+ ,3
+ ,88
+ ,28
+ ,108
+ ,105
+ ,95536
+ ,23567
+ ,1502
+ ,182999
+ ,88
+ ,430
+ ,39
+ ,6
+ ,73
+ ,34
+ ,129
+ ,104
+ ,225920
+ ,13065
+ ,1441
+ ,135649
+ ,46
+ ,451
+ ,46
+ ,2
+ ,99
+ ,30
+ ,116
+ ,108
+ ,62133
+ ,15358
+ ,1420
+ ,152299
+ ,53
+ ,537
+ ,63
+ ,0
+ ,62
+ ,33
+ ,122
+ ,98
+ ,61370
+ ,14587
+ ,1416
+ ,120221
+ ,37
+ ,519
+ ,35
+ ,1
+ ,53
+ ,22
+ ,85
+ ,69
+ ,43836
+ ,12770
+ ,2970
+ ,346485
+ ,90
+ ,1000
+ ,106
+ ,0
+ ,118
+ ,38
+ ,147
+ ,111
+ ,106117
+ ,24021
+ ,1317
+ ,145790
+ ,63
+ ,637
+ ,43
+ ,5
+ ,30
+ ,26
+ ,99
+ ,99
+ ,38692
+ ,9648
+ ,1644
+ ,193339
+ ,78
+ ,465
+ ,47
+ ,2
+ ,100
+ ,35
+ ,87
+ ,71
+ ,84651
+ ,20537
+ ,870
+ ,80953
+ ,25
+ ,437
+ ,31
+ ,0
+ ,49
+ ,8
+ ,28
+ ,27
+ ,56622
+ ,7905
+ ,1654
+ ,122774
+ ,45
+ ,711
+ ,162
+ ,0
+ ,24
+ ,24
+ ,90
+ ,69
+ ,15986
+ ,4527
+ ,1054
+ ,130585
+ ,46
+ ,299
+ ,57
+ ,5
+ ,67
+ ,29
+ ,109
+ ,107
+ ,95364
+ ,30495
+ ,937
+ ,112611
+ ,41
+ ,248
+ ,36
+ ,0
+ ,46
+ ,20
+ ,78
+ ,73
+ ,26706
+ ,7117
+ ,3004
+ ,286468
+ ,144
+ ,1162
+ ,263
+ ,1
+ ,57
+ ,29
+ ,111
+ ,107
+ ,89691
+ ,17719
+ ,2008
+ ,241066
+ ,82
+ ,714
+ ,78
+ ,0
+ ,75
+ ,45
+ ,158
+ ,93
+ ,67267
+ ,27056
+ ,2547
+ ,148446
+ ,91
+ ,905
+ ,63
+ ,1
+ ,135
+ ,37
+ ,141
+ ,129
+ ,126846
+ ,33473
+ ,1885
+ ,204713
+ ,71
+ ,649
+ ,54
+ ,1
+ ,68
+ ,33
+ ,122
+ ,69
+ ,41140
+ ,9758
+ ,1626
+ ,182079
+ ,63
+ ,512
+ ,63
+ ,2
+ ,124
+ ,33
+ ,124
+ ,118
+ ,102860
+ ,21115
+ ,1468
+ ,140344
+ ,53
+ ,472
+ ,77
+ ,6
+ ,33
+ ,25
+ ,93
+ ,73
+ ,51715
+ ,7236
+ ,2445
+ ,220516
+ ,62
+ ,905
+ ,79
+ ,1
+ ,98
+ ,32
+ ,124
+ ,119
+ ,55801
+ ,13790
+ ,1964
+ ,243060
+ ,63
+ ,786
+ ,110
+ ,4
+ ,58
+ ,29
+ ,112
+ ,104
+ ,111813
+ ,32902
+ ,1381
+ ,162765
+ ,32
+ ,489
+ ,56
+ ,2
+ ,68
+ ,28
+ ,108
+ ,107
+ ,120293
+ ,25131
+ ,1369
+ ,182613
+ ,39
+ ,479
+ ,56
+ ,3
+ ,81
+ ,28
+ ,99
+ ,99
+ ,138599
+ ,30910
+ ,1659
+ ,232138
+ ,62
+ ,617
+ ,43
+ ,0
+ ,131
+ ,31
+ ,117
+ ,90
+ ,161647
+ ,35947
+ ,2888
+ ,265318
+ ,117
+ ,925
+ ,111
+ ,10
+ ,110
+ ,52
+ ,199
+ ,197
+ ,115929
+ ,29848
+ ,1290
+ ,85574
+ ,34
+ ,351
+ ,71
+ ,0
+ ,37
+ ,21
+ ,78
+ ,36
+ ,24266
+ ,6943
+ ,2845
+ ,310839
+ ,92
+ ,1144
+ ,62
+ ,9
+ ,130
+ ,24
+ ,91
+ ,85
+ ,162901
+ ,42705
+ ,1982
+ ,225060
+ ,93
+ ,669
+ ,56
+ ,7
+ ,93
+ ,41
+ ,158
+ ,139
+ ,109825
+ ,31808
+ ,1904
+ ,232317
+ ,54
+ ,707
+ ,74
+ ,0
+ ,118
+ ,33
+ ,126
+ ,106
+ ,129838
+ ,26675
+ ,1391
+ ,144966
+ ,144
+ ,458
+ ,60
+ ,0
+ ,39
+ ,32
+ ,122
+ ,50
+ ,37510
+ ,8435
+ ,602
+ ,43287
+ ,14
+ ,214
+ ,43
+ ,4
+ ,13
+ ,19
+ ,71
+ ,64
+ ,43750
+ ,7409
+ ,1743
+ ,155754
+ ,61
+ ,599
+ ,68
+ ,4
+ ,74
+ ,20
+ ,75
+ ,31
+ ,40652
+ ,14993
+ ,1559
+ ,164709
+ ,109
+ ,572
+ ,53
+ ,0
+ ,81
+ ,31
+ ,115
+ ,63
+ ,87771
+ ,36867
+ ,2014
+ ,201940
+ ,38
+ ,897
+ ,87
+ ,0
+ ,109
+ ,31
+ ,119
+ ,92
+ ,85872
+ ,33835
+ ,2143
+ ,235454
+ ,73
+ ,819
+ ,46
+ ,0
+ ,151
+ ,32
+ ,124
+ ,106
+ ,89275
+ ,24164
+ ,2146
+ ,220801
+ ,75
+ ,720
+ ,105
+ ,1
+ ,51
+ ,18
+ ,72
+ ,63
+ ,44418
+ ,12607
+ ,874
+ ,99466
+ ,50
+ ,273
+ ,32
+ ,0
+ ,28
+ ,23
+ ,91
+ ,69
+ ,192565
+ ,22609
+ ,1590
+ ,92661
+ ,61
+ ,508
+ ,133
+ ,1
+ ,40
+ ,17
+ ,45
+ ,41
+ ,35232
+ ,5892
+ ,1590
+ ,133328
+ ,55
+ ,506
+ ,79
+ ,0
+ ,56
+ ,20
+ ,78
+ ,56
+ ,40909
+ ,17014
+ ,1210
+ ,61361
+ ,77
+ ,451
+ ,51
+ ,0
+ ,27
+ ,12
+ ,39
+ ,25
+ ,13294
+ ,5394
+ ,2072
+ ,125930
+ ,75
+ ,699
+ ,207
+ ,4
+ ,37
+ ,17
+ ,68
+ ,65
+ ,32387
+ ,9178
+ ,1281
+ ,100750
+ ,72
+ ,407
+ ,67
+ ,0
+ ,83
+ ,30
+ ,119
+ ,93
+ ,140867
+ ,6440
+ ,1401
+ ,224549
+ ,50
+ ,465
+ ,47
+ ,4
+ ,54
+ ,31
+ ,117
+ ,114
+ ,120662
+ ,21916
+ ,834
+ ,82316
+ ,32
+ ,245
+ ,34
+ ,4
+ ,27
+ ,10
+ ,39
+ ,38
+ ,21233
+ ,4011
+ ,1105
+ ,102010
+ ,53
+ ,370
+ ,66
+ ,3
+ ,28
+ ,13
+ ,50
+ ,44
+ ,44332
+ ,5818
+ ,1272
+ ,101523
+ ,42
+ ,316
+ ,76
+ ,0
+ ,59
+ ,22
+ ,88
+ ,87
+ ,61056
+ ,18647
+ ,1944
+ ,243511
+ ,71
+ ,603
+ ,65
+ ,0
+ ,133
+ ,42
+ ,155
+ ,110
+ ,101338
+ ,20556
+ ,391
+ ,22938
+ ,10
+ ,154
+ ,9
+ ,0
+ ,12
+ ,1
+ ,0
+ ,0
+ ,1168
+ ,238
+ ,761
+ ,41566
+ ,35
+ ,229
+ ,42
+ ,5
+ ,0
+ ,9
+ ,36
+ ,27
+ ,13497
+ ,70
+ ,1605
+ ,152474
+ ,65
+ ,577
+ ,45
+ ,0
+ ,106
+ ,32
+ ,123
+ ,83
+ ,65567
+ ,22392
+ ,530
+ ,61857
+ ,25
+ ,192
+ ,25
+ ,4
+ ,23
+ ,11
+ ,32
+ ,30
+ ,25162
+ ,3913
+ ,1988
+ ,99923
+ ,66
+ ,617
+ ,115
+ ,0
+ ,44
+ ,25
+ ,99
+ ,80
+ ,32334
+ ,12237
+ ,1386
+ ,132487
+ ,41
+ ,411
+ ,97
+ ,0
+ ,71
+ ,36
+ ,136
+ ,98
+ ,40735
+ ,8388
+ ,2395
+ ,317394
+ ,86
+ ,975
+ ,53
+ ,1
+ ,116
+ ,31
+ ,117
+ ,82
+ ,91413
+ ,22120
+ ,387
+ ,21054
+ ,16
+ ,146
+ ,2
+ ,0
+ ,4
+ ,0
+ ,0
+ ,0
+ ,855
+ ,338
+ ,1742
+ ,209641
+ ,42
+ ,705
+ ,52
+ ,5
+ ,62
+ ,24
+ ,88
+ ,60
+ ,97068
+ ,11727
+ ,620
+ ,22648
+ ,19
+ ,184
+ ,44
+ ,0
+ ,12
+ ,13
+ ,39
+ ,28
+ ,44339
+ ,3704
+ ,449
+ ,31414
+ ,19
+ ,200
+ ,22
+ ,0
+ ,18
+ ,8
+ ,25
+ ,9
+ ,14116
+ ,3988
+ ,800
+ ,46698
+ ,45
+ ,274
+ ,35
+ ,0
+ ,14
+ ,13
+ ,52
+ ,33
+ ,10288
+ ,3030
+ ,1684
+ ,131698
+ ,65
+ ,502
+ ,74
+ ,0
+ ,60
+ ,19
+ ,75
+ ,59
+ ,65622
+ ,13520
+ ,1050
+ ,91735
+ ,35
+ ,382
+ ,103
+ ,0
+ ,7
+ ,18
+ ,71
+ ,49
+ ,16563
+ ,1421
+ ,2699
+ ,244749
+ ,95
+ ,964
+ ,144
+ ,2
+ ,98
+ ,33
+ ,124
+ ,115
+ ,76643
+ ,20923
+ ,1606
+ ,184510
+ ,49
+ ,537
+ ,60
+ ,7
+ ,64
+ ,40
+ ,151
+ ,140
+ ,110681
+ ,20237
+ ,1502
+ ,79863
+ ,37
+ ,438
+ ,134
+ ,1
+ ,29
+ ,22
+ ,71
+ ,49
+ ,29011
+ ,3219
+ ,1204
+ ,128423
+ ,64
+ ,369
+ ,89
+ ,8
+ ,32
+ ,38
+ ,145
+ ,120
+ ,92696
+ ,3769
+ ,1138
+ ,97839
+ ,38
+ ,417
+ ,42
+ ,2
+ ,25
+ ,24
+ ,87
+ ,66
+ ,94785
+ ,12252
+ ,568
+ ,38214
+ ,34
+ ,276
+ ,52
+ ,0
+ ,16
+ ,8
+ ,27
+ ,21
+ ,8773
+ ,1888
+ ,1459
+ ,151101
+ ,32
+ ,514
+ ,98
+ ,2
+ ,48
+ ,35
+ ,131
+ ,124
+ ,83209
+ ,14497
+ ,2158
+ ,272458
+ ,65
+ ,822
+ ,99
+ ,0
+ ,100
+ ,43
+ ,162
+ ,152
+ ,93815
+ ,28864
+ ,1111
+ ,172494
+ ,52
+ ,389
+ ,52
+ ,0
+ ,46
+ ,43
+ ,165
+ ,139
+ ,86687
+ ,21721
+ ,1421
+ ,108043
+ ,62
+ ,466
+ ,29
+ ,1
+ ,45
+ ,14
+ ,54
+ ,38
+ ,34553
+ ,4821
+ ,2833
+ ,328107
+ ,65
+ ,1255
+ ,125
+ ,3
+ ,129
+ ,41
+ ,159
+ ,144
+ ,105547
+ ,33644
+ ,1955
+ ,250579
+ ,83
+ ,694
+ ,106
+ ,0
+ ,130
+ ,38
+ ,147
+ ,120
+ ,103487
+ ,15923
+ ,2922
+ ,351067
+ ,95
+ ,1024
+ ,95
+ ,3
+ ,136
+ ,45
+ ,170
+ ,160
+ ,213688
+ ,42935
+ ,1002
+ ,158015
+ ,29
+ ,400
+ ,40
+ ,0
+ ,59
+ ,31
+ ,119
+ ,114
+ ,71220
+ ,18864
+ ,1060
+ ,98866
+ ,18
+ ,397
+ ,140
+ ,0
+ ,25
+ ,13
+ ,49
+ ,39
+ ,23517
+ ,4977
+ ,956
+ ,85439
+ ,33
+ ,350
+ ,43
+ ,0
+ ,32
+ ,28
+ ,104
+ ,78
+ ,56926
+ ,7785
+ ,2186
+ ,229242
+ ,247
+ ,719
+ ,128
+ ,4
+ ,63
+ ,31
+ ,120
+ ,119
+ ,91721
+ ,17939
+ ,3604
+ ,351619
+ ,139
+ ,1277
+ ,142
+ ,4
+ ,95
+ ,40
+ ,150
+ ,141
+ ,115168
+ ,23436
+ ,1035
+ ,84207
+ ,29
+ ,356
+ ,73
+ ,11
+ ,14
+ ,30
+ ,112
+ ,101
+ ,111194
+ ,325
+ ,1417
+ ,120445
+ ,118
+ ,457
+ ,72
+ ,0
+ ,36
+ ,16
+ ,59
+ ,56
+ ,51009
+ ,13539
+ ,3261
+ ,324598
+ ,110
+ ,1402
+ ,128
+ ,0
+ ,113
+ ,37
+ ,136
+ ,133
+ ,135777
+ ,34538
+ ,1587
+ ,131069
+ ,67
+ ,600
+ ,61
+ ,4
+ ,47
+ ,30
+ ,107
+ ,83
+ ,51513
+ ,12198
+ ,1424
+ ,204271
+ ,42
+ ,480
+ ,73
+ ,0
+ ,92
+ ,35
+ ,130
+ ,116
+ ,74163
+ ,26924
+ ,1701
+ ,165543
+ ,65
+ ,595
+ ,148
+ ,1
+ ,70
+ ,32
+ ,115
+ ,90
+ ,51633
+ ,12716
+ ,1249
+ ,141722
+ ,94
+ ,436
+ ,64
+ ,0
+ ,19
+ ,27
+ ,107
+ ,36
+ ,75345
+ ,8172
+ ,946
+ ,116048
+ ,64
+ ,230
+ ,45
+ ,0
+ ,50
+ ,20
+ ,75
+ ,50
+ ,33416
+ ,10855
+ ,1926
+ ,250047
+ ,81
+ ,651
+ ,58
+ ,0
+ ,41
+ ,18
+ ,71
+ ,61
+ ,83305
+ ,11932
+ ,3352
+ ,299775
+ ,95
+ ,1367
+ ,97
+ ,9
+ ,91
+ ,31
+ ,120
+ ,97
+ ,98952
+ ,14300
+ ,1641
+ ,195838
+ ,67
+ ,564
+ ,50
+ ,1
+ ,111
+ ,31
+ ,116
+ ,98
+ ,102372
+ ,25515
+ ,2035
+ ,173260
+ ,63
+ ,716
+ ,37
+ ,3
+ ,41
+ ,21
+ ,79
+ ,78
+ ,37238
+ ,2805
+ ,2312
+ ,254488
+ ,83
+ ,747
+ ,50
+ ,10
+ ,120
+ ,39
+ ,150
+ ,117
+ ,103772
+ ,29402
+ ,1369
+ ,104389
+ ,45
+ ,467
+ ,105
+ ,5
+ ,135
+ ,41
+ ,156
+ ,148
+ ,123969
+ ,16440
+ ,1577
+ ,136084
+ ,30
+ ,671
+ ,69
+ ,0
+ ,27
+ ,13
+ ,51
+ ,41
+ ,27142
+ ,11221
+ ,2201
+ ,199476
+ ,70
+ ,861
+ ,46
+ ,2
+ ,87
+ ,32
+ ,118
+ ,105
+ ,135400
+ ,28732
+ ,961
+ ,92499
+ ,32
+ ,319
+ ,57
+ ,0
+ ,25
+ ,18
+ ,71
+ ,55
+ ,21399
+ ,5250
+ ,1900
+ ,224330
+ ,83
+ ,612
+ ,52
+ ,1
+ ,131
+ ,39
+ ,144
+ ,132
+ ,130115
+ ,28608
+ ,1254
+ ,135781
+ ,31
+ ,433
+ ,98
+ ,2
+ ,45
+ ,14
+ ,47
+ ,44
+ ,24874
+ ,8092
+ ,1335
+ ,74408
+ ,67
+ ,434
+ ,61
+ ,4
+ ,29
+ ,7
+ ,28
+ ,21
+ ,34988
+ ,4473
+ ,1597
+ ,81240
+ ,66
+ ,503
+ ,89
+ ,0
+ ,58
+ ,17
+ ,68
+ ,50
+ ,45549
+ ,1572
+ ,207
+ ,14688
+ ,10
+ ,85
+ ,0
+ ,0
+ ,4
+ ,0
+ ,0
+ ,0
+ ,6023
+ ,2065
+ ,1645
+ ,181633
+ ,70
+ ,564
+ ,48
+ ,2
+ ,47
+ ,30
+ ,110
+ ,73
+ ,64466
+ ,14817
+ ,2429
+ ,271856
+ ,103
+ ,824
+ ,91
+ ,1
+ ,109
+ ,37
+ ,147
+ ,86
+ ,54990
+ ,16714
+ ,151
+ ,7199
+ ,5
+ ,74
+ ,0
+ ,0
+ ,7
+ ,0
+ ,0
+ ,0
+ ,1644
+ ,556
+ ,474
+ ,46660
+ ,20
+ ,259
+ ,7
+ ,0
+ ,12
+ ,5
+ ,15
+ ,13
+ ,6179
+ ,2089
+ ,141
+ ,17547
+ ,5
+ ,69
+ ,3
+ ,0
+ ,0
+ ,1
+ ,4
+ ,4
+ ,3926
+ ,2658
+ ,1639
+ ,133368
+ ,36
+ ,535
+ ,54
+ ,1
+ ,37
+ ,16
+ ,64
+ ,57
+ ,32755
+ ,10695
+ ,872
+ ,95227
+ ,34
+ ,239
+ ,70
+ ,0
+ ,37
+ ,32
+ ,111
+ ,48
+ ,34777
+ ,1669
+ ,1318
+ ,152601
+ ,48
+ ,438
+ ,36
+ ,2
+ ,46
+ ,24
+ ,85
+ ,46
+ ,73224
+ ,16267
+ ,1018
+ ,98146
+ ,40
+ ,459
+ ,37
+ ,0
+ ,15
+ ,17
+ ,68
+ ,48
+ ,27114
+ ,7768
+ ,1383
+ ,79619
+ ,43
+ ,426
+ ,123
+ ,3
+ ,42
+ ,11
+ ,40
+ ,32
+ ,20760
+ ,7252
+ ,1314
+ ,59194
+ ,31
+ ,288
+ ,247
+ ,6
+ ,7
+ ,24
+ ,80
+ ,68
+ ,37636
+ ,6387
+ ,1335
+ ,139942
+ ,42
+ ,498
+ ,46
+ ,0
+ ,54
+ ,22
+ ,88
+ ,87
+ ,65461
+ ,18715
+ ,1403
+ ,118612
+ ,46
+ ,454
+ ,72
+ ,2
+ ,54
+ ,12
+ ,48
+ ,43
+ ,30080
+ ,7936
+ ,910
+ ,72880
+ ,33
+ ,376
+ ,41
+ ,0
+ ,14
+ ,19
+ ,76
+ ,67
+ ,24094
+ ,8643
+ ,616
+ ,65475
+ ,18
+ ,225
+ ,24
+ ,2
+ ,16
+ ,13
+ ,51
+ ,46
+ ,69008
+ ,7294
+ ,1407
+ ,99643
+ ,55
+ ,555
+ ,45
+ ,1
+ ,33
+ ,17
+ ,67
+ ,46
+ ,54968
+ ,4570
+ ,771
+ ,71965
+ ,35
+ ,252
+ ,33
+ ,1
+ ,32
+ ,15
+ ,59
+ ,56
+ ,46090
+ ,7185
+ ,766
+ ,77272
+ ,59
+ ,208
+ ,27
+ ,2
+ ,21
+ ,16
+ ,61
+ ,48
+ ,27507
+ ,10058
+ ,473
+ ,49289
+ ,19
+ ,130
+ ,36
+ ,1
+ ,15
+ ,24
+ ,76
+ ,44
+ ,10672
+ ,2342
+ ,1376
+ ,135131
+ ,66
+ ,481
+ ,87
+ ,0
+ ,38
+ ,15
+ ,60
+ ,60
+ ,34029
+ ,8509
+ ,1232
+ ,108446
+ ,60
+ ,389
+ ,90
+ ,1
+ ,22
+ ,17
+ ,68
+ ,65
+ ,46300
+ ,13275
+ ,1521
+ ,89746
+ ,36
+ ,565
+ ,114
+ ,3
+ ,28
+ ,18
+ ,71
+ ,55
+ ,24760
+ ,6816
+ ,572
+ ,44296
+ ,25
+ ,173
+ ,31
+ ,0
+ ,10
+ ,20
+ ,76
+ ,38
+ ,18779
+ ,1930
+ ,1059
+ ,77648
+ ,47
+ ,278
+ ,45
+ ,0
+ ,31
+ ,16
+ ,62
+ ,52
+ ,21280
+ ,8086
+ ,1544
+ ,181528
+ ,54
+ ,609
+ ,69
+ ,0
+ ,32
+ ,16
+ ,61
+ ,60
+ ,40662
+ ,10737
+ ,1230
+ ,134019
+ ,53
+ ,422
+ ,51
+ ,0
+ ,32
+ ,18
+ ,67
+ ,54
+ ,28987
+ ,8033
+ ,1206
+ ,124064
+ ,40
+ ,445
+ ,34
+ ,1
+ ,43
+ ,22
+ ,88
+ ,86
+ ,22827
+ ,7058
+ ,1205
+ ,92630
+ ,40
+ ,387
+ ,60
+ ,4
+ ,27
+ ,8
+ ,30
+ ,24
+ ,18513
+ ,6782
+ ,1255
+ ,121848
+ ,39
+ ,339
+ ,45
+ ,0
+ ,37
+ ,17
+ ,64
+ ,52
+ ,30594
+ ,5401
+ ,613
+ ,52915
+ ,14
+ ,181
+ ,54
+ ,0
+ ,20
+ ,18
+ ,68
+ ,49
+ ,24006
+ ,6521
+ ,721
+ ,81872
+ ,45
+ ,245
+ ,25
+ ,0
+ ,32
+ ,16
+ ,64
+ ,61
+ ,27913
+ ,10856
+ ,1109
+ ,58981
+ ,36
+ ,384
+ ,38
+ ,7
+ ,0
+ ,23
+ ,91
+ ,61
+ ,42744
+ ,2154
+ ,740
+ ,53515
+ ,28
+ ,212
+ ,52
+ ,2
+ ,5
+ ,22
+ ,88
+ ,81
+ ,12934
+ ,6117
+ ,1126
+ ,60812
+ ,44
+ ,399
+ ,67
+ ,0
+ ,26
+ ,13
+ ,52
+ ,43
+ ,22574
+ ,5238
+ ,728
+ ,56375
+ ,30
+ ,229
+ ,74
+ ,7
+ ,10
+ ,13
+ ,49
+ ,40
+ ,41385
+ ,4820
+ ,689
+ ,65490
+ ,22
+ ,224
+ ,38
+ ,3
+ ,27
+ ,16
+ ,62
+ ,40
+ ,18653
+ ,5615
+ ,592
+ ,80949
+ ,17
+ ,203
+ ,30
+ ,0
+ ,11
+ ,16
+ ,61
+ ,56
+ ,18472
+ ,4272
+ ,995
+ ,76302
+ ,31
+ ,333
+ ,26
+ ,0
+ ,29
+ ,20
+ ,76
+ ,68
+ ,30976
+ ,8702
+ ,1613
+ ,104011
+ ,55
+ ,384
+ ,67
+ ,6
+ ,25
+ ,22
+ ,88
+ ,79
+ ,63339
+ ,15340
+ ,2048
+ ,98104
+ ,54
+ ,636
+ ,132
+ ,2
+ ,55
+ ,17
+ ,66
+ ,47
+ ,25568
+ ,8030
+ ,705
+ ,67989
+ ,21
+ ,185
+ ,42
+ ,0
+ ,23
+ ,18
+ ,71
+ ,57
+ ,33747
+ ,9526
+ ,301
+ ,30989
+ ,14
+ ,93
+ ,35
+ ,0
+ ,5
+ ,17
+ ,68
+ ,41
+ ,4154
+ ,1278
+ ,1803
+ ,135458
+ ,81
+ ,581
+ ,118
+ ,3
+ ,43
+ ,12
+ ,48
+ ,29
+ ,19474
+ ,4236
+ ,799
+ ,73504
+ ,35
+ ,248
+ ,68
+ ,0
+ ,23
+ ,7
+ ,25
+ ,3
+ ,35130
+ ,3023
+ ,861
+ ,63123
+ ,43
+ ,304
+ ,43
+ ,1
+ ,34
+ ,17
+ ,68
+ ,60
+ ,39067
+ ,7196
+ ,1186
+ ,61254
+ ,46
+ ,344
+ ,76
+ ,1
+ ,36
+ ,14
+ ,41
+ ,30
+ ,13310
+ ,3394
+ ,1451
+ ,74914
+ ,30
+ ,407
+ ,64
+ ,0
+ ,35
+ ,23
+ ,90
+ ,79
+ ,65892
+ ,6371
+ ,628
+ ,31774
+ ,23
+ ,170
+ ,48
+ ,1
+ ,0
+ ,17
+ ,66
+ ,47
+ ,4143
+ ,1574
+ ,1161
+ ,81437
+ ,38
+ ,312
+ ,64
+ ,0
+ ,37
+ ,14
+ ,54
+ ,40
+ ,28579
+ ,9620
+ ,1463
+ ,87186
+ ,54
+ ,507
+ ,56
+ ,0
+ ,28
+ ,15
+ ,59
+ ,48
+ ,51776
+ ,6978
+ ,742
+ ,50090
+ ,20
+ ,224
+ ,71
+ ,0
+ ,16
+ ,17
+ ,60
+ ,36
+ ,21152
+ ,4911
+ ,979
+ ,65745
+ ,53
+ ,340
+ ,75
+ ,0
+ ,26
+ ,21
+ ,77
+ ,42
+ ,38084
+ ,8645
+ ,675
+ ,56653
+ ,45
+ ,168
+ ,39
+ ,0
+ ,38
+ ,18
+ ,68
+ ,49
+ ,27717
+ ,8987
+ ,1241
+ ,158399
+ ,39
+ ,443
+ ,42
+ ,0
+ ,23
+ ,18
+ ,72
+ ,57
+ ,32928
+ ,5544
+ ,676
+ ,46455
+ ,20
+ ,204
+ ,39
+ ,0
+ ,22
+ ,17
+ ,67
+ ,12
+ ,11342
+ ,3083
+ ,1049
+ ,73624
+ ,24
+ ,367
+ ,93
+ ,0
+ ,30
+ ,17
+ ,64
+ ,40
+ ,19499
+ ,6909
+ ,620
+ ,38395
+ ,31
+ ,210
+ ,38
+ ,0
+ ,16
+ ,16
+ ,63
+ ,43
+ ,16380
+ ,3189
+ ,1081
+ ,91899
+ ,35
+ ,335
+ ,60
+ ,0
+ ,18
+ ,15
+ ,59
+ ,33
+ ,36874
+ ,6745
+ ,1688
+ ,139526
+ ,151
+ ,364
+ ,71
+ ,0
+ ,28
+ ,21
+ ,84
+ ,77
+ ,48259
+ ,16724
+ ,736
+ ,52164
+ ,52
+ ,178
+ ,52
+ ,0
+ ,32
+ ,16
+ ,64
+ ,43
+ ,16734
+ ,4850
+ ,617
+ ,51567
+ ,30
+ ,206
+ ,27
+ ,2
+ ,21
+ ,14
+ ,56
+ ,45
+ ,28207
+ ,7025
+ ,812
+ ,70551
+ ,31
+ ,279
+ ,59
+ ,0
+ ,23
+ ,15
+ ,54
+ ,47
+ ,30143
+ ,6047
+ ,1051
+ ,84856
+ ,29
+ ,387
+ ,40
+ ,1
+ ,29
+ ,17
+ ,67
+ ,43
+ ,41369
+ ,7377
+ ,1656
+ ,102538
+ ,57
+ ,490
+ ,79
+ ,1
+ ,50
+ ,15
+ ,58
+ ,45
+ ,45833
+ ,9078
+ ,705
+ ,86678
+ ,40
+ ,238
+ ,44
+ ,0
+ ,12
+ ,15
+ ,59
+ ,50
+ ,29156
+ ,4605
+ ,945
+ ,85709
+ ,44
+ ,343
+ ,65
+ ,0
+ ,21
+ ,10
+ ,40
+ ,35
+ ,35944
+ ,3238
+ ,554
+ ,34662
+ ,25
+ ,232
+ ,10
+ ,0
+ ,18
+ ,6
+ ,22
+ ,7
+ ,36278
+ ,8100
+ ,1597
+ ,150580
+ ,77
+ ,530
+ ,124
+ ,0
+ ,27
+ ,22
+ ,83
+ ,71
+ ,45588
+ ,9653
+ ,982
+ ,99611
+ ,35
+ ,291
+ ,81
+ ,0
+ ,41
+ ,21
+ ,81
+ ,67
+ ,45097
+ ,8914
+ ,222
+ ,19349
+ ,11
+ ,67
+ ,15
+ ,0
+ ,13
+ ,1
+ ,2
+ ,0
+ ,3895
+ ,786
+ ,1212
+ ,99373
+ ,63
+ ,397
+ ,92
+ ,1
+ ,12
+ ,18
+ ,72
+ ,62
+ ,28394
+ ,6700
+ ,1143
+ ,86230
+ ,44
+ ,467
+ ,42
+ ,0
+ ,21
+ ,17
+ ,61
+ ,54
+ ,18632
+ ,5788
+ ,435
+ ,30837
+ ,19
+ ,178
+ ,10
+ ,0
+ ,8
+ ,4
+ ,15
+ ,4
+ ,2325
+ ,593
+ ,532
+ ,31706
+ ,13
+ ,175
+ ,24
+ ,0
+ ,26
+ ,10
+ ,32
+ ,25
+ ,25139
+ ,4506
+ ,882
+ ,89806
+ ,42
+ ,299
+ ,64
+ ,0
+ ,27
+ ,16
+ ,62
+ ,40
+ ,27975
+ ,6382
+ ,608
+ ,62088
+ ,38
+ ,154
+ ,45
+ ,1
+ ,13
+ ,16
+ ,58
+ ,38
+ ,14483
+ ,5621
+ ,459
+ ,40151
+ ,29
+ ,106
+ ,22
+ ,0
+ ,16
+ ,9
+ ,36
+ ,19
+ ,13127
+ ,3997
+ ,578
+ ,27634
+ ,20
+ ,189
+ ,56
+ ,0
+ ,2
+ ,16
+ ,59
+ ,17
+ ,5839
+ ,520
+ ,826
+ ,76990
+ ,27
+ ,194
+ ,94
+ ,0
+ ,42
+ ,17
+ ,68
+ ,67
+ ,24069
+ ,8891
+ ,509
+ ,37460
+ ,20
+ ,135
+ ,19
+ ,0
+ ,5
+ ,7
+ ,21
+ ,14
+ ,3738
+ ,999
+ ,717
+ ,54157
+ ,19
+ ,201
+ ,35
+ ,0
+ ,37
+ ,15
+ ,55
+ ,30
+ ,18625
+ ,7067
+ ,637
+ ,49862
+ ,37
+ ,207
+ ,32
+ ,0
+ ,17
+ ,14
+ ,54
+ ,54
+ ,36341
+ ,4639
+ ,857
+ ,84337
+ ,26
+ ,280
+ ,35
+ ,0
+ ,38
+ ,14
+ ,55
+ ,35
+ ,24548
+ ,5654
+ ,830
+ ,64175
+ ,42
+ ,260
+ ,48
+ ,0
+ ,37
+ ,18
+ ,72
+ ,59
+ ,21792
+ ,6928
+ ,652
+ ,59382
+ ,49
+ ,227
+ ,49
+ ,0
+ ,29
+ ,12
+ ,41
+ ,24
+ ,26263
+ ,1514
+ ,707
+ ,119308
+ ,30
+ ,239
+ ,48
+ ,0
+ ,32
+ ,16
+ ,61
+ ,58
+ ,23686
+ ,9238
+ ,954
+ ,76702
+ ,49
+ ,333
+ ,62
+ ,0
+ ,35
+ ,21
+ ,67
+ ,42
+ ,49303
+ ,8204
+ ,1461
+ ,103425
+ ,67
+ ,428
+ ,96
+ ,1
+ ,17
+ ,19
+ ,76
+ ,46
+ ,25659
+ ,5926
+ ,672
+ ,70344
+ ,28
+ ,230
+ ,45
+ ,0
+ ,20
+ ,16
+ ,64
+ ,61
+ ,28904
+ ,5785
+ ,778
+ ,43410
+ ,19
+ ,292
+ ,63
+ ,0
+ ,7
+ ,1
+ ,3
+ ,3
+ ,2781
+ ,4
+ ,1141
+ ,104838
+ ,49
+ ,350
+ ,71
+ ,1
+ ,46
+ ,16
+ ,63
+ ,52
+ ,29236
+ ,5930
+ ,680
+ ,62215
+ ,27
+ ,186
+ ,26
+ ,0
+ ,24
+ ,10
+ ,40
+ ,25
+ ,19546
+ ,3710
+ ,1090
+ ,69304
+ ,30
+ ,326
+ ,48
+ ,6
+ ,40
+ ,19
+ ,69
+ ,40
+ ,22818
+ ,705
+ ,616
+ ,53117
+ ,22
+ ,155
+ ,29
+ ,3
+ ,3
+ ,12
+ ,48
+ ,32
+ ,32689
+ ,443
+ ,285
+ ,19764
+ ,12
+ ,75
+ ,19
+ ,1
+ ,10
+ ,2
+ ,8
+ ,4
+ ,5752
+ ,2416
+ ,1145
+ ,86680
+ ,31
+ ,361
+ ,45
+ ,2
+ ,37
+ ,14
+ ,52
+ ,49
+ ,22197
+ ,7747
+ ,733
+ ,84105
+ ,20
+ ,261
+ ,45
+ ,0
+ ,17
+ ,17
+ ,66
+ ,63
+ ,20055
+ ,5432
+ ,888
+ ,77945
+ ,20
+ ,299
+ ,67
+ ,0
+ ,28
+ ,19
+ ,76
+ ,67
+ ,25272
+ ,4913
+ ,849
+ ,89113
+ ,39
+ ,300
+ ,30
+ ,0
+ ,19
+ ,14
+ ,43
+ ,32
+ ,82206
+ ,2650
+ ,1182
+ ,91005
+ ,29
+ ,450
+ ,36
+ ,3
+ ,29
+ ,11
+ ,39
+ ,23
+ ,32073
+ ,2370
+ ,528
+ ,40248
+ ,16
+ ,183
+ ,34
+ ,1
+ ,8
+ ,4
+ ,14
+ ,7
+ ,5444
+ ,775
+ ,642
+ ,64187
+ ,27
+ ,238
+ ,36
+ ,0
+ ,10
+ ,16
+ ,61
+ ,54
+ ,20154
+ ,5576
+ ,947
+ ,50857
+ ,21
+ ,165
+ ,34
+ ,0
+ ,15
+ ,20
+ ,71
+ ,37
+ ,36944
+ ,1352
+ ,819
+ ,56613
+ ,19
+ ,234
+ ,37
+ ,1
+ ,15
+ ,12
+ ,44
+ ,35
+ ,8019
+ ,3080
+ ,757
+ ,62792
+ ,35
+ ,176
+ ,46
+ ,0
+ ,28
+ ,15
+ ,60
+ ,51
+ ,30884
+ ,10205
+ ,894
+ ,72535
+ ,14
+ ,329
+ ,44
+ ,0
+ ,17
+ ,16
+ ,64
+ ,39
+ ,19540
+ ,6095)
+ ,dim=c(12
+ ,289)
+ ,dimnames=list(c('pageviews'
+ ,'time_in_rfc'
+ ,'logins'
+ ,'compendium_views_info'
+ ,'compendium_views_pr'
+ ,'shared_compendiums'
+ ,'blogged_computations'
+ ,'compendiums_reviewed'
+ ,'feedback_messages_p1'
+ ,'feedback_messages_p120'
+ ,'totsize'
+ ,'totrevisions
')
+ ,1:289))
> y <- array(NA,dim=c(12,289),dimnames=list(c('pageviews','time_in_rfc','logins','compendium_views_info','compendium_views_pr','shared_compendiums','blogged_computations','compendiums_reviewed','feedback_messages_p1','feedback_messages_p120','totsize','totrevisions
'),1:289))
> for (i in 1:dim(x)[1])
+ {
+ for (j in 1:dim(x)[2])
+ {
+ y[i,j] <- as.numeric(x[i,j])
+ }
+ }
> par3 = 'No Linear Trend'
> par2 = 'Do not include Seasonal Dummies'
> par1 = '11'
> par3 <- 'No Linear Trend'
> par2 <- 'Do not include Seasonal Dummies'
> par1 <- '11'
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Prof. Dr. P. Wessa
> #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #Technical description: Write here your technical program description (don't use hard returns!)
> library(lattice)
> library(lmtest)
Loading required package: zoo
Attaching package: 'zoo'
The following object(s) are masked from 'package:base':
as.Date, as.Date.numeric
> n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test
> par1 <- as.numeric(par1)
> x <- t(y)
> k <- length(x[1,])
> n <- length(x[,1])
> x1 <- cbind(x[,par1], x[,1:k!=par1])
> mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1])
> colnames(x1) <- mycolnames #colnames(x)[par1]
> x <- x1
> if (par3 == 'First Differences'){
+ x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep='')))
+ for (i in 1:n-1) {
+ for (j in 1:k) {
+ x2[i,j] <- x[i+1,j] - x[i,j]
+ }
+ }
+ x <- x2
+ }
> if (par2 == 'Include Monthly Dummies'){
+ x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep ='')))
+ for (i in 1:11){
+ x2[seq(i,n,12),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> if (par2 == 'Include Quarterly Dummies'){
+ x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep ='')))
+ for (i in 1:3){
+ x2[seq(i,n,4),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> k <- length(x[1,])
> if (par3 == 'Linear Trend'){
+ x <- cbind(x, c(1:n))
+ colnames(x)[k+1] <- 't'
+ }
> x
totsize pageviews time_in_rfc logins compendium_views_info
1 112285 1418 210907 56 396
2 84786 869 120982 56 297
3 83123 1530 176508 54 559
4 101193 2172 179321 89 967
5 38361 901 123185 40 270
6 68504 463 52746 25 143
7 119182 3201 385534 92 1562
8 22807 371 33170 18 109
9 17140 1192 101645 63 371
10 116174 1583 149061 44 656
11 57635 1439 165446 33 511
12 66198 1764 237213 84 655
13 71701 1495 173326 88 465
14 57793 1373 133131 55 525
15 80444 2187 258873 60 885
16 53855 1491 180083 66 497
17 97668 4041 324799 154 1436
18 133824 1706 230964 53 612
19 101481 2152 236785 119 865
20 99645 1036 135473 41 385
21 114789 1882 202925 61 567
22 99052 1929 215147 58 639
23 67654 2242 344297 75 963
24 65553 1220 153935 33 398
25 97500 1289 132943 40 410
26 69112 2515 174724 92 966
27 82753 2147 174415 100 801
28 85323 2352 225548 112 892
29 72654 1638 223632 73 513
30 30727 1222 124817 40 469
31 77873 1812 221698 45 683
32 117478 1677 210767 60 643
33 74007 1579 170266 62 535
34 90183 1731 260561 75 625
35 61542 807 84853 31 264
36 101494 2452 294424 77 992
37 27570 829 101011 34 238
38 55813 1940 215641 46 818
39 79215 2662 325107 99 937
40 1423 186 7176 17 70
41 55461 1499 167542 66 507
42 31081 865 106408 30 260
43 22996 1793 96560 76 503
44 83122 2527 265769 146 927
45 70106 2747 269651 67 1269
46 60578 1324 149112 56 537
47 39992 2702 175824 107 910
48 79892 1383 152871 58 532
49 49810 1179 111665 34 345
50 71570 2099 116408 61 918
51 100708 4308 362301 119 1635
52 33032 918 78800 42 330
53 82875 1831 183167 66 557
54 139077 3373 277965 89 1178
55 71595 1713 150629 44 740
56 72260 1438 168809 66 452
57 5950 496 24188 24 218
58 115762 2253 329267 259 764
59 32551 744 65029 17 255
60 31701 1161 101097 64 454
61 80670 2352 218946 41 866
62 143558 2144 244052 68 574
63 117105 4691 341570 168 1276
64 23789 1112 103597 43 379
65 120733 2694 233328 132 825
66 105195 1973 256462 105 798
67 73107 1769 206161 71 663
68 132068 3148 311473 112 1069
69 149193 2474 235800 94 921
70 46821 2084 177939 82 858
71 87011 1954 207176 70 711
72 95260 1226 196553 57 503
73 55183 1389 174184 53 382
74 106671 1496 143246 103 464
75 73511 2269 187559 121 717
76 92945 1833 187681 62 690
77 78664 1268 119016 52 462
78 70054 1943 182192 52 657
79 22618 893 73566 32 385
80 74011 1762 194979 62 577
81 83737 1403 167488 45 619
82 69094 1425 143756 46 479
83 93133 1857 275541 63 817
84 95536 1840 243199 75 752
85 225920 1502 182999 88 430
86 62133 1441 135649 46 451
87 61370 1420 152299 53 537
88 43836 1416 120221 37 519
89 106117 2970 346485 90 1000
90 38692 1317 145790 63 637
91 84651 1644 193339 78 465
92 56622 870 80953 25 437
93 15986 1654 122774 45 711
94 95364 1054 130585 46 299
95 26706 937 112611 41 248
96 89691 3004 286468 144 1162
97 67267 2008 241066 82 714
98 126846 2547 148446 91 905
99 41140 1885 204713 71 649
100 102860 1626 182079 63 512
101 51715 1468 140344 53 472
102 55801 2445 220516 62 905
103 111813 1964 243060 63 786
104 120293 1381 162765 32 489
105 138599 1369 182613 39 479
106 161647 1659 232138 62 617
107 115929 2888 265318 117 925
108 24266 1290 85574 34 351
109 162901 2845 310839 92 1144
110 109825 1982 225060 93 669
111 129838 1904 232317 54 707
112 37510 1391 144966 144 458
113 43750 602 43287 14 214
114 40652 1743 155754 61 599
115 87771 1559 164709 109 572
116 85872 2014 201940 38 897
117 89275 2143 235454 73 819
118 44418 2146 220801 75 720
119 192565 874 99466 50 273
120 35232 1590 92661 61 508
121 40909 1590 133328 55 506
122 13294 1210 61361 77 451
123 32387 2072 125930 75 699
124 140867 1281 100750 72 407
125 120662 1401 224549 50 465
126 21233 834 82316 32 245
127 44332 1105 102010 53 370
128 61056 1272 101523 42 316
129 101338 1944 243511 71 603
130 1168 391 22938 10 154
131 13497 761 41566 35 229
132 65567 1605 152474 65 577
133 25162 530 61857 25 192
134 32334 1988 99923 66 617
135 40735 1386 132487 41 411
136 91413 2395 317394 86 975
137 855 387 21054 16 146
138 97068 1742 209641 42 705
139 44339 620 22648 19 184
140 14116 449 31414 19 200
141 10288 800 46698 45 274
142 65622 1684 131698 65 502
143 16563 1050 91735 35 382
144 76643 2699 244749 95 964
145 110681 1606 184510 49 537
146 29011 1502 79863 37 438
147 92696 1204 128423 64 369
148 94785 1138 97839 38 417
149 8773 568 38214 34 276
150 83209 1459 151101 32 514
151 93815 2158 272458 65 822
152 86687 1111 172494 52 389
153 34553 1421 108043 62 466
154 105547 2833 328107 65 1255
155 103487 1955 250579 83 694
156 213688 2922 351067 95 1024
157 71220 1002 158015 29 400
158 23517 1060 98866 18 397
159 56926 956 85439 33 350
160 91721 2186 229242 247 719
161 115168 3604 351619 139 1277
162 111194 1035 84207 29 356
163 51009 1417 120445 118 457
164 135777 3261 324598 110 1402
165 51513 1587 131069 67 600
166 74163 1424 204271 42 480
167 51633 1701 165543 65 595
168 75345 1249 141722 94 436
169 33416 946 116048 64 230
170 83305 1926 250047 81 651
171 98952 3352 299775 95 1367
172 102372 1641 195838 67 564
173 37238 2035 173260 63 716
174 103772 2312 254488 83 747
175 123969 1369 104389 45 467
176 27142 1577 136084 30 671
177 135400 2201 199476 70 861
178 21399 961 92499 32 319
179 130115 1900 224330 83 612
180 24874 1254 135781 31 433
181 34988 1335 74408 67 434
182 45549 1597 81240 66 503
183 6023 207 14688 10 85
184 64466 1645 181633 70 564
185 54990 2429 271856 103 824
186 1644 151 7199 5 74
187 6179 474 46660 20 259
188 3926 141 17547 5 69
189 32755 1639 133368 36 535
190 34777 872 95227 34 239
191 73224 1318 152601 48 438
192 27114 1018 98146 40 459
193 20760 1383 79619 43 426
194 37636 1314 59194 31 288
195 65461 1335 139942 42 498
196 30080 1403 118612 46 454
197 24094 910 72880 33 376
198 69008 616 65475 18 225
199 54968 1407 99643 55 555
200 46090 771 71965 35 252
201 27507 766 77272 59 208
202 10672 473 49289 19 130
203 34029 1376 135131 66 481
204 46300 1232 108446 60 389
205 24760 1521 89746 36 565
206 18779 572 44296 25 173
207 21280 1059 77648 47 278
208 40662 1544 181528 54 609
209 28987 1230 134019 53 422
210 22827 1206 124064 40 445
211 18513 1205 92630 40 387
212 30594 1255 121848 39 339
213 24006 613 52915 14 181
214 27913 721 81872 45 245
215 42744 1109 58981 36 384
216 12934 740 53515 28 212
217 22574 1126 60812 44 399
218 41385 728 56375 30 229
219 18653 689 65490 22 224
220 18472 592 80949 17 203
221 30976 995 76302 31 333
222 63339 1613 104011 55 384
223 25568 2048 98104 54 636
224 33747 705 67989 21 185
225 4154 301 30989 14 93
226 19474 1803 135458 81 581
227 35130 799 73504 35 248
228 39067 861 63123 43 304
229 13310 1186 61254 46 344
230 65892 1451 74914 30 407
231 4143 628 31774 23 170
232 28579 1161 81437 38 312
233 51776 1463 87186 54 507
234 21152 742 50090 20 224
235 38084 979 65745 53 340
236 27717 675 56653 45 168
237 32928 1241 158399 39 443
238 11342 676 46455 20 204
239 19499 1049 73624 24 367
240 16380 620 38395 31 210
241 36874 1081 91899 35 335
242 48259 1688 139526 151 364
243 16734 736 52164 52 178
244 28207 617 51567 30 206
245 30143 812 70551 31 279
246 41369 1051 84856 29 387
247 45833 1656 102538 57 490
248 29156 705 86678 40 238
249 35944 945 85709 44 343
250 36278 554 34662 25 232
251 45588 1597 150580 77 530
252 45097 982 99611 35 291
253 3895 222 19349 11 67
254 28394 1212 99373 63 397
255 18632 1143 86230 44 467
256 2325 435 30837 19 178
257 25139 532 31706 13 175
258 27975 882 89806 42 299
259 14483 608 62088 38 154
260 13127 459 40151 29 106
261 5839 578 27634 20 189
262 24069 826 76990 27 194
263 3738 509 37460 20 135
264 18625 717 54157 19 201
265 36341 637 49862 37 207
266 24548 857 84337 26 280
267 21792 830 64175 42 260
268 26263 652 59382 49 227
269 23686 707 119308 30 239
270 49303 954 76702 49 333
271 25659 1461 103425 67 428
272 28904 672 70344 28 230
273 2781 778 43410 19 292
274 29236 1141 104838 49 350
275 19546 680 62215 27 186
276 22818 1090 69304 30 326
277 32689 616 53117 22 155
278 5752 285 19764 12 75
279 22197 1145 86680 31 361
280 20055 733 84105 20 261
281 25272 888 77945 20 299
282 82206 849 89113 39 300
283 32073 1182 91005 29 450
284 5444 528 40248 16 183
285 20154 642 64187 27 238
286 36944 947 50857 21 165
287 8019 819 56613 19 234
288 30884 757 62792 35 176
289 19540 894 72535 14 329
compendium_views_pr shared_compendiums blogged_computations
1 81 3 79
2 55 4 58
3 50 12 60
4 125 2 108
5 40 1 49
6 37 3 0
7 63 0 121
8 44 0 1
9 88 0 20
10 66 5 43
11 57 0 69
12 74 0 78
13 49 7 86
14 52 7 44
15 88 3 104
16 36 9 63
17 108 0 158
18 43 4 102
19 75 3 77
20 32 0 82
21 44 7 115
22 85 0 101
23 86 1 80
24 56 5 50
25 50 7 83
26 135 0 123
27 63 0 73
28 81 5 81
29 52 0 105
30 44 0 47
31 113 0 105
32 39 3 94
33 73 4 44
34 48 1 114
35 33 4 38
36 59 2 107
37 41 0 30
38 69 0 71
39 64 0 84
40 1 0 0
41 59 2 59
42 32 1 33
43 129 0 42
44 37 2 96
45 31 10 106
46 65 6 56
47 107 0 57
48 74 5 59
49 54 4 39
50 76 1 34
51 715 2 76
52 57 2 20
53 66 0 91
54 106 8 115
55 54 3 85
56 32 0 76
57 20 0 8
58 71 8 79
59 21 5 21
60 70 3 30
61 112 1 76
62 66 5 101
63 190 1 94
64 66 1 27
65 165 5 92
66 56 0 123
67 61 12 75
68 53 8 128
69 127 8 105
70 63 8 55
71 38 8 56
72 50 2 41
73 52 0 72
74 42 5 67
75 76 8 75
76 67 2 114
77 50 5 118
78 53 12 77
79 39 6 22
80 50 7 66
81 77 2 69
82 57 0 105
83 73 4 116
84 34 3 88
85 39 6 73
86 46 2 99
87 63 0 62
88 35 1 53
89 106 0 118
90 43 5 30
91 47 2 100
92 31 0 49
93 162 0 24
94 57 5 67
95 36 0 46
96 263 1 57
97 78 0 75
98 63 1 135
99 54 1 68
100 63 2 124
101 77 6 33
102 79 1 98
103 110 4 58
104 56 2 68
105 56 3 81
106 43 0 131
107 111 10 110
108 71 0 37
109 62 9 130
110 56 7 93
111 74 0 118
112 60 0 39
113 43 4 13
114 68 4 74
115 53 0 81
116 87 0 109
117 46 0 151
118 105 1 51
119 32 0 28
120 133 1 40
121 79 0 56
122 51 0 27
123 207 4 37
124 67 0 83
125 47 4 54
126 34 4 27
127 66 3 28
128 76 0 59
129 65 0 133
130 9 0 12
131 42 5 0
132 45 0 106
133 25 4 23
134 115 0 44
135 97 0 71
136 53 1 116
137 2 0 4
138 52 5 62
139 44 0 12
140 22 0 18
141 35 0 14
142 74 0 60
143 103 0 7
144 144 2 98
145 60 7 64
146 134 1 29
147 89 8 32
148 42 2 25
149 52 0 16
150 98 2 48
151 99 0 100
152 52 0 46
153 29 1 45
154 125 3 129
155 106 0 130
156 95 3 136
157 40 0 59
158 140 0 25
159 43 0 32
160 128 4 63
161 142 4 95
162 73 11 14
163 72 0 36
164 128 0 113
165 61 4 47
166 73 0 92
167 148 1 70
168 64 0 19
169 45 0 50
170 58 0 41
171 97 9 91
172 50 1 111
173 37 3 41
174 50 10 120
175 105 5 135
176 69 0 27
177 46 2 87
178 57 0 25
179 52 1 131
180 98 2 45
181 61 4 29
182 89 0 58
183 0 0 4
184 48 2 47
185 91 1 109
186 0 0 7
187 7 0 12
188 3 0 0
189 54 1 37
190 70 0 37
191 36 2 46
192 37 0 15
193 123 3 42
194 247 6 7
195 46 0 54
196 72 2 54
197 41 0 14
198 24 2 16
199 45 1 33
200 33 1 32
201 27 2 21
202 36 1 15
203 87 0 38
204 90 1 22
205 114 3 28
206 31 0 10
207 45 0 31
208 69 0 32
209 51 0 32
210 34 1 43
211 60 4 27
212 45 0 37
213 54 0 20
214 25 0 32
215 38 7 0
216 52 2 5
217 67 0 26
218 74 7 10
219 38 3 27
220 30 0 11
221 26 0 29
222 67 6 25
223 132 2 55
224 42 0 23
225 35 0 5
226 118 3 43
227 68 0 23
228 43 1 34
229 76 1 36
230 64 0 35
231 48 1 0
232 64 0 37
233 56 0 28
234 71 0 16
235 75 0 26
236 39 0 38
237 42 0 23
238 39 0 22
239 93 0 30
240 38 0 16
241 60 0 18
242 71 0 28
243 52 0 32
244 27 2 21
245 59 0 23
246 40 1 29
247 79 1 50
248 44 0 12
249 65 0 21
250 10 0 18
251 124 0 27
252 81 0 41
253 15 0 13
254 92 1 12
255 42 0 21
256 10 0 8
257 24 0 26
258 64 0 27
259 45 1 13
260 22 0 16
261 56 0 2
262 94 0 42
263 19 0 5
264 35 0 37
265 32 0 17
266 35 0 38
267 48 0 37
268 49 0 29
269 48 0 32
270 62 0 35
271 96 1 17
272 45 0 20
273 63 0 7
274 71 1 46
275 26 0 24
276 48 6 40
277 29 3 3
278 19 1 10
279 45 2 37
280 45 0 17
281 67 0 28
282 30 0 19
283 36 3 29
284 34 1 8
285 36 0 10
286 34 0 15
287 37 1 15
288 46 0 28
289 44 0 17
compendiums_reviewed feedback_messages_p1 feedback_messages_p120
1 30 115 94
2 28 109 103
3 38 146 93
4 30 116 103
5 22 68 51
6 26 101 70
7 25 96 91
8 18 67 22
9 11 44 38
10 26 100 93
11 25 93 60
12 38 140 123
13 44 166 148
14 30 99 90
15 40 139 124
16 34 130 70
17 47 181 168
18 30 116 115
19 31 116 71
20 23 88 66
21 36 139 134
22 36 135 117
23 30 108 108
24 25 89 84
25 39 156 156
26 34 129 120
27 31 118 114
28 31 118 94
29 33 125 120
30 25 95 81
31 33 126 110
32 35 135 133
33 42 154 122
34 43 165 158
35 30 113 109
36 33 127 124
37 13 52 39
38 32 121 92
39 36 136 126
40 0 0 0
41 28 108 70
42 14 46 37
43 17 54 38
44 32 124 120
45 30 115 93
46 35 128 95
47 20 80 77
48 28 97 90
49 28 104 80
50 39 59 31
51 34 125 110
52 26 82 66
53 39 149 138
54 39 149 133
55 33 122 113
56 28 118 100
57 4 12 7
58 39 144 140
59 18 67 61
60 14 52 41
61 29 108 96
62 44 166 164
63 21 80 78
64 16 60 49
65 28 107 102
66 35 127 124
67 28 107 99
68 38 146 129
69 23 84 62
70 36 141 73
71 32 123 114
72 29 111 99
73 25 98 70
74 27 105 104
75 36 135 116
76 28 107 91
77 23 85 74
78 40 155 138
79 23 88 67
80 40 155 151
81 28 104 72
82 34 132 120
83 33 127 115
84 28 108 105
85 34 129 104
86 30 116 108
87 33 122 98
88 22 85 69
89 38 147 111
90 26 99 99
91 35 87 71
92 8 28 27
93 24 90 69
94 29 109 107
95 20 78 73
96 29 111 107
97 45 158 93
98 37 141 129
99 33 122 69
100 33 124 118
101 25 93 73
102 32 124 119
103 29 112 104
104 28 108 107
105 28 99 99
106 31 117 90
107 52 199 197
108 21 78 36
109 24 91 85
110 41 158 139
111 33 126 106
112 32 122 50
113 19 71 64
114 20 75 31
115 31 115 63
116 31 119 92
117 32 124 106
118 18 72 63
119 23 91 69
120 17 45 41
121 20 78 56
122 12 39 25
123 17 68 65
124 30 119 93
125 31 117 114
126 10 39 38
127 13 50 44
128 22 88 87
129 42 155 110
130 1 0 0
131 9 36 27
132 32 123 83
133 11 32 30
134 25 99 80
135 36 136 98
136 31 117 82
137 0 0 0
138 24 88 60
139 13 39 28
140 8 25 9
141 13 52 33
142 19 75 59
143 18 71 49
144 33 124 115
145 40 151 140
146 22 71 49
147 38 145 120
148 24 87 66
149 8 27 21
150 35 131 124
151 43 162 152
152 43 165 139
153 14 54 38
154 41 159 144
155 38 147 120
156 45 170 160
157 31 119 114
158 13 49 39
159 28 104 78
160 31 120 119
161 40 150 141
162 30 112 101
163 16 59 56
164 37 136 133
165 30 107 83
166 35 130 116
167 32 115 90
168 27 107 36
169 20 75 50
170 18 71 61
171 31 120 97
172 31 116 98
173 21 79 78
174 39 150 117
175 41 156 148
176 13 51 41
177 32 118 105
178 18 71 55
179 39 144 132
180 14 47 44
181 7 28 21
182 17 68 50
183 0 0 0
184 30 110 73
185 37 147 86
186 0 0 0
187 5 15 13
188 1 4 4
189 16 64 57
190 32 111 48
191 24 85 46
192 17 68 48
193 11 40 32
194 24 80 68
195 22 88 87
196 12 48 43
197 19 76 67
198 13 51 46
199 17 67 46
200 15 59 56
201 16 61 48
202 24 76 44
203 15 60 60
204 17 68 65
205 18 71 55
206 20 76 38
207 16 62 52
208 16 61 60
209 18 67 54
210 22 88 86
211 8 30 24
212 17 64 52
213 18 68 49
214 16 64 61
215 23 91 61
216 22 88 81
217 13 52 43
218 13 49 40
219 16 62 40
220 16 61 56
221 20 76 68
222 22 88 79
223 17 66 47
224 18 71 57
225 17 68 41
226 12 48 29
227 7 25 3
228 17 68 60
229 14 41 30
230 23 90 79
231 17 66 47
232 14 54 40
233 15 59 48
234 17 60 36
235 21 77 42
236 18 68 49
237 18 72 57
238 17 67 12
239 17 64 40
240 16 63 43
241 15 59 33
242 21 84 77
243 16 64 43
244 14 56 45
245 15 54 47
246 17 67 43
247 15 58 45
248 15 59 50
249 10 40 35
250 6 22 7
251 22 83 71
252 21 81 67
253 1 2 0
254 18 72 62
255 17 61 54
256 4 15 4
257 10 32 25
258 16 62 40
259 16 58 38
260 9 36 19
261 16 59 17
262 17 68 67
263 7 21 14
264 15 55 30
265 14 54 54
266 14 55 35
267 18 72 59
268 12 41 24
269 16 61 58
270 21 67 42
271 19 76 46
272 16 64 61
273 1 3 3
274 16 63 52
275 10 40 25
276 19 69 40
277 12 48 32
278 2 8 4
279 14 52 49
280 17 66 63
281 19 76 67
282 14 43 32
283 11 39 23
284 4 14 7
285 16 61 54
286 20 71 37
287 12 44 35
288 15 60 51
289 16 64 39
totrevisions\r
1 24188
2 18273
3 14130
4 32287
5 8654
6 9245
7 33251
8 1271
9 5279
10 27101
11 16373
12 19716
13 17753
14 9028
15 18653
16 8828
17 29498
18 27563
19 18293
20 22530
21 15977
22 35082
23 16116
24 15849
25 16026
26 26569
27 24785
28 17569
29 23825
30 7869
31 14975
32 37791
33 9605
34 27295
35 2746
36 34461
37 8098
38 4787
39 24919
40 603
41 16329
42 12558
43 7784
44 28522
45 22265
46 14459
47 14526
48 22240
49 11802
50 7623
51 11912
52 7935
53 18220
54 19199
55 19918
56 21884
57 2694
58 15808
59 3597
60 5296
61 25239
62 29801
63 18450
64 7132
65 34861
66 35940
67 16688
68 24683
69 46230
70 10387
71 21436
72 30546
73 19746
74 15977
75 22583
76 17274
77 16469
78 14251
79 3007
80 16851
81 21113
82 17401
83 23958
84 23567
85 13065
86 15358
87 14587
88 12770
89 24021
90 9648
91 20537
92 7905
93 4527
94 30495
95 7117
96 17719
97 27056
98 33473
99 9758
100 21115
101 7236
102 13790
103 32902
104 25131
105 30910
106 35947
107 29848
108 6943
109 42705
110 31808
111 26675
112 8435
113 7409
114 14993
115 36867
116 33835
117 24164
118 12607
119 22609
120 5892
121 17014
122 5394
123 9178
124 6440
125 21916
126 4011
127 5818
128 18647
129 20556
130 238
131 70
132 22392
133 3913
134 12237
135 8388
136 22120
137 338
138 11727
139 3704
140 3988
141 3030
142 13520
143 1421
144 20923
145 20237
146 3219
147 3769
148 12252
149 1888
150 14497
151 28864
152 21721
153 4821
154 33644
155 15923
156 42935
157 18864
158 4977
159 7785
160 17939
161 23436
162 325
163 13539
164 34538
165 12198
166 26924
167 12716
168 8172
169 10855
170 11932
171 14300
172 25515
173 2805
174 29402
175 16440
176 11221
177 28732
178 5250
179 28608
180 8092
181 4473
182 1572
183 2065
184 14817
185 16714
186 556
187 2089
188 2658
189 10695
190 1669
191 16267
192 7768
193 7252
194 6387
195 18715
196 7936
197 8643
198 7294
199 4570
200 7185
201 10058
202 2342
203 8509
204 13275
205 6816
206 1930
207 8086
208 10737
209 8033
210 7058
211 6782
212 5401
213 6521
214 10856
215 2154
216 6117
217 5238
218 4820
219 5615
220 4272
221 8702
222 15340
223 8030
224 9526
225 1278
226 4236
227 3023
228 7196
229 3394
230 6371
231 1574
232 9620
233 6978
234 4911
235 8645
236 8987
237 5544
238 3083
239 6909
240 3189
241 6745
242 16724
243 4850
244 7025
245 6047
246 7377
247 9078
248 4605
249 3238
250 8100
251 9653
252 8914
253 786
254 6700
255 5788
256 593
257 4506
258 6382
259 5621
260 3997
261 520
262 8891
263 999
264 7067
265 4639
266 5654
267 6928
268 1514
269 9238
270 8204
271 5926
272 5785
273 4
274 5930
275 3710
276 705
277 443
278 2416
279 7747
280 5432
281 4913
282 2650
283 2370
284 775
285 5576
286 1352
287 3080
288 10205
289 6095
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) pageviews time_in_rfc
478.35220 2.32812 0.09643
logins compendium_views_info compendium_views_pr
36.50234 -26.92632 17.05420
shared_compendiums blogged_computations compendiums_reviewed
2235.70041 139.24423 1241.67399
feedback_messages_p1 feedback_messages_p120 `totrevisions\\r`
-274.84058 137.69466 1.88909
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-33042 -12610 -3314 6681 142656
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 478.35220 3409.99363 0.140 0.8885
pageviews 2.32812 7.84111 0.297 0.7668
time_in_rfc 0.09643 0.04599 2.097 0.0369 *
logins 36.50234 59.38488 0.615 0.5393
compendium_views_info -26.92632 17.62105 -1.528 0.1276
compendium_views_pr 17.05420 31.57822 0.540 0.5896
shared_compendiums 2235.70041 510.22856 4.382 1.67e-05 ***
blogged_computations 139.24423 77.18410 1.804 0.0723 .
compendiums_reviewed 1241.67399 666.26017 1.864 0.0634 .
feedback_messages_p1 -274.84058 203.15878 -1.353 0.1772
feedback_messages_p120 137.69466 100.04660 1.376 0.1698
`totrevisions\\r` 1.88909 0.24353 7.757 1.67e-13 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 20940 on 277 degrees of freedom
Multiple R-squared: 0.7419, Adjusted R-squared: 0.7317
F-statistic: 72.39 on 11 and 277 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.072785478 1.455710e-01 9.272145e-01
[2,] 0.027122822 5.424564e-02 9.728772e-01
[3,] 0.081633679 1.632674e-01 9.183663e-01
[4,] 0.051395084 1.027902e-01 9.486049e-01
[5,] 0.198758100 3.975162e-01 8.012419e-01
[6,] 0.133218362 2.664367e-01 8.667816e-01
[7,] 0.274895602 5.497912e-01 7.251044e-01
[8,] 0.297130659 5.942613e-01 7.028693e-01
[9,] 0.221911273 4.438225e-01 7.780887e-01
[10,] 0.184507056 3.690141e-01 8.154929e-01
[11,] 0.131128271 2.622565e-01 8.688717e-01
[12,] 0.096775372 1.935507e-01 9.032246e-01
[13,] 0.068933269 1.378665e-01 9.310667e-01
[14,] 0.048044139 9.608828e-02 9.519559e-01
[15,] 0.061992989 1.239860e-01 9.380070e-01
[16,] 0.050437076 1.008742e-01 9.495629e-01
[17,] 0.053353051 1.067061e-01 9.466469e-01
[18,] 0.053420660 1.068413e-01 9.465793e-01
[19,] 0.064057140 1.281143e-01 9.359429e-01
[20,] 0.056562838 1.131257e-01 9.434372e-01
[21,] 0.049947242 9.989448e-02 9.500528e-01
[22,] 0.051467095 1.029342e-01 9.485329e-01
[23,] 0.049858202 9.971640e-02 9.501418e-01
[24,] 0.038262220 7.652444e-02 9.617378e-01
[25,] 0.028657743 5.731549e-02 9.713423e-01
[26,] 0.026880973 5.376195e-02 9.731190e-01
[27,] 0.021640500 4.328100e-02 9.783595e-01
[28,] 0.017893415 3.578683e-02 9.821066e-01
[29,] 0.012939365 2.587873e-02 9.870606e-01
[30,] 0.009577722 1.915544e-02 9.904223e-01
[31,] 0.026702278 5.340456e-02 9.732977e-01
[32,] 0.023315468 4.663094e-02 9.766845e-01
[33,] 0.017429371 3.485874e-02 9.825706e-01
[34,] 0.012384614 2.476923e-02 9.876154e-01
[35,] 0.009361677 1.872335e-02 9.906383e-01
[36,] 0.020967566 4.193513e-02 9.790324e-01
[37,] 0.018354904 3.670981e-02 9.816451e-01
[38,] 0.016106900 3.221380e-02 9.838931e-01
[39,] 0.012466633 2.493327e-02 9.875334e-01
[40,] 0.037135071 7.427014e-02 9.628649e-01
[41,] 0.031374944 6.274989e-02 9.686251e-01
[42,] 0.023955942 4.791188e-02 9.760441e-01
[43,] 0.017917781 3.583556e-02 9.820822e-01
[44,] 0.023983248 4.796650e-02 9.760168e-01
[45,] 0.018238377 3.647675e-02 9.817616e-01
[46,] 0.013537078 2.707416e-02 9.864629e-01
[47,] 0.010049454 2.009891e-02 9.899505e-01
[48,] 0.010254757 2.050951e-02 9.897452e-01
[49,] 0.013918935 2.783787e-02 9.860811e-01
[50,] 0.011339879 2.267976e-02 9.886601e-01
[51,] 0.008626787 1.725357e-02 9.913732e-01
[52,] 0.006931788 1.386358e-02 9.930682e-01
[53,] 0.009097397 1.819479e-02 9.909026e-01
[54,] 0.007539744 1.507949e-02 9.924603e-01
[55,] 0.005629742 1.125948e-02 9.943703e-01
[56,] 0.005245794 1.049159e-02 9.947542e-01
[57,] 0.004193372 8.386744e-03 9.958066e-01
[58,] 0.003050171 6.100343e-03 9.969498e-01
[59,] 0.002725378 5.450756e-03 9.972746e-01
[60,] 0.004735715 9.471430e-03 9.952643e-01
[61,] 0.007761335 1.552267e-02 9.922387e-01
[62,] 0.007673986 1.534797e-02 9.923260e-01
[63,] 0.005717859 1.143572e-02 9.942821e-01
[64,] 0.007982824 1.596565e-02 9.920172e-01
[65,] 0.006776223 1.355245e-02 9.932238e-01
[66,] 0.006719144 1.343829e-02 9.932809e-01
[67,] 0.005593200 1.118640e-02 9.944068e-01
[68,] 0.004266800 8.533601e-03 9.957332e-01
[69,] 0.003426574 6.853148e-03 9.965734e-01
[70,] 0.002665397 5.330795e-03 9.973346e-01
[71,] 0.968592300 6.281540e-02 3.140770e-02
[72,] 0.964194714 7.161057e-02 3.580529e-02
[73,] 0.956636445 8.672711e-02 4.336356e-02
[74,] 0.948375686 1.032486e-01 5.162431e-02
[75,] 0.938134966 1.237301e-01 6.186503e-02
[76,] 0.936528429 1.269431e-01 6.347157e-02
[77,] 0.936513417 1.269732e-01 6.348658e-02
[78,] 0.944825096 1.103498e-01 5.517490e-02
[79,] 0.936192210 1.276156e-01 6.380779e-02
[80,] 0.928700349 1.425993e-01 7.129965e-02
[81,] 0.923001207 1.539976e-01 7.699879e-02
[82,] 0.918008949 1.639821e-01 8.199105e-02
[83,] 0.933320075 1.333598e-01 6.667992e-02
[84,] 0.931281948 1.374361e-01 6.871805e-02
[85,] 0.926790986 1.464180e-01 7.320901e-02
[86,] 0.914037573 1.719249e-01 8.596243e-02
[87,] 0.899577672 2.008447e-01 1.004223e-01
[88,] 0.892527669 2.149447e-01 1.074723e-01
[89,] 0.877000694 2.459986e-01 1.229993e-01
[90,] 0.900018651 1.999627e-01 9.998135e-02
[91,] 0.913972008 1.720560e-01 8.602799e-02
[92,] 0.951229593 9.754081e-02 4.877041e-02
[93,] 0.962667133 7.466573e-02 3.733287e-02
[94,] 0.956003769 8.799246e-02 4.399623e-02
[95,] 0.950924750 9.815050e-02 4.907525e-02
[96,] 0.951767257 9.646549e-02 4.823274e-02
[97,] 0.961280461 7.743908e-02 3.871954e-02
[98,] 0.956087286 8.782543e-02 4.391271e-02
[99,] 0.948018520 1.039630e-01 5.198148e-02
[100,] 0.951400083 9.719983e-02 4.859992e-02
[101,] 0.950102097 9.979581e-02 4.989790e-02
[102,] 0.945902929 1.081941e-01 5.409707e-02
[103,] 0.938283124 1.234338e-01 6.171688e-02
[104,] 0.930243372 1.395133e-01 6.975663e-02
[105,] 0.999998788 2.423879e-06 1.211940e-06
[106,] 0.999998196 3.607790e-06 1.803895e-06
[107,] 0.999997797 4.405931e-06 2.202965e-06
[108,] 0.999997320 5.359998e-06 2.679999e-06
[109,] 0.999996506 6.988249e-06 3.494124e-06
[110,] 1.000000000 5.288734e-10 2.644367e-10
[111,] 1.000000000 4.030383e-10 2.015191e-10
[112,] 1.000000000 5.615088e-10 2.807544e-10
[113,] 1.000000000 8.845432e-10 4.422716e-10
[114,] 0.999999999 1.387364e-09 6.936821e-10
[115,] 0.999999999 2.363156e-09 1.181578e-09
[116,] 0.999999998 3.917723e-09 1.958861e-09
[117,] 0.999999997 6.189351e-09 3.094676e-09
[118,] 0.999999996 7.169677e-09 3.584838e-09
[119,] 0.999999995 9.362253e-09 4.681127e-09
[120,] 0.999999994 1.262354e-08 6.311771e-09
[121,] 0.999999992 1.614392e-08 8.071960e-09
[122,] 0.999999987 2.552741e-08 1.276371e-08
[123,] 0.999999979 4.195859e-08 2.097929e-08
[124,] 0.999999988 2.395174e-08 1.197587e-08
[125,] 0.999999989 2.126805e-08 1.063403e-08
[126,] 0.999999983 3.486603e-08 1.743302e-08
[127,] 0.999999973 5.329261e-08 2.664630e-08
[128,] 0.999999969 6.156114e-08 3.078057e-08
[129,] 0.999999950 9.995044e-08 4.997522e-08
[130,] 0.999999939 1.217827e-07 6.089135e-08
[131,] 0.999999909 1.827996e-07 9.139978e-08
[132,] 0.999999855 2.895471e-07 1.447735e-07
[133,] 0.999999892 2.163055e-07 1.081528e-07
[134,] 0.999999979 4.223201e-08 2.111601e-08
[135,] 0.999999966 6.749587e-08 3.374794e-08
[136,] 0.999999953 9.318810e-08 4.659405e-08
[137,] 0.999999954 9.142758e-08 4.571379e-08
[138,] 0.999999926 1.478831e-07 7.394154e-08
[139,] 0.999999881 2.383914e-07 1.191957e-07
[140,] 0.999999924 1.514362e-07 7.571810e-08
[141,] 0.999999914 1.721387e-07 8.606936e-08
[142,] 0.999999999 1.524304e-09 7.621521e-10
[143,] 0.999999999 2.647189e-09 1.323594e-09
[144,] 0.999999998 4.102733e-09 2.051367e-09
[145,] 0.999999997 5.699237e-09 2.849619e-09
[146,] 0.999999996 8.405308e-09 4.202654e-09
[147,] 0.999999993 1.383299e-08 6.916493e-09
[148,] 1.000000000 9.088784e-11 4.544392e-11
[149,] 1.000000000 1.599688e-10 7.998441e-11
[150,] 1.000000000 2.232062e-10 1.116031e-10
[151,] 1.000000000 1.936360e-10 9.681798e-11
[152,] 1.000000000 1.714083e-10 8.570415e-11
[153,] 1.000000000 2.160208e-10 1.080104e-10
[154,] 1.000000000 1.995593e-11 9.977965e-12
[155,] 1.000000000 2.873622e-11 1.436811e-11
[156,] 1.000000000 2.467557e-12 1.233778e-12
[157,] 1.000000000 4.199165e-12 2.099582e-12
[158,] 1.000000000 7.851207e-12 3.925604e-12
[159,] 1.000000000 1.209926e-11 6.049632e-12
[160,] 1.000000000 7.126897e-12 3.563448e-12
[161,] 1.000000000 2.962964e-12 1.481482e-12
[162,] 1.000000000 3.676936e-12 1.838468e-12
[163,] 1.000000000 1.563482e-12 7.817409e-13
[164,] 1.000000000 2.956186e-12 1.478093e-12
[165,] 1.000000000 1.051092e-12 5.255458e-13
[166,] 1.000000000 1.450586e-12 7.252932e-13
[167,] 1.000000000 2.761008e-12 1.380504e-12
[168,] 1.000000000 1.460818e-12 7.304090e-13
[169,] 1.000000000 2.958271e-12 1.479135e-12
[170,] 1.000000000 5.937704e-12 2.968852e-12
[171,] 1.000000000 6.316221e-12 3.158110e-12
[172,] 1.000000000 1.267522e-11 6.337611e-12
[173,] 1.000000000 1.877206e-11 9.386029e-12
[174,] 1.000000000 3.555958e-11 1.777979e-11
[175,] 1.000000000 4.724507e-11 2.362253e-11
[176,] 1.000000000 8.960430e-11 4.480215e-11
[177,] 1.000000000 1.418515e-10 7.092573e-11
[178,] 1.000000000 2.420891e-10 1.210446e-10
[179,] 1.000000000 3.702810e-10 1.851405e-10
[180,] 1.000000000 4.451384e-10 2.225692e-10
[181,] 1.000000000 7.897290e-10 3.948645e-10
[182,] 0.999999999 1.409839e-09 7.049193e-10
[183,] 0.999999999 1.986819e-09 9.934095e-10
[184,] 1.000000000 8.253265e-11 4.126633e-11
[185,] 1.000000000 5.882540e-11 2.941270e-11
[186,] 1.000000000 5.296807e-11 2.648404e-11
[187,] 1.000000000 8.309498e-11 4.154749e-11
[188,] 1.000000000 4.298269e-11 2.149134e-11
[189,] 1.000000000 8.780055e-11 4.390027e-11
[190,] 1.000000000 1.536087e-10 7.680433e-11
[191,] 1.000000000 2.799029e-10 1.399514e-10
[192,] 1.000000000 5.641772e-10 2.820886e-10
[193,] 1.000000000 8.034508e-10 4.017254e-10
[194,] 0.999999999 1.493858e-09 7.469290e-10
[195,] 0.999999999 1.993789e-09 9.968946e-10
[196,] 0.999999999 1.604200e-09 8.021002e-10
[197,] 0.999999999 1.890236e-09 9.451180e-10
[198,] 0.999999998 3.302347e-09 1.651174e-09
[199,] 0.999999997 6.553624e-09 3.276812e-09
[200,] 0.999999995 1.087775e-08 5.438876e-09
[201,] 0.999999991 1.773268e-08 8.866338e-09
[202,] 0.999999990 1.927885e-08 9.639423e-09
[203,] 0.999999981 3.798017e-08 1.899008e-08
[204,] 0.999999989 2.211132e-08 1.105566e-08
[205,] 0.999999980 3.989768e-08 1.994884e-08
[206,] 0.999999969 6.209277e-08 3.104638e-08
[207,] 0.999999962 7.516272e-08 3.758136e-08
[208,] 0.999999966 6.720138e-08 3.360069e-08
[209,] 0.999999946 1.089803e-07 5.449016e-08
[210,] 0.999999901 1.988913e-07 9.944564e-08
[211,] 0.999999824 3.517803e-07 1.758902e-07
[212,] 0.999999738 5.240938e-07 2.620469e-07
[213,] 0.999999853 2.945846e-07 1.472923e-07
[214,] 0.999999753 4.932172e-07 2.466086e-07
[215,] 0.999999907 1.854361e-07 9.271803e-08
[216,] 0.999999971 5.857435e-08 2.928717e-08
[217,] 0.999999951 9.882226e-08 4.941113e-08
[218,] 0.999999899 2.026857e-07 1.013429e-07
[219,] 0.999999909 1.818344e-07 9.091718e-08
[220,] 0.999999814 3.722728e-07 1.861364e-07
[221,] 0.999999642 7.151508e-07 3.575754e-07
[222,] 0.999999297 1.406276e-06 7.031379e-07
[223,] 0.999999026 1.948608e-06 9.743040e-07
[224,] 0.999998187 3.626221e-06 1.813110e-06
[225,] 0.999996551 6.897638e-06 3.448819e-06
[226,] 0.999993200 1.359967e-05 6.799834e-06
[227,] 0.999989286 2.142785e-05 1.071393e-05
[228,] 0.999987450 2.509959e-05 1.254980e-05
[229,] 0.999979186 4.162792e-05 2.081396e-05
[230,] 0.999965221 6.955791e-05 3.477895e-05
[231,] 0.999935074 1.298512e-04 6.492562e-05
[232,] 0.999924051 1.518975e-04 7.594877e-05
[233,] 0.999878807 2.423852e-04 1.211926e-04
[234,] 0.999774814 4.503716e-04 2.251858e-04
[235,] 0.999747507 5.049861e-04 2.524930e-04
[236,] 0.999922853 1.542949e-04 7.714747e-05
[237,] 0.999849588 3.008231e-04 1.504115e-04
[238,] 0.999804679 3.906412e-04 1.953206e-04
[239,] 0.999624687 7.506259e-04 3.753129e-04
[240,] 0.999298199 1.403602e-03 7.018012e-04
[241,] 0.999541282 9.174361e-04 4.587180e-04
[242,] 0.999248792 1.502415e-03 7.512076e-04
[243,] 0.998637981 2.724037e-03 1.362019e-03
[244,] 0.997511413 4.977174e-03 2.488587e-03
[245,] 0.997014387 5.971225e-03 2.985613e-03
[246,] 0.994514540 1.097092e-02 5.485460e-03
[247,] 0.990576214 1.884757e-02 9.423786e-03
[248,] 0.990890864 1.821827e-02 9.109136e-03
[249,] 0.996734682 6.530635e-03 3.265318e-03
[250,] 0.993501139 1.299772e-02 6.498861e-03
[251,] 0.987557261 2.488548e-02 1.244274e-02
[252,] 0.976562000 4.687600e-02 2.343800e-02
[253,] 0.958400249 8.319950e-02 4.159975e-02
[254,] 0.955956877 8.808625e-02 4.404312e-02
[255,] 0.950589819 9.882036e-02 4.941018e-02
[256,] 0.930115324 1.397694e-01 6.988468e-02
[257,] 0.904578065 1.908439e-01 9.542194e-02
[258,] 0.851579833 2.968403e-01 1.484202e-01
[259,] 0.744692984 5.106140e-01 2.553070e-01
[260,] 0.720314971 5.593701e-01 2.796850e-01
> postscript(file="/var/wessaorg/rcomp/tmp/1iyzf1355259025.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/23r8r1355259025.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/3fnxm1355259025.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/4l98h1355259025.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/59cbk1355259025.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> qqnorm(mysum$resid, main='Residual Normal Q-Q Plot')
> qqline(mysum$resid)
> grid()
> dev.off()
null device
1
> (myerror <- as.ts(mysum$resid))
Time Series:
Start = 1
End = 289
Frequency = 1
1 2 3 4 5 6
13416.25568 5102.47318 -7449.36293 -1031.73135 -12023.45241 25833.05897
7 8 9 10 11 12
14848.53877 10288.33407 -9281.58093 26600.16800 -4846.92614 -21703.01642
13 14 15 16 17 18
-31050.59067 -8743.36636 -14915.46894 -18798.46310 -21807.54113 24493.46257
19 20 21 22 23 24
22453.82838 24575.03113 15856.73516 -17069.67360 -15778.86573 -10281.63742
25 26 27 28 29 30
5415.04800 -24324.94949 -1444.06853 1516.07614 -25713.83902 -11710.07795
31 32 33 34 35 36
2216.06244 -8626.79259 4492.54650 -25651.20409 15792.76531 -17806.80365
37 38 39 40 41 42
-6806.53227 11117.45739 -21730.00274 -72.27331 -12684.33025 -16702.04064
43 44 45 46 47 48
-14451.07196 -23764.81138 -29203.00115 -15941.86442 -13786.70180 -9509.77244
49 50 51 52 53 54
-10921.45782 18325.94966 22183.55148 -11882.81805 -1633.00009 34648.31272
55 56 57 58 59 60
-9172.86105 -6658.07555 -2198.50916 1318.13399 -3275.89066 -2161.24432
61 62 63 64 65 66
-6528.12040 13269.93435 32744.95029 -11350.38299 -3184.71888 -18558.20289
67 68 69 70 71 72
-24938.33995 10859.11826 5569.61179 -17799.00356 -9857.34533 -3434.25295
73 74 75 76 77 78
-18935.82560 27236.58067 -31303.48520 14407.01434 -661.82051 -28169.16422
79 80 81 82 83 84
-14256.40016 -21483.33819 7491.08683 -8272.57258 -12045.63592 1229.23678
85 86 87 88 89 90
142655.56121 -12610.30204 -2778.93266 -6700.86112 3269.37273 -17080.51837
91 92 93 94 95 96
-15943.30237 28925.19606 -11902.77347 -13661.36282 -15543.78652 12260.22554
97 98 99 100 101 102
-33042.15211 16875.45721 -16589.45048 6680.71027 -3465.79236 -14848.37766
103 104 105 106 107 108
1590.69189 30714.93498 30275.10341 43660.88234 -29174.94891 -8331.33909
109 110 111 112 113 114
16821.06988 -18618.62455 30563.88097 -8586.65882 4573.76565 -22220.86137
115 116 117 118 119 120
-18185.30615 -15035.98504 -7217.89005 -11880.16248 125781.63928 -2022.21295
121 122 123 124 125 126
-16904.24123 -9028.42275 -15308.04079 93805.89921 25265.57747 -11486.01195
127 128 129 130 131 132
8291.30755 -5046.19090 3339.09771 -2166.73755 -4898.49524 -15372.68665
133 134 135 136 137 138
-7216.28487 -14257.34057 -14378.70691 -840.24802 -437.06821 33039.56820
139 140 141 142 143 144
25797.11192 -462.28298 -5485.20066 13407.22697 -1197.24616 -14060.56644
145 146 147 148 149 150
10083.66237 -1412.84638 29795.46967 45272.50141 -2607.91980 12693.33942
151 152 153 154 155 156
-18138.07164 56.00017 4753.77212 -19143.03345 16654.85303 57100.36150
157 158 159 160 161 162
-3133.80949 -2245.96031 17382.93061 -1292.85801 4980.54410 59936.47814
163 164 165 166 167 168
435.44210 20099.48923 -10431.24259 -26557.16695 -14169.04186 38351.49242
169 170 171 172 173 174
-15939.85047 28327.13088 14773.35202 5079.18945 -1155.75109 -28317.82344
175 176 177 178 179 180
29932.53739 -7062.77012 37781.53758 -7595.97683 14598.57788 -19310.05542
181 182 183 184 185 186
7107.57548 22785.37996 1112.05219 -1621.22598 -26853.88421 -95.14192
187 188 189 190 191 192
-3271.08092 -2662.81701 -9947.06617 2945.29422 10009.59971 -722.77451
193 194 195 196 197 198
-16138.92387 -12829.75930 1487.18144 -10396.54191 -7516.93014 36836.33241
199 200 201 202 203 204
29281.95284 11436.72826 -15326.71249 -17187.28418 -5387.42725 -2501.48426
205 206 207 208 209 210
-9850.61177 1697.91263 -13727.90442 -3768.53773 -9722.75475 -19009.90735
211 212 213 214 215 216
-16273.43012 -3653.36288 -5072.30643 -13251.42331 10707.10774 -21635.70426
217 218 219 220 221 222
318.31532 3254.17359 -14572.96969 -7271.33387 -5571.32333 -3618.17458
223 224 225 226 227 228
-12976.00474 -3313.80549 -9792.58169 -13987.41142 18796.45530 5134.65422
229 230 231 232 233 234
-13454.98243 31974.01184 -12591.97555 -7908.15174 24109.53309 -2873.91257
235 236 237 238 239 240
4280.40653 -10256.82316 -12.20818 -4326.15657 -9357.12474 -2094.43091
241 242 243 244 245 246
9534.39607 -15608.97318 -10296.10549 -3544.40042 1463.98009 10115.83482
247 248 249 250 251 252
6149.31175 3211.07046 16237.89459 16155.59078 -60.36738 2316.04826
253 254 255 256 257 258
-1806.38407 -5330.17573 -8207.65406 -1840.10101 4995.98976 -1941.02370
259 260 261 262 263 264
-15229.06747 -4547.10217 -2498.71024 -17410.87549 -6388.05652 -10765.46144
265 266 267 268 269 270
12140.41666 -3123.89881 -11083.52325 8193.75376 -18417.12595 11511.68675
271 272 273 274 275 276
-5588.40214 91.43569 587.15568 -7139.23629 -151.72776 -10468.39404
277 278 279 280 281 282
14454.47639 -5065.98177 -15581.63864 -8982.68648 -3721.46545 59676.89022
283 284 285 286 287 288
9182.13330 -3277.76741 -5662.96636 17395.76129 -12618.52747 -7069.25255
289
-3945.50545
> postscript(file="/var/wessaorg/rcomp/tmp/6qzko1355259025.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 13416.25568 NA
1 5102.47318 13416.25568
2 -7449.36293 5102.47318
3 -1031.73135 -7449.36293
4 -12023.45241 -1031.73135
5 25833.05897 -12023.45241
6 14848.53877 25833.05897
7 10288.33407 14848.53877
8 -9281.58093 10288.33407
9 26600.16800 -9281.58093
10 -4846.92614 26600.16800
11 -21703.01642 -4846.92614
12 -31050.59067 -21703.01642
13 -8743.36636 -31050.59067
14 -14915.46894 -8743.36636
15 -18798.46310 -14915.46894
16 -21807.54113 -18798.46310
17 24493.46257 -21807.54113
18 22453.82838 24493.46257
19 24575.03113 22453.82838
20 15856.73516 24575.03113
21 -17069.67360 15856.73516
22 -15778.86573 -17069.67360
23 -10281.63742 -15778.86573
24 5415.04800 -10281.63742
25 -24324.94949 5415.04800
26 -1444.06853 -24324.94949
27 1516.07614 -1444.06853
28 -25713.83902 1516.07614
29 -11710.07795 -25713.83902
30 2216.06244 -11710.07795
31 -8626.79259 2216.06244
32 4492.54650 -8626.79259
33 -25651.20409 4492.54650
34 15792.76531 -25651.20409
35 -17806.80365 15792.76531
36 -6806.53227 -17806.80365
37 11117.45739 -6806.53227
38 -21730.00274 11117.45739
39 -72.27331 -21730.00274
40 -12684.33025 -72.27331
41 -16702.04064 -12684.33025
42 -14451.07196 -16702.04064
43 -23764.81138 -14451.07196
44 -29203.00115 -23764.81138
45 -15941.86442 -29203.00115
46 -13786.70180 -15941.86442
47 -9509.77244 -13786.70180
48 -10921.45782 -9509.77244
49 18325.94966 -10921.45782
50 22183.55148 18325.94966
51 -11882.81805 22183.55148
52 -1633.00009 -11882.81805
53 34648.31272 -1633.00009
54 -9172.86105 34648.31272
55 -6658.07555 -9172.86105
56 -2198.50916 -6658.07555
57 1318.13399 -2198.50916
58 -3275.89066 1318.13399
59 -2161.24432 -3275.89066
60 -6528.12040 -2161.24432
61 13269.93435 -6528.12040
62 32744.95029 13269.93435
63 -11350.38299 32744.95029
64 -3184.71888 -11350.38299
65 -18558.20289 -3184.71888
66 -24938.33995 -18558.20289
67 10859.11826 -24938.33995
68 5569.61179 10859.11826
69 -17799.00356 5569.61179
70 -9857.34533 -17799.00356
71 -3434.25295 -9857.34533
72 -18935.82560 -3434.25295
73 27236.58067 -18935.82560
74 -31303.48520 27236.58067
75 14407.01434 -31303.48520
76 -661.82051 14407.01434
77 -28169.16422 -661.82051
78 -14256.40016 -28169.16422
79 -21483.33819 -14256.40016
80 7491.08683 -21483.33819
81 -8272.57258 7491.08683
82 -12045.63592 -8272.57258
83 1229.23678 -12045.63592
84 142655.56121 1229.23678
85 -12610.30204 142655.56121
86 -2778.93266 -12610.30204
87 -6700.86112 -2778.93266
88 3269.37273 -6700.86112
89 -17080.51837 3269.37273
90 -15943.30237 -17080.51837
91 28925.19606 -15943.30237
92 -11902.77347 28925.19606
93 -13661.36282 -11902.77347
94 -15543.78652 -13661.36282
95 12260.22554 -15543.78652
96 -33042.15211 12260.22554
97 16875.45721 -33042.15211
98 -16589.45048 16875.45721
99 6680.71027 -16589.45048
100 -3465.79236 6680.71027
101 -14848.37766 -3465.79236
102 1590.69189 -14848.37766
103 30714.93498 1590.69189
104 30275.10341 30714.93498
105 43660.88234 30275.10341
106 -29174.94891 43660.88234
107 -8331.33909 -29174.94891
108 16821.06988 -8331.33909
109 -18618.62455 16821.06988
110 30563.88097 -18618.62455
111 -8586.65882 30563.88097
112 4573.76565 -8586.65882
113 -22220.86137 4573.76565
114 -18185.30615 -22220.86137
115 -15035.98504 -18185.30615
116 -7217.89005 -15035.98504
117 -11880.16248 -7217.89005
118 125781.63928 -11880.16248
119 -2022.21295 125781.63928
120 -16904.24123 -2022.21295
121 -9028.42275 -16904.24123
122 -15308.04079 -9028.42275
123 93805.89921 -15308.04079
124 25265.57747 93805.89921
125 -11486.01195 25265.57747
126 8291.30755 -11486.01195
127 -5046.19090 8291.30755
128 3339.09771 -5046.19090
129 -2166.73755 3339.09771
130 -4898.49524 -2166.73755
131 -15372.68665 -4898.49524
132 -7216.28487 -15372.68665
133 -14257.34057 -7216.28487
134 -14378.70691 -14257.34057
135 -840.24802 -14378.70691
136 -437.06821 -840.24802
137 33039.56820 -437.06821
138 25797.11192 33039.56820
139 -462.28298 25797.11192
140 -5485.20066 -462.28298
141 13407.22697 -5485.20066
142 -1197.24616 13407.22697
143 -14060.56644 -1197.24616
144 10083.66237 -14060.56644
145 -1412.84638 10083.66237
146 29795.46967 -1412.84638
147 45272.50141 29795.46967
148 -2607.91980 45272.50141
149 12693.33942 -2607.91980
150 -18138.07164 12693.33942
151 56.00017 -18138.07164
152 4753.77212 56.00017
153 -19143.03345 4753.77212
154 16654.85303 -19143.03345
155 57100.36150 16654.85303
156 -3133.80949 57100.36150
157 -2245.96031 -3133.80949
158 17382.93061 -2245.96031
159 -1292.85801 17382.93061
160 4980.54410 -1292.85801
161 59936.47814 4980.54410
162 435.44210 59936.47814
163 20099.48923 435.44210
164 -10431.24259 20099.48923
165 -26557.16695 -10431.24259
166 -14169.04186 -26557.16695
167 38351.49242 -14169.04186
168 -15939.85047 38351.49242
169 28327.13088 -15939.85047
170 14773.35202 28327.13088
171 5079.18945 14773.35202
172 -1155.75109 5079.18945
173 -28317.82344 -1155.75109
174 29932.53739 -28317.82344
175 -7062.77012 29932.53739
176 37781.53758 -7062.77012
177 -7595.97683 37781.53758
178 14598.57788 -7595.97683
179 -19310.05542 14598.57788
180 7107.57548 -19310.05542
181 22785.37996 7107.57548
182 1112.05219 22785.37996
183 -1621.22598 1112.05219
184 -26853.88421 -1621.22598
185 -95.14192 -26853.88421
186 -3271.08092 -95.14192
187 -2662.81701 -3271.08092
188 -9947.06617 -2662.81701
189 2945.29422 -9947.06617
190 10009.59971 2945.29422
191 -722.77451 10009.59971
192 -16138.92387 -722.77451
193 -12829.75930 -16138.92387
194 1487.18144 -12829.75930
195 -10396.54191 1487.18144
196 -7516.93014 -10396.54191
197 36836.33241 -7516.93014
198 29281.95284 36836.33241
199 11436.72826 29281.95284
200 -15326.71249 11436.72826
201 -17187.28418 -15326.71249
202 -5387.42725 -17187.28418
203 -2501.48426 -5387.42725
204 -9850.61177 -2501.48426
205 1697.91263 -9850.61177
206 -13727.90442 1697.91263
207 -3768.53773 -13727.90442
208 -9722.75475 -3768.53773
209 -19009.90735 -9722.75475
210 -16273.43012 -19009.90735
211 -3653.36288 -16273.43012
212 -5072.30643 -3653.36288
213 -13251.42331 -5072.30643
214 10707.10774 -13251.42331
215 -21635.70426 10707.10774
216 318.31532 -21635.70426
217 3254.17359 318.31532
218 -14572.96969 3254.17359
219 -7271.33387 -14572.96969
220 -5571.32333 -7271.33387
221 -3618.17458 -5571.32333
222 -12976.00474 -3618.17458
223 -3313.80549 -12976.00474
224 -9792.58169 -3313.80549
225 -13987.41142 -9792.58169
226 18796.45530 -13987.41142
227 5134.65422 18796.45530
228 -13454.98243 5134.65422
229 31974.01184 -13454.98243
230 -12591.97555 31974.01184
231 -7908.15174 -12591.97555
232 24109.53309 -7908.15174
233 -2873.91257 24109.53309
234 4280.40653 -2873.91257
235 -10256.82316 4280.40653
236 -12.20818 -10256.82316
237 -4326.15657 -12.20818
238 -9357.12474 -4326.15657
239 -2094.43091 -9357.12474
240 9534.39607 -2094.43091
241 -15608.97318 9534.39607
242 -10296.10549 -15608.97318
243 -3544.40042 -10296.10549
244 1463.98009 -3544.40042
245 10115.83482 1463.98009
246 6149.31175 10115.83482
247 3211.07046 6149.31175
248 16237.89459 3211.07046
249 16155.59078 16237.89459
250 -60.36738 16155.59078
251 2316.04826 -60.36738
252 -1806.38407 2316.04826
253 -5330.17573 -1806.38407
254 -8207.65406 -5330.17573
255 -1840.10101 -8207.65406
256 4995.98976 -1840.10101
257 -1941.02370 4995.98976
258 -15229.06747 -1941.02370
259 -4547.10217 -15229.06747
260 -2498.71024 -4547.10217
261 -17410.87549 -2498.71024
262 -6388.05652 -17410.87549
263 -10765.46144 -6388.05652
264 12140.41666 -10765.46144
265 -3123.89881 12140.41666
266 -11083.52325 -3123.89881
267 8193.75376 -11083.52325
268 -18417.12595 8193.75376
269 11511.68675 -18417.12595
270 -5588.40214 11511.68675
271 91.43569 -5588.40214
272 587.15568 91.43569
273 -7139.23629 587.15568
274 -151.72776 -7139.23629
275 -10468.39404 -151.72776
276 14454.47639 -10468.39404
277 -5065.98177 14454.47639
278 -15581.63864 -5065.98177
279 -8982.68648 -15581.63864
280 -3721.46545 -8982.68648
281 59676.89022 -3721.46545
282 9182.13330 59676.89022
283 -3277.76741 9182.13330
284 -5662.96636 -3277.76741
285 17395.76129 -5662.96636
286 -12618.52747 17395.76129
287 -7069.25255 -12618.52747
288 -3945.50545 -7069.25255
289 NA -3945.50545
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] 5102.47318 13416.25568
[2,] -7449.36293 5102.47318
[3,] -1031.73135 -7449.36293
[4,] -12023.45241 -1031.73135
[5,] 25833.05897 -12023.45241
[6,] 14848.53877 25833.05897
[7,] 10288.33407 14848.53877
[8,] -9281.58093 10288.33407
[9,] 26600.16800 -9281.58093
[10,] -4846.92614 26600.16800
[11,] -21703.01642 -4846.92614
[12,] -31050.59067 -21703.01642
[13,] -8743.36636 -31050.59067
[14,] -14915.46894 -8743.36636
[15,] -18798.46310 -14915.46894
[16,] -21807.54113 -18798.46310
[17,] 24493.46257 -21807.54113
[18,] 22453.82838 24493.46257
[19,] 24575.03113 22453.82838
[20,] 15856.73516 24575.03113
[21,] -17069.67360 15856.73516
[22,] -15778.86573 -17069.67360
[23,] -10281.63742 -15778.86573
[24,] 5415.04800 -10281.63742
[25,] -24324.94949 5415.04800
[26,] -1444.06853 -24324.94949
[27,] 1516.07614 -1444.06853
[28,] -25713.83902 1516.07614
[29,] -11710.07795 -25713.83902
[30,] 2216.06244 -11710.07795
[31,] -8626.79259 2216.06244
[32,] 4492.54650 -8626.79259
[33,] -25651.20409 4492.54650
[34,] 15792.76531 -25651.20409
[35,] -17806.80365 15792.76531
[36,] -6806.53227 -17806.80365
[37,] 11117.45739 -6806.53227
[38,] -21730.00274 11117.45739
[39,] -72.27331 -21730.00274
[40,] -12684.33025 -72.27331
[41,] -16702.04064 -12684.33025
[42,] -14451.07196 -16702.04064
[43,] -23764.81138 -14451.07196
[44,] -29203.00115 -23764.81138
[45,] -15941.86442 -29203.00115
[46,] -13786.70180 -15941.86442
[47,] -9509.77244 -13786.70180
[48,] -10921.45782 -9509.77244
[49,] 18325.94966 -10921.45782
[50,] 22183.55148 18325.94966
[51,] -11882.81805 22183.55148
[52,] -1633.00009 -11882.81805
[53,] 34648.31272 -1633.00009
[54,] -9172.86105 34648.31272
[55,] -6658.07555 -9172.86105
[56,] -2198.50916 -6658.07555
[57,] 1318.13399 -2198.50916
[58,] -3275.89066 1318.13399
[59,] -2161.24432 -3275.89066
[60,] -6528.12040 -2161.24432
[61,] 13269.93435 -6528.12040
[62,] 32744.95029 13269.93435
[63,] -11350.38299 32744.95029
[64,] -3184.71888 -11350.38299
[65,] -18558.20289 -3184.71888
[66,] -24938.33995 -18558.20289
[67,] 10859.11826 -24938.33995
[68,] 5569.61179 10859.11826
[69,] -17799.00356 5569.61179
[70,] -9857.34533 -17799.00356
[71,] -3434.25295 -9857.34533
[72,] -18935.82560 -3434.25295
[73,] 27236.58067 -18935.82560
[74,] -31303.48520 27236.58067
[75,] 14407.01434 -31303.48520
[76,] -661.82051 14407.01434
[77,] -28169.16422 -661.82051
[78,] -14256.40016 -28169.16422
[79,] -21483.33819 -14256.40016
[80,] 7491.08683 -21483.33819
[81,] -8272.57258 7491.08683
[82,] -12045.63592 -8272.57258
[83,] 1229.23678 -12045.63592
[84,] 142655.56121 1229.23678
[85,] -12610.30204 142655.56121
[86,] -2778.93266 -12610.30204
[87,] -6700.86112 -2778.93266
[88,] 3269.37273 -6700.86112
[89,] -17080.51837 3269.37273
[90,] -15943.30237 -17080.51837
[91,] 28925.19606 -15943.30237
[92,] -11902.77347 28925.19606
[93,] -13661.36282 -11902.77347
[94,] -15543.78652 -13661.36282
[95,] 12260.22554 -15543.78652
[96,] -33042.15211 12260.22554
[97,] 16875.45721 -33042.15211
[98,] -16589.45048 16875.45721
[99,] 6680.71027 -16589.45048
[100,] -3465.79236 6680.71027
[101,] -14848.37766 -3465.79236
[102,] 1590.69189 -14848.37766
[103,] 30714.93498 1590.69189
[104,] 30275.10341 30714.93498
[105,] 43660.88234 30275.10341
[106,] -29174.94891 43660.88234
[107,] -8331.33909 -29174.94891
[108,] 16821.06988 -8331.33909
[109,] -18618.62455 16821.06988
[110,] 30563.88097 -18618.62455
[111,] -8586.65882 30563.88097
[112,] 4573.76565 -8586.65882
[113,] -22220.86137 4573.76565
[114,] -18185.30615 -22220.86137
[115,] -15035.98504 -18185.30615
[116,] -7217.89005 -15035.98504
[117,] -11880.16248 -7217.89005
[118,] 125781.63928 -11880.16248
[119,] -2022.21295 125781.63928
[120,] -16904.24123 -2022.21295
[121,] -9028.42275 -16904.24123
[122,] -15308.04079 -9028.42275
[123,] 93805.89921 -15308.04079
[124,] 25265.57747 93805.89921
[125,] -11486.01195 25265.57747
[126,] 8291.30755 -11486.01195
[127,] -5046.19090 8291.30755
[128,] 3339.09771 -5046.19090
[129,] -2166.73755 3339.09771
[130,] -4898.49524 -2166.73755
[131,] -15372.68665 -4898.49524
[132,] -7216.28487 -15372.68665
[133,] -14257.34057 -7216.28487
[134,] -14378.70691 -14257.34057
[135,] -840.24802 -14378.70691
[136,] -437.06821 -840.24802
[137,] 33039.56820 -437.06821
[138,] 25797.11192 33039.56820
[139,] -462.28298 25797.11192
[140,] -5485.20066 -462.28298
[141,] 13407.22697 -5485.20066
[142,] -1197.24616 13407.22697
[143,] -14060.56644 -1197.24616
[144,] 10083.66237 -14060.56644
[145,] -1412.84638 10083.66237
[146,] 29795.46967 -1412.84638
[147,] 45272.50141 29795.46967
[148,] -2607.91980 45272.50141
[149,] 12693.33942 -2607.91980
[150,] -18138.07164 12693.33942
[151,] 56.00017 -18138.07164
[152,] 4753.77212 56.00017
[153,] -19143.03345 4753.77212
[154,] 16654.85303 -19143.03345
[155,] 57100.36150 16654.85303
[156,] -3133.80949 57100.36150
[157,] -2245.96031 -3133.80949
[158,] 17382.93061 -2245.96031
[159,] -1292.85801 17382.93061
[160,] 4980.54410 -1292.85801
[161,] 59936.47814 4980.54410
[162,] 435.44210 59936.47814
[163,] 20099.48923 435.44210
[164,] -10431.24259 20099.48923
[165,] -26557.16695 -10431.24259
[166,] -14169.04186 -26557.16695
[167,] 38351.49242 -14169.04186
[168,] -15939.85047 38351.49242
[169,] 28327.13088 -15939.85047
[170,] 14773.35202 28327.13088
[171,] 5079.18945 14773.35202
[172,] -1155.75109 5079.18945
[173,] -28317.82344 -1155.75109
[174,] 29932.53739 -28317.82344
[175,] -7062.77012 29932.53739
[176,] 37781.53758 -7062.77012
[177,] -7595.97683 37781.53758
[178,] 14598.57788 -7595.97683
[179,] -19310.05542 14598.57788
[180,] 7107.57548 -19310.05542
[181,] 22785.37996 7107.57548
[182,] 1112.05219 22785.37996
[183,] -1621.22598 1112.05219
[184,] -26853.88421 -1621.22598
[185,] -95.14192 -26853.88421
[186,] -3271.08092 -95.14192
[187,] -2662.81701 -3271.08092
[188,] -9947.06617 -2662.81701
[189,] 2945.29422 -9947.06617
[190,] 10009.59971 2945.29422
[191,] -722.77451 10009.59971
[192,] -16138.92387 -722.77451
[193,] -12829.75930 -16138.92387
[194,] 1487.18144 -12829.75930
[195,] -10396.54191 1487.18144
[196,] -7516.93014 -10396.54191
[197,] 36836.33241 -7516.93014
[198,] 29281.95284 36836.33241
[199,] 11436.72826 29281.95284
[200,] -15326.71249 11436.72826
[201,] -17187.28418 -15326.71249
[202,] -5387.42725 -17187.28418
[203,] -2501.48426 -5387.42725
[204,] -9850.61177 -2501.48426
[205,] 1697.91263 -9850.61177
[206,] -13727.90442 1697.91263
[207,] -3768.53773 -13727.90442
[208,] -9722.75475 -3768.53773
[209,] -19009.90735 -9722.75475
[210,] -16273.43012 -19009.90735
[211,] -3653.36288 -16273.43012
[212,] -5072.30643 -3653.36288
[213,] -13251.42331 -5072.30643
[214,] 10707.10774 -13251.42331
[215,] -21635.70426 10707.10774
[216,] 318.31532 -21635.70426
[217,] 3254.17359 318.31532
[218,] -14572.96969 3254.17359
[219,] -7271.33387 -14572.96969
[220,] -5571.32333 -7271.33387
[221,] -3618.17458 -5571.32333
[222,] -12976.00474 -3618.17458
[223,] -3313.80549 -12976.00474
[224,] -9792.58169 -3313.80549
[225,] -13987.41142 -9792.58169
[226,] 18796.45530 -13987.41142
[227,] 5134.65422 18796.45530
[228,] -13454.98243 5134.65422
[229,] 31974.01184 -13454.98243
[230,] -12591.97555 31974.01184
[231,] -7908.15174 -12591.97555
[232,] 24109.53309 -7908.15174
[233,] -2873.91257 24109.53309
[234,] 4280.40653 -2873.91257
[235,] -10256.82316 4280.40653
[236,] -12.20818 -10256.82316
[237,] -4326.15657 -12.20818
[238,] -9357.12474 -4326.15657
[239,] -2094.43091 -9357.12474
[240,] 9534.39607 -2094.43091
[241,] -15608.97318 9534.39607
[242,] -10296.10549 -15608.97318
[243,] -3544.40042 -10296.10549
[244,] 1463.98009 -3544.40042
[245,] 10115.83482 1463.98009
[246,] 6149.31175 10115.83482
[247,] 3211.07046 6149.31175
[248,] 16237.89459 3211.07046
[249,] 16155.59078 16237.89459
[250,] -60.36738 16155.59078
[251,] 2316.04826 -60.36738
[252,] -1806.38407 2316.04826
[253,] -5330.17573 -1806.38407
[254,] -8207.65406 -5330.17573
[255,] -1840.10101 -8207.65406
[256,] 4995.98976 -1840.10101
[257,] -1941.02370 4995.98976
[258,] -15229.06747 -1941.02370
[259,] -4547.10217 -15229.06747
[260,] -2498.71024 -4547.10217
[261,] -17410.87549 -2498.71024
[262,] -6388.05652 -17410.87549
[263,] -10765.46144 -6388.05652
[264,] 12140.41666 -10765.46144
[265,] -3123.89881 12140.41666
[266,] -11083.52325 -3123.89881
[267,] 8193.75376 -11083.52325
[268,] -18417.12595 8193.75376
[269,] 11511.68675 -18417.12595
[270,] -5588.40214 11511.68675
[271,] 91.43569 -5588.40214
[272,] 587.15568 91.43569
[273,] -7139.23629 587.15568
[274,] -151.72776 -7139.23629
[275,] -10468.39404 -151.72776
[276,] 14454.47639 -10468.39404
[277,] -5065.98177 14454.47639
[278,] -15581.63864 -5065.98177
[279,] -8982.68648 -15581.63864
[280,] -3721.46545 -8982.68648
[281,] 59676.89022 -3721.46545
[282,] 9182.13330 59676.89022
[283,] -3277.76741 9182.13330
[284,] -5662.96636 -3277.76741
[285,] 17395.76129 -5662.96636
[286,] -12618.52747 17395.76129
[287,] -7069.25255 -12618.52747
[288,] -3945.50545 -7069.25255
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 5102.47318 13416.25568
2 -7449.36293 5102.47318
3 -1031.73135 -7449.36293
4 -12023.45241 -1031.73135
5 25833.05897 -12023.45241
6 14848.53877 25833.05897
7 10288.33407 14848.53877
8 -9281.58093 10288.33407
9 26600.16800 -9281.58093
10 -4846.92614 26600.16800
11 -21703.01642 -4846.92614
12 -31050.59067 -21703.01642
13 -8743.36636 -31050.59067
14 -14915.46894 -8743.36636
15 -18798.46310 -14915.46894
16 -21807.54113 -18798.46310
17 24493.46257 -21807.54113
18 22453.82838 24493.46257
19 24575.03113 22453.82838
20 15856.73516 24575.03113
21 -17069.67360 15856.73516
22 -15778.86573 -17069.67360
23 -10281.63742 -15778.86573
24 5415.04800 -10281.63742
25 -24324.94949 5415.04800
26 -1444.06853 -24324.94949
27 1516.07614 -1444.06853
28 -25713.83902 1516.07614
29 -11710.07795 -25713.83902
30 2216.06244 -11710.07795
31 -8626.79259 2216.06244
32 4492.54650 -8626.79259
33 -25651.20409 4492.54650
34 15792.76531 -25651.20409
35 -17806.80365 15792.76531
36 -6806.53227 -17806.80365
37 11117.45739 -6806.53227
38 -21730.00274 11117.45739
39 -72.27331 -21730.00274
40 -12684.33025 -72.27331
41 -16702.04064 -12684.33025
42 -14451.07196 -16702.04064
43 -23764.81138 -14451.07196
44 -29203.00115 -23764.81138
45 -15941.86442 -29203.00115
46 -13786.70180 -15941.86442
47 -9509.77244 -13786.70180
48 -10921.45782 -9509.77244
49 18325.94966 -10921.45782
50 22183.55148 18325.94966
51 -11882.81805 22183.55148
52 -1633.00009 -11882.81805
53 34648.31272 -1633.00009
54 -9172.86105 34648.31272
55 -6658.07555 -9172.86105
56 -2198.50916 -6658.07555
57 1318.13399 -2198.50916
58 -3275.89066 1318.13399
59 -2161.24432 -3275.89066
60 -6528.12040 -2161.24432
61 13269.93435 -6528.12040
62 32744.95029 13269.93435
63 -11350.38299 32744.95029
64 -3184.71888 -11350.38299
65 -18558.20289 -3184.71888
66 -24938.33995 -18558.20289
67 10859.11826 -24938.33995
68 5569.61179 10859.11826
69 -17799.00356 5569.61179
70 -9857.34533 -17799.00356
71 -3434.25295 -9857.34533
72 -18935.82560 -3434.25295
73 27236.58067 -18935.82560
74 -31303.48520 27236.58067
75 14407.01434 -31303.48520
76 -661.82051 14407.01434
77 -28169.16422 -661.82051
78 -14256.40016 -28169.16422
79 -21483.33819 -14256.40016
80 7491.08683 -21483.33819
81 -8272.57258 7491.08683
82 -12045.63592 -8272.57258
83 1229.23678 -12045.63592
84 142655.56121 1229.23678
85 -12610.30204 142655.56121
86 -2778.93266 -12610.30204
87 -6700.86112 -2778.93266
88 3269.37273 -6700.86112
89 -17080.51837 3269.37273
90 -15943.30237 -17080.51837
91 28925.19606 -15943.30237
92 -11902.77347 28925.19606
93 -13661.36282 -11902.77347
94 -15543.78652 -13661.36282
95 12260.22554 -15543.78652
96 -33042.15211 12260.22554
97 16875.45721 -33042.15211
98 -16589.45048 16875.45721
99 6680.71027 -16589.45048
100 -3465.79236 6680.71027
101 -14848.37766 -3465.79236
102 1590.69189 -14848.37766
103 30714.93498 1590.69189
104 30275.10341 30714.93498
105 43660.88234 30275.10341
106 -29174.94891 43660.88234
107 -8331.33909 -29174.94891
108 16821.06988 -8331.33909
109 -18618.62455 16821.06988
110 30563.88097 -18618.62455
111 -8586.65882 30563.88097
112 4573.76565 -8586.65882
113 -22220.86137 4573.76565
114 -18185.30615 -22220.86137
115 -15035.98504 -18185.30615
116 -7217.89005 -15035.98504
117 -11880.16248 -7217.89005
118 125781.63928 -11880.16248
119 -2022.21295 125781.63928
120 -16904.24123 -2022.21295
121 -9028.42275 -16904.24123
122 -15308.04079 -9028.42275
123 93805.89921 -15308.04079
124 25265.57747 93805.89921
125 -11486.01195 25265.57747
126 8291.30755 -11486.01195
127 -5046.19090 8291.30755
128 3339.09771 -5046.19090
129 -2166.73755 3339.09771
130 -4898.49524 -2166.73755
131 -15372.68665 -4898.49524
132 -7216.28487 -15372.68665
133 -14257.34057 -7216.28487
134 -14378.70691 -14257.34057
135 -840.24802 -14378.70691
136 -437.06821 -840.24802
137 33039.56820 -437.06821
138 25797.11192 33039.56820
139 -462.28298 25797.11192
140 -5485.20066 -462.28298
141 13407.22697 -5485.20066
142 -1197.24616 13407.22697
143 -14060.56644 -1197.24616
144 10083.66237 -14060.56644
145 -1412.84638 10083.66237
146 29795.46967 -1412.84638
147 45272.50141 29795.46967
148 -2607.91980 45272.50141
149 12693.33942 -2607.91980
150 -18138.07164 12693.33942
151 56.00017 -18138.07164
152 4753.77212 56.00017
153 -19143.03345 4753.77212
154 16654.85303 -19143.03345
155 57100.36150 16654.85303
156 -3133.80949 57100.36150
157 -2245.96031 -3133.80949
158 17382.93061 -2245.96031
159 -1292.85801 17382.93061
160 4980.54410 -1292.85801
161 59936.47814 4980.54410
162 435.44210 59936.47814
163 20099.48923 435.44210
164 -10431.24259 20099.48923
165 -26557.16695 -10431.24259
166 -14169.04186 -26557.16695
167 38351.49242 -14169.04186
168 -15939.85047 38351.49242
169 28327.13088 -15939.85047
170 14773.35202 28327.13088
171 5079.18945 14773.35202
172 -1155.75109 5079.18945
173 -28317.82344 -1155.75109
174 29932.53739 -28317.82344
175 -7062.77012 29932.53739
176 37781.53758 -7062.77012
177 -7595.97683 37781.53758
178 14598.57788 -7595.97683
179 -19310.05542 14598.57788
180 7107.57548 -19310.05542
181 22785.37996 7107.57548
182 1112.05219 22785.37996
183 -1621.22598 1112.05219
184 -26853.88421 -1621.22598
185 -95.14192 -26853.88421
186 -3271.08092 -95.14192
187 -2662.81701 -3271.08092
188 -9947.06617 -2662.81701
189 2945.29422 -9947.06617
190 10009.59971 2945.29422
191 -722.77451 10009.59971
192 -16138.92387 -722.77451
193 -12829.75930 -16138.92387
194 1487.18144 -12829.75930
195 -10396.54191 1487.18144
196 -7516.93014 -10396.54191
197 36836.33241 -7516.93014
198 29281.95284 36836.33241
199 11436.72826 29281.95284
200 -15326.71249 11436.72826
201 -17187.28418 -15326.71249
202 -5387.42725 -17187.28418
203 -2501.48426 -5387.42725
204 -9850.61177 -2501.48426
205 1697.91263 -9850.61177
206 -13727.90442 1697.91263
207 -3768.53773 -13727.90442
208 -9722.75475 -3768.53773
209 -19009.90735 -9722.75475
210 -16273.43012 -19009.90735
211 -3653.36288 -16273.43012
212 -5072.30643 -3653.36288
213 -13251.42331 -5072.30643
214 10707.10774 -13251.42331
215 -21635.70426 10707.10774
216 318.31532 -21635.70426
217 3254.17359 318.31532
218 -14572.96969 3254.17359
219 -7271.33387 -14572.96969
220 -5571.32333 -7271.33387
221 -3618.17458 -5571.32333
222 -12976.00474 -3618.17458
223 -3313.80549 -12976.00474
224 -9792.58169 -3313.80549
225 -13987.41142 -9792.58169
226 18796.45530 -13987.41142
227 5134.65422 18796.45530
228 -13454.98243 5134.65422
229 31974.01184 -13454.98243
230 -12591.97555 31974.01184
231 -7908.15174 -12591.97555
232 24109.53309 -7908.15174
233 -2873.91257 24109.53309
234 4280.40653 -2873.91257
235 -10256.82316 4280.40653
236 -12.20818 -10256.82316
237 -4326.15657 -12.20818
238 -9357.12474 -4326.15657
239 -2094.43091 -9357.12474
240 9534.39607 -2094.43091
241 -15608.97318 9534.39607
242 -10296.10549 -15608.97318
243 -3544.40042 -10296.10549
244 1463.98009 -3544.40042
245 10115.83482 1463.98009
246 6149.31175 10115.83482
247 3211.07046 6149.31175
248 16237.89459 3211.07046
249 16155.59078 16237.89459
250 -60.36738 16155.59078
251 2316.04826 -60.36738
252 -1806.38407 2316.04826
253 -5330.17573 -1806.38407
254 -8207.65406 -5330.17573
255 -1840.10101 -8207.65406
256 4995.98976 -1840.10101
257 -1941.02370 4995.98976
258 -15229.06747 -1941.02370
259 -4547.10217 -15229.06747
260 -2498.71024 -4547.10217
261 -17410.87549 -2498.71024
262 -6388.05652 -17410.87549
263 -10765.46144 -6388.05652
264 12140.41666 -10765.46144
265 -3123.89881 12140.41666
266 -11083.52325 -3123.89881
267 8193.75376 -11083.52325
268 -18417.12595 8193.75376
269 11511.68675 -18417.12595
270 -5588.40214 11511.68675
271 91.43569 -5588.40214
272 587.15568 91.43569
273 -7139.23629 587.15568
274 -151.72776 -7139.23629
275 -10468.39404 -151.72776
276 14454.47639 -10468.39404
277 -5065.98177 14454.47639
278 -15581.63864 -5065.98177
279 -8982.68648 -15581.63864
280 -3721.46545 -8982.68648
281 59676.89022 -3721.46545
282 9182.13330 59676.89022
283 -3277.76741 9182.13330
284 -5662.96636 -3277.76741
285 17395.76129 -5662.96636
286 -12618.52747 17395.76129
287 -7069.25255 -12618.52747
288 -3945.50545 -7069.25255
> 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/7uvj11355259025.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/8t0tc1355259025.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/9zf7x1355259025.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/10mxqq1355259025.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/11pm5j1355259025.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/120tyf1355259025.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/13nc5j1355259025.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/14jyyn1355259025.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/15sbsb1355259025.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/165puk1355259025.tab")
+ }
>
> try(system("convert tmp/1iyzf1355259025.ps tmp/1iyzf1355259025.png",intern=TRUE))
character(0)
> try(system("convert tmp/23r8r1355259025.ps tmp/23r8r1355259025.png",intern=TRUE))
character(0)
> try(system("convert tmp/3fnxm1355259025.ps tmp/3fnxm1355259025.png",intern=TRUE))
character(0)
> try(system("convert tmp/4l98h1355259025.ps tmp/4l98h1355259025.png",intern=TRUE))
character(0)
> try(system("convert tmp/59cbk1355259025.ps tmp/59cbk1355259025.png",intern=TRUE))
character(0)
> try(system("convert tmp/6qzko1355259025.ps tmp/6qzko1355259025.png",intern=TRUE))
character(0)
> try(system("convert tmp/7uvj11355259025.ps tmp/7uvj11355259025.png",intern=TRUE))
character(0)
> try(system("convert tmp/8t0tc1355259025.ps tmp/8t0tc1355259025.png",intern=TRUE))
character(0)
> try(system("convert tmp/9zf7x1355259025.ps tmp/9zf7x1355259025.png",intern=TRUE))
character(0)
> try(system("convert tmp/10mxqq1355259025.ps tmp/10mxqq1355259025.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
14.699 1.343 16.133