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(210907
+ ,56
+ ,396
+ ,81
+ ,3
+ ,79
+ ,30
+ ,115
+ ,94
+ ,112285
+ ,24188
+ ,146283
+ ,144
+ ,145
+ ,120982
+ ,56
+ ,297
+ ,55
+ ,4
+ ,58
+ ,28
+ ,109
+ ,103
+ ,84786
+ ,18273
+ ,98364
+ ,103
+ ,101
+ ,176508
+ ,54
+ ,559
+ ,50
+ ,12
+ ,60
+ ,38
+ ,146
+ ,93
+ ,83123
+ ,14130
+ ,86146
+ ,98
+ ,98
+ ,179321
+ ,89
+ ,967
+ ,125
+ ,2
+ ,108
+ ,30
+ ,116
+ ,103
+ ,101193
+ ,32287
+ ,96933
+ ,135
+ ,132
+ ,123185
+ ,40
+ ,270
+ ,40
+ ,1
+ ,49
+ ,22
+ ,68
+ ,51
+ ,38361
+ ,8654
+ ,79234
+ ,61
+ ,60
+ ,52746
+ ,25
+ ,143
+ ,37
+ ,3
+ ,0
+ ,26
+ ,101
+ ,70
+ ,68504
+ ,9245
+ ,42551
+ ,39
+ ,38
+ ,385534
+ ,92
+ ,1562
+ ,63
+ ,0
+ ,121
+ ,25
+ ,96
+ ,91
+ ,119182
+ ,33251
+ ,195663
+ ,150
+ ,144
+ ,33170
+ ,18
+ ,109
+ ,44
+ ,0
+ ,1
+ ,18
+ ,67
+ ,22
+ ,22807
+ ,1271
+ ,6853
+ ,5
+ ,5
+ ,101645
+ ,63
+ ,371
+ ,88
+ ,0
+ ,20
+ ,11
+ ,44
+ ,38
+ ,17140
+ ,5279
+ ,21529
+ ,28
+ ,28
+ ,149061
+ ,44
+ ,656
+ ,66
+ ,5
+ ,43
+ ,26
+ ,100
+ ,93
+ ,116174
+ ,27101
+ ,95757
+ ,84
+ ,84
+ ,165446
+ ,33
+ ,511
+ ,57
+ ,0
+ ,69
+ ,25
+ ,93
+ ,60
+ ,57635
+ ,16373
+ ,85584
+ ,80
+ ,79
+ ,237213
+ ,84
+ ,655
+ ,74
+ ,0
+ ,78
+ ,38
+ ,140
+ ,123
+ ,66198
+ ,19716
+ ,143983
+ ,130
+ ,127
+ ,173326
+ ,88
+ ,465
+ ,49
+ ,7
+ ,86
+ ,44
+ ,166
+ ,148
+ ,71701
+ ,17753
+ ,75851
+ ,82
+ ,78
+ ,133131
+ ,55
+ ,525
+ ,52
+ ,7
+ ,44
+ ,30
+ ,99
+ ,90
+ ,57793
+ ,9028
+ ,59238
+ ,60
+ ,60
+ ,258873
+ ,60
+ ,885
+ ,88
+ ,3
+ ,104
+ ,40
+ ,139
+ ,124
+ ,80444
+ ,18653
+ ,93163
+ ,131
+ ,131
+ ,180083
+ ,66
+ ,497
+ ,36
+ ,9
+ ,63
+ ,34
+ ,130
+ ,70
+ ,53855
+ ,8828
+ ,96037
+ ,84
+ ,84
+ ,324799
+ ,154
+ ,1436
+ ,108
+ ,0
+ ,158
+ ,47
+ ,181
+ ,168
+ ,97668
+ ,29498
+ ,151511
+ ,140
+ ,133
+ ,230964
+ ,53
+ ,612
+ ,43
+ ,4
+ ,102
+ ,30
+ ,116
+ ,115
+ ,133824
+ ,27563
+ ,136368
+ ,151
+ ,150
+ ,236785
+ ,119
+ ,865
+ ,75
+ ,3
+ ,77
+ ,31
+ ,116
+ ,71
+ ,101481
+ ,18293
+ ,112642
+ ,91
+ ,91
+ ,135473
+ ,41
+ ,385
+ ,32
+ ,0
+ ,82
+ ,23
+ ,88
+ ,66
+ ,99645
+ ,22530
+ ,94728
+ ,138
+ ,132
+ ,202925
+ ,61
+ ,567
+ ,44
+ ,7
+ ,115
+ ,36
+ ,139
+ ,134
+ ,114789
+ ,15977
+ ,105499
+ ,150
+ ,136
+ ,215147
+ ,58
+ ,639
+ ,85
+ ,0
+ ,101
+ ,36
+ ,135
+ ,117
+ ,99052
+ ,35082
+ ,121527
+ ,124
+ ,124
+ ,344297
+ ,75
+ ,963
+ ,86
+ ,1
+ ,80
+ ,30
+ ,108
+ ,108
+ ,67654
+ ,16116
+ ,127766
+ ,119
+ ,118
+ ,153935
+ ,33
+ ,398
+ ,56
+ ,5
+ ,50
+ ,25
+ ,89
+ ,84
+ ,65553
+ ,15849
+ ,98958
+ ,73
+ ,70
+ ,132943
+ ,40
+ ,410
+ ,50
+ ,7
+ ,83
+ ,39
+ ,156
+ ,156
+ ,97500
+ ,16026
+ ,77900
+ ,110
+ ,107
+ ,174724
+ ,92
+ ,966
+ ,135
+ ,0
+ ,123
+ ,34
+ ,129
+ ,120
+ ,69112
+ ,26569
+ ,85646
+ ,123
+ ,119
+ ,174415
+ ,100
+ ,801
+ ,63
+ ,0
+ ,73
+ ,31
+ ,118
+ ,114
+ ,82753
+ ,24785
+ ,98579
+ ,90
+ ,89
+ ,225548
+ ,112
+ ,892
+ ,81
+ ,5
+ ,81
+ ,31
+ ,118
+ ,94
+ ,85323
+ ,17569
+ ,130767
+ ,116
+ ,112
+ ,223632
+ ,73
+ ,513
+ ,52
+ ,0
+ ,105
+ ,33
+ ,125
+ ,120
+ ,72654
+ ,23825
+ ,131741
+ ,113
+ ,108
+ ,124817
+ ,40
+ ,469
+ ,44
+ ,0
+ ,47
+ ,25
+ ,95
+ ,81
+ ,30727
+ ,7869
+ ,53907
+ ,56
+ ,52
+ ,221698
+ ,45
+ ,683
+ ,113
+ ,0
+ ,105
+ ,33
+ ,126
+ ,110
+ ,77873
+ ,14975
+ ,178812
+ ,115
+ ,112
+ ,210767
+ ,60
+ ,643
+ ,39
+ ,3
+ ,94
+ ,35
+ ,135
+ ,133
+ ,117478
+ ,37791
+ ,146761
+ ,119
+ ,116
+ ,170266
+ ,62
+ ,535
+ ,73
+ ,4
+ ,44
+ ,42
+ ,154
+ ,122
+ ,74007
+ ,9605
+ ,82036
+ ,129
+ ,123
+ ,260561
+ ,75
+ ,625
+ ,48
+ ,1
+ ,114
+ ,43
+ ,165
+ ,158
+ ,90183
+ ,27295
+ ,163253
+ ,127
+ ,125
+ ,84853
+ ,31
+ ,264
+ ,33
+ ,4
+ ,38
+ ,30
+ ,113
+ ,109
+ ,61542
+ ,2746
+ ,27032
+ ,27
+ ,27
+ ,294424
+ ,77
+ ,992
+ ,59
+ ,2
+ ,107
+ ,33
+ ,127
+ ,124
+ ,101494
+ ,34461
+ ,171975
+ ,175
+ ,162
+ ,101011
+ ,34
+ ,238
+ ,41
+ ,0
+ ,30
+ ,13
+ ,52
+ ,39
+ ,27570
+ ,8098
+ ,65990
+ ,35
+ ,32
+ ,215641
+ ,46
+ ,818
+ ,69
+ ,0
+ ,71
+ ,32
+ ,121
+ ,92
+ ,55813
+ ,4787
+ ,86572
+ ,64
+ ,64
+ ,325107
+ ,99
+ ,937
+ ,64
+ ,0
+ ,84
+ ,36
+ ,136
+ ,126
+ ,79215
+ ,24919
+ ,159676
+ ,96
+ ,92
+ ,7176
+ ,17
+ ,70
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,1423
+ ,603
+ ,1929
+ ,0
+ ,0
+ ,167542
+ ,66
+ ,507
+ ,59
+ ,2
+ ,59
+ ,28
+ ,108
+ ,70
+ ,55461
+ ,16329
+ ,85371
+ ,84
+ ,83
+ ,106408
+ ,30
+ ,260
+ ,32
+ ,1
+ ,33
+ ,14
+ ,46
+ ,37
+ ,31081
+ ,12558
+ ,58391
+ ,41
+ ,41
+ ,96560
+ ,76
+ ,503
+ ,129
+ ,0
+ ,42
+ ,17
+ ,54
+ ,38
+ ,22996
+ ,7784
+ ,31580
+ ,47
+ ,47
+ ,265769
+ ,146
+ ,927
+ ,37
+ ,2
+ ,96
+ ,32
+ ,124
+ ,120
+ ,83122
+ ,28522
+ ,136815
+ ,126
+ ,120
+ ,269651
+ ,67
+ ,1269
+ ,31
+ ,10
+ ,106
+ ,30
+ ,115
+ ,93
+ ,70106
+ ,22265
+ ,120642
+ ,105
+ ,105
+ ,149112
+ ,56
+ ,537
+ ,65
+ ,6
+ ,56
+ ,35
+ ,128
+ ,95
+ ,60578
+ ,14459
+ ,69107
+ ,80
+ ,79
+ ,175824
+ ,107
+ ,910
+ ,107
+ ,0
+ ,57
+ ,20
+ ,80
+ ,77
+ ,39992
+ ,14526
+ ,50495
+ ,70
+ ,65
+ ,152871
+ ,58
+ ,532
+ ,74
+ ,5
+ ,59
+ ,28
+ ,97
+ ,90
+ ,79892
+ ,22240
+ ,108016
+ ,73
+ ,70
+ ,111665
+ ,34
+ ,345
+ ,54
+ ,4
+ ,39
+ ,28
+ ,104
+ ,80
+ ,49810
+ ,11802
+ ,46341
+ ,57
+ ,55
+ ,116408
+ ,61
+ ,918
+ ,76
+ ,1
+ ,34
+ ,39
+ ,59
+ ,31
+ ,71570
+ ,7623
+ ,78348
+ ,40
+ ,39
+ ,362301
+ ,119
+ ,1635
+ ,715
+ ,2
+ ,76
+ ,34
+ ,125
+ ,110
+ ,100708
+ ,11912
+ ,79336
+ ,68
+ ,67
+ ,78800
+ ,42
+ ,330
+ ,57
+ ,2
+ ,20
+ ,26
+ ,82
+ ,66
+ ,33032
+ ,7935
+ ,56968
+ ,21
+ ,21
+ ,183167
+ ,66
+ ,557
+ ,66
+ ,0
+ ,91
+ ,39
+ ,149
+ ,138
+ ,82875
+ ,18220
+ ,93176
+ ,127
+ ,127
+ ,277965
+ ,89
+ ,1178
+ ,106
+ ,8
+ ,115
+ ,39
+ ,149
+ ,133
+ ,139077
+ ,19199
+ ,161632
+ ,154
+ ,152
+ ,150629
+ ,44
+ ,740
+ ,54
+ ,3
+ ,85
+ ,33
+ ,122
+ ,113
+ ,71595
+ ,19918
+ ,87850
+ ,116
+ ,113
+ ,168809
+ ,66
+ ,452
+ ,32
+ ,0
+ ,76
+ ,28
+ ,118
+ ,100
+ ,72260
+ ,21884
+ ,127969
+ ,102
+ ,99
+ ,24188
+ ,24
+ ,218
+ ,20
+ ,0
+ ,8
+ ,4
+ ,12
+ ,7
+ ,5950
+ ,2694
+ ,15049
+ ,7
+ ,7
+ ,329267
+ ,259
+ ,764
+ ,71
+ ,8
+ ,79
+ ,39
+ ,144
+ ,140
+ ,115762
+ ,15808
+ ,155135
+ ,148
+ ,141
+ ,65029
+ ,17
+ ,255
+ ,21
+ ,5
+ ,21
+ ,18
+ ,67
+ ,61
+ ,32551
+ ,3597
+ ,25109
+ ,21
+ ,21
+ ,101097
+ ,64
+ ,454
+ ,70
+ ,3
+ ,30
+ ,14
+ ,52
+ ,41
+ ,31701
+ ,5296
+ ,45824
+ ,35
+ ,35
+ ,218946
+ ,41
+ ,866
+ ,112
+ ,1
+ ,76
+ ,29
+ ,108
+ ,96
+ ,80670
+ ,25239
+ ,102996
+ ,112
+ ,109
+ ,244052
+ ,68
+ ,574
+ ,66
+ ,5
+ ,101
+ ,44
+ ,166
+ ,164
+ ,143558
+ ,29801
+ ,160604
+ ,137
+ ,133
+ ,341570
+ ,168
+ ,1276
+ ,190
+ ,1
+ ,94
+ ,21
+ ,80
+ ,78
+ ,117105
+ ,18450
+ ,158051
+ ,135
+ ,123
+ ,103597
+ ,43
+ ,379
+ ,66
+ ,1
+ ,27
+ ,16
+ ,60
+ ,49
+ ,23789
+ ,7132
+ ,44547
+ ,26
+ ,26
+ ,233328
+ ,132
+ ,825
+ ,165
+ ,5
+ ,92
+ ,28
+ ,107
+ ,102
+ ,120733
+ ,34861
+ ,162647
+ ,230
+ ,230
+ ,256462
+ ,105
+ ,798
+ ,56
+ ,0
+ ,123
+ ,35
+ ,127
+ ,124
+ ,105195
+ ,35940
+ ,174141
+ ,181
+ ,166
+ ,206161
+ ,71
+ ,663
+ ,61
+ ,12
+ ,75
+ ,28
+ ,107
+ ,99
+ ,73107
+ ,16688
+ ,60622
+ ,71
+ ,68
+ ,311473
+ ,112
+ ,1069
+ ,53
+ ,8
+ ,128
+ ,38
+ ,146
+ ,129
+ ,132068
+ ,24683
+ ,179566
+ ,147
+ ,147
+ ,235800
+ ,94
+ ,921
+ ,127
+ ,8
+ ,105
+ ,23
+ ,84
+ ,62
+ ,149193
+ ,46230
+ ,184301
+ ,190
+ ,179
+ ,177939
+ ,82
+ ,858
+ ,63
+ ,8
+ ,55
+ ,36
+ ,141
+ ,73
+ ,46821
+ ,10387
+ ,75661
+ ,64
+ ,61
+ ,207176
+ ,70
+ ,711
+ ,38
+ ,8
+ ,56
+ ,32
+ ,123
+ ,114
+ ,87011
+ ,21436
+ ,96144
+ ,105
+ ,101
+ ,196553
+ ,57
+ ,503
+ ,50
+ ,2
+ ,41
+ ,29
+ ,111
+ ,99
+ ,95260
+ ,30546
+ ,129847
+ ,107
+ ,108
+ ,174184
+ ,53
+ ,382
+ ,52
+ ,0
+ ,72
+ ,25
+ ,98
+ ,70
+ ,55183
+ ,19746
+ ,117286
+ ,94
+ ,90
+ ,143246
+ ,103
+ ,464
+ ,42
+ ,5
+ ,67
+ ,27
+ ,105
+ ,104
+ ,106671
+ ,15977
+ ,71180
+ ,116
+ ,114
+ ,187559
+ ,121
+ ,717
+ ,76
+ ,8
+ ,75
+ ,36
+ ,135
+ ,116
+ ,73511
+ ,22583
+ ,109377
+ ,106
+ ,103
+ ,187681
+ ,62
+ ,690
+ ,67
+ ,2
+ ,114
+ ,28
+ ,107
+ ,91
+ ,92945
+ ,17274
+ ,85298
+ ,143
+ ,142
+ ,119016
+ ,52
+ ,462
+ ,50
+ ,5
+ ,118
+ ,23
+ ,85
+ ,74
+ ,78664
+ ,16469
+ ,73631
+ ,81
+ ,79
+ ,182192
+ ,52
+ ,657
+ ,53
+ ,12
+ ,77
+ ,40
+ ,155
+ ,138
+ ,70054
+ ,14251
+ ,86767
+ ,89
+ ,88
+ ,73566
+ ,32
+ ,385
+ ,39
+ ,6
+ ,22
+ ,23
+ ,88
+ ,67
+ ,22618
+ ,3007
+ ,23824
+ ,26
+ ,25
+ ,194979
+ ,62
+ ,577
+ ,50
+ ,7
+ ,66
+ ,40
+ ,155
+ ,151
+ ,74011
+ ,16851
+ ,93487
+ ,84
+ ,83
+ ,167488
+ ,45
+ ,619
+ ,77
+ ,2
+ ,69
+ ,28
+ ,104
+ ,72
+ ,83737
+ ,21113
+ ,82981
+ ,113
+ ,113
+ ,143756
+ ,46
+ ,479
+ ,57
+ ,0
+ ,105
+ ,34
+ ,132
+ ,120
+ ,69094
+ ,17401
+ ,73815
+ ,120
+ ,118
+ ,275541
+ ,63
+ ,817
+ ,73
+ ,4
+ ,116
+ ,33
+ ,127
+ ,115
+ ,93133
+ ,23958
+ ,94552
+ ,110
+ ,110
+ ,243199
+ ,75
+ ,752
+ ,34
+ ,3
+ ,88
+ ,28
+ ,108
+ ,105
+ ,95536
+ ,23567
+ ,132190
+ ,134
+ ,129
+ ,182999
+ ,88
+ ,430
+ ,39
+ ,6
+ ,73
+ ,34
+ ,129
+ ,104
+ ,225920
+ ,13065
+ ,128754
+ ,54
+ ,51
+ ,135649
+ ,46
+ ,451
+ ,46
+ ,2
+ ,99
+ ,30
+ ,116
+ ,108
+ ,62133
+ ,15358
+ ,66363
+ ,96
+ ,93
+ ,152299
+ ,53
+ ,537
+ ,63
+ ,0
+ ,62
+ ,33
+ ,122
+ ,98
+ ,61370
+ ,14587
+ ,67808
+ ,78
+ ,76
+ ,120221
+ ,37
+ ,519
+ ,35
+ ,1
+ ,53
+ ,22
+ ,85
+ ,69
+ ,43836
+ ,12770
+ ,61724
+ ,51
+ ,49
+ ,346485
+ ,90
+ ,1000
+ ,106
+ ,0
+ ,118
+ ,38
+ ,147
+ ,111
+ ,106117
+ ,24021
+ ,131722
+ ,121
+ ,118
+ ,145790
+ ,63
+ ,637
+ ,43
+ ,5
+ ,30
+ ,26
+ ,99
+ ,99
+ ,38692
+ ,9648
+ ,68580
+ ,38
+ ,38
+ ,193339
+ ,78
+ ,465
+ ,47
+ ,2
+ ,100
+ ,35
+ ,87
+ ,71
+ ,84651
+ ,20537
+ ,106175
+ ,145
+ ,141
+ ,80953
+ ,25
+ ,437
+ ,31
+ ,0
+ ,49
+ ,8
+ ,28
+ ,27
+ ,56622
+ ,7905
+ ,55792
+ ,59
+ ,58
+ ,122774
+ ,45
+ ,711
+ ,162
+ ,0
+ ,24
+ ,24
+ ,90
+ ,69
+ ,15986
+ ,4527
+ ,25157
+ ,27
+ ,27
+ ,130585
+ ,46
+ ,299
+ ,57
+ ,5
+ ,67
+ ,29
+ ,109
+ ,107
+ ,95364
+ ,30495
+ ,76669
+ ,91
+ ,91
+ ,112611
+ ,41
+ ,248
+ ,36
+ ,0
+ ,46
+ ,20
+ ,78
+ ,73
+ ,26706
+ ,7117
+ ,57283
+ ,48
+ ,48
+ ,286468
+ ,144
+ ,1162
+ ,263
+ ,1
+ ,57
+ ,29
+ ,111
+ ,107
+ ,89691
+ ,17719
+ ,105805
+ ,68
+ ,63
+ ,241066
+ ,82
+ ,714
+ ,78
+ ,0
+ ,75
+ ,45
+ ,158
+ ,93
+ ,67267
+ ,27056
+ ,129484
+ ,58
+ ,56
+ ,148446
+ ,91
+ ,905
+ ,63
+ ,1
+ ,135
+ ,37
+ ,141
+ ,129
+ ,126846
+ ,33473
+ ,72413
+ ,150
+ ,144
+ ,204713
+ ,71
+ ,649
+ ,54
+ ,1
+ ,68
+ ,33
+ ,122
+ ,69
+ ,41140
+ ,9758
+ ,87831
+ ,74
+ ,73
+ ,182079
+ ,63
+ ,512
+ ,63
+ ,2
+ ,124
+ ,33
+ ,124
+ ,118
+ ,102860
+ ,21115
+ ,96971
+ ,181
+ ,168
+ ,140344
+ ,53
+ ,472
+ ,77
+ ,6
+ ,33
+ ,25
+ ,93
+ ,73
+ ,51715
+ ,7236
+ ,71299
+ ,65
+ ,64
+ ,220516
+ ,62
+ ,905
+ ,79
+ ,1
+ ,98
+ ,32
+ ,124
+ ,119
+ ,55801
+ ,13790
+ ,77494
+ ,97
+ ,97
+ ,243060
+ ,63
+ ,786
+ ,110
+ ,4
+ ,58
+ ,29
+ ,112
+ ,104
+ ,111813
+ ,32902
+ ,120336
+ ,121
+ ,117
+ ,162765
+ ,32
+ ,489
+ ,56
+ ,2
+ ,68
+ ,28
+ ,108
+ ,107
+ ,120293
+ ,25131
+ ,93913
+ ,99
+ ,100
+ ,182613
+ ,39
+ ,479
+ ,56
+ ,3
+ ,81
+ ,28
+ ,99
+ ,99
+ ,138599
+ ,30910
+ ,136048
+ ,152
+ ,149
+ ,232138
+ ,62
+ ,617
+ ,43
+ ,0
+ ,131
+ ,31
+ ,117
+ ,90
+ ,161647
+ ,35947
+ ,181248
+ ,188
+ ,187
+ ,265318
+ ,117
+ ,925
+ ,111
+ ,10
+ ,110
+ ,52
+ ,199
+ ,197
+ ,115929
+ ,29848
+ ,146123
+ ,138
+ ,127
+ ,85574
+ ,34
+ ,351
+ ,71
+ ,0
+ ,37
+ ,21
+ ,78
+ ,36
+ ,24266
+ ,6943
+ ,32036
+ ,40
+ ,37
+ ,310839
+ ,92
+ ,1144
+ ,62
+ ,9
+ ,130
+ ,24
+ ,91
+ ,85
+ ,162901
+ ,42705
+ ,186646
+ ,254
+ ,245
+ ,225060
+ ,93
+ ,669
+ ,56
+ ,7
+ ,93
+ ,41
+ ,158
+ ,139
+ ,109825
+ ,31808
+ ,102255
+ ,87
+ ,87
+ ,232317
+ ,54
+ ,707
+ ,74
+ ,0
+ ,118
+ ,33
+ ,126
+ ,106
+ ,129838
+ ,26675
+ ,168237
+ ,178
+ ,177
+ ,144966
+ ,144
+ ,458
+ ,60
+ ,0
+ ,39
+ ,32
+ ,122
+ ,50
+ ,37510
+ ,8435
+ ,64219
+ ,51
+ ,49
+ ,43287
+ ,14
+ ,214
+ ,43
+ ,4
+ ,13
+ ,19
+ ,71
+ ,64
+ ,43750
+ ,7409
+ ,19630
+ ,49
+ ,49
+ ,155754
+ ,61
+ ,599
+ ,68
+ ,4
+ ,74
+ ,20
+ ,75
+ ,31
+ ,40652
+ ,14993
+ ,76825
+ ,73
+ ,73
+ ,164709
+ ,109
+ ,572
+ ,53
+ ,0
+ ,81
+ ,31
+ ,115
+ ,63
+ ,87771
+ ,36867
+ ,115338
+ ,176
+ ,177
+ ,201940
+ ,38
+ ,897
+ ,87
+ ,0
+ ,109
+ ,31
+ ,119
+ ,92
+ ,85872
+ ,33835
+ ,109427
+ ,94
+ ,94
+ ,235454
+ ,73
+ ,819
+ ,46
+ ,0
+ ,151
+ ,32
+ ,124
+ ,106
+ ,89275
+ ,24164
+ ,118168
+ ,120
+ ,117
+ ,220801
+ ,75
+ ,720
+ ,105
+ ,1
+ ,51
+ ,18
+ ,72
+ ,63
+ ,44418
+ ,12607
+ ,84845
+ ,66
+ ,60
+ ,99466
+ ,50
+ ,273
+ ,32
+ ,0
+ ,28
+ ,23
+ ,91
+ ,69
+ ,192565
+ ,22609
+ ,153197
+ ,56
+ ,55
+ ,92661
+ ,61
+ ,508
+ ,133
+ ,1
+ ,40
+ ,17
+ ,45
+ ,41
+ ,35232
+ ,5892
+ ,29877
+ ,39
+ ,39
+ ,133328
+ ,55
+ ,506
+ ,79
+ ,0
+ ,56
+ ,20
+ ,78
+ ,56
+ ,40909
+ ,17014
+ ,63506
+ ,66
+ ,64
+ ,61361
+ ,77
+ ,451
+ ,51
+ ,0
+ ,27
+ ,12
+ ,39
+ ,25
+ ,13294
+ ,5394
+ ,22445
+ ,27
+ ,26
+ ,125930
+ ,75
+ ,699
+ ,207
+ ,4
+ ,37
+ ,17
+ ,68
+ ,65
+ ,32387
+ ,9178
+ ,47695
+ ,65
+ ,64
+ ,100750
+ ,72
+ ,407
+ ,67
+ ,0
+ ,83
+ ,30
+ ,119
+ ,93
+ ,140867
+ ,6440
+ ,68370
+ ,58
+ ,58
+ ,224549
+ ,50
+ ,465
+ ,47
+ ,4
+ ,54
+ ,31
+ ,117
+ ,114
+ ,120662
+ ,21916
+ ,146304
+ ,98
+ ,95
+ ,82316
+ ,32
+ ,245
+ ,34
+ ,4
+ ,27
+ ,10
+ ,39
+ ,38
+ ,21233
+ ,4011
+ ,38233
+ ,25
+ ,25
+ ,102010
+ ,53
+ ,370
+ ,66
+ ,3
+ ,28
+ ,13
+ ,50
+ ,44
+ ,44332
+ ,5818
+ ,42071
+ ,26
+ ,26
+ ,101523
+ ,42
+ ,316
+ ,76
+ ,0
+ ,59
+ ,22
+ ,88
+ ,87
+ ,61056
+ ,18647
+ ,50517
+ ,77
+ ,76
+ ,243511
+ ,71
+ ,603
+ ,65
+ ,0
+ ,133
+ ,42
+ ,155
+ ,110
+ ,101338
+ ,20556
+ ,103950
+ ,130
+ ,129
+ ,22938
+ ,10
+ ,154
+ ,9
+ ,0
+ ,12
+ ,1
+ ,0
+ ,0
+ ,1168
+ ,238
+ ,5841
+ ,11
+ ,11
+ ,41566
+ ,35
+ ,229
+ ,42
+ ,5
+ ,0
+ ,9
+ ,36
+ ,27
+ ,13497
+ ,70
+ ,2341
+ ,2
+ ,2
+ ,152474
+ ,65
+ ,577
+ ,45
+ ,0
+ ,106
+ ,32
+ ,123
+ ,83
+ ,65567
+ ,22392
+ ,84396
+ ,101
+ ,101
+ ,61857
+ ,25
+ ,192
+ ,25
+ ,4
+ ,23
+ ,11
+ ,32
+ ,30
+ ,25162
+ ,3913
+ ,24610
+ ,31
+ ,28
+ ,99923
+ ,66
+ ,617
+ ,115
+ ,0
+ ,44
+ ,25
+ ,99
+ ,80
+ ,32334
+ ,12237
+ ,35753
+ ,36
+ ,36
+ ,132487
+ ,41
+ ,411
+ ,97
+ ,0
+ ,71
+ ,36
+ ,136
+ ,98
+ ,40735
+ ,8388
+ ,55515
+ ,120
+ ,89
+ ,317394
+ ,86
+ ,975
+ ,53
+ ,1
+ ,116
+ ,31
+ ,117
+ ,82
+ ,91413
+ ,22120
+ ,209056
+ ,195
+ ,193
+ ,21054
+ ,16
+ ,146
+ ,2
+ ,0
+ ,4
+ ,0
+ ,0
+ ,0
+ ,855
+ ,338
+ ,6622
+ ,4
+ ,4
+ ,209641
+ ,42
+ ,705
+ ,52
+ ,5
+ ,62
+ ,24
+ ,88
+ ,60
+ ,97068
+ ,11727
+ ,115814
+ ,89
+ ,84
+ ,22648
+ ,19
+ ,184
+ ,44
+ ,0
+ ,12
+ ,13
+ ,39
+ ,28
+ ,44339
+ ,3704
+ ,11609
+ ,24
+ ,23
+ ,31414
+ ,19
+ ,200
+ ,22
+ ,0
+ ,18
+ ,8
+ ,25
+ ,9
+ ,14116
+ ,3988
+ ,13155
+ ,39
+ ,39
+ ,46698
+ ,45
+ ,274
+ ,35
+ ,0
+ ,14
+ ,13
+ ,52
+ ,33
+ ,10288
+ ,3030
+ ,18274
+ ,14
+ ,14
+ ,131698
+ ,65
+ ,502
+ ,74
+ ,0
+ ,60
+ ,19
+ ,75
+ ,59
+ ,65622
+ ,13520
+ ,72875
+ ,78
+ ,78
+ ,91735
+ ,35
+ ,382
+ ,103
+ ,0
+ ,7
+ ,18
+ ,71
+ ,49
+ ,16563
+ ,1421
+ ,10112
+ ,15
+ ,14
+ ,244749
+ ,95
+ ,964
+ ,144
+ ,2
+ ,98
+ ,33
+ ,124
+ ,115
+ ,76643
+ ,20923
+ ,142775
+ ,106
+ ,101
+ ,184510
+ ,49
+ ,537
+ ,60
+ ,7
+ ,64
+ ,40
+ ,151
+ ,140
+ ,110681
+ ,20237
+ ,68847
+ ,83
+ ,82
+ ,79863
+ ,37
+ ,438
+ ,134
+ ,1
+ ,29
+ ,22
+ ,71
+ ,49
+ ,29011
+ ,3219
+ ,17659
+ ,24
+ ,24
+ ,128423
+ ,64
+ ,369
+ ,89
+ ,8
+ ,32
+ ,38
+ ,145
+ ,120
+ ,92696
+ ,3769
+ ,20112
+ ,37
+ ,36
+ ,97839
+ ,38
+ ,417
+ ,42
+ ,2
+ ,25
+ ,24
+ ,87
+ ,66
+ ,94785
+ ,12252
+ ,61023
+ ,77
+ ,75
+ ,38214
+ ,34
+ ,276
+ ,52
+ ,0
+ ,16
+ ,8
+ ,27
+ ,21
+ ,8773
+ ,1888
+ ,13983
+ ,16
+ ,16
+ ,151101
+ ,32
+ ,514
+ ,98
+ ,2
+ ,48
+ ,35
+ ,131
+ ,124
+ ,83209
+ ,14497
+ ,65176
+ ,56
+ ,55
+ ,272458
+ ,65
+ ,822
+ ,99
+ ,0
+ ,100
+ ,43
+ ,162
+ ,152
+ ,93815
+ ,28864
+ ,132432
+ ,132
+ ,131
+ ,172494
+ ,52
+ ,389
+ ,52
+ ,0
+ ,46
+ ,43
+ ,165
+ ,139
+ ,86687
+ ,21721
+ ,112494
+ ,144
+ ,131
+ ,108043
+ ,62
+ ,466
+ ,29
+ ,1
+ ,45
+ ,14
+ ,54
+ ,38
+ ,34553
+ ,4821
+ ,45109
+ ,40
+ ,39
+ ,328107
+ ,65
+ ,1255
+ ,125
+ ,3
+ ,129
+ ,41
+ ,159
+ ,144
+ ,105547
+ ,33644
+ ,170875
+ ,153
+ ,144
+ ,250579
+ ,83
+ ,694
+ ,106
+ ,0
+ ,130
+ ,38
+ ,147
+ ,120
+ ,103487
+ ,15923
+ ,180759
+ ,143
+ ,139
+ ,351067
+ ,95
+ ,1024
+ ,95
+ ,3
+ ,136
+ ,45
+ ,170
+ ,160
+ ,213688
+ ,42935
+ ,214921
+ ,220
+ ,211
+ ,158015
+ ,29
+ ,400
+ ,40
+ ,0
+ ,59
+ ,31
+ ,119
+ ,114
+ ,71220
+ ,18864
+ ,100226
+ ,79
+ ,78
+ ,98866
+ ,18
+ ,397
+ ,140
+ ,0
+ ,25
+ ,13
+ ,49
+ ,39
+ ,23517
+ ,4977
+ ,32043
+ ,50
+ ,50
+ ,85439
+ ,33
+ ,350
+ ,43
+ ,0
+ ,32
+ ,28
+ ,104
+ ,78
+ ,56926
+ ,7785
+ ,54454
+ ,39
+ ,39
+ ,229242
+ ,247
+ ,719
+ ,128
+ ,4
+ ,63
+ ,31
+ ,120
+ ,119
+ ,91721
+ ,17939
+ ,78876
+ ,95
+ ,90
+ ,351619
+ ,139
+ ,1277
+ ,142
+ ,4
+ ,95
+ ,40
+ ,150
+ ,141
+ ,115168
+ ,23436
+ ,170745
+ ,169
+ ,166
+ ,84207
+ ,29
+ ,356
+ ,73
+ ,11
+ ,14
+ ,30
+ ,112
+ ,101
+ ,111194
+ ,325
+ ,6940
+ ,12
+ ,12
+ ,120445
+ ,118
+ ,457
+ ,72
+ ,0
+ ,36
+ ,16
+ ,59
+ ,56
+ ,51009
+ ,13539
+ ,49025
+ ,63
+ ,57
+ ,324598
+ ,110
+ ,1402
+ ,128
+ ,0
+ ,113
+ ,37
+ ,136
+ ,133
+ ,135777
+ ,34538
+ ,122037
+ ,134
+ ,133
+ ,131069
+ ,67
+ ,600
+ ,61
+ ,4
+ ,47
+ ,30
+ ,107
+ ,83
+ ,51513
+ ,12198
+ ,53782
+ ,69
+ ,69
+ ,204271
+ ,42
+ ,480
+ ,73
+ ,0
+ ,92
+ ,35
+ ,130
+ ,116
+ ,74163
+ ,26924
+ ,127748
+ ,119
+ ,119
+ ,165543
+ ,65
+ ,595
+ ,148
+ ,1
+ ,70
+ ,32
+ ,115
+ ,90
+ ,51633
+ ,12716
+ ,86839
+ ,119
+ ,119
+ ,141722
+ ,94
+ ,436
+ ,64
+ ,0
+ ,19
+ ,27
+ ,107
+ ,36
+ ,75345
+ ,8172
+ ,44830
+ ,75
+ ,65
+ ,116048
+ ,64
+ ,230
+ ,45
+ ,0
+ ,50
+ ,20
+ ,75
+ ,50
+ ,33416
+ ,10855
+ ,77395
+ ,63
+ ,61
+ ,250047
+ ,81
+ ,651
+ ,58
+ ,0
+ ,41
+ ,18
+ ,71
+ ,61
+ ,83305
+ ,11932
+ ,89324
+ ,55
+ ,49
+ ,299775
+ ,95
+ ,1367
+ ,97
+ ,9
+ ,91
+ ,31
+ ,120
+ ,97
+ ,98952
+ ,14300
+ ,103300
+ ,103
+ ,101
+ ,195838
+ ,67
+ ,564
+ ,50
+ ,1
+ ,111
+ ,31
+ ,116
+ ,98
+ ,102372
+ ,25515
+ ,112283
+ ,197
+ ,196
+ ,173260
+ ,63
+ ,716
+ ,37
+ ,3
+ ,41
+ ,21
+ ,79
+ ,78
+ ,37238
+ ,2805
+ ,10901
+ ,16
+ ,15
+ ,254488
+ ,83
+ ,747
+ ,50
+ ,10
+ ,120
+ ,39
+ ,150
+ ,117
+ ,103772
+ ,29402
+ ,120691
+ ,140
+ ,136
+ ,104389
+ ,45
+ ,467
+ ,105
+ ,5
+ ,135
+ ,41
+ ,156
+ ,148
+ ,123969
+ ,16440
+ ,58106
+ ,89
+ ,89
+ ,136084
+ ,30
+ ,671
+ ,69
+ ,0
+ ,27
+ ,13
+ ,51
+ ,41
+ ,27142
+ ,11221
+ ,57140
+ ,40
+ ,40
+ ,199476
+ ,70
+ ,861
+ ,46
+ ,2
+ ,87
+ ,32
+ ,118
+ ,105
+ ,135400
+ ,28732
+ ,122422
+ ,125
+ ,123
+ ,92499
+ ,32
+ ,319
+ ,57
+ ,0
+ ,25
+ ,18
+ ,71
+ ,55
+ ,21399
+ ,5250
+ ,25899
+ ,21
+ ,21
+ ,224330
+ ,83
+ ,612
+ ,52
+ ,1
+ ,131
+ ,39
+ ,144
+ ,132
+ ,130115
+ ,28608
+ ,139296
+ ,167
+ ,163
+ ,135781
+ ,31
+ ,433
+ ,98
+ ,2
+ ,45
+ ,14
+ ,47
+ ,44
+ ,24874
+ ,8092
+ ,52678
+ ,32
+ ,29
+ ,74408
+ ,67
+ ,434
+ ,61
+ ,4
+ ,29
+ ,7
+ ,28
+ ,21
+ ,34988
+ ,4473
+ ,23853
+ ,36
+ ,35
+ ,81240
+ ,66
+ ,503
+ ,89
+ ,0
+ ,58
+ ,17
+ ,68
+ ,50
+ ,45549
+ ,1572
+ ,17306
+ ,13
+ ,13
+ ,14688
+ ,10
+ ,85
+ ,0
+ ,0
+ ,4
+ ,0
+ ,0
+ ,0
+ ,6023
+ ,2065
+ ,7953
+ ,5
+ ,5
+ ,181633
+ ,70
+ ,564
+ ,48
+ ,2
+ ,47
+ ,30
+ ,110
+ ,73
+ ,64466
+ ,14817
+ ,89455
+ ,96
+ ,96
+ ,271856
+ ,103
+ ,824
+ ,91
+ ,1
+ ,109
+ ,37
+ ,147
+ ,86
+ ,54990
+ ,16714
+ ,147866
+ ,151
+ ,151
+ ,7199
+ ,5
+ ,74
+ ,0
+ ,0
+ ,7
+ ,0
+ ,0
+ ,0
+ ,1644
+ ,556
+ ,4245
+ ,6
+ ,6
+ ,46660
+ ,20
+ ,259
+ ,7
+ ,0
+ ,12
+ ,5
+ ,15
+ ,13
+ ,6179
+ ,2089
+ ,21509
+ ,13
+ ,13
+ ,17547
+ ,5
+ ,69
+ ,3
+ ,0
+ ,0
+ ,1
+ ,4
+ ,4
+ ,3926
+ ,2658
+ ,7670
+ ,3
+ ,3
+ ,133368
+ ,36
+ ,535
+ ,54
+ ,1
+ ,37
+ ,16
+ ,64
+ ,57
+ ,32755
+ ,10695
+ ,66675
+ ,57
+ ,56
+ ,95227
+ ,34
+ ,239
+ ,70
+ ,0
+ ,37
+ ,32
+ ,111
+ ,48
+ ,34777
+ ,1669
+ ,14336
+ ,23
+ ,23
+ ,152601
+ ,48
+ ,438
+ ,36
+ ,2
+ ,46
+ ,24
+ ,85
+ ,46
+ ,73224
+ ,16267
+ ,53608
+ ,61
+ ,57
+ ,98146
+ ,40
+ ,459
+ ,37
+ ,0
+ ,15
+ ,17
+ ,68
+ ,48
+ ,27114
+ ,7768
+ ,30059
+ ,21
+ ,14
+ ,79619
+ ,43
+ ,426
+ ,123
+ ,3
+ ,42
+ ,11
+ ,40
+ ,32
+ ,20760
+ ,7252
+ ,29668
+ ,43
+ ,43
+ ,59194
+ ,31
+ ,288
+ ,247
+ ,6
+ ,7
+ ,24
+ ,80
+ ,68
+ ,37636
+ ,6387
+ ,22097
+ ,20
+ ,20
+ ,139942
+ ,42
+ ,498
+ ,46
+ ,0
+ ,54
+ ,22
+ ,88
+ ,87
+ ,65461
+ ,18715
+ ,96841
+ ,82
+ ,72
+ ,118612
+ ,46
+ ,454
+ ,72
+ ,2
+ ,54
+ ,12
+ ,48
+ ,43
+ ,30080
+ ,7936
+ ,41907
+ ,90
+ ,87
+ ,72880
+ ,33
+ ,376
+ ,41
+ ,0
+ ,14
+ ,19
+ ,76
+ ,67
+ ,24094
+ ,8643
+ ,27080
+ ,25
+ ,21
+ ,65475
+ ,18
+ ,225
+ ,24
+ ,2
+ ,16
+ ,13
+ ,51
+ ,46
+ ,69008
+ ,7294
+ ,35885
+ ,60
+ ,56
+ ,99643
+ ,55
+ ,555
+ ,45
+ ,1
+ ,33
+ ,17
+ ,67
+ ,46
+ ,54968
+ ,4570
+ ,41247
+ ,61
+ ,59
+ ,71965
+ ,35
+ ,252
+ ,33
+ ,1
+ ,32
+ ,15
+ ,59
+ ,56
+ ,46090
+ ,7185
+ ,28313
+ ,85
+ ,82
+ ,77272
+ ,59
+ ,208
+ ,27
+ ,2
+ ,21
+ ,16
+ ,61
+ ,48
+ ,27507
+ ,10058
+ ,36845
+ ,43
+ ,43
+ ,49289
+ ,19
+ ,130
+ ,36
+ ,1
+ ,15
+ ,24
+ ,76
+ ,44
+ ,10672
+ ,2342
+ ,16548
+ ,25
+ ,25
+ ,135131
+ ,66
+ ,481
+ ,87
+ ,0
+ ,38
+ ,15
+ ,60
+ ,60
+ ,34029
+ ,8509
+ ,36134
+ ,41
+ ,38
+ ,108446
+ ,60
+ ,389
+ ,90
+ ,1
+ ,22
+ ,17
+ ,68
+ ,65
+ ,46300
+ ,13275
+ ,55764
+ ,26
+ ,25
+ ,89746
+ ,36
+ ,565
+ ,114
+ ,3
+ ,28
+ ,18
+ ,71
+ ,55
+ ,24760
+ ,6816
+ ,28910
+ ,38
+ ,38
+ ,44296
+ ,25
+ ,173
+ ,31
+ ,0
+ ,10
+ ,20
+ ,76
+ ,38
+ ,18779
+ ,1930
+ ,13339
+ ,12
+ ,12
+ ,77648
+ ,47
+ ,278
+ ,45
+ ,0
+ ,31
+ ,16
+ ,62
+ ,52
+ ,21280
+ ,8086
+ ,25319
+ ,29
+ ,29
+ ,181528
+ ,54
+ ,609
+ ,69
+ ,0
+ ,32
+ ,16
+ ,61
+ ,60
+ ,40662
+ ,10737
+ ,66956
+ ,49
+ ,47
+ ,134019
+ ,53
+ ,422
+ ,51
+ ,0
+ ,32
+ ,18
+ ,67
+ ,54
+ ,28987
+ ,8033
+ ,47487
+ ,46
+ ,45
+ ,124064
+ ,40
+ ,445
+ ,34
+ ,1
+ ,43
+ ,22
+ ,88
+ ,86
+ ,22827
+ ,7058
+ ,52785
+ ,41
+ ,40
+ ,92630
+ ,40
+ ,387
+ ,60
+ ,4
+ ,27
+ ,8
+ ,30
+ ,24
+ ,18513
+ ,6782
+ ,44683
+ ,31
+ ,30
+ ,121848
+ ,39
+ ,339
+ ,45
+ ,0
+ ,37
+ ,17
+ ,64
+ ,52
+ ,30594
+ ,5401
+ ,35619
+ ,41
+ ,41
+ ,52915
+ ,14
+ ,181
+ ,54
+ ,0
+ ,20
+ ,18
+ ,68
+ ,49
+ ,24006
+ ,6521
+ ,21920
+ ,26
+ ,25
+ ,81872
+ ,45
+ ,245
+ ,25
+ ,0
+ ,32
+ ,16
+ ,64
+ ,61
+ ,27913
+ ,10856
+ ,45608
+ ,23
+ ,23
+ ,58981
+ ,36
+ ,384
+ ,38
+ ,7
+ ,0
+ ,23
+ ,91
+ ,61
+ ,42744
+ ,2154
+ ,7721
+ ,14
+ ,14
+ ,53515
+ ,28
+ ,212
+ ,52
+ ,2
+ ,5
+ ,22
+ ,88
+ ,81
+ ,12934
+ ,6117
+ ,20634
+ ,16
+ ,16
+ ,60812
+ ,44
+ ,399
+ ,67
+ ,0
+ ,26
+ ,13
+ ,52
+ ,43
+ ,22574
+ ,5238
+ ,29788
+ ,25
+ ,26
+ ,56375
+ ,30
+ ,229
+ ,74
+ ,7
+ ,10
+ ,13
+ ,49
+ ,40
+ ,41385
+ ,4820
+ ,31931
+ ,21
+ ,21
+ ,65490
+ ,22
+ ,224
+ ,38
+ ,3
+ ,27
+ ,16
+ ,62
+ ,40
+ ,18653
+ ,5615
+ ,37754
+ ,32
+ ,27
+ ,80949
+ ,17
+ ,203
+ ,30
+ ,0
+ ,11
+ ,16
+ ,61
+ ,56
+ ,18472
+ ,4272
+ ,32505
+ ,9
+ ,9
+ ,76302
+ ,31
+ ,333
+ ,26
+ ,0
+ ,29
+ ,20
+ ,76
+ ,68
+ ,30976
+ ,8702
+ ,40557
+ ,35
+ ,33
+ ,104011
+ ,55
+ ,384
+ ,67
+ ,6
+ ,25
+ ,22
+ ,88
+ ,79
+ ,63339
+ ,15340
+ ,94238
+ ,42
+ ,42
+ ,98104
+ ,54
+ ,636
+ ,132
+ ,2
+ ,55
+ ,17
+ ,66
+ ,47
+ ,25568
+ ,8030
+ ,44197
+ ,68
+ ,68
+ ,67989
+ ,21
+ ,185
+ ,42
+ ,0
+ ,23
+ ,18
+ ,71
+ ,57
+ ,33747
+ ,9526
+ ,43228
+ ,32
+ ,32
+ ,30989
+ ,14
+ ,93
+ ,35
+ ,0
+ ,5
+ ,17
+ ,68
+ ,41
+ ,4154
+ ,1278
+ ,4103
+ ,6
+ ,6
+ ,135458
+ ,81
+ ,581
+ ,118
+ ,3
+ ,43
+ ,12
+ ,48
+ ,29
+ ,19474
+ ,4236
+ ,44144
+ ,68
+ ,67
+ ,73504
+ ,35
+ ,248
+ ,68
+ ,0
+ ,23
+ ,7
+ ,25
+ ,3
+ ,35130
+ ,3023
+ ,32868
+ ,33
+ ,33
+ ,63123
+ ,43
+ ,304
+ ,43
+ ,1
+ ,34
+ ,17
+ ,68
+ ,60
+ ,39067
+ ,7196
+ ,27640
+ ,84
+ ,77
+ ,61254
+ ,46
+ ,344
+ ,76
+ ,1
+ ,36
+ ,14
+ ,41
+ ,30
+ ,13310
+ ,3394
+ ,14063
+ ,46
+ ,46
+ ,74914
+ ,30
+ ,407
+ ,64
+ ,0
+ ,35
+ ,23
+ ,90
+ ,79
+ ,65892
+ ,6371
+ ,28990
+ ,30
+ ,30
+ ,31774
+ ,23
+ ,170
+ ,48
+ ,1
+ ,0
+ ,17
+ ,66
+ ,47
+ ,4143
+ ,1574
+ ,4694
+ ,0
+ ,0
+ ,81437
+ ,38
+ ,312
+ ,64
+ ,0
+ ,37
+ ,14
+ ,54
+ ,40
+ ,28579
+ ,9620
+ ,42648
+ ,36
+ ,36
+ ,87186
+ ,54
+ ,507
+ ,56
+ ,0
+ ,28
+ ,15
+ ,59
+ ,48
+ ,51776
+ ,6978
+ ,64329
+ ,47
+ ,46
+ ,50090
+ ,20
+ ,224
+ ,71
+ ,0
+ ,16
+ ,17
+ ,60
+ ,36
+ ,21152
+ ,4911
+ ,21928
+ ,20
+ ,18
+ ,65745
+ ,53
+ ,340
+ ,75
+ ,0
+ ,26
+ ,21
+ ,77
+ ,42
+ ,38084
+ ,8645
+ ,25836
+ ,50
+ ,48
+ ,56653
+ ,45
+ ,168
+ ,39
+ ,0
+ ,38
+ ,18
+ ,68
+ ,49
+ ,27717
+ ,8987
+ ,22779
+ ,30
+ ,29
+ ,158399
+ ,39
+ ,443
+ ,42
+ ,0
+ ,23
+ ,18
+ ,72
+ ,57
+ ,32928
+ ,5544
+ ,40820
+ ,30
+ ,28
+ ,46455
+ ,20
+ ,204
+ ,39
+ ,0
+ ,22
+ ,17
+ ,67
+ ,12
+ ,11342
+ ,3083
+ ,27530
+ ,34
+ ,34
+ ,73624
+ ,24
+ ,367
+ ,93
+ ,0
+ ,30
+ ,17
+ ,64
+ ,40
+ ,19499
+ ,6909
+ ,32378
+ ,33
+ ,33
+ ,38395
+ ,31
+ ,210
+ ,38
+ ,0
+ ,16
+ ,16
+ ,63
+ ,43
+ ,16380
+ ,3189
+ ,10824
+ ,34
+ ,34
+ ,91899
+ ,35
+ ,335
+ ,60
+ ,0
+ ,18
+ ,15
+ ,59
+ ,33
+ ,36874
+ ,6745
+ ,39613
+ ,37
+ ,33
+ ,139526
+ ,151
+ ,364
+ ,71
+ ,0
+ ,28
+ ,21
+ ,84
+ ,77
+ ,48259
+ ,16724
+ ,60865
+ ,83
+ ,80
+ ,52164
+ ,52
+ ,178
+ ,52
+ ,0
+ ,32
+ ,16
+ ,64
+ ,43
+ ,16734
+ ,4850
+ ,19787
+ ,32
+ ,32
+ ,51567
+ ,30
+ ,206
+ ,27
+ ,2
+ ,21
+ ,14
+ ,56
+ ,45
+ ,28207
+ ,7025
+ ,20107
+ ,30
+ ,30
+ ,70551
+ ,31
+ ,279
+ ,59
+ ,0
+ ,23
+ ,15
+ ,54
+ ,47
+ ,30143
+ ,6047
+ ,36605
+ ,43
+ ,41
+ ,84856
+ ,29
+ ,387
+ ,40
+ ,1
+ ,29
+ ,17
+ ,67
+ ,43
+ ,41369
+ ,7377
+ ,40961
+ ,41
+ ,41
+ ,102538
+ ,57
+ ,490
+ ,79
+ ,1
+ ,50
+ ,15
+ ,58
+ ,45
+ ,45833
+ ,9078
+ ,48231
+ ,51
+ ,51
+ ,86678
+ ,40
+ ,238
+ ,44
+ ,0
+ ,12
+ ,15
+ ,59
+ ,50
+ ,29156
+ ,4605
+ ,39725
+ ,19
+ ,18
+ ,85709
+ ,44
+ ,343
+ ,65
+ ,0
+ ,21
+ ,10
+ ,40
+ ,35
+ ,35944
+ ,3238
+ ,21455
+ ,37
+ ,34
+ ,34662
+ ,25
+ ,232
+ ,10
+ ,0
+ ,18
+ ,6
+ ,22
+ ,7
+ ,36278
+ ,8100
+ ,23430
+ ,33
+ ,31
+ ,150580
+ ,77
+ ,530
+ ,124
+ ,0
+ ,27
+ ,22
+ ,83
+ ,71
+ ,45588
+ ,9653
+ ,62991
+ ,41
+ ,39
+ ,99611
+ ,35
+ ,291
+ ,81
+ ,0
+ ,41
+ ,21
+ ,81
+ ,67
+ ,45097
+ ,8914
+ ,49363
+ ,54
+ ,54
+ ,19349
+ ,11
+ ,67
+ ,15
+ ,0
+ ,13
+ ,1
+ ,2
+ ,0
+ ,3895
+ ,786
+ ,9604
+ ,14
+ ,14
+ ,99373
+ ,63
+ ,397
+ ,92
+ ,1
+ ,12
+ ,18
+ ,72
+ ,62
+ ,28394
+ ,6700
+ ,24552
+ ,25
+ ,24
+ ,86230
+ ,44
+ ,467
+ ,42
+ ,0
+ ,21
+ ,17
+ ,61
+ ,54
+ ,18632
+ ,5788
+ ,31493
+ ,25
+ ,24
+ ,30837
+ ,19
+ ,178
+ ,10
+ ,0
+ ,8
+ ,4
+ ,15
+ ,4
+ ,2325
+ ,593
+ ,3439
+ ,8
+ ,8
+ ,31706
+ ,13
+ ,175
+ ,24
+ ,0
+ ,26
+ ,10
+ ,32
+ ,25
+ ,25139
+ ,4506
+ ,19555
+ ,26
+ ,26
+ ,89806
+ ,42
+ ,299
+ ,64
+ ,0
+ ,27
+ ,16
+ ,62
+ ,40
+ ,27975
+ ,6382
+ ,21228
+ ,20
+ ,19
+ ,62088
+ ,38
+ ,154
+ ,45
+ ,1
+ ,13
+ ,16
+ ,58
+ ,38
+ ,14483
+ ,5621
+ ,23177
+ ,11
+ ,11
+ ,40151
+ ,29
+ ,106
+ ,22
+ ,0
+ ,16
+ ,9
+ ,36
+ ,19
+ ,13127
+ ,3997
+ ,22094
+ ,14
+ ,14
+ ,27634
+ ,20
+ ,189
+ ,56
+ ,0
+ ,2
+ ,16
+ ,59
+ ,17
+ ,5839
+ ,520
+ ,2342
+ ,3
+ ,1
+ ,76990
+ ,27
+ ,194
+ ,94
+ ,0
+ ,42
+ ,17
+ ,68
+ ,67
+ ,24069
+ ,8891
+ ,38798
+ ,40
+ ,39
+ ,37460
+ ,20
+ ,135
+ ,19
+ ,0
+ ,5
+ ,7
+ ,21
+ ,14
+ ,3738
+ ,999
+ ,3255
+ ,5
+ ,5
+ ,54157
+ ,19
+ ,201
+ ,35
+ ,0
+ ,37
+ ,15
+ ,55
+ ,30
+ ,18625
+ ,7067
+ ,24261
+ ,38
+ ,37
+ ,49862
+ ,37
+ ,207
+ ,32
+ ,0
+ ,17
+ ,14
+ ,54
+ ,54
+ ,36341
+ ,4639
+ ,18511
+ ,32
+ ,32
+ ,84337
+ ,26
+ ,280
+ ,35
+ ,0
+ ,38
+ ,14
+ ,55
+ ,35
+ ,24548
+ ,5654
+ ,40798
+ ,41
+ ,38
+ ,64175
+ ,42
+ ,260
+ ,48
+ ,0
+ ,37
+ ,18
+ ,72
+ ,59
+ ,21792
+ ,6928
+ ,28893
+ ,46
+ ,47
+ ,59382
+ ,49
+ ,227
+ ,49
+ ,0
+ ,29
+ ,12
+ ,41
+ ,24
+ ,26263
+ ,1514
+ ,21425
+ ,47
+ ,47
+ ,119308
+ ,30
+ ,239
+ ,48
+ ,0
+ ,32
+ ,16
+ ,61
+ ,58
+ ,23686
+ ,9238
+ ,50276
+ ,37
+ ,37
+ ,76702
+ ,49
+ ,333
+ ,62
+ ,0
+ ,35
+ ,21
+ ,67
+ ,42
+ ,49303
+ ,8204
+ ,37643
+ ,51
+ ,51
+ ,103425
+ ,67
+ ,428
+ ,96
+ ,1
+ ,17
+ ,19
+ ,76
+ ,46
+ ,25659
+ ,5926
+ ,30377
+ ,49
+ ,45
+ ,70344
+ ,28
+ ,230
+ ,45
+ ,0
+ ,20
+ ,16
+ ,64
+ ,61
+ ,28904
+ ,5785
+ ,27126
+ ,21
+ ,21
+ ,43410
+ ,19
+ ,292
+ ,63
+ ,0
+ ,7
+ ,1
+ ,3
+ ,3
+ ,2781
+ ,4
+ ,13
+ ,1
+ ,1
+ ,104838
+ ,49
+ ,350
+ ,71
+ ,1
+ ,46
+ ,16
+ ,63
+ ,52
+ ,29236
+ ,5930
+ ,42097
+ ,44
+ ,42
+ ,62215
+ ,27
+ ,186
+ ,26
+ ,0
+ ,24
+ ,10
+ ,40
+ ,25
+ ,19546
+ ,3710
+ ,24451
+ ,26
+ ,26
+ ,69304
+ ,30
+ ,326
+ ,48
+ ,6
+ ,40
+ ,19
+ ,69
+ ,40
+ ,22818
+ ,705
+ ,14335
+ ,21
+ ,21
+ ,53117
+ ,22
+ ,155
+ ,29
+ ,3
+ ,3
+ ,12
+ ,48
+ ,32
+ ,32689
+ ,443
+ ,5084
+ ,4
+ ,4
+ ,19764
+ ,12
+ ,75
+ ,19
+ ,1
+ ,10
+ ,2
+ ,8
+ ,4
+ ,5752
+ ,2416
+ ,9927
+ ,10
+ ,10
+ ,86680
+ ,31
+ ,361
+ ,45
+ ,2
+ ,37
+ ,14
+ ,52
+ ,49
+ ,22197
+ ,7747
+ ,43527
+ ,43
+ ,43
+ ,84105
+ ,20
+ ,261
+ ,45
+ ,0
+ ,17
+ ,17
+ ,66
+ ,63
+ ,20055
+ ,5432
+ ,27184
+ ,34
+ ,34
+ ,77945
+ ,20
+ ,299
+ ,67
+ ,0
+ ,28
+ ,19
+ ,76
+ ,67
+ ,25272
+ ,4913
+ ,21610
+ ,32
+ ,31
+ ,89113
+ ,39
+ ,300
+ ,30
+ ,0
+ ,19
+ ,14
+ ,43
+ ,32
+ ,82206
+ ,2650
+ ,20484
+ ,20
+ ,19
+ ,91005
+ ,29
+ ,450
+ ,36
+ ,3
+ ,29
+ ,11
+ ,39
+ ,23
+ ,32073
+ ,2370
+ ,20156
+ ,34
+ ,34
+ ,40248
+ ,16
+ ,183
+ ,34
+ ,1
+ ,8
+ ,4
+ ,14
+ ,7
+ ,5444
+ ,775
+ ,6012
+ ,6
+ ,6
+ ,64187
+ ,27
+ ,238
+ ,36
+ ,0
+ ,10
+ ,16
+ ,61
+ ,54
+ ,20154
+ ,5576
+ ,18475
+ ,12
+ ,11
+ ,50857
+ ,21
+ ,165
+ ,34
+ ,0
+ ,15
+ ,20
+ ,71
+ ,37
+ ,36944
+ ,1352
+ ,12645
+ ,24
+ ,24
+ ,56613
+ ,19
+ ,234
+ ,37
+ ,1
+ ,15
+ ,12
+ ,44
+ ,35
+ ,8019
+ ,3080
+ ,11017
+ ,16
+ ,16
+ ,62792
+ ,35
+ ,176
+ ,46
+ ,0
+ ,28
+ ,15
+ ,60
+ ,51
+ ,30884
+ ,10205
+ ,37623
+ ,72
+ ,72
+ ,72535
+ ,14
+ ,329
+ ,44
+ ,0
+ ,17
+ ,16
+ ,64
+ ,39
+ ,19540
+ ,6095
+ ,35873
+ ,27
+ ,21)
+ ,dim=c(14
+ ,289)
+ ,dimnames=list(c('time_in_rfc'
+ ,'logins'
+ ,'compendium_views_info'
+ ,'compendium_views_pr'
+ ,'shared_compendiums'
+ ,'blogged_computations'
+ ,'compendiums_reviewed'
+ ,'feedback_messages_p1'
+ ,'feedback_messages_p120'
+ ,'totsize'
+ ,'totrevisions'
+ ,'totseconds'
+ ,'tothyperlinks'
+ ,'totblogs')
+ ,1:289))
> y <- array(NA,dim=c(14,289),dimnames=list(c('time_in_rfc','logins','compendium_views_info','compendium_views_pr','shared_compendiums','blogged_computations','compendiums_reviewed','feedback_messages_p1','feedback_messages_p120','totsize','totrevisions','totseconds','tothyperlinks','totblogs'),1:289))
> for (i in 1:dim(x)[1])
+ {
+ for (j in 1:dim(x)[2])
+ {
+ y[i,j] <- as.numeric(x[i,j])
+ }
+ }
> par3 = 'No Linear Trend'
> par2 = 'Do not include Seasonal Dummies'
> par1 = '1'
> par3 <- 'No Linear Trend'
> par2 <- 'Do not include Seasonal Dummies'
> par1 <- '1'
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Prof. Dr. P. Wessa
> #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #Technical description: Write here your technical program description (don't use hard returns!)
> library(lattice)
> library(lmtest)
Loading required package: zoo
Attaching package: 'zoo'
The following object(s) are masked from 'package:base':
as.Date, as.Date.numeric
> n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test
> par1 <- as.numeric(par1)
> x <- t(y)
> k <- length(x[1,])
> n <- length(x[,1])
> x1 <- cbind(x[,par1], x[,1:k!=par1])
> mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1])
> colnames(x1) <- mycolnames #colnames(x)[par1]
> x <- x1
> if (par3 == 'First Differences'){
+ x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep='')))
+ for (i in 1:n-1) {
+ for (j in 1:k) {
+ x2[i,j] <- x[i+1,j] - x[i,j]
+ }
+ }
+ x <- x2
+ }
> if (par2 == 'Include Monthly Dummies'){
+ x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep ='')))
+ for (i in 1:11){
+ x2[seq(i,n,12),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> if (par2 == 'Include Quarterly Dummies'){
+ x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep ='')))
+ for (i in 1:3){
+ x2[seq(i,n,4),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> k <- length(x[1,])
> if (par3 == 'Linear Trend'){
+ x <- cbind(x, c(1:n))
+ colnames(x)[k+1] <- 't'
+ }
> x
time_in_rfc logins compendium_views_info compendium_views_pr
1 210907 56 396 81
2 120982 56 297 55
3 176508 54 559 50
4 179321 89 967 125
5 123185 40 270 40
6 52746 25 143 37
7 385534 92 1562 63
8 33170 18 109 44
9 101645 63 371 88
10 149061 44 656 66
11 165446 33 511 57
12 237213 84 655 74
13 173326 88 465 49
14 133131 55 525 52
15 258873 60 885 88
16 180083 66 497 36
17 324799 154 1436 108
18 230964 53 612 43
19 236785 119 865 75
20 135473 41 385 32
21 202925 61 567 44
22 215147 58 639 85
23 344297 75 963 86
24 153935 33 398 56
25 132943 40 410 50
26 174724 92 966 135
27 174415 100 801 63
28 225548 112 892 81
29 223632 73 513 52
30 124817 40 469 44
31 221698 45 683 113
32 210767 60 643 39
33 170266 62 535 73
34 260561 75 625 48
35 84853 31 264 33
36 294424 77 992 59
37 101011 34 238 41
38 215641 46 818 69
39 325107 99 937 64
40 7176 17 70 1
41 167542 66 507 59
42 106408 30 260 32
43 96560 76 503 129
44 265769 146 927 37
45 269651 67 1269 31
46 149112 56 537 65
47 175824 107 910 107
48 152871 58 532 74
49 111665 34 345 54
50 116408 61 918 76
51 362301 119 1635 715
52 78800 42 330 57
53 183167 66 557 66
54 277965 89 1178 106
55 150629 44 740 54
56 168809 66 452 32
57 24188 24 218 20
58 329267 259 764 71
59 65029 17 255 21
60 101097 64 454 70
61 218946 41 866 112
62 244052 68 574 66
63 341570 168 1276 190
64 103597 43 379 66
65 233328 132 825 165
66 256462 105 798 56
67 206161 71 663 61
68 311473 112 1069 53
69 235800 94 921 127
70 177939 82 858 63
71 207176 70 711 38
72 196553 57 503 50
73 174184 53 382 52
74 143246 103 464 42
75 187559 121 717 76
76 187681 62 690 67
77 119016 52 462 50
78 182192 52 657 53
79 73566 32 385 39
80 194979 62 577 50
81 167488 45 619 77
82 143756 46 479 57
83 275541 63 817 73
84 243199 75 752 34
85 182999 88 430 39
86 135649 46 451 46
87 152299 53 537 63
88 120221 37 519 35
89 346485 90 1000 106
90 145790 63 637 43
91 193339 78 465 47
92 80953 25 437 31
93 122774 45 711 162
94 130585 46 299 57
95 112611 41 248 36
96 286468 144 1162 263
97 241066 82 714 78
98 148446 91 905 63
99 204713 71 649 54
100 182079 63 512 63
101 140344 53 472 77
102 220516 62 905 79
103 243060 63 786 110
104 162765 32 489 56
105 182613 39 479 56
106 232138 62 617 43
107 265318 117 925 111
108 85574 34 351 71
109 310839 92 1144 62
110 225060 93 669 56
111 232317 54 707 74
112 144966 144 458 60
113 43287 14 214 43
114 155754 61 599 68
115 164709 109 572 53
116 201940 38 897 87
117 235454 73 819 46
118 220801 75 720 105
119 99466 50 273 32
120 92661 61 508 133
121 133328 55 506 79
122 61361 77 451 51
123 125930 75 699 207
124 100750 72 407 67
125 224549 50 465 47
126 82316 32 245 34
127 102010 53 370 66
128 101523 42 316 76
129 243511 71 603 65
130 22938 10 154 9
131 41566 35 229 42
132 152474 65 577 45
133 61857 25 192 25
134 99923 66 617 115
135 132487 41 411 97
136 317394 86 975 53
137 21054 16 146 2
138 209641 42 705 52
139 22648 19 184 44
140 31414 19 200 22
141 46698 45 274 35
142 131698 65 502 74
143 91735 35 382 103
144 244749 95 964 144
145 184510 49 537 60
146 79863 37 438 134
147 128423 64 369 89
148 97839 38 417 42
149 38214 34 276 52
150 151101 32 514 98
151 272458 65 822 99
152 172494 52 389 52
153 108043 62 466 29
154 328107 65 1255 125
155 250579 83 694 106
156 351067 95 1024 95
157 158015 29 400 40
158 98866 18 397 140
159 85439 33 350 43
160 229242 247 719 128
161 351619 139 1277 142
162 84207 29 356 73
163 120445 118 457 72
164 324598 110 1402 128
165 131069 67 600 61
166 204271 42 480 73
167 165543 65 595 148
168 141722 94 436 64
169 116048 64 230 45
170 250047 81 651 58
171 299775 95 1367 97
172 195838 67 564 50
173 173260 63 716 37
174 254488 83 747 50
175 104389 45 467 105
176 136084 30 671 69
177 199476 70 861 46
178 92499 32 319 57
179 224330 83 612 52
180 135781 31 433 98
181 74408 67 434 61
182 81240 66 503 89
183 14688 10 85 0
184 181633 70 564 48
185 271856 103 824 91
186 7199 5 74 0
187 46660 20 259 7
188 17547 5 69 3
189 133368 36 535 54
190 95227 34 239 70
191 152601 48 438 36
192 98146 40 459 37
193 79619 43 426 123
194 59194 31 288 247
195 139942 42 498 46
196 118612 46 454 72
197 72880 33 376 41
198 65475 18 225 24
199 99643 55 555 45
200 71965 35 252 33
201 77272 59 208 27
202 49289 19 130 36
203 135131 66 481 87
204 108446 60 389 90
205 89746 36 565 114
206 44296 25 173 31
207 77648 47 278 45
208 181528 54 609 69
209 134019 53 422 51
210 124064 40 445 34
211 92630 40 387 60
212 121848 39 339 45
213 52915 14 181 54
214 81872 45 245 25
215 58981 36 384 38
216 53515 28 212 52
217 60812 44 399 67
218 56375 30 229 74
219 65490 22 224 38
220 80949 17 203 30
221 76302 31 333 26
222 104011 55 384 67
223 98104 54 636 132
224 67989 21 185 42
225 30989 14 93 35
226 135458 81 581 118
227 73504 35 248 68
228 63123 43 304 43
229 61254 46 344 76
230 74914 30 407 64
231 31774 23 170 48
232 81437 38 312 64
233 87186 54 507 56
234 50090 20 224 71
235 65745 53 340 75
236 56653 45 168 39
237 158399 39 443 42
238 46455 20 204 39
239 73624 24 367 93
240 38395 31 210 38
241 91899 35 335 60
242 139526 151 364 71
243 52164 52 178 52
244 51567 30 206 27
245 70551 31 279 59
246 84856 29 387 40
247 102538 57 490 79
248 86678 40 238 44
249 85709 44 343 65
250 34662 25 232 10
251 150580 77 530 124
252 99611 35 291 81
253 19349 11 67 15
254 99373 63 397 92
255 86230 44 467 42
256 30837 19 178 10
257 31706 13 175 24
258 89806 42 299 64
259 62088 38 154 45
260 40151 29 106 22
261 27634 20 189 56
262 76990 27 194 94
263 37460 20 135 19
264 54157 19 201 35
265 49862 37 207 32
266 84337 26 280 35
267 64175 42 260 48
268 59382 49 227 49
269 119308 30 239 48
270 76702 49 333 62
271 103425 67 428 96
272 70344 28 230 45
273 43410 19 292 63
274 104838 49 350 71
275 62215 27 186 26
276 69304 30 326 48
277 53117 22 155 29
278 19764 12 75 19
279 86680 31 361 45
280 84105 20 261 45
281 77945 20 299 67
282 89113 39 300 30
283 91005 29 450 36
284 40248 16 183 34
285 64187 27 238 36
286 50857 21 165 34
287 56613 19 234 37
288 62792 35 176 46
289 72535 14 329 44
shared_compendiums blogged_computations compendiums_reviewed
1 3 79 30
2 4 58 28
3 12 60 38
4 2 108 30
5 1 49 22
6 3 0 26
7 0 121 25
8 0 1 18
9 0 20 11
10 5 43 26
11 0 69 25
12 0 78 38
13 7 86 44
14 7 44 30
15 3 104 40
16 9 63 34
17 0 158 47
18 4 102 30
19 3 77 31
20 0 82 23
21 7 115 36
22 0 101 36
23 1 80 30
24 5 50 25
25 7 83 39
26 0 123 34
27 0 73 31
28 5 81 31
29 0 105 33
30 0 47 25
31 0 105 33
32 3 94 35
33 4 44 42
34 1 114 43
35 4 38 30
36 2 107 33
37 0 30 13
38 0 71 32
39 0 84 36
40 0 0 0
41 2 59 28
42 1 33 14
43 0 42 17
44 2 96 32
45 10 106 30
46 6 56 35
47 0 57 20
48 5 59 28
49 4 39 28
50 1 34 39
51 2 76 34
52 2 20 26
53 0 91 39
54 8 115 39
55 3 85 33
56 0 76 28
57 0 8 4
58 8 79 39
59 5 21 18
60 3 30 14
61 1 76 29
62 5 101 44
63 1 94 21
64 1 27 16
65 5 92 28
66 0 123 35
67 12 75 28
68 8 128 38
69 8 105 23
70 8 55 36
71 8 56 32
72 2 41 29
73 0 72 25
74 5 67 27
75 8 75 36
76 2 114 28
77 5 118 23
78 12 77 40
79 6 22 23
80 7 66 40
81 2 69 28
82 0 105 34
83 4 116 33
84 3 88 28
85 6 73 34
86 2 99 30
87 0 62 33
88 1 53 22
89 0 118 38
90 5 30 26
91 2 100 35
92 0 49 8
93 0 24 24
94 5 67 29
95 0 46 20
96 1 57 29
97 0 75 45
98 1 135 37
99 1 68 33
100 2 124 33
101 6 33 25
102 1 98 32
103 4 58 29
104 2 68 28
105 3 81 28
106 0 131 31
107 10 110 52
108 0 37 21
109 9 130 24
110 7 93 41
111 0 118 33
112 0 39 32
113 4 13 19
114 4 74 20
115 0 81 31
116 0 109 31
117 0 151 32
118 1 51 18
119 0 28 23
120 1 40 17
121 0 56 20
122 0 27 12
123 4 37 17
124 0 83 30
125 4 54 31
126 4 27 10
127 3 28 13
128 0 59 22
129 0 133 42
130 0 12 1
131 5 0 9
132 0 106 32
133 4 23 11
134 0 44 25
135 0 71 36
136 1 116 31
137 0 4 0
138 5 62 24
139 0 12 13
140 0 18 8
141 0 14 13
142 0 60 19
143 0 7 18
144 2 98 33
145 7 64 40
146 1 29 22
147 8 32 38
148 2 25 24
149 0 16 8
150 2 48 35
151 0 100 43
152 0 46 43
153 1 45 14
154 3 129 41
155 0 130 38
156 3 136 45
157 0 59 31
158 0 25 13
159 0 32 28
160 4 63 31
161 4 95 40
162 11 14 30
163 0 36 16
164 0 113 37
165 4 47 30
166 0 92 35
167 1 70 32
168 0 19 27
169 0 50 20
170 0 41 18
171 9 91 31
172 1 111 31
173 3 41 21
174 10 120 39
175 5 135 41
176 0 27 13
177 2 87 32
178 0 25 18
179 1 131 39
180 2 45 14
181 4 29 7
182 0 58 17
183 0 4 0
184 2 47 30
185 1 109 37
186 0 7 0
187 0 12 5
188 0 0 1
189 1 37 16
190 0 37 32
191 2 46 24
192 0 15 17
193 3 42 11
194 6 7 24
195 0 54 22
196 2 54 12
197 0 14 19
198 2 16 13
199 1 33 17
200 1 32 15
201 2 21 16
202 1 15 24
203 0 38 15
204 1 22 17
205 3 28 18
206 0 10 20
207 0 31 16
208 0 32 16
209 0 32 18
210 1 43 22
211 4 27 8
212 0 37 17
213 0 20 18
214 0 32 16
215 7 0 23
216 2 5 22
217 0 26 13
218 7 10 13
219 3 27 16
220 0 11 16
221 0 29 20
222 6 25 22
223 2 55 17
224 0 23 18
225 0 5 17
226 3 43 12
227 0 23 7
228 1 34 17
229 1 36 14
230 0 35 23
231 1 0 17
232 0 37 14
233 0 28 15
234 0 16 17
235 0 26 21
236 0 38 18
237 0 23 18
238 0 22 17
239 0 30 17
240 0 16 16
241 0 18 15
242 0 28 21
243 0 32 16
244 2 21 14
245 0 23 15
246 1 29 17
247 1 50 15
248 0 12 15
249 0 21 10
250 0 18 6
251 0 27 22
252 0 41 21
253 0 13 1
254 1 12 18
255 0 21 17
256 0 8 4
257 0 26 10
258 0 27 16
259 1 13 16
260 0 16 9
261 0 2 16
262 0 42 17
263 0 5 7
264 0 37 15
265 0 17 14
266 0 38 14
267 0 37 18
268 0 29 12
269 0 32 16
270 0 35 21
271 1 17 19
272 0 20 16
273 0 7 1
274 1 46 16
275 0 24 10
276 6 40 19
277 3 3 12
278 1 10 2
279 2 37 14
280 0 17 17
281 0 28 19
282 0 19 14
283 3 29 11
284 1 8 4
285 0 10 16
286 0 15 20
287 1 15 12
288 0 28 15
289 0 17 16
feedback_messages_p1 feedback_messages_p120 totsize totrevisions totseconds
1 115 94 112285 24188 146283
2 109 103 84786 18273 98364
3 146 93 83123 14130 86146
4 116 103 101193 32287 96933
5 68 51 38361 8654 79234
6 101 70 68504 9245 42551
7 96 91 119182 33251 195663
8 67 22 22807 1271 6853
9 44 38 17140 5279 21529
10 100 93 116174 27101 95757
11 93 60 57635 16373 85584
12 140 123 66198 19716 143983
13 166 148 71701 17753 75851
14 99 90 57793 9028 59238
15 139 124 80444 18653 93163
16 130 70 53855 8828 96037
17 181 168 97668 29498 151511
18 116 115 133824 27563 136368
19 116 71 101481 18293 112642
20 88 66 99645 22530 94728
21 139 134 114789 15977 105499
22 135 117 99052 35082 121527
23 108 108 67654 16116 127766
24 89 84 65553 15849 98958
25 156 156 97500 16026 77900
26 129 120 69112 26569 85646
27 118 114 82753 24785 98579
28 118 94 85323 17569 130767
29 125 120 72654 23825 131741
30 95 81 30727 7869 53907
31 126 110 77873 14975 178812
32 135 133 117478 37791 146761
33 154 122 74007 9605 82036
34 165 158 90183 27295 163253
35 113 109 61542 2746 27032
36 127 124 101494 34461 171975
37 52 39 27570 8098 65990
38 121 92 55813 4787 86572
39 136 126 79215 24919 159676
40 0 0 1423 603 1929
41 108 70 55461 16329 85371
42 46 37 31081 12558 58391
43 54 38 22996 7784 31580
44 124 120 83122 28522 136815
45 115 93 70106 22265 120642
46 128 95 60578 14459 69107
47 80 77 39992 14526 50495
48 97 90 79892 22240 108016
49 104 80 49810 11802 46341
50 59 31 71570 7623 78348
51 125 110 100708 11912 79336
52 82 66 33032 7935 56968
53 149 138 82875 18220 93176
54 149 133 139077 19199 161632
55 122 113 71595 19918 87850
56 118 100 72260 21884 127969
57 12 7 5950 2694 15049
58 144 140 115762 15808 155135
59 67 61 32551 3597 25109
60 52 41 31701 5296 45824
61 108 96 80670 25239 102996
62 166 164 143558 29801 160604
63 80 78 117105 18450 158051
64 60 49 23789 7132 44547
65 107 102 120733 34861 162647
66 127 124 105195 35940 174141
67 107 99 73107 16688 60622
68 146 129 132068 24683 179566
69 84 62 149193 46230 184301
70 141 73 46821 10387 75661
71 123 114 87011 21436 96144
72 111 99 95260 30546 129847
73 98 70 55183 19746 117286
74 105 104 106671 15977 71180
75 135 116 73511 22583 109377
76 107 91 92945 17274 85298
77 85 74 78664 16469 73631
78 155 138 70054 14251 86767
79 88 67 22618 3007 23824
80 155 151 74011 16851 93487
81 104 72 83737 21113 82981
82 132 120 69094 17401 73815
83 127 115 93133 23958 94552
84 108 105 95536 23567 132190
85 129 104 225920 13065 128754
86 116 108 62133 15358 66363
87 122 98 61370 14587 67808
88 85 69 43836 12770 61724
89 147 111 106117 24021 131722
90 99 99 38692 9648 68580
91 87 71 84651 20537 106175
92 28 27 56622 7905 55792
93 90 69 15986 4527 25157
94 109 107 95364 30495 76669
95 78 73 26706 7117 57283
96 111 107 89691 17719 105805
97 158 93 67267 27056 129484
98 141 129 126846 33473 72413
99 122 69 41140 9758 87831
100 124 118 102860 21115 96971
101 93 73 51715 7236 71299
102 124 119 55801 13790 77494
103 112 104 111813 32902 120336
104 108 107 120293 25131 93913
105 99 99 138599 30910 136048
106 117 90 161647 35947 181248
107 199 197 115929 29848 146123
108 78 36 24266 6943 32036
109 91 85 162901 42705 186646
110 158 139 109825 31808 102255
111 126 106 129838 26675 168237
112 122 50 37510 8435 64219
113 71 64 43750 7409 19630
114 75 31 40652 14993 76825
115 115 63 87771 36867 115338
116 119 92 85872 33835 109427
117 124 106 89275 24164 118168
118 72 63 44418 12607 84845
119 91 69 192565 22609 153197
120 45 41 35232 5892 29877
121 78 56 40909 17014 63506
122 39 25 13294 5394 22445
123 68 65 32387 9178 47695
124 119 93 140867 6440 68370
125 117 114 120662 21916 146304
126 39 38 21233 4011 38233
127 50 44 44332 5818 42071
128 88 87 61056 18647 50517
129 155 110 101338 20556 103950
130 0 0 1168 238 5841
131 36 27 13497 70 2341
132 123 83 65567 22392 84396
133 32 30 25162 3913 24610
134 99 80 32334 12237 35753
135 136 98 40735 8388 55515
136 117 82 91413 22120 209056
137 0 0 855 338 6622
138 88 60 97068 11727 115814
139 39 28 44339 3704 11609
140 25 9 14116 3988 13155
141 52 33 10288 3030 18274
142 75 59 65622 13520 72875
143 71 49 16563 1421 10112
144 124 115 76643 20923 142775
145 151 140 110681 20237 68847
146 71 49 29011 3219 17659
147 145 120 92696 3769 20112
148 87 66 94785 12252 61023
149 27 21 8773 1888 13983
150 131 124 83209 14497 65176
151 162 152 93815 28864 132432
152 165 139 86687 21721 112494
153 54 38 34553 4821 45109
154 159 144 105547 33644 170875
155 147 120 103487 15923 180759
156 170 160 213688 42935 214921
157 119 114 71220 18864 100226
158 49 39 23517 4977 32043
159 104 78 56926 7785 54454
160 120 119 91721 17939 78876
161 150 141 115168 23436 170745
162 112 101 111194 325 6940
163 59 56 51009 13539 49025
164 136 133 135777 34538 122037
165 107 83 51513 12198 53782
166 130 116 74163 26924 127748
167 115 90 51633 12716 86839
168 107 36 75345 8172 44830
169 75 50 33416 10855 77395
170 71 61 83305 11932 89324
171 120 97 98952 14300 103300
172 116 98 102372 25515 112283
173 79 78 37238 2805 10901
174 150 117 103772 29402 120691
175 156 148 123969 16440 58106
176 51 41 27142 11221 57140
177 118 105 135400 28732 122422
178 71 55 21399 5250 25899
179 144 132 130115 28608 139296
180 47 44 24874 8092 52678
181 28 21 34988 4473 23853
182 68 50 45549 1572 17306
183 0 0 6023 2065 7953
184 110 73 64466 14817 89455
185 147 86 54990 16714 147866
186 0 0 1644 556 4245
187 15 13 6179 2089 21509
188 4 4 3926 2658 7670
189 64 57 32755 10695 66675
190 111 48 34777 1669 14336
191 85 46 73224 16267 53608
192 68 48 27114 7768 30059
193 40 32 20760 7252 29668
194 80 68 37636 6387 22097
195 88 87 65461 18715 96841
196 48 43 30080 7936 41907
197 76 67 24094 8643 27080
198 51 46 69008 7294 35885
199 67 46 54968 4570 41247
200 59 56 46090 7185 28313
201 61 48 27507 10058 36845
202 76 44 10672 2342 16548
203 60 60 34029 8509 36134
204 68 65 46300 13275 55764
205 71 55 24760 6816 28910
206 76 38 18779 1930 13339
207 62 52 21280 8086 25319
208 61 60 40662 10737 66956
209 67 54 28987 8033 47487
210 88 86 22827 7058 52785
211 30 24 18513 6782 44683
212 64 52 30594 5401 35619
213 68 49 24006 6521 21920
214 64 61 27913 10856 45608
215 91 61 42744 2154 7721
216 88 81 12934 6117 20634
217 52 43 22574 5238 29788
218 49 40 41385 4820 31931
219 62 40 18653 5615 37754
220 61 56 18472 4272 32505
221 76 68 30976 8702 40557
222 88 79 63339 15340 94238
223 66 47 25568 8030 44197
224 71 57 33747 9526 43228
225 68 41 4154 1278 4103
226 48 29 19474 4236 44144
227 25 3 35130 3023 32868
228 68 60 39067 7196 27640
229 41 30 13310 3394 14063
230 90 79 65892 6371 28990
231 66 47 4143 1574 4694
232 54 40 28579 9620 42648
233 59 48 51776 6978 64329
234 60 36 21152 4911 21928
235 77 42 38084 8645 25836
236 68 49 27717 8987 22779
237 72 57 32928 5544 40820
238 67 12 11342 3083 27530
239 64 40 19499 6909 32378
240 63 43 16380 3189 10824
241 59 33 36874 6745 39613
242 84 77 48259 16724 60865
243 64 43 16734 4850 19787
244 56 45 28207 7025 20107
245 54 47 30143 6047 36605
246 67 43 41369 7377 40961
247 58 45 45833 9078 48231
248 59 50 29156 4605 39725
249 40 35 35944 3238 21455
250 22 7 36278 8100 23430
251 83 71 45588 9653 62991
252 81 67 45097 8914 49363
253 2 0 3895 786 9604
254 72 62 28394 6700 24552
255 61 54 18632 5788 31493
256 15 4 2325 593 3439
257 32 25 25139 4506 19555
258 62 40 27975 6382 21228
259 58 38 14483 5621 23177
260 36 19 13127 3997 22094
261 59 17 5839 520 2342
262 68 67 24069 8891 38798
263 21 14 3738 999 3255
264 55 30 18625 7067 24261
265 54 54 36341 4639 18511
266 55 35 24548 5654 40798
267 72 59 21792 6928 28893
268 41 24 26263 1514 21425
269 61 58 23686 9238 50276
270 67 42 49303 8204 37643
271 76 46 25659 5926 30377
272 64 61 28904 5785 27126
273 3 3 2781 4 13
274 63 52 29236 5930 42097
275 40 25 19546 3710 24451
276 69 40 22818 705 14335
277 48 32 32689 443 5084
278 8 4 5752 2416 9927
279 52 49 22197 7747 43527
280 66 63 20055 5432 27184
281 76 67 25272 4913 21610
282 43 32 82206 2650 20484
283 39 23 32073 2370 20156
284 14 7 5444 775 6012
285 61 54 20154 5576 18475
286 71 37 36944 1352 12645
287 44 35 8019 3080 11017
288 60 51 30884 10205 37623
289 64 39 19540 6095 35873
tothyperlinks totblogs
1 144 145
2 103 101
3 98 98
4 135 132
5 61 60
6 39 38
7 150 144
8 5 5
9 28 28
10 84 84
11 80 79
12 130 127
13 82 78
14 60 60
15 131 131
16 84 84
17 140 133
18 151 150
19 91 91
20 138 132
21 150 136
22 124 124
23 119 118
24 73 70
25 110 107
26 123 119
27 90 89
28 116 112
29 113 108
30 56 52
31 115 112
32 119 116
33 129 123
34 127 125
35 27 27
36 175 162
37 35 32
38 64 64
39 96 92
40 0 0
41 84 83
42 41 41
43 47 47
44 126 120
45 105 105
46 80 79
47 70 65
48 73 70
49 57 55
50 40 39
51 68 67
52 21 21
53 127 127
54 154 152
55 116 113
56 102 99
57 7 7
58 148 141
59 21 21
60 35 35
61 112 109
62 137 133
63 135 123
64 26 26
65 230 230
66 181 166
67 71 68
68 147 147
69 190 179
70 64 61
71 105 101
72 107 108
73 94 90
74 116 114
75 106 103
76 143 142
77 81 79
78 89 88
79 26 25
80 84 83
81 113 113
82 120 118
83 110 110
84 134 129
85 54 51
86 96 93
87 78 76
88 51 49
89 121 118
90 38 38
91 145 141
92 59 58
93 27 27
94 91 91
95 48 48
96 68 63
97 58 56
98 150 144
99 74 73
100 181 168
101 65 64
102 97 97
103 121 117
104 99 100
105 152 149
106 188 187
107 138 127
108 40 37
109 254 245
110 87 87
111 178 177
112 51 49
113 49 49
114 73 73
115 176 177
116 94 94
117 120 117
118 66 60
119 56 55
120 39 39
121 66 64
122 27 26
123 65 64
124 58 58
125 98 95
126 25 25
127 26 26
128 77 76
129 130 129
130 11 11
131 2 2
132 101 101
133 31 28
134 36 36
135 120 89
136 195 193
137 4 4
138 89 84
139 24 23
140 39 39
141 14 14
142 78 78
143 15 14
144 106 101
145 83 82
146 24 24
147 37 36
148 77 75
149 16 16
150 56 55
151 132 131
152 144 131
153 40 39
154 153 144
155 143 139
156 220 211
157 79 78
158 50 50
159 39 39
160 95 90
161 169 166
162 12 12
163 63 57
164 134 133
165 69 69
166 119 119
167 119 119
168 75 65
169 63 61
170 55 49
171 103 101
172 197 196
173 16 15
174 140 136
175 89 89
176 40 40
177 125 123
178 21 21
179 167 163
180 32 29
181 36 35
182 13 13
183 5 5
184 96 96
185 151 151
186 6 6
187 13 13
188 3 3
189 57 56
190 23 23
191 61 57
192 21 14
193 43 43
194 20 20
195 82 72
196 90 87
197 25 21
198 60 56
199 61 59
200 85 82
201 43 43
202 25 25
203 41 38
204 26 25
205 38 38
206 12 12
207 29 29
208 49 47
209 46 45
210 41 40
211 31 30
212 41 41
213 26 25
214 23 23
215 14 14
216 16 16
217 25 26
218 21 21
219 32 27
220 9 9
221 35 33
222 42 42
223 68 68
224 32 32
225 6 6
226 68 67
227 33 33
228 84 77
229 46 46
230 30 30
231 0 0
232 36 36
233 47 46
234 20 18
235 50 48
236 30 29
237 30 28
238 34 34
239 33 33
240 34 34
241 37 33
242 83 80
243 32 32
244 30 30
245 43 41
246 41 41
247 51 51
248 19 18
249 37 34
250 33 31
251 41 39
252 54 54
253 14 14
254 25 24
255 25 24
256 8 8
257 26 26
258 20 19
259 11 11
260 14 14
261 3 1
262 40 39
263 5 5
264 38 37
265 32 32
266 41 38
267 46 47
268 47 47
269 37 37
270 51 51
271 49 45
272 21 21
273 1 1
274 44 42
275 26 26
276 21 21
277 4 4
278 10 10
279 43 43
280 34 34
281 32 31
282 20 19
283 34 34
284 6 6
285 12 11
286 24 24
287 16 16
288 72 72
289 27 21
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) logins compendium_views_info
-4.826e+03 1.924e+02 1.227e+02
compendium_views_pr shared_compendiums blogged_computations
2.390e+01 -1.811e+02 4.993e+01
compendiums_reviewed feedback_messages_p1 feedback_messages_p120
-9.581e+02 3.873e+02 1.889e+02
totsize totrevisions totseconds
-9.793e-02 -8.406e-01 8.623e-01
tothyperlinks totblogs
1.235e+02 -1.131e+02
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-62155 -11318 710 10219 85846
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -4.826e+03 3.414e+03 -1.414 0.15860
logins 1.924e+02 5.470e+01 3.517 0.00051 ***
compendium_views_info 1.227e+02 9.270e+00 13.235 < 2e-16 ***
compendium_views_pr 2.390e+01 3.073e+01 0.778 0.43738
shared_compendiums -1.811e+02 5.468e+02 -0.331 0.74068
blogged_computations 4.993e+01 8.966e+01 0.557 0.57805
compendiums_reviewed -9.581e+02 6.878e+02 -1.393 0.16474
feedback_messages_p1 3.873e+02 2.082e+02 1.860 0.06396 .
feedback_messages_p120 1.889e+02 1.038e+02 1.820 0.06982 .
totsize -9.793e-02 6.461e-02 -1.516 0.13074
totrevisions -8.406e-01 3.162e-01 -2.658 0.00832 **
totseconds 8.623e-01 7.050e-02 12.231 < 2e-16 ***
tothyperlinks 1.235e+02 4.355e+02 0.284 0.77700
totblogs -1.131e+02 4.513e+02 -0.251 0.80230
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 21520 on 275 degrees of freedom
Multiple R-squared: 0.9348, Adjusted R-squared: 0.9317
F-statistic: 303.2 on 13 and 275 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.7667839 4.664322e-01 2.332161e-01
[2,] 0.6349092 7.301817e-01 3.650908e-01
[3,] 0.5892813 8.214375e-01 4.107187e-01
[4,] 0.4889577 9.779154e-01 5.110423e-01
[5,] 0.4839913 9.679826e-01 5.160087e-01
[6,] 0.5786982 8.426036e-01 4.213018e-01
[7,] 0.8957680 2.084640e-01 1.042320e-01
[8,] 0.8516667 2.966666e-01 1.483333e-01
[9,] 0.9220175 1.559649e-01 7.798247e-02
[10,] 0.9509139 9.817223e-02 4.908612e-02
[11,] 0.9485322 1.029357e-01 5.146783e-02
[12,] 0.9489899 1.020202e-01 5.101009e-02
[13,] 0.9406649 1.186701e-01 5.933507e-02
[14,] 0.9245430 1.509139e-01 7.545695e-02
[15,] 0.9793799 4.124025e-02 2.062013e-02
[16,] 0.9705139 5.897223e-02 2.948612e-02
[17,] 0.9583115 8.337691e-02 4.168846e-02
[18,] 0.9422974 1.154052e-01 5.770261e-02
[19,] 0.9231254 1.537491e-01 7.687457e-02
[20,] 0.9055201 1.889597e-01 9.447985e-02
[21,] 0.8784042 2.431915e-01 1.215958e-01
[22,] 0.8464252 3.071495e-01 1.535748e-01
[23,] 0.9305808 1.388384e-01 6.941921e-02
[24,] 0.9424826 1.150349e-01 5.751744e-02
[25,] 0.9280082 1.439837e-01 7.199185e-02
[26,] 0.9112272 1.775455e-01 8.877276e-02
[27,] 0.8909376 2.181248e-01 1.090624e-01
[28,] 0.8677070 2.645859e-01 1.322930e-01
[29,] 0.8832127 2.335745e-01 1.167873e-01
[30,] 0.8555791 2.888418e-01 1.444209e-01
[31,] 0.8431325 3.137350e-01 1.568675e-01
[32,] 0.8175200 3.649600e-01 1.824800e-01
[33,] 0.7983199 4.033602e-01 2.016801e-01
[34,] 0.9193071 1.613858e-01 8.069290e-02
[35,] 0.9861993 2.760134e-02 1.380067e-02
[36,] 0.9870142 2.597151e-02 1.298575e-02
[37,] 0.9835936 3.281286e-02 1.640643e-02
[38,] 0.9931406 1.371880e-02 6.859398e-03
[39,] 0.9968147 6.370626e-03 3.185313e-03
[40,] 0.9971143 5.771386e-03 2.885693e-03
[41,] 0.9968451 6.309885e-03 3.154943e-03
[42,] 0.9970892 5.821543e-03 2.910771e-03
[43,] 0.9959695 8.060933e-03 4.030467e-03
[44,] 0.9951746 9.650876e-03 4.825438e-03
[45,] 0.9938811 1.223782e-02 6.118911e-03
[46,] 0.9922587 1.548260e-02 7.741301e-03
[47,] 0.9902958 1.940846e-02 9.704231e-03
[48,] 0.9871818 2.563643e-02 1.281821e-02
[49,] 0.9947045 1.059103e-02 5.295516e-03
[50,] 0.9937816 1.243684e-02 6.218421e-03
[51,] 0.9981800 3.639910e-03 1.819955e-03
[52,] 0.9976838 4.632399e-03 2.316199e-03
[53,] 0.9976270 4.746036e-03 2.373018e-03
[54,] 0.9978098 4.380335e-03 2.190168e-03
[55,] 0.9972640 5.471914e-03 2.735957e-03
[56,] 0.9966160 6.767975e-03 3.383988e-03
[57,] 0.9957444 8.511132e-03 4.255566e-03
[58,] 0.9943340 1.133198e-02 5.665990e-03
[59,] 0.9960093 7.981413e-03 3.990706e-03
[60,] 0.9950748 9.850419e-03 4.925210e-03
[61,] 0.9941130 1.177402e-02 5.887011e-03
[62,] 0.9930492 1.390164e-02 6.950818e-03
[63,] 0.9925289 1.494223e-02 7.471114e-03
[64,] 0.9902626 1.947470e-02 9.737351e-03
[65,] 0.9886929 2.261415e-02 1.130707e-02
[66,] 0.9858015 2.839707e-02 1.419853e-02
[67,] 0.9989609 2.078161e-03 1.039080e-03
[68,] 0.9987340 2.531991e-03 1.265995e-03
[69,] 0.9982800 3.440087e-03 1.720044e-03
[70,] 0.9977040 4.592044e-03 2.296022e-03
[71,] 0.9969769 6.046186e-03 3.023093e-03
[72,] 0.9964039 7.192121e-03 3.596060e-03
[73,] 0.9999183 1.633259e-04 8.166297e-05
[74,] 0.9999298 1.403344e-04 7.016722e-05
[75,] 0.9999630 7.398051e-05 3.699026e-05
[76,] 0.9999632 7.368486e-05 3.684243e-05
[77,] 0.9999512 9.759962e-05 4.879981e-05
[78,] 0.9999468 1.064850e-04 5.324248e-05
[79,] 0.9999269 1.461218e-04 7.306091e-05
[80,] 0.9999108 1.784607e-04 8.923034e-05
[81,] 0.9998889 2.222180e-04 1.111090e-04
[82,] 0.9999782 4.357638e-05 2.178819e-05
[83,] 0.9999734 5.320193e-05 2.660096e-05
[84,] 0.9999642 7.151925e-05 3.575962e-05
[85,] 0.9999482 1.036834e-04 5.184171e-05
[86,] 0.9999261 1.478195e-04 7.390974e-05
[87,] 0.9999602 7.955131e-05 3.977565e-05
[88,] 0.9999501 9.981199e-05 4.990600e-05
[89,] 0.9999292 1.416457e-04 7.082285e-05
[90,] 0.9999059 1.882445e-04 9.412223e-05
[91,] 0.9999460 1.080818e-04 5.404092e-05
[92,] 0.9999224 1.551520e-04 7.757600e-05
[93,] 0.9999005 1.989266e-04 9.946331e-05
[94,] 0.9999021 1.957567e-04 9.787834e-05
[95,] 0.9998940 2.120724e-04 1.060362e-04
[96,] 0.9998617 2.765173e-04 1.382586e-04
[97,] 0.9998201 3.598962e-04 1.799481e-04
[98,] 0.9997503 4.993598e-04 2.496799e-04
[99,] 0.9997191 5.618759e-04 2.809380e-04
[100,] 0.9996449 7.101353e-04 3.550676e-04
[101,] 0.9995135 9.729205e-04 4.864602e-04
[102,] 0.9997597 4.805815e-04 2.402908e-04
[103,] 0.9999722 5.552612e-05 2.776306e-05
[104,] 0.9999641 7.188732e-05 3.594366e-05
[105,] 0.9999483 1.033617e-04 5.168086e-05
[106,] 0.9999704 5.917252e-05 2.958626e-05
[107,] 0.9999730 5.408896e-05 2.704448e-05
[108,] 0.9999861 2.779344e-05 1.389672e-05
[109,] 0.9999865 2.704352e-05 1.352176e-05
[110,] 0.9999814 3.712264e-05 1.856132e-05
[111,] 0.9999739 5.210535e-05 2.605268e-05
[112,] 0.9999633 7.340007e-05 3.670004e-05
[113,] 0.9999910 1.800998e-05 9.004989e-06
[114,] 0.9999868 2.637962e-05 1.318981e-05
[115,] 0.9999808 3.843941e-05 1.921971e-05
[116,] 0.9999768 4.631023e-05 2.315512e-05
[117,] 0.9999693 6.132770e-05 3.066385e-05
[118,] 0.9999825 3.501987e-05 1.750993e-05
[119,] 0.9999766 4.679976e-05 2.339988e-05
[120,] 0.9999692 6.164087e-05 3.082043e-05
[121,] 0.9999557 8.852969e-05 4.426484e-05
[122,] 0.9999422 1.156214e-04 5.781071e-05
[123,] 0.9999334 1.331192e-04 6.655961e-05
[124,] 0.9999088 1.823384e-04 9.116921e-05
[125,] 0.9999008 1.983348e-04 9.916739e-05
[126,] 0.9998663 2.673773e-04 1.336886e-04
[127,] 0.9998610 2.779100e-04 1.389550e-04
[128,] 0.9999235 1.530303e-04 7.651514e-05
[129,] 0.9999499 1.002121e-04 5.010606e-05
[130,] 0.9999327 1.345473e-04 6.727363e-05
[131,] 0.9999542 9.157706e-05 4.578853e-05
[132,] 0.9999495 1.009860e-04 5.049298e-05
[133,] 0.9999452 1.096135e-04 5.480675e-05
[134,] 0.9999249 1.502737e-04 7.513686e-05
[135,] 0.9999374 1.251297e-04 6.256487e-05
[136,] 0.9999156 1.688330e-04 8.441650e-05
[137,] 0.9998902 2.196453e-04 1.098227e-04
[138,] 0.9998504 2.992738e-04 1.496369e-04
[139,] 0.9999205 1.589930e-04 7.949652e-05
[140,] 0.9999033 1.933186e-04 9.665930e-05
[141,] 0.9998635 2.729443e-04 1.364721e-04
[142,] 0.9998707 2.585687e-04 1.292843e-04
[143,] 0.9999054 1.892907e-04 9.464533e-05
[144,] 0.9998793 2.414076e-04 1.207038e-04
[145,] 0.9998348 3.303177e-04 1.651589e-04
[146,] 0.9998134 3.731964e-04 1.865982e-04
[147,] 0.9998118 3.763631e-04 1.881816e-04
[148,] 0.9998040 3.919178e-04 1.959589e-04
[149,] 0.9997992 4.016876e-04 2.008438e-04
[150,] 0.9997948 4.104205e-04 2.052103e-04
[151,] 0.9997352 5.295739e-04 2.647869e-04
[152,] 0.9997233 5.533838e-04 2.766919e-04
[153,] 0.9996413 7.174924e-04 3.587462e-04
[154,] 0.9999916 1.688024e-05 8.440120e-06
[155,] 0.9999884 2.315714e-05 1.157857e-05
[156,] 0.9999876 2.479861e-05 1.239930e-05
[157,] 0.9999971 5.848313e-06 2.924157e-06
[158,] 0.9999995 1.047213e-06 5.236063e-07
[159,] 0.9999996 7.870130e-07 3.935065e-07
[160,] 0.9999994 1.176356e-06 5.881782e-07
[161,] 0.9999994 1.273432e-06 6.367162e-07
[162,] 0.9999992 1.534558e-06 7.672792e-07
[163,] 0.9999988 2.317113e-06 1.158557e-06
[164,] 0.9999991 1.809290e-06 9.046452e-07
[165,] 0.9999987 2.638984e-06 1.319492e-06
[166,] 0.9999992 1.565006e-06 7.825030e-07
[167,] 0.9999987 2.512135e-06 1.256067e-06
[168,] 0.9999986 2.704025e-06 1.352012e-06
[169,] 0.9999984 3.198305e-06 1.599152e-06
[170,] 0.9999976 4.748236e-06 2.374118e-06
[171,] 0.9999968 6.314311e-06 3.157155e-06
[172,] 0.9999951 9.793713e-06 4.896856e-06
[173,] 0.9999931 1.386618e-05 6.933090e-06
[174,] 0.9999947 1.056725e-05 5.283627e-06
[175,] 0.9999999 2.656020e-07 1.328010e-07
[176,] 0.9999998 4.528517e-07 2.264259e-07
[177,] 0.9999996 7.471969e-07 3.735984e-07
[178,] 0.9999994 1.182239e-06 5.911193e-07
[179,] 0.9999994 1.235141e-06 6.175706e-07
[180,] 0.9999994 1.166619e-06 5.833097e-07
[181,] 0.9999993 1.435801e-06 7.179007e-07
[182,] 0.9999989 2.237351e-06 1.118676e-06
[183,] 0.9999989 2.150353e-06 1.075176e-06
[184,] 0.9999984 3.133049e-06 1.566524e-06
[185,] 0.9999978 4.418465e-06 2.209233e-06
[186,] 0.9999965 7.026223e-06 3.513112e-06
[187,] 0.9999962 7.676510e-06 3.838255e-06
[188,] 0.9999937 1.255340e-05 6.276699e-06
[189,] 0.9999909 1.819342e-05 9.096712e-06
[190,] 0.9999858 2.843947e-05 1.421973e-05
[191,] 0.9999778 4.433776e-05 2.216888e-05
[192,] 0.9999950 1.006416e-05 5.032081e-06
[193,] 0.9999972 5.687352e-06 2.843676e-06
[194,] 0.9999952 9.511580e-06 4.755790e-06
[195,] 0.9999928 1.435522e-05 7.177610e-06
[196,] 0.9999978 4.393610e-06 2.196805e-06
[197,] 0.9999963 7.319827e-06 3.659914e-06
[198,] 0.9999939 1.212615e-05 6.063074e-06
[199,] 0.9999903 1.945965e-05 9.729826e-06
[200,] 0.9999858 2.832426e-05 1.416213e-05
[201,] 0.9999939 1.219623e-05 6.098117e-06
[202,] 0.9999897 2.057357e-05 1.028678e-05
[203,] 0.9999834 3.319659e-05 1.659829e-05
[204,] 0.9999759 4.828694e-05 2.414347e-05
[205,] 0.9999704 5.929517e-05 2.964758e-05
[206,] 0.9999925 1.505194e-05 7.525972e-06
[207,] 0.9999939 1.222493e-05 6.112463e-06
[208,] 0.9999895 2.108392e-05 1.054196e-05
[209,] 0.9999818 3.641330e-05 1.820665e-05
[210,] 0.9999739 5.227525e-05 2.613762e-05
[211,] 0.9999632 7.369870e-05 3.684935e-05
[212,] 0.9999569 8.610214e-05 4.305107e-05
[213,] 0.9999275 1.449722e-04 7.248611e-05
[214,] 0.9999569 8.624007e-05 4.312004e-05
[215,] 0.9999583 8.330281e-05 4.165140e-05
[216,] 0.9999283 1.434773e-04 7.173867e-05
[217,] 0.9999997 6.456191e-07 3.228095e-07
[218,] 0.9999994 1.255483e-06 6.277417e-07
[219,] 0.9999988 2.491338e-06 1.245669e-06
[220,] 0.9999976 4.881911e-06 2.440955e-06
[221,] 1.0000000 1.335167e-08 6.675834e-09
[222,] 1.0000000 3.178479e-08 1.589240e-08
[223,] 1.0000000 7.033220e-08 3.516610e-08
[224,] 0.9999999 1.230386e-07 6.151930e-08
[225,] 0.9999999 2.666599e-07 1.333299e-07
[226,] 0.9999998 4.277507e-07 2.138753e-07
[227,] 0.9999996 8.366667e-07 4.183333e-07
[228,] 0.9999991 1.808693e-06 9.043467e-07
[229,] 0.9999987 2.641786e-06 1.320893e-06
[230,] 0.9999975 5.033529e-06 2.516764e-06
[231,] 0.9999977 4.537374e-06 2.268687e-06
[232,] 0.9999959 8.220721e-06 4.110361e-06
[233,] 0.9999926 1.485258e-05 7.426292e-06
[234,] 0.9999908 1.838725e-05 9.193624e-06
[235,] 0.9999845 3.092819e-05 1.546409e-05
[236,] 0.9999750 5.003463e-05 2.501732e-05
[237,] 0.9999446 1.108555e-04 5.542777e-05
[238,] 0.9998819 2.361580e-04 1.180790e-04
[239,] 0.9998907 2.186932e-04 1.093466e-04
[240,] 0.9997600 4.799648e-04 2.399824e-04
[241,] 0.9996815 6.370872e-04 3.185436e-04
[242,] 0.9997874 4.251645e-04 2.125822e-04
[243,] 0.9995431 9.138324e-04 4.569162e-04
[244,] 0.9991772 1.645660e-03 8.228302e-04
[245,] 0.9987628 2.474343e-03 1.237172e-03
[246,] 0.9973090 5.382085e-03 2.691043e-03
[247,] 0.9965743 6.851479e-03 3.425740e-03
[248,] 0.9970672 5.865539e-03 2.932770e-03
[249,] 0.9954792 9.041569e-03 4.520784e-03
[250,] 0.9906574 1.868522e-02 9.342608e-03
[251,] 0.9822134 3.557326e-02 1.778663e-02
[252,] 0.9686631 6.267370e-02 3.133685e-02
[253,] 0.9999392 1.215578e-04 6.077889e-05
[254,] 0.9999548 9.031436e-05 4.515718e-05
[255,] 0.9996148 7.703845e-04 3.851923e-04
[256,] 0.9981820 3.636010e-03 1.818005e-03
> postscript(file="/var/wessaorg/rcomp/tmp/1empx1353091884.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/2rfmd1353091884.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/3mvt51353091884.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/45vy51353091884.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/5qdf11353091884.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
21298.6550 -22181.4982 7359.3814 -43546.7872 11040.8391 -15280.2836
7 8 9 10 11 12
11725.2811 4539.3801 19311.5600 -18549.8468 17688.8664 -3977.4460
13 14 15 16 17 18
5120.9825 -4261.0696 39757.0455 6875.8629 -41988.2923 24215.5644
19 20 21 22 23 24
6100.3998 710.1252 6835.4833 16562.0500 85845.5017 7586.2831
25 26 27 28 29 30
-22449.9309 -52566.1802 -38067.9395 -30732.4510 17335.4532 -5025.3756
31 32 33 34 35 36
-47335.7030 -8645.6132 -5571.9448 4533.0930 -1500.7788 2678.3949
37 38 39 40 41 42
4470.0794 6657.3966 37792.5963 -898.4406 9808.4745 22765.9485
43 44 45 46 47 48
-10857.1302 -4417.9117 -11227.5702 -2413.8504 -12060.2031 -17999.6860
49 50 51 52 53 54
10567.6587 -52474.4635 37143.4506 -25690.3421 -3799.8365 -42586.5580
55 56 57 58 59 60
-39307.8527 -23110.2705 -15567.1212 30015.7503 -1229.0147 -11509.8189
61 62 63 64 65 66
9971.4910 5296.9777 10806.3679 3385.2601 -31476.0034 -16592.0433
67 68 69 70 71 72
47478.7699 -9795.3726 -28154.6978 -27917.5267 12631.2775 13960.4966
73 74 75 76 77 78
9325.9635 -5486.8923 -32249.7755 5665.2571 -16897.0447 -11318.0673
79 80 81 82 83 84
-16828.8326 3450.7093 8959.2775 -10580.7625 68493.5417 15062.5583
85 86 87 88 89 90
-1197.2187 -5534.6098 1065.6950 -12876.0611 76797.3749 -21040.1444
91 92 93 94 95 96
38827.4220 -20756.5940 -14858.8141 19353.2699 9835.4733 7829.1549
97 98 99 100 101 102
18111.6052 -50976.5154 18697.7657 6797.0358 -2312.1926 5164.8627
103 104 105 106 107 108
32430.8009 12868.8608 6443.3629 -3145.3865 -28879.6859 194.7747
109 110 111 112 113 114
10219.0056 25451.5954 -16738.8462 -8889.2839 -10093.2839 4597.8522
115 116 117 118 119 120
-15243.3883 -10365.6745 5303.1562 35870.0186 -62154.7766 -8219.6818
121 122 123 124 125 126
1819.2371 -28769.2914 -28857.7477 -39215.7262 25193.8727 9013.1120
127 128 129 130 131 132
6276.1051 2482.0355 48041.0277 2253.1194 384.7113 -13598.7697
133 134 135 136 137 138
13646.3440 -35845.8617 -5134.5911 -6892.6140 -742.6953 12047.4618
139 140 141 142 143 144
-11262.1119 -4108.9658 -18552.5637 -10180.4599 14468.1463 -32760.5293
145 146 147 148 149 150
31571.9291 -6242.3953 25377.1822 -14926.6128 -15941.9516 5148.1761
151 152 153 154 155 156
24089.5488 -6126.3172 -5659.0726 -3260.2994 -32204.9828 18913.7914
157 158 159 160 161 162
1841.1070 11375.9878 -25076.4781 8096.9822 868.3404 10665.9752
163 164 165 166 167 168
-2347.3877 20798.5843 -13857.7541 15569.1531 -13235.3157 23801.4107
169 170 171 172 173 174
2222.3731 74126.3810 9619.4674 10783.8362 46845.8243 32396.4975
175 176 177 178 179 180
-39240.9942 -2821.2590 -24430.2853 12736.9168 1656.0763 27334.1971
181 182 183 184 185 186
-11058.5201 -22137.9425 2377.4647 13651.7721 973.0693 -1460.0007
187 188 189 190 191 192
-4700.4073 7500.8093 -3551.0166 31478.3504 46578.1554 468.3511
193 194 195 196 197 198
-9483.9658 -13485.4028 -20987.0986 11825.5130 -14777.3598 3109.0055
199 200 201 202 203 204
-22299.0808 2796.8683 5465.3280 6815.8902 22631.3482 -4166.8058
205 206 207 208 209 210
-23063.5822 -3769.5000 4943.1791 32908.9598 23207.0064 -3314.6020
211 212 213 214 215 216
283.8547 32731.6549 753.9489 -2899.9777 -15409.7773 -14039.9061
217 218 219 220 221 222
-29705.1333 -7391.4721 -6255.7399 14737.1290 -16090.5277 -41143.6391
223 224 225 226 227 228
-38742.5708 -3387.8759 712.7598 3006.9550 12150.3207 -17168.5340
229 230 231 232 233 234
-4987.8028 -20487.1022 -10119.3559 -4023.2716 -45997.4772 -5802.8252
235 236 237 238 239 240
-14093.0493 1236.3612 49960.4051 -11943.7461 -11514.5374 -12808.9109
241 242 243 244 245 246
6131.7396 5690.5282 -6951.9362 -1741.3961 -6997.1318 -8877.8388
247 248 249 250 251 252
-14514.4075 7195.7562 11881.2251 -9428.5606 4386.7814 2896.0774
253 254 255 256 257 258
5629.9033 5271.6057 -15133.9913 4474.7394 -7728.3003 20229.5153
259 260 261 262 263 264
10897.6682 1607.7132 -8015.6267 1620.6585 15447.0399 1723.7035
265 266 267 268 269 270
-5976.6153 3787.5628 -13077.4974 278.4139 33393.1050 -7010.0034
271 272 273 274 275 276
306.9863 2596.0749 6382.5311 5456.2244 10252.4853 476.3422
277 278 279 280 281 282
20404.0745 4299.9040 -6583.8814 12700.5066 2879.2526 30972.0198
283 284 285 286 287 288
11765.7406 11538.2798 5255.3484 8113.9317 9175.3048 -3290.9784
289
-9276.6030
> postscript(file="/var/wessaorg/rcomp/tmp/68qo91353091884.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 21298.6550 NA
1 -22181.4982 21298.6550
2 7359.3814 -22181.4982
3 -43546.7872 7359.3814
4 11040.8391 -43546.7872
5 -15280.2836 11040.8391
6 11725.2811 -15280.2836
7 4539.3801 11725.2811
8 19311.5600 4539.3801
9 -18549.8468 19311.5600
10 17688.8664 -18549.8468
11 -3977.4460 17688.8664
12 5120.9825 -3977.4460
13 -4261.0696 5120.9825
14 39757.0455 -4261.0696
15 6875.8629 39757.0455
16 -41988.2923 6875.8629
17 24215.5644 -41988.2923
18 6100.3998 24215.5644
19 710.1252 6100.3998
20 6835.4833 710.1252
21 16562.0500 6835.4833
22 85845.5017 16562.0500
23 7586.2831 85845.5017
24 -22449.9309 7586.2831
25 -52566.1802 -22449.9309
26 -38067.9395 -52566.1802
27 -30732.4510 -38067.9395
28 17335.4532 -30732.4510
29 -5025.3756 17335.4532
30 -47335.7030 -5025.3756
31 -8645.6132 -47335.7030
32 -5571.9448 -8645.6132
33 4533.0930 -5571.9448
34 -1500.7788 4533.0930
35 2678.3949 -1500.7788
36 4470.0794 2678.3949
37 6657.3966 4470.0794
38 37792.5963 6657.3966
39 -898.4406 37792.5963
40 9808.4745 -898.4406
41 22765.9485 9808.4745
42 -10857.1302 22765.9485
43 -4417.9117 -10857.1302
44 -11227.5702 -4417.9117
45 -2413.8504 -11227.5702
46 -12060.2031 -2413.8504
47 -17999.6860 -12060.2031
48 10567.6587 -17999.6860
49 -52474.4635 10567.6587
50 37143.4506 -52474.4635
51 -25690.3421 37143.4506
52 -3799.8365 -25690.3421
53 -42586.5580 -3799.8365
54 -39307.8527 -42586.5580
55 -23110.2705 -39307.8527
56 -15567.1212 -23110.2705
57 30015.7503 -15567.1212
58 -1229.0147 30015.7503
59 -11509.8189 -1229.0147
60 9971.4910 -11509.8189
61 5296.9777 9971.4910
62 10806.3679 5296.9777
63 3385.2601 10806.3679
64 -31476.0034 3385.2601
65 -16592.0433 -31476.0034
66 47478.7699 -16592.0433
67 -9795.3726 47478.7699
68 -28154.6978 -9795.3726
69 -27917.5267 -28154.6978
70 12631.2775 -27917.5267
71 13960.4966 12631.2775
72 9325.9635 13960.4966
73 -5486.8923 9325.9635
74 -32249.7755 -5486.8923
75 5665.2571 -32249.7755
76 -16897.0447 5665.2571
77 -11318.0673 -16897.0447
78 -16828.8326 -11318.0673
79 3450.7093 -16828.8326
80 8959.2775 3450.7093
81 -10580.7625 8959.2775
82 68493.5417 -10580.7625
83 15062.5583 68493.5417
84 -1197.2187 15062.5583
85 -5534.6098 -1197.2187
86 1065.6950 -5534.6098
87 -12876.0611 1065.6950
88 76797.3749 -12876.0611
89 -21040.1444 76797.3749
90 38827.4220 -21040.1444
91 -20756.5940 38827.4220
92 -14858.8141 -20756.5940
93 19353.2699 -14858.8141
94 9835.4733 19353.2699
95 7829.1549 9835.4733
96 18111.6052 7829.1549
97 -50976.5154 18111.6052
98 18697.7657 -50976.5154
99 6797.0358 18697.7657
100 -2312.1926 6797.0358
101 5164.8627 -2312.1926
102 32430.8009 5164.8627
103 12868.8608 32430.8009
104 6443.3629 12868.8608
105 -3145.3865 6443.3629
106 -28879.6859 -3145.3865
107 194.7747 -28879.6859
108 10219.0056 194.7747
109 25451.5954 10219.0056
110 -16738.8462 25451.5954
111 -8889.2839 -16738.8462
112 -10093.2839 -8889.2839
113 4597.8522 -10093.2839
114 -15243.3883 4597.8522
115 -10365.6745 -15243.3883
116 5303.1562 -10365.6745
117 35870.0186 5303.1562
118 -62154.7766 35870.0186
119 -8219.6818 -62154.7766
120 1819.2371 -8219.6818
121 -28769.2914 1819.2371
122 -28857.7477 -28769.2914
123 -39215.7262 -28857.7477
124 25193.8727 -39215.7262
125 9013.1120 25193.8727
126 6276.1051 9013.1120
127 2482.0355 6276.1051
128 48041.0277 2482.0355
129 2253.1194 48041.0277
130 384.7113 2253.1194
131 -13598.7697 384.7113
132 13646.3440 -13598.7697
133 -35845.8617 13646.3440
134 -5134.5911 -35845.8617
135 -6892.6140 -5134.5911
136 -742.6953 -6892.6140
137 12047.4618 -742.6953
138 -11262.1119 12047.4618
139 -4108.9658 -11262.1119
140 -18552.5637 -4108.9658
141 -10180.4599 -18552.5637
142 14468.1463 -10180.4599
143 -32760.5293 14468.1463
144 31571.9291 -32760.5293
145 -6242.3953 31571.9291
146 25377.1822 -6242.3953
147 -14926.6128 25377.1822
148 -15941.9516 -14926.6128
149 5148.1761 -15941.9516
150 24089.5488 5148.1761
151 -6126.3172 24089.5488
152 -5659.0726 -6126.3172
153 -3260.2994 -5659.0726
154 -32204.9828 -3260.2994
155 18913.7914 -32204.9828
156 1841.1070 18913.7914
157 11375.9878 1841.1070
158 -25076.4781 11375.9878
159 8096.9822 -25076.4781
160 868.3404 8096.9822
161 10665.9752 868.3404
162 -2347.3877 10665.9752
163 20798.5843 -2347.3877
164 -13857.7541 20798.5843
165 15569.1531 -13857.7541
166 -13235.3157 15569.1531
167 23801.4107 -13235.3157
168 2222.3731 23801.4107
169 74126.3810 2222.3731
170 9619.4674 74126.3810
171 10783.8362 9619.4674
172 46845.8243 10783.8362
173 32396.4975 46845.8243
174 -39240.9942 32396.4975
175 -2821.2590 -39240.9942
176 -24430.2853 -2821.2590
177 12736.9168 -24430.2853
178 1656.0763 12736.9168
179 27334.1971 1656.0763
180 -11058.5201 27334.1971
181 -22137.9425 -11058.5201
182 2377.4647 -22137.9425
183 13651.7721 2377.4647
184 973.0693 13651.7721
185 -1460.0007 973.0693
186 -4700.4073 -1460.0007
187 7500.8093 -4700.4073
188 -3551.0166 7500.8093
189 31478.3504 -3551.0166
190 46578.1554 31478.3504
191 468.3511 46578.1554
192 -9483.9658 468.3511
193 -13485.4028 -9483.9658
194 -20987.0986 -13485.4028
195 11825.5130 -20987.0986
196 -14777.3598 11825.5130
197 3109.0055 -14777.3598
198 -22299.0808 3109.0055
199 2796.8683 -22299.0808
200 5465.3280 2796.8683
201 6815.8902 5465.3280
202 22631.3482 6815.8902
203 -4166.8058 22631.3482
204 -23063.5822 -4166.8058
205 -3769.5000 -23063.5822
206 4943.1791 -3769.5000
207 32908.9598 4943.1791
208 23207.0064 32908.9598
209 -3314.6020 23207.0064
210 283.8547 -3314.6020
211 32731.6549 283.8547
212 753.9489 32731.6549
213 -2899.9777 753.9489
214 -15409.7773 -2899.9777
215 -14039.9061 -15409.7773
216 -29705.1333 -14039.9061
217 -7391.4721 -29705.1333
218 -6255.7399 -7391.4721
219 14737.1290 -6255.7399
220 -16090.5277 14737.1290
221 -41143.6391 -16090.5277
222 -38742.5708 -41143.6391
223 -3387.8759 -38742.5708
224 712.7598 -3387.8759
225 3006.9550 712.7598
226 12150.3207 3006.9550
227 -17168.5340 12150.3207
228 -4987.8028 -17168.5340
229 -20487.1022 -4987.8028
230 -10119.3559 -20487.1022
231 -4023.2716 -10119.3559
232 -45997.4772 -4023.2716
233 -5802.8252 -45997.4772
234 -14093.0493 -5802.8252
235 1236.3612 -14093.0493
236 49960.4051 1236.3612
237 -11943.7461 49960.4051
238 -11514.5374 -11943.7461
239 -12808.9109 -11514.5374
240 6131.7396 -12808.9109
241 5690.5282 6131.7396
242 -6951.9362 5690.5282
243 -1741.3961 -6951.9362
244 -6997.1318 -1741.3961
245 -8877.8388 -6997.1318
246 -14514.4075 -8877.8388
247 7195.7562 -14514.4075
248 11881.2251 7195.7562
249 -9428.5606 11881.2251
250 4386.7814 -9428.5606
251 2896.0774 4386.7814
252 5629.9033 2896.0774
253 5271.6057 5629.9033
254 -15133.9913 5271.6057
255 4474.7394 -15133.9913
256 -7728.3003 4474.7394
257 20229.5153 -7728.3003
258 10897.6682 20229.5153
259 1607.7132 10897.6682
260 -8015.6267 1607.7132
261 1620.6585 -8015.6267
262 15447.0399 1620.6585
263 1723.7035 15447.0399
264 -5976.6153 1723.7035
265 3787.5628 -5976.6153
266 -13077.4974 3787.5628
267 278.4139 -13077.4974
268 33393.1050 278.4139
269 -7010.0034 33393.1050
270 306.9863 -7010.0034
271 2596.0749 306.9863
272 6382.5311 2596.0749
273 5456.2244 6382.5311
274 10252.4853 5456.2244
275 476.3422 10252.4853
276 20404.0745 476.3422
277 4299.9040 20404.0745
278 -6583.8814 4299.9040
279 12700.5066 -6583.8814
280 2879.2526 12700.5066
281 30972.0198 2879.2526
282 11765.7406 30972.0198
283 11538.2798 11765.7406
284 5255.3484 11538.2798
285 8113.9317 5255.3484
286 9175.3048 8113.9317
287 -3290.9784 9175.3048
288 -9276.6030 -3290.9784
289 NA -9276.6030
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] -22181.4982 21298.6550
[2,] 7359.3814 -22181.4982
[3,] -43546.7872 7359.3814
[4,] 11040.8391 -43546.7872
[5,] -15280.2836 11040.8391
[6,] 11725.2811 -15280.2836
[7,] 4539.3801 11725.2811
[8,] 19311.5600 4539.3801
[9,] -18549.8468 19311.5600
[10,] 17688.8664 -18549.8468
[11,] -3977.4460 17688.8664
[12,] 5120.9825 -3977.4460
[13,] -4261.0696 5120.9825
[14,] 39757.0455 -4261.0696
[15,] 6875.8629 39757.0455
[16,] -41988.2923 6875.8629
[17,] 24215.5644 -41988.2923
[18,] 6100.3998 24215.5644
[19,] 710.1252 6100.3998
[20,] 6835.4833 710.1252
[21,] 16562.0500 6835.4833
[22,] 85845.5017 16562.0500
[23,] 7586.2831 85845.5017
[24,] -22449.9309 7586.2831
[25,] -52566.1802 -22449.9309
[26,] -38067.9395 -52566.1802
[27,] -30732.4510 -38067.9395
[28,] 17335.4532 -30732.4510
[29,] -5025.3756 17335.4532
[30,] -47335.7030 -5025.3756
[31,] -8645.6132 -47335.7030
[32,] -5571.9448 -8645.6132
[33,] 4533.0930 -5571.9448
[34,] -1500.7788 4533.0930
[35,] 2678.3949 -1500.7788
[36,] 4470.0794 2678.3949
[37,] 6657.3966 4470.0794
[38,] 37792.5963 6657.3966
[39,] -898.4406 37792.5963
[40,] 9808.4745 -898.4406
[41,] 22765.9485 9808.4745
[42,] -10857.1302 22765.9485
[43,] -4417.9117 -10857.1302
[44,] -11227.5702 -4417.9117
[45,] -2413.8504 -11227.5702
[46,] -12060.2031 -2413.8504
[47,] -17999.6860 -12060.2031
[48,] 10567.6587 -17999.6860
[49,] -52474.4635 10567.6587
[50,] 37143.4506 -52474.4635
[51,] -25690.3421 37143.4506
[52,] -3799.8365 -25690.3421
[53,] -42586.5580 -3799.8365
[54,] -39307.8527 -42586.5580
[55,] -23110.2705 -39307.8527
[56,] -15567.1212 -23110.2705
[57,] 30015.7503 -15567.1212
[58,] -1229.0147 30015.7503
[59,] -11509.8189 -1229.0147
[60,] 9971.4910 -11509.8189
[61,] 5296.9777 9971.4910
[62,] 10806.3679 5296.9777
[63,] 3385.2601 10806.3679
[64,] -31476.0034 3385.2601
[65,] -16592.0433 -31476.0034
[66,] 47478.7699 -16592.0433
[67,] -9795.3726 47478.7699
[68,] -28154.6978 -9795.3726
[69,] -27917.5267 -28154.6978
[70,] 12631.2775 -27917.5267
[71,] 13960.4966 12631.2775
[72,] 9325.9635 13960.4966
[73,] -5486.8923 9325.9635
[74,] -32249.7755 -5486.8923
[75,] 5665.2571 -32249.7755
[76,] -16897.0447 5665.2571
[77,] -11318.0673 -16897.0447
[78,] -16828.8326 -11318.0673
[79,] 3450.7093 -16828.8326
[80,] 8959.2775 3450.7093
[81,] -10580.7625 8959.2775
[82,] 68493.5417 -10580.7625
[83,] 15062.5583 68493.5417
[84,] -1197.2187 15062.5583
[85,] -5534.6098 -1197.2187
[86,] 1065.6950 -5534.6098
[87,] -12876.0611 1065.6950
[88,] 76797.3749 -12876.0611
[89,] -21040.1444 76797.3749
[90,] 38827.4220 -21040.1444
[91,] -20756.5940 38827.4220
[92,] -14858.8141 -20756.5940
[93,] 19353.2699 -14858.8141
[94,] 9835.4733 19353.2699
[95,] 7829.1549 9835.4733
[96,] 18111.6052 7829.1549
[97,] -50976.5154 18111.6052
[98,] 18697.7657 -50976.5154
[99,] 6797.0358 18697.7657
[100,] -2312.1926 6797.0358
[101,] 5164.8627 -2312.1926
[102,] 32430.8009 5164.8627
[103,] 12868.8608 32430.8009
[104,] 6443.3629 12868.8608
[105,] -3145.3865 6443.3629
[106,] -28879.6859 -3145.3865
[107,] 194.7747 -28879.6859
[108,] 10219.0056 194.7747
[109,] 25451.5954 10219.0056
[110,] -16738.8462 25451.5954
[111,] -8889.2839 -16738.8462
[112,] -10093.2839 -8889.2839
[113,] 4597.8522 -10093.2839
[114,] -15243.3883 4597.8522
[115,] -10365.6745 -15243.3883
[116,] 5303.1562 -10365.6745
[117,] 35870.0186 5303.1562
[118,] -62154.7766 35870.0186
[119,] -8219.6818 -62154.7766
[120,] 1819.2371 -8219.6818
[121,] -28769.2914 1819.2371
[122,] -28857.7477 -28769.2914
[123,] -39215.7262 -28857.7477
[124,] 25193.8727 -39215.7262
[125,] 9013.1120 25193.8727
[126,] 6276.1051 9013.1120
[127,] 2482.0355 6276.1051
[128,] 48041.0277 2482.0355
[129,] 2253.1194 48041.0277
[130,] 384.7113 2253.1194
[131,] -13598.7697 384.7113
[132,] 13646.3440 -13598.7697
[133,] -35845.8617 13646.3440
[134,] -5134.5911 -35845.8617
[135,] -6892.6140 -5134.5911
[136,] -742.6953 -6892.6140
[137,] 12047.4618 -742.6953
[138,] -11262.1119 12047.4618
[139,] -4108.9658 -11262.1119
[140,] -18552.5637 -4108.9658
[141,] -10180.4599 -18552.5637
[142,] 14468.1463 -10180.4599
[143,] -32760.5293 14468.1463
[144,] 31571.9291 -32760.5293
[145,] -6242.3953 31571.9291
[146,] 25377.1822 -6242.3953
[147,] -14926.6128 25377.1822
[148,] -15941.9516 -14926.6128
[149,] 5148.1761 -15941.9516
[150,] 24089.5488 5148.1761
[151,] -6126.3172 24089.5488
[152,] -5659.0726 -6126.3172
[153,] -3260.2994 -5659.0726
[154,] -32204.9828 -3260.2994
[155,] 18913.7914 -32204.9828
[156,] 1841.1070 18913.7914
[157,] 11375.9878 1841.1070
[158,] -25076.4781 11375.9878
[159,] 8096.9822 -25076.4781
[160,] 868.3404 8096.9822
[161,] 10665.9752 868.3404
[162,] -2347.3877 10665.9752
[163,] 20798.5843 -2347.3877
[164,] -13857.7541 20798.5843
[165,] 15569.1531 -13857.7541
[166,] -13235.3157 15569.1531
[167,] 23801.4107 -13235.3157
[168,] 2222.3731 23801.4107
[169,] 74126.3810 2222.3731
[170,] 9619.4674 74126.3810
[171,] 10783.8362 9619.4674
[172,] 46845.8243 10783.8362
[173,] 32396.4975 46845.8243
[174,] -39240.9942 32396.4975
[175,] -2821.2590 -39240.9942
[176,] -24430.2853 -2821.2590
[177,] 12736.9168 -24430.2853
[178,] 1656.0763 12736.9168
[179,] 27334.1971 1656.0763
[180,] -11058.5201 27334.1971
[181,] -22137.9425 -11058.5201
[182,] 2377.4647 -22137.9425
[183,] 13651.7721 2377.4647
[184,] 973.0693 13651.7721
[185,] -1460.0007 973.0693
[186,] -4700.4073 -1460.0007
[187,] 7500.8093 -4700.4073
[188,] -3551.0166 7500.8093
[189,] 31478.3504 -3551.0166
[190,] 46578.1554 31478.3504
[191,] 468.3511 46578.1554
[192,] -9483.9658 468.3511
[193,] -13485.4028 -9483.9658
[194,] -20987.0986 -13485.4028
[195,] 11825.5130 -20987.0986
[196,] -14777.3598 11825.5130
[197,] 3109.0055 -14777.3598
[198,] -22299.0808 3109.0055
[199,] 2796.8683 -22299.0808
[200,] 5465.3280 2796.8683
[201,] 6815.8902 5465.3280
[202,] 22631.3482 6815.8902
[203,] -4166.8058 22631.3482
[204,] -23063.5822 -4166.8058
[205,] -3769.5000 -23063.5822
[206,] 4943.1791 -3769.5000
[207,] 32908.9598 4943.1791
[208,] 23207.0064 32908.9598
[209,] -3314.6020 23207.0064
[210,] 283.8547 -3314.6020
[211,] 32731.6549 283.8547
[212,] 753.9489 32731.6549
[213,] -2899.9777 753.9489
[214,] -15409.7773 -2899.9777
[215,] -14039.9061 -15409.7773
[216,] -29705.1333 -14039.9061
[217,] -7391.4721 -29705.1333
[218,] -6255.7399 -7391.4721
[219,] 14737.1290 -6255.7399
[220,] -16090.5277 14737.1290
[221,] -41143.6391 -16090.5277
[222,] -38742.5708 -41143.6391
[223,] -3387.8759 -38742.5708
[224,] 712.7598 -3387.8759
[225,] 3006.9550 712.7598
[226,] 12150.3207 3006.9550
[227,] -17168.5340 12150.3207
[228,] -4987.8028 -17168.5340
[229,] -20487.1022 -4987.8028
[230,] -10119.3559 -20487.1022
[231,] -4023.2716 -10119.3559
[232,] -45997.4772 -4023.2716
[233,] -5802.8252 -45997.4772
[234,] -14093.0493 -5802.8252
[235,] 1236.3612 -14093.0493
[236,] 49960.4051 1236.3612
[237,] -11943.7461 49960.4051
[238,] -11514.5374 -11943.7461
[239,] -12808.9109 -11514.5374
[240,] 6131.7396 -12808.9109
[241,] 5690.5282 6131.7396
[242,] -6951.9362 5690.5282
[243,] -1741.3961 -6951.9362
[244,] -6997.1318 -1741.3961
[245,] -8877.8388 -6997.1318
[246,] -14514.4075 -8877.8388
[247,] 7195.7562 -14514.4075
[248,] 11881.2251 7195.7562
[249,] -9428.5606 11881.2251
[250,] 4386.7814 -9428.5606
[251,] 2896.0774 4386.7814
[252,] 5629.9033 2896.0774
[253,] 5271.6057 5629.9033
[254,] -15133.9913 5271.6057
[255,] 4474.7394 -15133.9913
[256,] -7728.3003 4474.7394
[257,] 20229.5153 -7728.3003
[258,] 10897.6682 20229.5153
[259,] 1607.7132 10897.6682
[260,] -8015.6267 1607.7132
[261,] 1620.6585 -8015.6267
[262,] 15447.0399 1620.6585
[263,] 1723.7035 15447.0399
[264,] -5976.6153 1723.7035
[265,] 3787.5628 -5976.6153
[266,] -13077.4974 3787.5628
[267,] 278.4139 -13077.4974
[268,] 33393.1050 278.4139
[269,] -7010.0034 33393.1050
[270,] 306.9863 -7010.0034
[271,] 2596.0749 306.9863
[272,] 6382.5311 2596.0749
[273,] 5456.2244 6382.5311
[274,] 10252.4853 5456.2244
[275,] 476.3422 10252.4853
[276,] 20404.0745 476.3422
[277,] 4299.9040 20404.0745
[278,] -6583.8814 4299.9040
[279,] 12700.5066 -6583.8814
[280,] 2879.2526 12700.5066
[281,] 30972.0198 2879.2526
[282,] 11765.7406 30972.0198
[283,] 11538.2798 11765.7406
[284,] 5255.3484 11538.2798
[285,] 8113.9317 5255.3484
[286,] 9175.3048 8113.9317
[287,] -3290.9784 9175.3048
[288,] -9276.6030 -3290.9784
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 -22181.4982 21298.6550
2 7359.3814 -22181.4982
3 -43546.7872 7359.3814
4 11040.8391 -43546.7872
5 -15280.2836 11040.8391
6 11725.2811 -15280.2836
7 4539.3801 11725.2811
8 19311.5600 4539.3801
9 -18549.8468 19311.5600
10 17688.8664 -18549.8468
11 -3977.4460 17688.8664
12 5120.9825 -3977.4460
13 -4261.0696 5120.9825
14 39757.0455 -4261.0696
15 6875.8629 39757.0455
16 -41988.2923 6875.8629
17 24215.5644 -41988.2923
18 6100.3998 24215.5644
19 710.1252 6100.3998
20 6835.4833 710.1252
21 16562.0500 6835.4833
22 85845.5017 16562.0500
23 7586.2831 85845.5017
24 -22449.9309 7586.2831
25 -52566.1802 -22449.9309
26 -38067.9395 -52566.1802
27 -30732.4510 -38067.9395
28 17335.4532 -30732.4510
29 -5025.3756 17335.4532
30 -47335.7030 -5025.3756
31 -8645.6132 -47335.7030
32 -5571.9448 -8645.6132
33 4533.0930 -5571.9448
34 -1500.7788 4533.0930
35 2678.3949 -1500.7788
36 4470.0794 2678.3949
37 6657.3966 4470.0794
38 37792.5963 6657.3966
39 -898.4406 37792.5963
40 9808.4745 -898.4406
41 22765.9485 9808.4745
42 -10857.1302 22765.9485
43 -4417.9117 -10857.1302
44 -11227.5702 -4417.9117
45 -2413.8504 -11227.5702
46 -12060.2031 -2413.8504
47 -17999.6860 -12060.2031
48 10567.6587 -17999.6860
49 -52474.4635 10567.6587
50 37143.4506 -52474.4635
51 -25690.3421 37143.4506
52 -3799.8365 -25690.3421
53 -42586.5580 -3799.8365
54 -39307.8527 -42586.5580
55 -23110.2705 -39307.8527
56 -15567.1212 -23110.2705
57 30015.7503 -15567.1212
58 -1229.0147 30015.7503
59 -11509.8189 -1229.0147
60 9971.4910 -11509.8189
61 5296.9777 9971.4910
62 10806.3679 5296.9777
63 3385.2601 10806.3679
64 -31476.0034 3385.2601
65 -16592.0433 -31476.0034
66 47478.7699 -16592.0433
67 -9795.3726 47478.7699
68 -28154.6978 -9795.3726
69 -27917.5267 -28154.6978
70 12631.2775 -27917.5267
71 13960.4966 12631.2775
72 9325.9635 13960.4966
73 -5486.8923 9325.9635
74 -32249.7755 -5486.8923
75 5665.2571 -32249.7755
76 -16897.0447 5665.2571
77 -11318.0673 -16897.0447
78 -16828.8326 -11318.0673
79 3450.7093 -16828.8326
80 8959.2775 3450.7093
81 -10580.7625 8959.2775
82 68493.5417 -10580.7625
83 15062.5583 68493.5417
84 -1197.2187 15062.5583
85 -5534.6098 -1197.2187
86 1065.6950 -5534.6098
87 -12876.0611 1065.6950
88 76797.3749 -12876.0611
89 -21040.1444 76797.3749
90 38827.4220 -21040.1444
91 -20756.5940 38827.4220
92 -14858.8141 -20756.5940
93 19353.2699 -14858.8141
94 9835.4733 19353.2699
95 7829.1549 9835.4733
96 18111.6052 7829.1549
97 -50976.5154 18111.6052
98 18697.7657 -50976.5154
99 6797.0358 18697.7657
100 -2312.1926 6797.0358
101 5164.8627 -2312.1926
102 32430.8009 5164.8627
103 12868.8608 32430.8009
104 6443.3629 12868.8608
105 -3145.3865 6443.3629
106 -28879.6859 -3145.3865
107 194.7747 -28879.6859
108 10219.0056 194.7747
109 25451.5954 10219.0056
110 -16738.8462 25451.5954
111 -8889.2839 -16738.8462
112 -10093.2839 -8889.2839
113 4597.8522 -10093.2839
114 -15243.3883 4597.8522
115 -10365.6745 -15243.3883
116 5303.1562 -10365.6745
117 35870.0186 5303.1562
118 -62154.7766 35870.0186
119 -8219.6818 -62154.7766
120 1819.2371 -8219.6818
121 -28769.2914 1819.2371
122 -28857.7477 -28769.2914
123 -39215.7262 -28857.7477
124 25193.8727 -39215.7262
125 9013.1120 25193.8727
126 6276.1051 9013.1120
127 2482.0355 6276.1051
128 48041.0277 2482.0355
129 2253.1194 48041.0277
130 384.7113 2253.1194
131 -13598.7697 384.7113
132 13646.3440 -13598.7697
133 -35845.8617 13646.3440
134 -5134.5911 -35845.8617
135 -6892.6140 -5134.5911
136 -742.6953 -6892.6140
137 12047.4618 -742.6953
138 -11262.1119 12047.4618
139 -4108.9658 -11262.1119
140 -18552.5637 -4108.9658
141 -10180.4599 -18552.5637
142 14468.1463 -10180.4599
143 -32760.5293 14468.1463
144 31571.9291 -32760.5293
145 -6242.3953 31571.9291
146 25377.1822 -6242.3953
147 -14926.6128 25377.1822
148 -15941.9516 -14926.6128
149 5148.1761 -15941.9516
150 24089.5488 5148.1761
151 -6126.3172 24089.5488
152 -5659.0726 -6126.3172
153 -3260.2994 -5659.0726
154 -32204.9828 -3260.2994
155 18913.7914 -32204.9828
156 1841.1070 18913.7914
157 11375.9878 1841.1070
158 -25076.4781 11375.9878
159 8096.9822 -25076.4781
160 868.3404 8096.9822
161 10665.9752 868.3404
162 -2347.3877 10665.9752
163 20798.5843 -2347.3877
164 -13857.7541 20798.5843
165 15569.1531 -13857.7541
166 -13235.3157 15569.1531
167 23801.4107 -13235.3157
168 2222.3731 23801.4107
169 74126.3810 2222.3731
170 9619.4674 74126.3810
171 10783.8362 9619.4674
172 46845.8243 10783.8362
173 32396.4975 46845.8243
174 -39240.9942 32396.4975
175 -2821.2590 -39240.9942
176 -24430.2853 -2821.2590
177 12736.9168 -24430.2853
178 1656.0763 12736.9168
179 27334.1971 1656.0763
180 -11058.5201 27334.1971
181 -22137.9425 -11058.5201
182 2377.4647 -22137.9425
183 13651.7721 2377.4647
184 973.0693 13651.7721
185 -1460.0007 973.0693
186 -4700.4073 -1460.0007
187 7500.8093 -4700.4073
188 -3551.0166 7500.8093
189 31478.3504 -3551.0166
190 46578.1554 31478.3504
191 468.3511 46578.1554
192 -9483.9658 468.3511
193 -13485.4028 -9483.9658
194 -20987.0986 -13485.4028
195 11825.5130 -20987.0986
196 -14777.3598 11825.5130
197 3109.0055 -14777.3598
198 -22299.0808 3109.0055
199 2796.8683 -22299.0808
200 5465.3280 2796.8683
201 6815.8902 5465.3280
202 22631.3482 6815.8902
203 -4166.8058 22631.3482
204 -23063.5822 -4166.8058
205 -3769.5000 -23063.5822
206 4943.1791 -3769.5000
207 32908.9598 4943.1791
208 23207.0064 32908.9598
209 -3314.6020 23207.0064
210 283.8547 -3314.6020
211 32731.6549 283.8547
212 753.9489 32731.6549
213 -2899.9777 753.9489
214 -15409.7773 -2899.9777
215 -14039.9061 -15409.7773
216 -29705.1333 -14039.9061
217 -7391.4721 -29705.1333
218 -6255.7399 -7391.4721
219 14737.1290 -6255.7399
220 -16090.5277 14737.1290
221 -41143.6391 -16090.5277
222 -38742.5708 -41143.6391
223 -3387.8759 -38742.5708
224 712.7598 -3387.8759
225 3006.9550 712.7598
226 12150.3207 3006.9550
227 -17168.5340 12150.3207
228 -4987.8028 -17168.5340
229 -20487.1022 -4987.8028
230 -10119.3559 -20487.1022
231 -4023.2716 -10119.3559
232 -45997.4772 -4023.2716
233 -5802.8252 -45997.4772
234 -14093.0493 -5802.8252
235 1236.3612 -14093.0493
236 49960.4051 1236.3612
237 -11943.7461 49960.4051
238 -11514.5374 -11943.7461
239 -12808.9109 -11514.5374
240 6131.7396 -12808.9109
241 5690.5282 6131.7396
242 -6951.9362 5690.5282
243 -1741.3961 -6951.9362
244 -6997.1318 -1741.3961
245 -8877.8388 -6997.1318
246 -14514.4075 -8877.8388
247 7195.7562 -14514.4075
248 11881.2251 7195.7562
249 -9428.5606 11881.2251
250 4386.7814 -9428.5606
251 2896.0774 4386.7814
252 5629.9033 2896.0774
253 5271.6057 5629.9033
254 -15133.9913 5271.6057
255 4474.7394 -15133.9913
256 -7728.3003 4474.7394
257 20229.5153 -7728.3003
258 10897.6682 20229.5153
259 1607.7132 10897.6682
260 -8015.6267 1607.7132
261 1620.6585 -8015.6267
262 15447.0399 1620.6585
263 1723.7035 15447.0399
264 -5976.6153 1723.7035
265 3787.5628 -5976.6153
266 -13077.4974 3787.5628
267 278.4139 -13077.4974
268 33393.1050 278.4139
269 -7010.0034 33393.1050
270 306.9863 -7010.0034
271 2596.0749 306.9863
272 6382.5311 2596.0749
273 5456.2244 6382.5311
274 10252.4853 5456.2244
275 476.3422 10252.4853
276 20404.0745 476.3422
277 4299.9040 20404.0745
278 -6583.8814 4299.9040
279 12700.5066 -6583.8814
280 2879.2526 12700.5066
281 30972.0198 2879.2526
282 11765.7406 30972.0198
283 11538.2798 11765.7406
284 5255.3484 11538.2798
285 8113.9317 5255.3484
286 9175.3048 8113.9317
287 -3290.9784 9175.3048
288 -9276.6030 -3290.9784
> 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/717v21353091884.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/8xeut1353091884.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/9fmws1353091884.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/1029011353091884.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/11lh5o1353091884.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/120qqm1353091884.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/13bl1h1353091884.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/14tnpx1353091884.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/15fpgy1353091884.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/166o881353091885.tab")
+ }
>
> try(system("convert tmp/1empx1353091884.ps tmp/1empx1353091884.png",intern=TRUE))
character(0)
> try(system("convert tmp/2rfmd1353091884.ps tmp/2rfmd1353091884.png",intern=TRUE))
character(0)
> try(system("convert tmp/3mvt51353091884.ps tmp/3mvt51353091884.png",intern=TRUE))
character(0)
> try(system("convert tmp/45vy51353091884.ps tmp/45vy51353091884.png",intern=TRUE))
character(0)
> try(system("convert tmp/5qdf11353091884.ps tmp/5qdf11353091884.png",intern=TRUE))
character(0)
> try(system("convert tmp/68qo91353091884.ps tmp/68qo91353091884.png",intern=TRUE))
character(0)
> try(system("convert tmp/717v21353091884.ps tmp/717v21353091884.png",intern=TRUE))
character(0)
> try(system("convert tmp/8xeut1353091884.ps tmp/8xeut1353091884.png",intern=TRUE))
character(0)
> try(system("convert tmp/9fmws1353091884.ps tmp/9fmws1353091884.png",intern=TRUE))
character(0)
> try(system("convert tmp/1029011353091884.ps tmp/1029011353091884.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
19.239 1.373 20.638