R version 2.13.0 (2011-04-13)
Copyright (C) 2011 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i486-pc-linux-gnu (32-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> x <- array(list(10
+ ,14688
+ ,10
+ ,4
+ ,0
+ ,0
+ ,6023
+ ,10
+ ,46660
+ ,20
+ ,12
+ ,5
+ ,13
+ ,6179
+ ,10
+ ,7199
+ ,5
+ ,7
+ ,0
+ ,0
+ ,1644
+ ,10
+ ,21054
+ ,16
+ ,4
+ ,0
+ ,0
+ ,855
+ ,10
+ ,80953
+ ,25
+ ,49
+ ,8
+ ,27
+ ,56622
+ ,10
+ ,19349
+ ,11
+ ,13
+ ,1
+ ,0
+ ,3895
+ ,11
+ ,173260
+ ,63
+ ,41
+ ,21
+ ,78
+ ,37238
+ ,11
+ ,38214
+ ,34
+ ,16
+ ,8
+ ,21
+ ,8773
+ ,11
+ ,17547
+ ,5
+ ,0
+ ,1
+ ,4
+ ,3926
+ ,11
+ ,43287
+ ,14
+ ,13
+ ,19
+ ,64
+ ,43750
+ ,11
+ ,78800
+ ,42
+ ,20
+ ,26
+ ,66
+ ,33032
+ ,11
+ ,65029
+ ,17
+ ,21
+ ,18
+ ,61
+ ,32551
+ ,11
+ ,31414
+ ,19
+ ,18
+ ,8
+ ,9
+ ,14116
+ ,11
+ ,7176
+ ,17
+ ,0
+ ,0
+ ,0
+ ,1423
+ ,11
+ ,24188
+ ,24
+ ,8
+ ,4
+ ,7
+ ,5950
+ ,11
+ ,61857
+ ,25
+ ,23
+ ,11
+ ,30
+ ,25162
+ ,11
+ ,22938
+ ,10
+ ,12
+ ,1
+ ,0
+ ,1168
+ ,11
+ ,43410
+ ,19
+ ,7
+ ,1
+ ,3
+ ,2781
+ ,11
+ ,19764
+ ,12
+ ,10
+ ,2
+ ,4
+ ,5752
+ ,11
+ ,53117
+ ,22
+ ,3
+ ,12
+ ,32
+ ,32689
+ ,11
+ ,33170
+ ,18
+ ,1
+ ,18
+ ,22
+ ,22807
+ ,11
+ ,232138
+ ,62
+ ,131
+ ,31
+ ,90
+ ,161647
+ ,11
+ ,111665
+ ,34
+ ,39
+ ,28
+ ,80
+ ,49810
+ ,11
+ ,152871
+ ,58
+ ,59
+ ,28
+ ,90
+ ,79892
+ ,11
+ ,100750
+ ,72
+ ,83
+ ,30
+ ,93
+ ,140867
+ ,11
+ ,221698
+ ,45
+ ,105
+ ,33
+ ,110
+ ,77873
+ ,11
+ ,209641
+ ,42
+ ,62
+ ,24
+ ,60
+ ,97068
+ ,11
+ ,52746
+ ,25
+ ,0
+ ,26
+ ,70
+ ,68504
+ ,11
+ ,215641
+ ,46
+ ,71
+ ,32
+ ,92
+ ,55813
+ ,11
+ ,85439
+ ,33
+ ,32
+ ,28
+ ,78
+ ,56926
+ ,11
+ ,145790
+ ,63
+ ,30
+ ,26
+ ,99
+ ,38692
+ ,11
+ ,132943
+ ,40
+ ,83
+ ,39
+ ,156
+ ,97500
+ ,11
+ ,132487
+ ,41
+ ,71
+ ,36
+ ,98
+ ,40735
+ ,11
+ ,135473
+ ,41
+ ,82
+ ,23
+ ,66
+ ,99645
+ ,11
+ ,84853
+ ,31
+ ,38
+ ,30
+ ,109
+ ,61542
+ ,11
+ ,210767
+ ,60
+ ,94
+ ,35
+ ,133
+ ,117478
+ ,11
+ ,97839
+ ,38
+ ,25
+ ,24
+ ,66
+ ,94785
+ ,11
+ ,120221
+ ,37
+ ,53
+ ,22
+ ,69
+ ,43836
+ ,11
+ ,224549
+ ,50
+ ,54
+ ,31
+ ,114
+ ,120662
+ ,11
+ ,116408
+ ,61
+ ,34
+ ,39
+ ,31
+ ,71570
+ ,11
+ ,158015
+ ,29
+ ,59
+ ,31
+ ,114
+ ,71220
+ ,12
+ ,150629
+ ,44
+ ,85
+ ,33
+ ,113
+ ,71595
+ ,12
+ ,123185
+ ,40
+ ,49
+ ,22
+ ,51
+ ,38361
+ ,12
+ ,124817
+ ,40
+ ,47
+ ,25
+ ,81
+ ,30727
+ ,12
+ ,104389
+ ,45
+ ,135
+ ,41
+ ,148
+ ,123969
+ ,12
+ ,162765
+ ,32
+ ,68
+ ,28
+ ,107
+ ,120293
+ ,12
+ ,149061
+ ,44
+ ,43
+ ,26
+ ,93
+ ,116174
+ ,12
+ ,84207
+ ,29
+ ,14
+ ,30
+ ,101
+ ,111194
+ ,12
+ ,201940
+ ,38
+ ,109
+ ,31
+ ,92
+ ,85872
+ ,12
+ ,168809
+ ,66
+ ,76
+ ,28
+ ,100
+ ,72260
+ ,12
+ ,176508
+ ,54
+ ,60
+ ,38
+ ,93
+ ,83123
+ ,12
+ ,99466
+ ,50
+ ,28
+ ,23
+ ,69
+ ,192565
+ ,12
+ ,106408
+ ,30
+ ,33
+ ,14
+ ,37
+ ,31081
+ ,12
+ ,344297
+ ,75
+ ,80
+ ,30
+ ,108
+ ,67654
+ ,12
+ ,218946
+ ,41
+ ,76
+ ,29
+ ,96
+ ,80670
+ ,12
+ ,122774
+ ,45
+ ,24
+ ,24
+ ,69
+ ,15986
+ ,12
+ ,153935
+ ,33
+ ,50
+ ,25
+ ,84
+ ,65553
+ ,12
+ ,206161
+ ,71
+ ,75
+ ,28
+ ,99
+ ,73107
+ ,12
+ ,183167
+ ,66
+ ,91
+ ,39
+ ,138
+ ,82875
+ ,12
+ ,144966
+ ,144
+ ,39
+ ,32
+ ,50
+ ,37510
+ ,12
+ ,34662
+ ,25
+ ,18
+ ,6
+ ,7
+ ,36278
+ ,12
+ ,116048
+ ,64
+ ,50
+ ,20
+ ,50
+ ,33416
+ ,12
+ ,77945
+ ,20
+ ,28
+ ,19
+ ,67
+ ,25272
+ ,12
+ ,65490
+ ,22
+ ,27
+ ,16
+ ,40
+ ,18653
+ ,12
+ ,81872
+ ,45
+ ,32
+ ,16
+ ,61
+ ,27913
+ ,12
+ ,101011
+ ,34
+ ,30
+ ,13
+ ,39
+ ,27570
+ ,12
+ ,71965
+ ,35
+ ,32
+ ,15
+ ,56
+ ,46090
+ ,12
+ ,72880
+ ,33
+ ,14
+ ,19
+ ,67
+ ,24094
+ ,12
+ ,62792
+ ,35
+ ,28
+ ,15
+ ,51
+ ,30884
+ ,12
+ ,31774
+ ,23
+ ,0
+ ,17
+ ,47
+ ,4143
+ ,12
+ ,182192
+ ,52
+ ,77
+ ,40
+ ,138
+ ,70054
+ ,12
+ ,67989
+ ,21
+ ,23
+ ,18
+ ,57
+ ,33747
+ ,12
+ ,52915
+ ,14
+ ,20
+ ,18
+ ,49
+ ,24006
+ ,12
+ ,56375
+ ,30
+ ,10
+ ,13
+ ,40
+ ,41385
+ ,12
+ ,30989
+ ,14
+ ,5
+ ,17
+ ,41
+ ,4154
+ ,12
+ ,135131
+ ,66
+ ,38
+ ,15
+ ,60
+ ,34029
+ ,12
+ ,76990
+ ,27
+ ,42
+ ,17
+ ,67
+ ,24069
+ ,12
+ ,64175
+ ,42
+ ,37
+ ,18
+ ,59
+ ,21792
+ ,12
+ ,59194
+ ,31
+ ,7
+ ,24
+ ,68
+ ,37636
+ ,12
+ ,53515
+ ,28
+ ,5
+ ,22
+ ,81
+ ,12934
+ ,12
+ ,89746
+ ,36
+ ,28
+ ,18
+ ,55
+ ,24760
+ ,12
+ ,49862
+ ,37
+ ,17
+ ,14
+ ,54
+ ,36341
+ ,12
+ ,85574
+ ,34
+ ,37
+ ,21
+ ,36
+ ,24266
+ ,12
+ ,220801
+ ,75
+ ,51
+ ,18
+ ,63
+ ,44418
+ ,12
+ ,49289
+ ,19
+ ,15
+ ,24
+ ,44
+ ,10672
+ ,12
+ ,135781
+ ,31
+ ,45
+ ,14
+ ,44
+ ,24874
+ ,12
+ ,82316
+ ,32
+ ,27
+ ,10
+ ,38
+ ,21233
+ ,12
+ ,133368
+ ,36
+ ,37
+ ,16
+ ,57
+ ,32755
+ ,12
+ ,92499
+ ,32
+ ,25
+ ,18
+ ,55
+ ,21399
+ ,12
+ ,37460
+ ,20
+ ,5
+ ,7
+ ,14
+ ,3738
+ ,12
+ ,91735
+ ,35
+ ,7
+ ,18
+ ,49
+ ,16563
+ ,12
+ ,187681
+ ,62
+ ,114
+ ,28
+ ,91
+ ,92945
+ ,12
+ ,79619
+ ,43
+ ,42
+ ,11
+ ,32
+ ,20760
+ ,12
+ ,98104
+ ,54
+ ,55
+ ,17
+ ,47
+ ,25568
+ ,12
+ ,80949
+ ,17
+ ,11
+ ,16
+ ,56
+ ,18472
+ ,12
+ ,31706
+ ,13
+ ,26
+ ,10
+ ,25
+ ,25139
+ ,12
+ ,70344
+ ,28
+ ,20
+ ,16
+ ,61
+ ,28904
+ ,12
+ ,99923
+ ,66
+ ,44
+ ,25
+ ,80
+ ,32334
+ ,12
+ ,91899
+ ,35
+ ,18
+ ,15
+ ,33
+ ,36874
+ ,12
+ ,235454
+ ,73
+ ,151
+ ,32
+ ,106
+ ,89275
+ ,12
+ ,46455
+ ,20
+ ,22
+ ,17
+ ,12
+ ,11342
+ ,12
+ ,60812
+ ,44
+ ,26
+ ,13
+ ,43
+ ,22574
+ ,12
+ ,77648
+ ,47
+ ,31
+ ,16
+ ,52
+ ,21280
+ ,12
+ ,101523
+ ,42
+ ,59
+ ,22
+ ,87
+ ,61056
+ ,12
+ ,254488
+ ,83
+ ,120
+ ,39
+ ,117
+ ,103772
+ ,12
+ ,103597
+ ,43
+ ,27
+ ,16
+ ,49
+ ,23789
+ ,12
+ ,136084
+ ,30
+ ,27
+ ,13
+ ,41
+ ,27142
+ ,12
+ ,62088
+ ,38
+ ,13
+ ,16
+ ,38
+ ,14483
+ ,12
+ ,86680
+ ,31
+ ,37
+ ,14
+ ,49
+ ,22197
+ ,12
+ ,38395
+ ,31
+ ,16
+ ,16
+ ,43
+ ,16380
+ ,12
+ ,179321
+ ,89
+ ,108
+ ,30
+ ,103
+ ,101193
+ ,12
+ ,40248
+ ,16
+ ,8
+ ,4
+ ,7
+ ,5444
+ ,12
+ ,70551
+ ,31
+ ,23
+ ,15
+ ,47
+ ,30143
+ ,12
+ ,84105
+ ,20
+ ,17
+ ,17
+ ,63
+ ,20055
+ ,12
+ ,112611
+ ,41
+ ,46
+ ,20
+ ,73
+ ,26706
+ ,12
+ ,73624
+ ,24
+ ,30
+ ,17
+ ,40
+ ,19499
+ ,12
+ ,89806
+ ,42
+ ,27
+ ,16
+ ,40
+ ,27975
+ ,12
+ ,98866
+ ,18
+ ,25
+ ,13
+ ,39
+ ,23517
+ ,12
+ ,61254
+ ,46
+ ,36
+ ,14
+ ,30
+ ,13310
+ ,12
+ ,99643
+ ,55
+ ,33
+ ,17
+ ,46
+ ,54968
+ ,12
+ ,310839
+ ,92
+ ,130
+ ,24
+ ,85
+ ,162901
+ ,12
+ ,181633
+ ,70
+ ,47
+ ,30
+ ,73
+ ,64466
+ ,12
+ ,50090
+ ,20
+ ,16
+ ,17
+ ,36
+ ,21152
+ ,12
+ ,92661
+ ,61
+ ,40
+ ,17
+ ,41
+ ,35232
+ ,12
+ ,46698
+ ,45
+ ,14
+ ,13
+ ,33
+ ,10288
+ ,12
+ ,61361
+ ,77
+ ,27
+ ,12
+ ,25
+ ,13294
+ ,12
+ ,84337
+ ,26
+ ,38
+ ,14
+ ,35
+ ,24548
+ ,12
+ ,72535
+ ,14
+ ,17
+ ,16
+ ,39
+ ,19540
+ ,12
+ ,84856
+ ,29
+ ,29
+ ,17
+ ,43
+ ,41369
+ ,12
+ ,124064
+ ,40
+ ,43
+ ,22
+ ,86
+ ,22827
+ ,12
+ ,54157
+ ,19
+ ,37
+ ,15
+ ,30
+ ,18625
+ ,12
+ ,76302
+ ,31
+ ,29
+ ,20
+ ,68
+ ,30976
+ ,12
+ ,59382
+ ,49
+ ,29
+ ,12
+ ,24
+ ,26263
+ ,12
+ ,65745
+ ,53
+ ,26
+ ,21
+ ,42
+ ,38084
+ ,12
+ ,74914
+ ,30
+ ,35
+ ,23
+ ,79
+ ,65892
+ ,12
+ ,101645
+ ,63
+ ,20
+ ,11
+ ,38
+ ,17140
+ ,12
+ ,73504
+ ,35
+ ,23
+ ,7
+ ,3
+ ,35130
+ ,12
+ ,99373
+ ,63
+ ,12
+ ,18
+ ,62
+ ,28394
+ ,12
+ ,135458
+ ,81
+ ,43
+ ,12
+ ,29
+ ,19474
+ ,12
+ ,184510
+ ,49
+ ,64
+ ,40
+ ,140
+ ,110681
+ ,12
+ ,385534
+ ,92
+ ,121
+ ,25
+ ,91
+ ,119182
+ ,12
+ ,108043
+ ,62
+ ,45
+ ,14
+ ,38
+ ,34553
+ ,12
+ ,328107
+ ,65
+ ,129
+ ,41
+ ,144
+ ,105547
+ ,12
+ ,362301
+ ,119
+ ,76
+ ,34
+ ,110
+ ,100708
+ ,12
+ ,130585
+ ,46
+ ,67
+ ,29
+ ,107
+ ,95364
+ ,12
+ ,85709
+ ,44
+ ,21
+ ,10
+ ,35
+ ,35944
+ ,12
+ ,174184
+ ,53
+ ,72
+ ,25
+ ,70
+ ,55183
+ ,12
+ ,81437
+ ,38
+ ,37
+ ,14
+ ,40
+ ,28579
+ ,12
+ ,120982
+ ,56
+ ,58
+ ,28
+ ,103
+ ,84786
+ ,12
+ ,187559
+ ,121
+ ,75
+ ,36
+ ,116
+ ,73511
+ ,12
+ ,235800
+ ,94
+ ,105
+ ,23
+ ,62
+ ,149193
+ ,12
+ ,220516
+ ,62
+ ,98
+ ,32
+ ,119
+ ,55801
+ ,12
+ ,210907
+ ,56
+ ,79
+ ,30
+ ,94
+ ,112285
+ ,12
+ ,155754
+ ,61
+ ,74
+ ,20
+ ,31
+ ,40652
+ ,12
+ ,152299
+ ,53
+ ,62
+ ,33
+ ,98
+ ,61370
+ ,12
+ ,79863
+ ,37
+ ,29
+ ,22
+ ,49
+ ,29011
+ ,12
+ ,346485
+ ,90
+ ,118
+ ,38
+ ,111
+ ,106117
+ ,12
+ ,135649
+ ,46
+ ,99
+ ,30
+ ,108
+ ,62133
+ ,12
+ ,134019
+ ,53
+ ,32
+ ,18
+ ,54
+ ,28987
+ ,12
+ ,167488
+ ,45
+ ,69
+ ,28
+ ,72
+ ,83737
+ ,12
+ ,108446
+ ,60
+ ,22
+ ,17
+ ,65
+ ,46300
+ ,12
+ ,27634
+ ,20
+ ,2
+ ,16
+ ,17
+ ,5839
+ ,12
+ ,232317
+ ,54
+ ,118
+ ,33
+ ,106
+ ,129838
+ ,12
+ ,22648
+ ,19
+ ,12
+ ,13
+ ,28
+ ,44339
+ ,12
+ ,119308
+ ,30
+ ,32
+ ,16
+ ,58
+ ,23686
+ ,12
+ ,194979
+ ,62
+ ,66
+ ,40
+ ,151
+ ,74011
+ ,12
+ ,81240
+ ,66
+ ,58
+ ,17
+ ,50
+ ,45549
+ ,12
+ ,121848
+ ,39
+ ,37
+ ,17
+ ,52
+ ,30594
+ ,12
+ ,86678
+ ,40
+ ,12
+ ,15
+ ,50
+ ,29156
+ ,12
+ ,73566
+ ,32
+ ,22
+ ,23
+ ,67
+ ,22618
+ ,12
+ ,102010
+ ,53
+ ,28
+ ,13
+ ,44
+ ,44332
+ ,12
+ ,151101
+ ,32
+ ,48
+ ,35
+ ,124
+ ,83209
+ ,12
+ ,207176
+ ,70
+ ,56
+ ,32
+ ,114
+ ,87011
+ ,12
+ ,165543
+ ,65
+ ,70
+ ,32
+ ,90
+ ,51633
+ ,12
+ ,182999
+ ,88
+ ,73
+ ,34
+ ,104
+ ,225920
+ ,12
+ ,256462
+ ,105
+ ,123
+ ,35
+ ,124
+ ,105195
+ ,12
+ ,244052
+ ,68
+ ,101
+ ,44
+ ,164
+ ,143558
+ ,12
+ ,63123
+ ,43
+ ,34
+ ,17
+ ,60
+ ,39067
+ ,12
+ ,87186
+ ,54
+ ,28
+ ,15
+ ,48
+ ,51776
+ ,12
+ ,182079
+ ,63
+ ,124
+ ,33
+ ,118
+ ,102860
+ ,12
+ ,95227
+ ,34
+ ,37
+ ,32
+ ,48
+ ,34777
+ ,12
+ ,165446
+ ,33
+ ,69
+ ,25
+ ,60
+ ,57635
+ ,12
+ ,229242
+ ,247
+ ,63
+ ,31
+ ,119
+ ,91721
+ ,12
+ ,167542
+ ,66
+ ,59
+ ,28
+ ,70
+ ,55461
+ ,12
+ ,225060
+ ,93
+ ,93
+ ,41
+ ,139
+ ,109825
+ ,12
+ ,164709
+ ,109
+ ,81
+ ,31
+ ,63
+ ,87771
+ ,12
+ ,131698
+ ,65
+ ,60
+ ,19
+ ,59
+ ,65622
+ ,12
+ ,56613
+ ,19
+ ,15
+ ,12
+ ,35
+ ,8019
+ ,12
+ ,250579
+ ,83
+ ,130
+ ,38
+ ,120
+ ,103487
+ ,12
+ ,44296
+ ,25
+ ,10
+ ,20
+ ,38
+ ,18779
+ ,12
+ ,152474
+ ,65
+ ,106
+ ,32
+ ,83
+ ,65567
+ ,12
+ ,86230
+ ,44
+ ,21
+ ,17
+ ,54
+ ,18632
+ ,12
+ ,258873
+ ,60
+ ,104
+ ,40
+ ,124
+ ,80444
+ ,12
+ ,41566
+ ,35
+ ,0
+ ,9
+ ,27
+ ,13497
+ ,12
+ ,250047
+ ,81
+ ,41
+ ,18
+ ,61
+ ,83305
+ ,12
+ ,243511
+ ,71
+ ,133
+ ,42
+ ,110
+ ,101338
+ ,12
+ ,62215
+ ,27
+ ,24
+ ,10
+ ,25
+ ,19546
+ ,12
+ ,104838
+ ,49
+ ,46
+ ,16
+ ,52
+ ,29236
+ ,12
+ ,77272
+ ,59
+ ,21
+ ,16
+ ,48
+ ,27507
+ ,12
+ ,148446
+ ,91
+ ,135
+ ,37
+ ,129
+ ,126846
+ ,12
+ ,69304
+ ,30
+ ,40
+ ,19
+ ,40
+ ,22818
+ ,12
+ ,102538
+ ,57
+ ,50
+ ,15
+ ,45
+ ,45833
+ ,12
+ ,101097
+ ,64
+ ,30
+ ,14
+ ,41
+ ,31701
+ ,12
+ ,223632
+ ,73
+ ,105
+ ,33
+ ,120
+ ,72654
+ ,12
+ ,141722
+ ,94
+ ,19
+ ,27
+ ,36
+ ,75345
+ ,12
+ ,351067
+ ,95
+ ,136
+ ,45
+ ,160
+ ,213688
+ ,12
+ ,311473
+ ,112
+ ,128
+ ,38
+ ,129
+ ,132068
+ ,12
+ ,152601
+ ,48
+ ,46
+ ,24
+ ,46
+ ,73224
+ ,12
+ ,215147
+ ,58
+ ,101
+ ,36
+ ,117
+ ,99052
+ ,12
+ ,52164
+ ,52
+ ,32
+ ,16
+ ,43
+ ,16734
+ ,12
+ ,91005
+ ,29
+ ,29
+ ,11
+ ,23
+ ,32073
+ ,12
+ ,265318
+ ,117
+ ,110
+ ,52
+ ,197
+ ,115929
+ ,12
+ ,299775
+ ,95
+ ,91
+ ,31
+ ,97
+ ,98952
+ ,12
+ ,241066
+ ,82
+ ,75
+ ,45
+ ,93
+ ,67267
+ ,12
+ ,118612
+ ,46
+ ,54
+ ,12
+ ,43
+ ,30080
+ ,12
+ ,173326
+ ,88
+ ,86
+ ,44
+ ,148
+ ,71701
+ ,12
+ ,225548
+ ,112
+ ,81
+ ,31
+ ,94
+ ,85323
+ ,12
+ ,277965
+ ,89
+ ,115
+ ,39
+ ,133
+ ,139077
+ ,12
+ ,317394
+ ,86
+ ,116
+ ,31
+ ,82
+ ,91413
+ ,12
+ ,177939
+ ,82
+ ,55
+ ,36
+ ,73
+ ,46821
+ ,12
+ ,92630
+ ,40
+ ,27
+ ,8
+ ,24
+ ,18513
+ ,12
+ ,140344
+ ,53
+ ,33
+ ,25
+ ,73
+ ,51715
+ ,12
+ ,230964
+ ,53
+ ,102
+ ,30
+ ,115
+ ,133824
+ ,12
+ ,199476
+ ,70
+ ,87
+ ,32
+ ,105
+ ,135400
+ ,12
+ ,174724
+ ,92
+ ,123
+ ,34
+ ,120
+ ,69112
+ ,12
+ ,128423
+ ,64
+ ,32
+ ,38
+ ,120
+ ,92696
+ ,12
+ ,341570
+ ,168
+ ,94
+ ,21
+ ,78
+ ,117105
+ ,12
+ ,64187
+ ,27
+ ,10
+ ,16
+ ,54
+ ,20154
+ ,12
+ ,150580
+ ,77
+ ,27
+ ,22
+ ,71
+ ,45588
+ ,12
+ ,204713
+ ,71
+ ,68
+ ,33
+ ,69
+ ,41140
+ ,12
+ ,244749
+ ,95
+ ,98
+ ,33
+ ,115
+ ,76643
+ ,12
+ ,98146
+ ,40
+ ,15
+ ,17
+ ,48
+ ,27114
+ ,12
+ ,351619
+ ,139
+ ,95
+ ,40
+ ,141
+ ,115168
+ ,12
+ ,58981
+ ,36
+ ,0
+ ,23
+ ,61
+ ,42744
+ ,12
+ ,233328
+ ,132
+ ,92
+ ,28
+ ,102
+ ,120733
+ ,12
+ ,182613
+ ,39
+ ,81
+ ,28
+ ,99
+ ,138599
+ ,12
+ ,89113
+ ,39
+ ,19
+ ,14
+ ,32
+ ,82206
+ ,12
+ ,324799
+ ,154
+ ,158
+ ,47
+ ,168
+ ,97668
+ ,12
+ ,143246
+ ,103
+ ,67
+ ,27
+ ,104
+ ,106671
+ ,12
+ ,40151
+ ,29
+ ,16
+ ,9
+ ,19
+ ,13127
+ ,12
+ ,158399
+ ,39
+ ,23
+ ,18
+ ,57
+ ,32928
+ ,12
+ ,195838
+ ,67
+ ,111
+ ,31
+ ,98
+ ,102372
+ ,12
+ ,286468
+ ,144
+ ,57
+ ,29
+ ,107
+ ,89691
+ ,12
+ ,175824
+ ,107
+ ,57
+ ,20
+ ,77
+ ,39992
+ ,12
+ ,143756
+ ,46
+ ,105
+ ,34
+ ,120
+ ,69094
+ ,12
+ ,139942
+ ,42
+ ,54
+ ,22
+ ,87
+ ,65461
+ ,12
+ ,104011
+ ,55
+ ,25
+ ,22
+ ,79
+ ,63339
+ ,12
+ ,196553
+ ,57
+ ,41
+ ,29
+ ,99
+ ,95260
+ ,12
+ ,180083
+ ,66
+ ,63
+ ,34
+ ,70
+ ,53855
+ ,12
+ ,260561
+ ,75
+ ,114
+ ,43
+ ,158
+ ,90183
+ ,12
+ ,294424
+ ,77
+ ,107
+ ,33
+ ,124
+ ,101494
+ ,12
+ ,99611
+ ,35
+ ,41
+ ,21
+ ,67
+ ,45097
+ ,12
+ ,131069
+ ,67
+ ,47
+ ,30
+ ,83
+ ,51513
+ ,12
+ ,65475
+ ,18
+ ,16
+ ,13
+ ,46
+ ,69008
+ ,12
+ ,237213
+ ,84
+ ,78
+ ,38
+ ,123
+ ,66198
+ ,12
+ ,324598
+ ,110
+ ,113
+ ,37
+ ,133
+ ,135777
+ ,12
+ ,170266
+ ,62
+ ,44
+ ,42
+ ,122
+ ,74007
+ ,12
+ ,269651
+ ,67
+ ,106
+ ,30
+ ,93
+ ,70106
+ ,12
+ ,243060
+ ,63
+ ,58
+ ,29
+ ,104
+ ,111813
+ ,12
+ ,149112
+ ,56
+ ,56
+ ,35
+ ,95
+ ,60578
+ ,12
+ ,174415
+ ,100
+ ,73
+ ,31
+ ,114
+ ,82753
+ ,12
+ ,133131
+ ,55
+ ,44
+ ,30
+ ,90
+ ,57793
+ ,12
+ ,133328
+ ,55
+ ,56
+ ,20
+ ,56
+ ,40909
+ ,12
+ ,56653
+ ,45
+ ,38
+ ,18
+ ,49
+ ,27717
+ ,12
+ ,50857
+ ,21
+ ,15
+ ,20
+ ,37
+ ,36944
+ ,12
+ ,74408
+ ,67
+ ,29
+ ,7
+ ,21
+ ,34988
+ ,12
+ ,193339
+ ,78
+ ,100
+ ,35
+ ,71
+ ,84651
+ ,12
+ ,275541
+ ,63
+ ,116
+ ,33
+ ,115
+ ,93133
+ ,12
+ ,96560
+ ,76
+ ,42
+ ,17
+ ,38
+ ,22996
+ ,12
+ ,243199
+ ,75
+ ,88
+ ,28
+ ,105
+ ,95536
+ ,12
+ ,76702
+ ,49
+ ,35
+ ,21
+ ,42
+ ,49303
+ ,12
+ ,272458
+ ,65
+ ,100
+ ,43
+ ,152
+ ,93815
+ ,12
+ ,103425
+ ,67
+ ,17
+ ,19
+ ,46
+ ,25659
+ ,12
+ ,271856
+ ,103
+ ,109
+ ,37
+ ,86
+ ,54990
+ ,12
+ ,139526
+ ,151
+ ,28
+ ,21
+ ,77
+ ,48259
+ ,12
+ ,172494
+ ,52
+ ,46
+ ,43
+ ,139
+ ,86687
+ ,12
+ ,120445
+ ,118
+ ,36
+ ,16
+ ,56
+ ,51009
+ ,12
+ ,181528
+ ,54
+ ,32
+ ,16
+ ,60
+ ,40662
+ ,12
+ ,236785
+ ,119
+ ,77
+ ,31
+ ,71
+ ,101481
+ ,12
+ ,224330
+ ,83
+ ,131
+ ,39
+ ,132
+ ,130115
+ ,12
+ ,202925
+ ,61
+ ,115
+ ,36
+ ,134
+ ,114789
+ ,12
+ ,30837
+ ,19
+ ,8
+ ,4
+ ,4
+ ,2325
+ ,12
+ ,329267
+ ,259
+ ,79
+ ,39
+ ,140
+ ,115762
+ ,12
+ ,125930
+ ,75
+ ,37
+ ,17
+ ,65
+ ,32387
+ ,12
+ ,204271
+ ,42
+ ,92
+ ,35
+ ,116
+ ,74163
+ ,12
+ ,51567
+ ,30
+ ,21
+ ,14
+ ,45
+ ,28207
+ ,12
+ ,325107
+ ,99
+ ,84
+ ,36
+ ,126
+ ,79215
+ ,12
+ ,119016
+ ,52
+ ,118
+ ,23
+ ,74
+ ,78664
+ ,12
+ ,265769
+ ,146
+ ,96
+ ,32
+ ,120
+ ,83122)
+ ,dim=c(7
+ ,289)
+ ,dimnames=list(c('Maand'
+ ,'Time_in_rfc'
+ ,'Logins'
+ ,'Blogged_computations'
+ ,'Compendiums_reviewed'
+ ,'Feedback_messages_p120'
+ ,'Total_number_of_characters')
+ ,1:289))
> y <- array(NA,dim=c(7,289),dimnames=list(c('Maand','Time_in_rfc','Logins','Blogged_computations','Compendiums_reviewed','Feedback_messages_p120','Total_number_of_characters'),1:289))
> for (i in 1:dim(x)[1])
+ {
+ for (j in 1:dim(x)[2])
+ {
+ y[i,j] <- as.numeric(x[i,j])
+ }
+ }
> par3 = 'Linear Trend'
> par2 = 'Do not include Seasonal Dummies'
> par1 = '2'
> library(lattice)
> library(lmtest)
Loading required package: zoo
> n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test
> par1 <- as.numeric(par1)
> x <- t(y)
> k <- length(x[1,])
> n <- length(x[,1])
> x1 <- cbind(x[,par1], x[,1:k!=par1])
> mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1])
> colnames(x1) <- mycolnames #colnames(x)[par1]
> x <- x1
> if (par3 == 'First Differences'){
+ x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep='')))
+ for (i in 1:n-1) {
+ for (j in 1:k) {
+ x2[i,j] <- x[i+1,j] - x[i,j]
+ }
+ }
+ x <- x2
+ }
> if (par2 == 'Include Monthly Dummies'){
+ x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep ='')))
+ for (i in 1:11){
+ x2[seq(i,n,12),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> if (par2 == 'Include Quarterly Dummies'){
+ x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep ='')))
+ for (i in 1:3){
+ x2[seq(i,n,4),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> k <- length(x[1,])
> if (par3 == 'Linear Trend'){
+ x <- cbind(x, c(1:n))
+ colnames(x)[k+1] <- 't'
+ }
> x
Time_in_rfc Maand Logins Blogged_computations Compendiums_reviewed
1 14688 10 10 4 0
2 46660 10 20 12 5
3 7199 10 5 7 0
4 21054 10 16 4 0
5 80953 10 25 49 8
6 19349 10 11 13 1
7 173260 11 63 41 21
8 38214 11 34 16 8
9 17547 11 5 0 1
10 43287 11 14 13 19
11 78800 11 42 20 26
12 65029 11 17 21 18
13 31414 11 19 18 8
14 7176 11 17 0 0
15 24188 11 24 8 4
16 61857 11 25 23 11
17 22938 11 10 12 1
18 43410 11 19 7 1
19 19764 11 12 10 2
20 53117 11 22 3 12
21 33170 11 18 1 18
22 232138 11 62 131 31
23 111665 11 34 39 28
24 152871 11 58 59 28
25 100750 11 72 83 30
26 221698 11 45 105 33
27 209641 11 42 62 24
28 52746 11 25 0 26
29 215641 11 46 71 32
30 85439 11 33 32 28
31 145790 11 63 30 26
32 132943 11 40 83 39
33 132487 11 41 71 36
34 135473 11 41 82 23
35 84853 11 31 38 30
36 210767 11 60 94 35
37 97839 11 38 25 24
38 120221 11 37 53 22
39 224549 11 50 54 31
40 116408 11 61 34 39
41 158015 11 29 59 31
42 150629 12 44 85 33
43 123185 12 40 49 22
44 124817 12 40 47 25
45 104389 12 45 135 41
46 162765 12 32 68 28
47 149061 12 44 43 26
48 84207 12 29 14 30
49 201940 12 38 109 31
50 168809 12 66 76 28
51 176508 12 54 60 38
52 99466 12 50 28 23
53 106408 12 30 33 14
54 344297 12 75 80 30
55 218946 12 41 76 29
56 122774 12 45 24 24
57 153935 12 33 50 25
58 206161 12 71 75 28
59 183167 12 66 91 39
60 144966 12 144 39 32
61 34662 12 25 18 6
62 116048 12 64 50 20
63 77945 12 20 28 19
64 65490 12 22 27 16
65 81872 12 45 32 16
66 101011 12 34 30 13
67 71965 12 35 32 15
68 72880 12 33 14 19
69 62792 12 35 28 15
70 31774 12 23 0 17
71 182192 12 52 77 40
72 67989 12 21 23 18
73 52915 12 14 20 18
74 56375 12 30 10 13
75 30989 12 14 5 17
76 135131 12 66 38 15
77 76990 12 27 42 17
78 64175 12 42 37 18
79 59194 12 31 7 24
80 53515 12 28 5 22
81 89746 12 36 28 18
82 49862 12 37 17 14
83 85574 12 34 37 21
84 220801 12 75 51 18
85 49289 12 19 15 24
86 135781 12 31 45 14
87 82316 12 32 27 10
88 133368 12 36 37 16
89 92499 12 32 25 18
90 37460 12 20 5 7
91 91735 12 35 7 18
92 187681 12 62 114 28
93 79619 12 43 42 11
94 98104 12 54 55 17
95 80949 12 17 11 16
96 31706 12 13 26 10
97 70344 12 28 20 16
98 99923 12 66 44 25
99 91899 12 35 18 15
100 235454 12 73 151 32
101 46455 12 20 22 17
102 60812 12 44 26 13
103 77648 12 47 31 16
104 101523 12 42 59 22
105 254488 12 83 120 39
106 103597 12 43 27 16
107 136084 12 30 27 13
108 62088 12 38 13 16
109 86680 12 31 37 14
110 38395 12 31 16 16
111 179321 12 89 108 30
112 40248 12 16 8 4
113 70551 12 31 23 15
114 84105 12 20 17 17
115 112611 12 41 46 20
116 73624 12 24 30 17
117 89806 12 42 27 16
118 98866 12 18 25 13
119 61254 12 46 36 14
120 99643 12 55 33 17
121 310839 12 92 130 24
122 181633 12 70 47 30
123 50090 12 20 16 17
124 92661 12 61 40 17
125 46698 12 45 14 13
126 61361 12 77 27 12
127 84337 12 26 38 14
128 72535 12 14 17 16
129 84856 12 29 29 17
130 124064 12 40 43 22
131 54157 12 19 37 15
132 76302 12 31 29 20
133 59382 12 49 29 12
134 65745 12 53 26 21
135 74914 12 30 35 23
136 101645 12 63 20 11
137 73504 12 35 23 7
138 99373 12 63 12 18
139 135458 12 81 43 12
140 184510 12 49 64 40
141 385534 12 92 121 25
142 108043 12 62 45 14
143 328107 12 65 129 41
144 362301 12 119 76 34
145 130585 12 46 67 29
146 85709 12 44 21 10
147 174184 12 53 72 25
148 81437 12 38 37 14
149 120982 12 56 58 28
150 187559 12 121 75 36
151 235800 12 94 105 23
152 220516 12 62 98 32
153 210907 12 56 79 30
154 155754 12 61 74 20
155 152299 12 53 62 33
156 79863 12 37 29 22
157 346485 12 90 118 38
158 135649 12 46 99 30
159 134019 12 53 32 18
160 167488 12 45 69 28
161 108446 12 60 22 17
162 27634 12 20 2 16
163 232317 12 54 118 33
164 22648 12 19 12 13
165 119308 12 30 32 16
166 194979 12 62 66 40
167 81240 12 66 58 17
168 121848 12 39 37 17
169 86678 12 40 12 15
170 73566 12 32 22 23
171 102010 12 53 28 13
172 151101 12 32 48 35
173 207176 12 70 56 32
174 165543 12 65 70 32
175 182999 12 88 73 34
176 256462 12 105 123 35
177 244052 12 68 101 44
178 63123 12 43 34 17
179 87186 12 54 28 15
180 182079 12 63 124 33
181 95227 12 34 37 32
182 165446 12 33 69 25
183 229242 12 247 63 31
184 167542 12 66 59 28
185 225060 12 93 93 41
186 164709 12 109 81 31
187 131698 12 65 60 19
188 56613 12 19 15 12
189 250579 12 83 130 38
190 44296 12 25 10 20
191 152474 12 65 106 32
192 86230 12 44 21 17
193 258873 12 60 104 40
194 41566 12 35 0 9
195 250047 12 81 41 18
196 243511 12 71 133 42
197 62215 12 27 24 10
198 104838 12 49 46 16
199 77272 12 59 21 16
200 148446 12 91 135 37
201 69304 12 30 40 19
202 102538 12 57 50 15
203 101097 12 64 30 14
204 223632 12 73 105 33
205 141722 12 94 19 27
206 351067 12 95 136 45
207 311473 12 112 128 38
208 152601 12 48 46 24
209 215147 12 58 101 36
210 52164 12 52 32 16
211 91005 12 29 29 11
212 265318 12 117 110 52
213 299775 12 95 91 31
214 241066 12 82 75 45
215 118612 12 46 54 12
216 173326 12 88 86 44
217 225548 12 112 81 31
218 277965 12 89 115 39
219 317394 12 86 116 31
220 177939 12 82 55 36
221 92630 12 40 27 8
222 140344 12 53 33 25
223 230964 12 53 102 30
224 199476 12 70 87 32
225 174724 12 92 123 34
226 128423 12 64 32 38
227 341570 12 168 94 21
228 64187 12 27 10 16
229 150580 12 77 27 22
230 204713 12 71 68 33
231 244749 12 95 98 33
232 98146 12 40 15 17
233 351619 12 139 95 40
234 58981 12 36 0 23
235 233328 12 132 92 28
236 182613 12 39 81 28
237 89113 12 39 19 14
238 324799 12 154 158 47
239 143246 12 103 67 27
240 40151 12 29 16 9
241 158399 12 39 23 18
242 195838 12 67 111 31
243 286468 12 144 57 29
244 175824 12 107 57 20
245 143756 12 46 105 34
246 139942 12 42 54 22
247 104011 12 55 25 22
248 196553 12 57 41 29
249 180083 12 66 63 34
250 260561 12 75 114 43
251 294424 12 77 107 33
252 99611 12 35 41 21
253 131069 12 67 47 30
254 65475 12 18 16 13
255 237213 12 84 78 38
256 324598 12 110 113 37
257 170266 12 62 44 42
258 269651 12 67 106 30
259 243060 12 63 58 29
260 149112 12 56 56 35
261 174415 12 100 73 31
262 133131 12 55 44 30
263 133328 12 55 56 20
264 56653 12 45 38 18
265 50857 12 21 15 20
266 74408 12 67 29 7
267 193339 12 78 100 35
268 275541 12 63 116 33
269 96560 12 76 42 17
270 243199 12 75 88 28
271 76702 12 49 35 21
272 272458 12 65 100 43
273 103425 12 67 17 19
274 271856 12 103 109 37
275 139526 12 151 28 21
276 172494 12 52 46 43
277 120445 12 118 36 16
278 181528 12 54 32 16
279 236785 12 119 77 31
280 224330 12 83 131 39
281 202925 12 61 115 36
282 30837 12 19 8 4
283 329267 12 259 79 39
284 125930 12 75 37 17
285 204271 12 42 92 35
286 51567 12 30 21 14
287 325107 12 99 84 36
288 119016 12 52 118 23
289 265769 12 146 96 32
Feedback_messages_p120 Total_number_of_characters t
1 0 6023 1
2 13 6179 2
3 0 1644 3
4 0 855 4
5 27 56622 5
6 0 3895 6
7 78 37238 7
8 21 8773 8
9 4 3926 9
10 64 43750 10
11 66 33032 11
12 61 32551 12
13 9 14116 13
14 0 1423 14
15 7 5950 15
16 30 25162 16
17 0 1168 17
18 3 2781 18
19 4 5752 19
20 32 32689 20
21 22 22807 21
22 90 161647 22
23 80 49810 23
24 90 79892 24
25 93 140867 25
26 110 77873 26
27 60 97068 27
28 70 68504 28
29 92 55813 29
30 78 56926 30
31 99 38692 31
32 156 97500 32
33 98 40735 33
34 66 99645 34
35 109 61542 35
36 133 117478 36
37 66 94785 37
38 69 43836 38
39 114 120662 39
40 31 71570 40
41 114 71220 41
42 113 71595 42
43 51 38361 43
44 81 30727 44
45 148 123969 45
46 107 120293 46
47 93 116174 47
48 101 111194 48
49 92 85872 49
50 100 72260 50
51 93 83123 51
52 69 192565 52
53 37 31081 53
54 108 67654 54
55 96 80670 55
56 69 15986 56
57 84 65553 57
58 99 73107 58
59 138 82875 59
60 50 37510 60
61 7 36278 61
62 50 33416 62
63 67 25272 63
64 40 18653 64
65 61 27913 65
66 39 27570 66
67 56 46090 67
68 67 24094 68
69 51 30884 69
70 47 4143 70
71 138 70054 71
72 57 33747 72
73 49 24006 73
74 40 41385 74
75 41 4154 75
76 60 34029 76
77 67 24069 77
78 59 21792 78
79 68 37636 79
80 81 12934 80
81 55 24760 81
82 54 36341 82
83 36 24266 83
84 63 44418 84
85 44 10672 85
86 44 24874 86
87 38 21233 87
88 57 32755 88
89 55 21399 89
90 14 3738 90
91 49 16563 91
92 91 92945 92
93 32 20760 93
94 47 25568 94
95 56 18472 95
96 25 25139 96
97 61 28904 97
98 80 32334 98
99 33 36874 99
100 106 89275 100
101 12 11342 101
102 43 22574 102
103 52 21280 103
104 87 61056 104
105 117 103772 105
106 49 23789 106
107 41 27142 107
108 38 14483 108
109 49 22197 109
110 43 16380 110
111 103 101193 111
112 7 5444 112
113 47 30143 113
114 63 20055 114
115 73 26706 115
116 40 19499 116
117 40 27975 117
118 39 23517 118
119 30 13310 119
120 46 54968 120
121 85 162901 121
122 73 64466 122
123 36 21152 123
124 41 35232 124
125 33 10288 125
126 25 13294 126
127 35 24548 127
128 39 19540 128
129 43 41369 129
130 86 22827 130
131 30 18625 131
132 68 30976 132
133 24 26263 133
134 42 38084 134
135 79 65892 135
136 38 17140 136
137 3 35130 137
138 62 28394 138
139 29 19474 139
140 140 110681 140
141 91 119182 141
142 38 34553 142
143 144 105547 143
144 110 100708 144
145 107 95364 145
146 35 35944 146
147 70 55183 147
148 40 28579 148
149 103 84786 149
150 116 73511 150
151 62 149193 151
152 119 55801 152
153 94 112285 153
154 31 40652 154
155 98 61370 155
156 49 29011 156
157 111 106117 157
158 108 62133 158
159 54 28987 159
160 72 83737 160
161 65 46300 161
162 17 5839 162
163 106 129838 163
164 28 44339 164
165 58 23686 165
166 151 74011 166
167 50 45549 167
168 52 30594 168
169 50 29156 169
170 67 22618 170
171 44 44332 171
172 124 83209 172
173 114 87011 173
174 90 51633 174
175 104 225920 175
176 124 105195 176
177 164 143558 177
178 60 39067 178
179 48 51776 179
180 118 102860 180
181 48 34777 181
182 60 57635 182
183 119 91721 183
184 70 55461 184
185 139 109825 185
186 63 87771 186
187 59 65622 187
188 35 8019 188
189 120 103487 189
190 38 18779 190
191 83 65567 191
192 54 18632 192
193 124 80444 193
194 27 13497 194
195 61 83305 195
196 110 101338 196
197 25 19546 197
198 52 29236 198
199 48 27507 199
200 129 126846 200
201 40 22818 201
202 45 45833 202
203 41 31701 203
204 120 72654 204
205 36 75345 205
206 160 213688 206
207 129 132068 207
208 46 73224 208
209 117 99052 209
210 43 16734 210
211 23 32073 211
212 197 115929 212
213 97 98952 213
214 93 67267 214
215 43 30080 215
216 148 71701 216
217 94 85323 217
218 133 139077 218
219 82 91413 219
220 73 46821 220
221 24 18513 221
222 73 51715 222
223 115 133824 223
224 105 135400 224
225 120 69112 225
226 120 92696 226
227 78 117105 227
228 54 20154 228
229 71 45588 229
230 69 41140 230
231 115 76643 231
232 48 27114 232
233 141 115168 233
234 61 42744 234
235 102 120733 235
236 99 138599 236
237 32 82206 237
238 168 97668 238
239 104 106671 239
240 19 13127 240
241 57 32928 241
242 98 102372 242
243 107 89691 243
244 77 39992 244
245 120 69094 245
246 87 65461 246
247 79 63339 247
248 99 95260 248
249 70 53855 249
250 158 90183 250
251 124 101494 251
252 67 45097 252
253 83 51513 253
254 46 69008 254
255 123 66198 255
256 133 135777 256
257 122 74007 257
258 93 70106 258
259 104 111813 259
260 95 60578 260
261 114 82753 261
262 90 57793 262
263 56 40909 263
264 49 27717 264
265 37 36944 265
266 21 34988 266
267 71 84651 267
268 115 93133 268
269 38 22996 269
270 105 95536 270
271 42 49303 271
272 152 93815 272
273 46 25659 273
274 86 54990 274
275 77 48259 275
276 139 86687 276
277 56 51009 277
278 60 40662 278
279 71 101481 279
280 132 130115 280
281 134 114789 281
282 4 2325 282
283 140 115762 283
284 65 32387 284
285 116 74163 285
286 45 28207 286
287 126 79215 287
288 74 78664 288
289 120 83122 289
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) Maand
-5997.0124 585.5924
Logins Blogged_computations
725.8121 923.2352
Compendiums_reviewed Feedback_messages_p120
498.1369 265.2097
Total_number_of_characters t
0.2482 28.0229
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-145646 -18544 -2553 15235 153080
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -5.997e+03 7.273e+04 -0.082 0.93435
Maand 5.856e+02 6.386e+03 0.092 0.92700
Logins 7.258e+02 8.598e+01 8.442 1.68e-15 ***
Blogged_computations 9.232e+02 1.056e+02 8.745 < 2e-16 ***
Compendiums_reviewed 4.981e+02 5.338e+02 0.933 0.35149
Feedback_messages_p120 2.652e+02 1.498e+02 1.771 0.07772 .
Total_number_of_characters 2.482e-01 9.195e-02 2.699 0.00737 **
t 2.802e+01 3.473e+01 0.807 0.42047
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 36780 on 281 degrees of freedom
Multiple R-squared: 0.8053, Adjusted R-squared: 0.8005
F-statistic: 166 on 7 and 281 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,] 1.070793e-01 2.141586e-01 0.8929207
[2,] 3.888741e-02 7.777483e-02 0.9611126
[3,] 1.836929e-02 3.673859e-02 0.9816307
[4,] 6.382970e-03 1.276594e-02 0.9936170
[5,] 1.968840e-03 3.937679e-03 0.9980312
[6,] 7.178486e-04 1.435697e-03 0.9992822
[7,] 2.070586e-04 4.141173e-04 0.9997929
[8,] 1.251743e-04 2.503486e-04 0.9998748
[9,] 3.872526e-05 7.745051e-05 0.9999613
[10,] 3.368473e-05 6.736946e-05 0.9999663
[11,] 3.747388e-05 7.494777e-05 0.9999625
[12,] 1.195510e-05 2.391021e-05 0.9999880
[13,] 4.298353e-06 8.596706e-06 0.9999957
[14,] 4.635105e-06 9.270211e-06 0.9999954
[15,] 1.677985e-03 3.355971e-03 0.9983220
[16,] 1.167861e-03 2.335723e-03 0.9988321
[17,] 8.926185e-02 1.785237e-01 0.9107381
[18,] 6.451885e-02 1.290377e-01 0.9354812
[19,] 5.577483e-02 1.115497e-01 0.9442252
[20,] 5.236767e-02 1.047353e-01 0.9476323
[21,] 3.716812e-02 7.433624e-02 0.9628319
[22,] 6.425261e-02 1.285052e-01 0.9357474
[23,] 1.081028e-01 2.162056e-01 0.8918972
[24,] 8.608288e-02 1.721658e-01 0.9139171
[25,] 6.852569e-02 1.370514e-01 0.9314743
[26,] 5.730937e-02 1.146187e-01 0.9426906
[27,] 4.796031e-02 9.592063e-02 0.9520397
[28,] 3.539095e-02 7.078190e-02 0.9646091
[29,] 1.356792e-01 2.713584e-01 0.8643208
[30,] 1.121769e-01 2.243539e-01 0.8878231
[31,] 9.322993e-02 1.864599e-01 0.9067701
[32,] 8.427442e-02 1.685488e-01 0.9157256
[33,] 6.583602e-02 1.316720e-01 0.9341640
[34,] 5.029258e-02 1.005852e-01 0.9497074
[35,] 3.398813e-01 6.797626e-01 0.6601187
[36,] 3.423093e-01 6.846186e-01 0.6576907
[37,] 3.062270e-01 6.124540e-01 0.6937730
[38,] 2.700745e-01 5.401490e-01 0.7299255
[39,] 2.665830e-01 5.331660e-01 0.7334170
[40,] 2.405487e-01 4.810974e-01 0.7594513
[41,] 2.124502e-01 4.249005e-01 0.7875498
[42,] 2.038256e-01 4.076513e-01 0.7961744
[43,] 1.762348e-01 3.524697e-01 0.8237652
[44,] 6.118894e-01 7.762212e-01 0.3881106
[45,] 6.431736e-01 7.136527e-01 0.3568264
[46,] 6.244410e-01 7.511179e-01 0.3755590
[47,] 5.899290e-01 8.201419e-01 0.4100710
[48,] 5.716608e-01 8.566784e-01 0.4283392
[49,] 6.118212e-01 7.763575e-01 0.3881788
[50,] 8.072425e-01 3.855151e-01 0.1927575
[51,] 8.135053e-01 3.729895e-01 0.1864947
[52,] 8.052088e-01 3.895825e-01 0.1947912
[53,] 7.833569e-01 4.332863e-01 0.2166431
[54,] 7.586754e-01 4.826493e-01 0.2413246
[55,] 7.503413e-01 4.993173e-01 0.2496587
[56,] 7.198922e-01 5.602156e-01 0.2801078
[57,] 7.067634e-01 5.864733e-01 0.2932366
[58,] 6.754956e-01 6.490088e-01 0.3245044
[59,] 6.602218e-01 6.795564e-01 0.3397782
[60,] 6.322475e-01 7.355049e-01 0.3677525
[61,] 5.945214e-01 8.109573e-01 0.4054786
[62,] 5.560044e-01 8.879913e-01 0.4439956
[63,] 5.178327e-01 9.643345e-01 0.4821673
[64,] 4.790074e-01 9.580149e-01 0.5209926
[65,] 4.423345e-01 8.846691e-01 0.5576655
[66,] 4.061451e-01 8.122902e-01 0.5938549
[67,] 3.797922e-01 7.595843e-01 0.6202078
[68,] 3.807523e-01 7.615045e-01 0.6192477
[69,] 3.470522e-01 6.941043e-01 0.6529478
[70,] 3.155069e-01 6.310139e-01 0.6844931
[71,] 2.821825e-01 5.643651e-01 0.7178175
[72,] 2.649815e-01 5.299631e-01 0.7350185
[73,] 2.344489e-01 4.688978e-01 0.7655511
[74,] 3.330874e-01 6.661748e-01 0.6669126
[75,] 2.995845e-01 5.991690e-01 0.7004155
[76,] 3.091381e-01 6.182762e-01 0.6908619
[77,] 2.774900e-01 5.549800e-01 0.7225100
[78,] 2.734173e-01 5.468346e-01 0.7265827
[79,] 2.450575e-01 4.901150e-01 0.7549425
[80,] 2.168716e-01 4.337433e-01 0.7831284
[81,] 2.029484e-01 4.058967e-01 0.7970516
[82,] 1.945915e-01 3.891830e-01 0.8054085
[83,] 1.780888e-01 3.561776e-01 0.8219112
[84,] 1.698197e-01 3.396393e-01 0.8301803
[85,] 1.588891e-01 3.177782e-01 0.8411109
[86,] 1.443607e-01 2.887213e-01 0.8556393
[87,] 1.248573e-01 2.497147e-01 0.8751427
[88,] 1.300063e-01 2.600125e-01 0.8699937
[89,] 1.182519e-01 2.365038e-01 0.8817481
[90,] 1.072219e-01 2.144438e-01 0.8927781
[91,] 9.140201e-02 1.828040e-01 0.9085980
[92,] 8.384396e-02 1.676879e-01 0.9161560
[93,] 7.370535e-02 1.474107e-01 0.9262946
[94,] 7.116335e-02 1.423267e-01 0.9288367
[95,] 6.126315e-02 1.225263e-01 0.9387368
[96,] 5.311610e-02 1.062322e-01 0.9468839
[97,] 7.624082e-02 1.524816e-01 0.9237592
[98,] 6.436560e-02 1.287312e-01 0.9356344
[99,] 5.380371e-02 1.076074e-01 0.9461963
[100,] 4.934258e-02 9.868516e-02 0.9506574
[101,] 6.027359e-02 1.205472e-01 0.9397264
[102,] 5.166252e-02 1.033250e-01 0.9483375
[103,] 4.289873e-02 8.579746e-02 0.9571013
[104,] 3.785777e-02 7.571555e-02 0.9621422
[105,] 3.108215e-02 6.216430e-02 0.9689178
[106,] 2.539989e-02 5.079978e-02 0.9746001
[107,] 2.069971e-02 4.139942e-02 0.9793003
[108,] 2.139061e-02 4.278122e-02 0.9786094
[109,] 1.964565e-02 3.929130e-02 0.9803544
[110,] 1.592699e-02 3.185398e-02 0.9840730
[111,] 2.124947e-02 4.249893e-02 0.9787505
[112,] 2.059682e-02 4.119364e-02 0.9794032
[113,] 1.667578e-02 3.335156e-02 0.9833242
[114,] 1.465350e-02 2.930699e-02 0.9853465
[115,] 1.289100e-02 2.578199e-02 0.9871090
[116,] 1.426649e-02 2.853298e-02 0.9857335
[117,] 1.145836e-02 2.291673e-02 0.9885416
[118,] 9.867756e-03 1.973551e-02 0.9901322
[119,] 7.826102e-03 1.565220e-02 0.9921739
[120,] 6.299745e-03 1.259949e-02 0.9937003
[121,] 5.154514e-03 1.030903e-02 0.9948455
[122,] 4.106555e-03 8.213110e-03 0.9958934
[123,] 3.608119e-03 7.216237e-03 0.9963919
[124,] 3.337765e-03 6.675530e-03 0.9966622
[125,] 3.033791e-03 6.067581e-03 0.9969662
[126,] 2.394575e-03 4.789149e-03 0.9976054
[127,] 1.885167e-03 3.770335e-03 0.9981148
[128,] 1.432483e-03 2.864967e-03 0.9985675
[129,] 1.115557e-03 2.231114e-03 0.9988844
[130,] 8.662463e-04 1.732493e-03 0.9991338
[131,] 2.252983e-02 4.505966e-02 0.9774702
[132,] 1.873482e-02 3.746963e-02 0.9812652
[133,] 3.316767e-02 6.633534e-02 0.9668323
[134,] 1.616717e-01 3.233433e-01 0.8383283
[135,] 1.630641e-01 3.261282e-01 0.8369359
[136,] 1.442090e-01 2.884180e-01 0.8557910
[137,] 1.310761e-01 2.621523e-01 0.8689239
[138,] 1.153573e-01 2.307147e-01 0.8846427
[139,] 1.217306e-01 2.434613e-01 0.8782694
[140,] 1.391801e-01 2.783603e-01 0.8608199
[141,] 1.214366e-01 2.428731e-01 0.8785634
[142,] 1.098975e-01 2.197950e-01 0.8901025
[143,] 1.014260e-01 2.028520e-01 0.8985740
[144,] 8.876944e-02 1.775389e-01 0.9112306
[145,] 7.604981e-02 1.520996e-01 0.9239502
[146,] 6.502752e-02 1.300550e-01 0.9349725
[147,] 1.421331e-01 2.842662e-01 0.8578669
[148,] 1.619315e-01 3.238630e-01 0.8380685
[149,] 1.564370e-01 3.128740e-01 0.8435630
[150,] 1.396665e-01 2.793330e-01 0.8603335
[151,] 1.223884e-01 2.447768e-01 0.8776116
[152,] 1.059947e-01 2.119895e-01 0.8940053
[153,] 9.152422e-02 1.830484e-01 0.9084758
[154,] 8.624493e-02 1.724899e-01 0.9137551
[155,] 8.613403e-02 1.722681e-01 0.9138660
[156,] 7.365054e-02 1.473011e-01 0.9263495
[157,] 9.371298e-02 1.874260e-01 0.9062870
[158,] 8.745704e-02 1.749141e-01 0.9125430
[159,] 7.675593e-02 1.535119e-01 0.9232441
[160,] 6.547927e-02 1.309585e-01 0.9345207
[161,] 5.583106e-02 1.116621e-01 0.9441689
[162,] 4.723120e-02 9.446240e-02 0.9527688
[163,] 4.487575e-02 8.975150e-02 0.9551242
[164,] 3.756475e-02 7.512949e-02 0.9624353
[165,] 5.456187e-02 1.091237e-01 0.9454381
[166,] 4.851862e-02 9.703723e-02 0.9514814
[167,] 4.063376e-02 8.126752e-02 0.9593662
[168,] 4.011481e-02 8.022962e-02 0.9598852
[169,] 3.426253e-02 6.852507e-02 0.9657375
[170,] 4.218296e-02 8.436592e-02 0.9578170
[171,] 3.538789e-02 7.077579e-02 0.9646121
[172,] 3.390583e-02 6.781167e-02 0.9660942
[173,] 7.307196e-02 1.461439e-01 0.9269280
[174,] 6.307378e-02 1.261476e-01 0.9369262
[175,] 5.571741e-02 1.114348e-01 0.9442826
[176,] 6.555364e-02 1.311073e-01 0.9344464
[177,] 5.719280e-02 1.143856e-01 0.9428072
[178,] 4.828963e-02 9.657927e-02 0.9517104
[179,] 4.066869e-02 8.133739e-02 0.9593313
[180,] 3.451793e-02 6.903586e-02 0.9654821
[181,] 4.098342e-02 8.196684e-02 0.9590166
[182,] 3.374205e-02 6.748410e-02 0.9662580
[183,] 3.511881e-02 7.023762e-02 0.9648812
[184,] 2.871116e-02 5.742232e-02 0.9712888
[185,] 8.677542e-02 1.735508e-01 0.9132246
[186,] 7.497472e-02 1.499494e-01 0.9250253
[187,] 6.306305e-02 1.261261e-01 0.9369369
[188,] 5.278232e-02 1.055646e-01 0.9472177
[189,] 4.502221e-02 9.004443e-02 0.9549778
[190,] 2.748935e-01 5.497871e-01 0.7251065
[191,] 2.544590e-01 5.089181e-01 0.7455410
[192,] 2.362317e-01 4.724635e-01 0.7637683
[193,] 2.090479e-01 4.180958e-01 0.7909521
[194,] 1.831434e-01 3.662868e-01 0.8168566
[195,] 1.664815e-01 3.329630e-01 0.8335185
[196,] 1.538211e-01 3.076422e-01 0.8461789
[197,] 1.355704e-01 2.711408e-01 0.8644296
[198,] 1.208314e-01 2.416627e-01 0.8791686
[199,] 1.037097e-01 2.074195e-01 0.8962903
[200,] 1.122263e-01 2.244526e-01 0.8877737
[201,] 9.683632e-02 1.936726e-01 0.9031637
[202,] 1.000790e-01 2.001580e-01 0.8999210
[203,] 1.404193e-01 2.808387e-01 0.8595807
[204,] 1.346254e-01 2.692508e-01 0.8653746
[205,] 1.148658e-01 2.297316e-01 0.8851342
[206,] 1.548090e-01 3.096179e-01 0.8451910
[207,] 1.340137e-01 2.680274e-01 0.8659863
[208,] 1.150074e-01 2.300148e-01 0.8849926
[209,] 1.907916e-01 3.815833e-01 0.8092084
[210,] 1.654880e-01 3.309759e-01 0.8345120
[211,] 1.491809e-01 2.983617e-01 0.8508191
[212,] 1.303293e-01 2.606587e-01 0.8696707
[213,] 1.109385e-01 2.218770e-01 0.8890615
[214,] 9.936180e-02 1.987236e-01 0.9006382
[215,] 1.800852e-01 3.601704e-01 0.8199148
[216,] 2.025702e-01 4.051403e-01 0.7974298
[217,] 2.847234e-01 5.694467e-01 0.7152766
[218,] 2.513867e-01 5.027733e-01 0.7486133
[219,] 2.233600e-01 4.467200e-01 0.7766400
[220,] 2.190205e-01 4.380411e-01 0.7809795
[221,] 1.897837e-01 3.795673e-01 0.8102163
[222,] 1.682612e-01 3.365224e-01 0.8317388
[223,] 2.203229e-01 4.406458e-01 0.7796771
[224,] 1.994521e-01 3.989042e-01 0.8005479
[225,] 1.769698e-01 3.539396e-01 0.8230302
[226,] 1.505386e-01 3.010771e-01 0.8494614
[227,] 1.257698e-01 2.515397e-01 0.8742302
[228,] 1.236103e-01 2.472206e-01 0.8763897
[229,] 2.115503e-01 4.231005e-01 0.7884497
[230,] 1.819087e-01 3.638173e-01 0.8180913
[231,] 2.494903e-01 4.989805e-01 0.7505097
[232,] 2.443885e-01 4.887770e-01 0.7556115
[233,] 2.671497e-01 5.342995e-01 0.7328503
[234,] 2.294951e-01 4.589901e-01 0.7705049
[235,] 3.564464e-01 7.128929e-01 0.6435536
[236,] 3.161209e-01 6.322418e-01 0.6838791
[237,] 2.933231e-01 5.866463e-01 0.7066769
[238,] 2.778773e-01 5.557546e-01 0.7221227
[239,] 2.432295e-01 4.864589e-01 0.7567705
[240,] 2.337302e-01 4.674603e-01 0.7662698
[241,] 2.298445e-01 4.596891e-01 0.7701555
[242,] 1.994500e-01 3.989001e-01 0.8005500
[243,] 1.829561e-01 3.659122e-01 0.8170439
[244,] 1.490389e-01 2.980778e-01 0.8509611
[245,] 1.207819e-01 2.415638e-01 0.8792181
[246,] 1.147468e-01 2.294936e-01 0.8852532
[247,] 9.650634e-02 1.930127e-01 0.9034937
[248,] 1.248189e-01 2.496378e-01 0.8751811
[249,] 2.474367e-01 4.948735e-01 0.7525633
[250,] 2.112191e-01 4.224383e-01 0.7887809
[251,] 2.080722e-01 4.161443e-01 0.7919278
[252,] 1.772923e-01 3.545846e-01 0.8227077
[253,] 1.379090e-01 2.758180e-01 0.8620910
[254,] 1.727641e-01 3.455281e-01 0.8272359
[255,] 1.375973e-01 2.751946e-01 0.8624027
[256,] 1.038914e-01 2.077828e-01 0.8961086
[257,] 7.716385e-02 1.543277e-01 0.9228361
[258,] 7.840829e-02 1.568166e-01 0.9215917
[259,] 6.726738e-02 1.345348e-01 0.9327326
[260,] 1.139384e-01 2.278768e-01 0.8860616
[261,] 8.823702e-02 1.764740e-01 0.9117630
[262,] 7.260190e-02 1.452038e-01 0.9273981
[263,] 4.745579e-02 9.491159e-02 0.9525442
[264,] 2.850618e-02 5.701236e-02 0.9714938
[265,] 3.454739e-02 6.909477e-02 0.9654526
[266,] 1.456245e-01 2.912490e-01 0.8543755
[267,] 1.425141e-01 2.850282e-01 0.8574859
[268,] 2.431916e-01 4.863831e-01 0.7568084
> postscript(file="/var/wessaorg/rcomp/tmp/1hfmo1323954524.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/2xx4e1323954524.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/3h5d91323954524.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/4okrc1323954524.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/5lz3f1323954524.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
2355.0288 13677.8715 -3243.7456 5564.8420 -7629.8150 -2128.9579
7 8 9 10 11 12
48650.4142 -13636.1166 10687.7774 -16898.4159 -9555.8765 -702.2978
13 14 15 16 17 18
-9679.1943 -6352.8345 -6808.0907 1903.1345 2892.1168 20223.9669
19 20 21 22 23 24
-2639.8582 10796.3416 -2312.8652 -14300.8541 2364.1832 -2460.8157
25 26 27 28 29 30
-103855.4234 25983.1382 68753.9708 -15148.1371 61253.7002 -18105.3690
31 32 33 34 35 36
22242.4642 -59059.8719 -18224.6106 -25081.8417 -33282.7875 -2886.1817
37 38 39 40 41 42
-7289.9013 2786.0526 61240.5056 -6235.2594 17548.3892 -26166.6653
43 44 45 46 47 48
12672.5805 8567.1686 -145645.5446 2256.5480 8627.2496 -21471.9062
49 50 51 52 53 54
10167.3725 -10096.2483 15235.0193 -42715.8559 27150.2432 153079.7188
55 56 57 58 59 60
56521.1877 31132.7744 30192.0237 24380.7181 -28031.2173 -36779.8506
61 62 63 64 65 66
-16690.9177 -10850.6935 1276.1932 -1437.2367 -14260.9528 22094.5891
67 68 69 70 71 72
-19653.4238 -146.7942 -20089.1833 -9872.8997 -3571.6530 -3995.0035
73 74 75 76 77 78
-6707.1345 -5091.8876 -7293.3526 17153.6177 -16782.2326 -33707.5443
79 80 81 82 83 84
-12343.9437 -10347.2033 4767.5248 -26331.4982 -2650.8228 79196.3640
85 86 87 88 89 90
-8035.4208 43478.0319 10364.9721 38365.6553 13803.5330 6647.7187
91 92 93 94 95 96
30215.8906 -27328.2547 -13122.1008 -22811.4801 27355.4477 -23305.3541
97 98 99 100 101 102
-3514.0341 -34075.1426 20696.4727 -26982.8310 -6698.8771 -22499.1919
103 104 105 106 107 108
-16044.9402 -36563.4909 3270.6423 16589.0480 61267.4197 -3194.4860
109 110 111 112 113 114
456.7992 -26429.4573 -56504.2413 11880.1891 -4799.0026 19514.6426
115 116 117 118 119 120
179.3858 310.1202 4563.4845 35727.5684 -27968.5402 -9447.9087
121 122 123 124 125 126
44690.9095 32679.6206 -6940.9532 -21134.7774 -21203.0905 -39922.4663
127 128 129 130 131 132
3444.6048 18898.1939 2247.9926 11226.3705 -18545.5010 -13386.6968
133 134 135 136 137 138
-26575.1964 -32564.8722 -32750.6187 12801.1116 8994.6943 5213.5884
139 140 141 142 143 144
13540.4420 377.7720 135896.7115 -9140.6020 71982.7638 129591.5312
145 146 147 148 149 150
-36246.3361 6078.4838 19378.4327 -10156.9736 -40725.8804 -41683.9317
151 152 153 154 155 156
440.9151 18398.2658 24264.1481 9539.9235 -6445.0481 -10322.6732
157 158 159 160 161 162
92083.5781 -53605.2041 30039.0146 11782.0894 1845.1877 -8226.5762
163 164 165 166 167 168
1805.5966 -32754.1827 33104.7581 5020.7409 -58955.7357 23790.6840
169 170 171 172 173 174
12831.3486 -10605.3856 2720.6640 6735.5023 31018.5467 -4792.4641
175 176 177 178 179 180
-54797.0392 -15699.6796 -5585.2214 -39572.8579 -16958.1269 -57466.7917
181 182 183 184 185 186
-7015.0344 28989.0636 -84124.0614 12702.6719 -19063.2697 -49365.1105
187 188 189 190 191 192
-18544.3531 5425.2565 -12451.3836 -14137.9616 -53176.4233 1081.5176
193 194 195 196 197 198
40090.7323 -5297.9592 101087.5393 -12582.5485 -2553.3206 -8792.0328
199 200 201 202 203 204
-19073.3023 -133002.0205 -21799.1355 -22468.5157 -5486.8261 664.0490
205 206 207 208 209 210
7480.6969 31862.6186 19255.1157 26104.4063 -631.9995 -45564.4296
211 212 213 214 215 216
16699.5018 -35053.1265 74080.9117 41502.6981 3467.0784 -55992.7728
217 218 219 220 221 222
813.8630 10835.9131 80832.6136 11534.6135 16501.7526 19508.0542
223 224 225 226 227 228
12386.9522 -16354.7068 -78860.0683 -28698.3091 65244.1790 644.4465
229 230 231 232 233 234
21213.4329 37975.5976 11855.1450 19805.1940 69558.3294 -12980.1623
235 236 237 238 239 240
-25998.7983 -2724.6155 -271.7253 -32756.7095 -68605.3980 -16205.3857
241 242 243 244 245 246
68817.7955 -29924.8068 56401.2851 -2640.3245 -60378.4673 1399.0658
247 248 249 250 251 252
-14573.3137 45002.9792 17139.0493 7133.0257 57170.0093 -11160.4286
253 254 255 256 257 258
-18814.9604 -6313.2639 28075.4344 44823.8293 4764.6329 57888.9659
259 260 261 262 263 264
65717.1710 -9216.5379 -40122.8092 -8940.9487 -1661.4986 -48361.0254
265 266 267 268 269 270
-15635.1489 -27220.1859 -21386.0629 44125.0537 -30200.0084 33414.1774
271 272 273 274 275 276
-33637.5549 39286.6916 2387.1506 32868.1253 -47519.0458 3717.8309
277 278 279 280 281 282
-42712.6641 69994.6309 11014.9295 -52463.7184 -38389.0381 -2902.3382
283 284 285 286 287 288
-25904.5990 -5399.9761 13225.4431 -24549.5861 95615.7907 -87376.6075
289
-6356.2017
> postscript(file="/var/wessaorg/rcomp/tmp/6spt01323954524.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 2355.0288 NA
1 13677.8715 2355.0288
2 -3243.7456 13677.8715
3 5564.8420 -3243.7456
4 -7629.8150 5564.8420
5 -2128.9579 -7629.8150
6 48650.4142 -2128.9579
7 -13636.1166 48650.4142
8 10687.7774 -13636.1166
9 -16898.4159 10687.7774
10 -9555.8765 -16898.4159
11 -702.2978 -9555.8765
12 -9679.1943 -702.2978
13 -6352.8345 -9679.1943
14 -6808.0907 -6352.8345
15 1903.1345 -6808.0907
16 2892.1168 1903.1345
17 20223.9669 2892.1168
18 -2639.8582 20223.9669
19 10796.3416 -2639.8582
20 -2312.8652 10796.3416
21 -14300.8541 -2312.8652
22 2364.1832 -14300.8541
23 -2460.8157 2364.1832
24 -103855.4234 -2460.8157
25 25983.1382 -103855.4234
26 68753.9708 25983.1382
27 -15148.1371 68753.9708
28 61253.7002 -15148.1371
29 -18105.3690 61253.7002
30 22242.4642 -18105.3690
31 -59059.8719 22242.4642
32 -18224.6106 -59059.8719
33 -25081.8417 -18224.6106
34 -33282.7875 -25081.8417
35 -2886.1817 -33282.7875
36 -7289.9013 -2886.1817
37 2786.0526 -7289.9013
38 61240.5056 2786.0526
39 -6235.2594 61240.5056
40 17548.3892 -6235.2594
41 -26166.6653 17548.3892
42 12672.5805 -26166.6653
43 8567.1686 12672.5805
44 -145645.5446 8567.1686
45 2256.5480 -145645.5446
46 8627.2496 2256.5480
47 -21471.9062 8627.2496
48 10167.3725 -21471.9062
49 -10096.2483 10167.3725
50 15235.0193 -10096.2483
51 -42715.8559 15235.0193
52 27150.2432 -42715.8559
53 153079.7188 27150.2432
54 56521.1877 153079.7188
55 31132.7744 56521.1877
56 30192.0237 31132.7744
57 24380.7181 30192.0237
58 -28031.2173 24380.7181
59 -36779.8506 -28031.2173
60 -16690.9177 -36779.8506
61 -10850.6935 -16690.9177
62 1276.1932 -10850.6935
63 -1437.2367 1276.1932
64 -14260.9528 -1437.2367
65 22094.5891 -14260.9528
66 -19653.4238 22094.5891
67 -146.7942 -19653.4238
68 -20089.1833 -146.7942
69 -9872.8997 -20089.1833
70 -3571.6530 -9872.8997
71 -3995.0035 -3571.6530
72 -6707.1345 -3995.0035
73 -5091.8876 -6707.1345
74 -7293.3526 -5091.8876
75 17153.6177 -7293.3526
76 -16782.2326 17153.6177
77 -33707.5443 -16782.2326
78 -12343.9437 -33707.5443
79 -10347.2033 -12343.9437
80 4767.5248 -10347.2033
81 -26331.4982 4767.5248
82 -2650.8228 -26331.4982
83 79196.3640 -2650.8228
84 -8035.4208 79196.3640
85 43478.0319 -8035.4208
86 10364.9721 43478.0319
87 38365.6553 10364.9721
88 13803.5330 38365.6553
89 6647.7187 13803.5330
90 30215.8906 6647.7187
91 -27328.2547 30215.8906
92 -13122.1008 -27328.2547
93 -22811.4801 -13122.1008
94 27355.4477 -22811.4801
95 -23305.3541 27355.4477
96 -3514.0341 -23305.3541
97 -34075.1426 -3514.0341
98 20696.4727 -34075.1426
99 -26982.8310 20696.4727
100 -6698.8771 -26982.8310
101 -22499.1919 -6698.8771
102 -16044.9402 -22499.1919
103 -36563.4909 -16044.9402
104 3270.6423 -36563.4909
105 16589.0480 3270.6423
106 61267.4197 16589.0480
107 -3194.4860 61267.4197
108 456.7992 -3194.4860
109 -26429.4573 456.7992
110 -56504.2413 -26429.4573
111 11880.1891 -56504.2413
112 -4799.0026 11880.1891
113 19514.6426 -4799.0026
114 179.3858 19514.6426
115 310.1202 179.3858
116 4563.4845 310.1202
117 35727.5684 4563.4845
118 -27968.5402 35727.5684
119 -9447.9087 -27968.5402
120 44690.9095 -9447.9087
121 32679.6206 44690.9095
122 -6940.9532 32679.6206
123 -21134.7774 -6940.9532
124 -21203.0905 -21134.7774
125 -39922.4663 -21203.0905
126 3444.6048 -39922.4663
127 18898.1939 3444.6048
128 2247.9926 18898.1939
129 11226.3705 2247.9926
130 -18545.5010 11226.3705
131 -13386.6968 -18545.5010
132 -26575.1964 -13386.6968
133 -32564.8722 -26575.1964
134 -32750.6187 -32564.8722
135 12801.1116 -32750.6187
136 8994.6943 12801.1116
137 5213.5884 8994.6943
138 13540.4420 5213.5884
139 377.7720 13540.4420
140 135896.7115 377.7720
141 -9140.6020 135896.7115
142 71982.7638 -9140.6020
143 129591.5312 71982.7638
144 -36246.3361 129591.5312
145 6078.4838 -36246.3361
146 19378.4327 6078.4838
147 -10156.9736 19378.4327
148 -40725.8804 -10156.9736
149 -41683.9317 -40725.8804
150 440.9151 -41683.9317
151 18398.2658 440.9151
152 24264.1481 18398.2658
153 9539.9235 24264.1481
154 -6445.0481 9539.9235
155 -10322.6732 -6445.0481
156 92083.5781 -10322.6732
157 -53605.2041 92083.5781
158 30039.0146 -53605.2041
159 11782.0894 30039.0146
160 1845.1877 11782.0894
161 -8226.5762 1845.1877
162 1805.5966 -8226.5762
163 -32754.1827 1805.5966
164 33104.7581 -32754.1827
165 5020.7409 33104.7581
166 -58955.7357 5020.7409
167 23790.6840 -58955.7357
168 12831.3486 23790.6840
169 -10605.3856 12831.3486
170 2720.6640 -10605.3856
171 6735.5023 2720.6640
172 31018.5467 6735.5023
173 -4792.4641 31018.5467
174 -54797.0392 -4792.4641
175 -15699.6796 -54797.0392
176 -5585.2214 -15699.6796
177 -39572.8579 -5585.2214
178 -16958.1269 -39572.8579
179 -57466.7917 -16958.1269
180 -7015.0344 -57466.7917
181 28989.0636 -7015.0344
182 -84124.0614 28989.0636
183 12702.6719 -84124.0614
184 -19063.2697 12702.6719
185 -49365.1105 -19063.2697
186 -18544.3531 -49365.1105
187 5425.2565 -18544.3531
188 -12451.3836 5425.2565
189 -14137.9616 -12451.3836
190 -53176.4233 -14137.9616
191 1081.5176 -53176.4233
192 40090.7323 1081.5176
193 -5297.9592 40090.7323
194 101087.5393 -5297.9592
195 -12582.5485 101087.5393
196 -2553.3206 -12582.5485
197 -8792.0328 -2553.3206
198 -19073.3023 -8792.0328
199 -133002.0205 -19073.3023
200 -21799.1355 -133002.0205
201 -22468.5157 -21799.1355
202 -5486.8261 -22468.5157
203 664.0490 -5486.8261
204 7480.6969 664.0490
205 31862.6186 7480.6969
206 19255.1157 31862.6186
207 26104.4063 19255.1157
208 -631.9995 26104.4063
209 -45564.4296 -631.9995
210 16699.5018 -45564.4296
211 -35053.1265 16699.5018
212 74080.9117 -35053.1265
213 41502.6981 74080.9117
214 3467.0784 41502.6981
215 -55992.7728 3467.0784
216 813.8630 -55992.7728
217 10835.9131 813.8630
218 80832.6136 10835.9131
219 11534.6135 80832.6136
220 16501.7526 11534.6135
221 19508.0542 16501.7526
222 12386.9522 19508.0542
223 -16354.7068 12386.9522
224 -78860.0683 -16354.7068
225 -28698.3091 -78860.0683
226 65244.1790 -28698.3091
227 644.4465 65244.1790
228 21213.4329 644.4465
229 37975.5976 21213.4329
230 11855.1450 37975.5976
231 19805.1940 11855.1450
232 69558.3294 19805.1940
233 -12980.1623 69558.3294
234 -25998.7983 -12980.1623
235 -2724.6155 -25998.7983
236 -271.7253 -2724.6155
237 -32756.7095 -271.7253
238 -68605.3980 -32756.7095
239 -16205.3857 -68605.3980
240 68817.7955 -16205.3857
241 -29924.8068 68817.7955
242 56401.2851 -29924.8068
243 -2640.3245 56401.2851
244 -60378.4673 -2640.3245
245 1399.0658 -60378.4673
246 -14573.3137 1399.0658
247 45002.9792 -14573.3137
248 17139.0493 45002.9792
249 7133.0257 17139.0493
250 57170.0093 7133.0257
251 -11160.4286 57170.0093
252 -18814.9604 -11160.4286
253 -6313.2639 -18814.9604
254 28075.4344 -6313.2639
255 44823.8293 28075.4344
256 4764.6329 44823.8293
257 57888.9659 4764.6329
258 65717.1710 57888.9659
259 -9216.5379 65717.1710
260 -40122.8092 -9216.5379
261 -8940.9487 -40122.8092
262 -1661.4986 -8940.9487
263 -48361.0254 -1661.4986
264 -15635.1489 -48361.0254
265 -27220.1859 -15635.1489
266 -21386.0629 -27220.1859
267 44125.0537 -21386.0629
268 -30200.0084 44125.0537
269 33414.1774 -30200.0084
270 -33637.5549 33414.1774
271 39286.6916 -33637.5549
272 2387.1506 39286.6916
273 32868.1253 2387.1506
274 -47519.0458 32868.1253
275 3717.8309 -47519.0458
276 -42712.6641 3717.8309
277 69994.6309 -42712.6641
278 11014.9295 69994.6309
279 -52463.7184 11014.9295
280 -38389.0381 -52463.7184
281 -2902.3382 -38389.0381
282 -25904.5990 -2902.3382
283 -5399.9761 -25904.5990
284 13225.4431 -5399.9761
285 -24549.5861 13225.4431
286 95615.7907 -24549.5861
287 -87376.6075 95615.7907
288 -6356.2017 -87376.6075
289 NA -6356.2017
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] 13677.8715 2355.0288
[2,] -3243.7456 13677.8715
[3,] 5564.8420 -3243.7456
[4,] -7629.8150 5564.8420
[5,] -2128.9579 -7629.8150
[6,] 48650.4142 -2128.9579
[7,] -13636.1166 48650.4142
[8,] 10687.7774 -13636.1166
[9,] -16898.4159 10687.7774
[10,] -9555.8765 -16898.4159
[11,] -702.2978 -9555.8765
[12,] -9679.1943 -702.2978
[13,] -6352.8345 -9679.1943
[14,] -6808.0907 -6352.8345
[15,] 1903.1345 -6808.0907
[16,] 2892.1168 1903.1345
[17,] 20223.9669 2892.1168
[18,] -2639.8582 20223.9669
[19,] 10796.3416 -2639.8582
[20,] -2312.8652 10796.3416
[21,] -14300.8541 -2312.8652
[22,] 2364.1832 -14300.8541
[23,] -2460.8157 2364.1832
[24,] -103855.4234 -2460.8157
[25,] 25983.1382 -103855.4234
[26,] 68753.9708 25983.1382
[27,] -15148.1371 68753.9708
[28,] 61253.7002 -15148.1371
[29,] -18105.3690 61253.7002
[30,] 22242.4642 -18105.3690
[31,] -59059.8719 22242.4642
[32,] -18224.6106 -59059.8719
[33,] -25081.8417 -18224.6106
[34,] -33282.7875 -25081.8417
[35,] -2886.1817 -33282.7875
[36,] -7289.9013 -2886.1817
[37,] 2786.0526 -7289.9013
[38,] 61240.5056 2786.0526
[39,] -6235.2594 61240.5056
[40,] 17548.3892 -6235.2594
[41,] -26166.6653 17548.3892
[42,] 12672.5805 -26166.6653
[43,] 8567.1686 12672.5805
[44,] -145645.5446 8567.1686
[45,] 2256.5480 -145645.5446
[46,] 8627.2496 2256.5480
[47,] -21471.9062 8627.2496
[48,] 10167.3725 -21471.9062
[49,] -10096.2483 10167.3725
[50,] 15235.0193 -10096.2483
[51,] -42715.8559 15235.0193
[52,] 27150.2432 -42715.8559
[53,] 153079.7188 27150.2432
[54,] 56521.1877 153079.7188
[55,] 31132.7744 56521.1877
[56,] 30192.0237 31132.7744
[57,] 24380.7181 30192.0237
[58,] -28031.2173 24380.7181
[59,] -36779.8506 -28031.2173
[60,] -16690.9177 -36779.8506
[61,] -10850.6935 -16690.9177
[62,] 1276.1932 -10850.6935
[63,] -1437.2367 1276.1932
[64,] -14260.9528 -1437.2367
[65,] 22094.5891 -14260.9528
[66,] -19653.4238 22094.5891
[67,] -146.7942 -19653.4238
[68,] -20089.1833 -146.7942
[69,] -9872.8997 -20089.1833
[70,] -3571.6530 -9872.8997
[71,] -3995.0035 -3571.6530
[72,] -6707.1345 -3995.0035
[73,] -5091.8876 -6707.1345
[74,] -7293.3526 -5091.8876
[75,] 17153.6177 -7293.3526
[76,] -16782.2326 17153.6177
[77,] -33707.5443 -16782.2326
[78,] -12343.9437 -33707.5443
[79,] -10347.2033 -12343.9437
[80,] 4767.5248 -10347.2033
[81,] -26331.4982 4767.5248
[82,] -2650.8228 -26331.4982
[83,] 79196.3640 -2650.8228
[84,] -8035.4208 79196.3640
[85,] 43478.0319 -8035.4208
[86,] 10364.9721 43478.0319
[87,] 38365.6553 10364.9721
[88,] 13803.5330 38365.6553
[89,] 6647.7187 13803.5330
[90,] 30215.8906 6647.7187
[91,] -27328.2547 30215.8906
[92,] -13122.1008 -27328.2547
[93,] -22811.4801 -13122.1008
[94,] 27355.4477 -22811.4801
[95,] -23305.3541 27355.4477
[96,] -3514.0341 -23305.3541
[97,] -34075.1426 -3514.0341
[98,] 20696.4727 -34075.1426
[99,] -26982.8310 20696.4727
[100,] -6698.8771 -26982.8310
[101,] -22499.1919 -6698.8771
[102,] -16044.9402 -22499.1919
[103,] -36563.4909 -16044.9402
[104,] 3270.6423 -36563.4909
[105,] 16589.0480 3270.6423
[106,] 61267.4197 16589.0480
[107,] -3194.4860 61267.4197
[108,] 456.7992 -3194.4860
[109,] -26429.4573 456.7992
[110,] -56504.2413 -26429.4573
[111,] 11880.1891 -56504.2413
[112,] -4799.0026 11880.1891
[113,] 19514.6426 -4799.0026
[114,] 179.3858 19514.6426
[115,] 310.1202 179.3858
[116,] 4563.4845 310.1202
[117,] 35727.5684 4563.4845
[118,] -27968.5402 35727.5684
[119,] -9447.9087 -27968.5402
[120,] 44690.9095 -9447.9087
[121,] 32679.6206 44690.9095
[122,] -6940.9532 32679.6206
[123,] -21134.7774 -6940.9532
[124,] -21203.0905 -21134.7774
[125,] -39922.4663 -21203.0905
[126,] 3444.6048 -39922.4663
[127,] 18898.1939 3444.6048
[128,] 2247.9926 18898.1939
[129,] 11226.3705 2247.9926
[130,] -18545.5010 11226.3705
[131,] -13386.6968 -18545.5010
[132,] -26575.1964 -13386.6968
[133,] -32564.8722 -26575.1964
[134,] -32750.6187 -32564.8722
[135,] 12801.1116 -32750.6187
[136,] 8994.6943 12801.1116
[137,] 5213.5884 8994.6943
[138,] 13540.4420 5213.5884
[139,] 377.7720 13540.4420
[140,] 135896.7115 377.7720
[141,] -9140.6020 135896.7115
[142,] 71982.7638 -9140.6020
[143,] 129591.5312 71982.7638
[144,] -36246.3361 129591.5312
[145,] 6078.4838 -36246.3361
[146,] 19378.4327 6078.4838
[147,] -10156.9736 19378.4327
[148,] -40725.8804 -10156.9736
[149,] -41683.9317 -40725.8804
[150,] 440.9151 -41683.9317
[151,] 18398.2658 440.9151
[152,] 24264.1481 18398.2658
[153,] 9539.9235 24264.1481
[154,] -6445.0481 9539.9235
[155,] -10322.6732 -6445.0481
[156,] 92083.5781 -10322.6732
[157,] -53605.2041 92083.5781
[158,] 30039.0146 -53605.2041
[159,] 11782.0894 30039.0146
[160,] 1845.1877 11782.0894
[161,] -8226.5762 1845.1877
[162,] 1805.5966 -8226.5762
[163,] -32754.1827 1805.5966
[164,] 33104.7581 -32754.1827
[165,] 5020.7409 33104.7581
[166,] -58955.7357 5020.7409
[167,] 23790.6840 -58955.7357
[168,] 12831.3486 23790.6840
[169,] -10605.3856 12831.3486
[170,] 2720.6640 -10605.3856
[171,] 6735.5023 2720.6640
[172,] 31018.5467 6735.5023
[173,] -4792.4641 31018.5467
[174,] -54797.0392 -4792.4641
[175,] -15699.6796 -54797.0392
[176,] -5585.2214 -15699.6796
[177,] -39572.8579 -5585.2214
[178,] -16958.1269 -39572.8579
[179,] -57466.7917 -16958.1269
[180,] -7015.0344 -57466.7917
[181,] 28989.0636 -7015.0344
[182,] -84124.0614 28989.0636
[183,] 12702.6719 -84124.0614
[184,] -19063.2697 12702.6719
[185,] -49365.1105 -19063.2697
[186,] -18544.3531 -49365.1105
[187,] 5425.2565 -18544.3531
[188,] -12451.3836 5425.2565
[189,] -14137.9616 -12451.3836
[190,] -53176.4233 -14137.9616
[191,] 1081.5176 -53176.4233
[192,] 40090.7323 1081.5176
[193,] -5297.9592 40090.7323
[194,] 101087.5393 -5297.9592
[195,] -12582.5485 101087.5393
[196,] -2553.3206 -12582.5485
[197,] -8792.0328 -2553.3206
[198,] -19073.3023 -8792.0328
[199,] -133002.0205 -19073.3023
[200,] -21799.1355 -133002.0205
[201,] -22468.5157 -21799.1355
[202,] -5486.8261 -22468.5157
[203,] 664.0490 -5486.8261
[204,] 7480.6969 664.0490
[205,] 31862.6186 7480.6969
[206,] 19255.1157 31862.6186
[207,] 26104.4063 19255.1157
[208,] -631.9995 26104.4063
[209,] -45564.4296 -631.9995
[210,] 16699.5018 -45564.4296
[211,] -35053.1265 16699.5018
[212,] 74080.9117 -35053.1265
[213,] 41502.6981 74080.9117
[214,] 3467.0784 41502.6981
[215,] -55992.7728 3467.0784
[216,] 813.8630 -55992.7728
[217,] 10835.9131 813.8630
[218,] 80832.6136 10835.9131
[219,] 11534.6135 80832.6136
[220,] 16501.7526 11534.6135
[221,] 19508.0542 16501.7526
[222,] 12386.9522 19508.0542
[223,] -16354.7068 12386.9522
[224,] -78860.0683 -16354.7068
[225,] -28698.3091 -78860.0683
[226,] 65244.1790 -28698.3091
[227,] 644.4465 65244.1790
[228,] 21213.4329 644.4465
[229,] 37975.5976 21213.4329
[230,] 11855.1450 37975.5976
[231,] 19805.1940 11855.1450
[232,] 69558.3294 19805.1940
[233,] -12980.1623 69558.3294
[234,] -25998.7983 -12980.1623
[235,] -2724.6155 -25998.7983
[236,] -271.7253 -2724.6155
[237,] -32756.7095 -271.7253
[238,] -68605.3980 -32756.7095
[239,] -16205.3857 -68605.3980
[240,] 68817.7955 -16205.3857
[241,] -29924.8068 68817.7955
[242,] 56401.2851 -29924.8068
[243,] -2640.3245 56401.2851
[244,] -60378.4673 -2640.3245
[245,] 1399.0658 -60378.4673
[246,] -14573.3137 1399.0658
[247,] 45002.9792 -14573.3137
[248,] 17139.0493 45002.9792
[249,] 7133.0257 17139.0493
[250,] 57170.0093 7133.0257
[251,] -11160.4286 57170.0093
[252,] -18814.9604 -11160.4286
[253,] -6313.2639 -18814.9604
[254,] 28075.4344 -6313.2639
[255,] 44823.8293 28075.4344
[256,] 4764.6329 44823.8293
[257,] 57888.9659 4764.6329
[258,] 65717.1710 57888.9659
[259,] -9216.5379 65717.1710
[260,] -40122.8092 -9216.5379
[261,] -8940.9487 -40122.8092
[262,] -1661.4986 -8940.9487
[263,] -48361.0254 -1661.4986
[264,] -15635.1489 -48361.0254
[265,] -27220.1859 -15635.1489
[266,] -21386.0629 -27220.1859
[267,] 44125.0537 -21386.0629
[268,] -30200.0084 44125.0537
[269,] 33414.1774 -30200.0084
[270,] -33637.5549 33414.1774
[271,] 39286.6916 -33637.5549
[272,] 2387.1506 39286.6916
[273,] 32868.1253 2387.1506
[274,] -47519.0458 32868.1253
[275,] 3717.8309 -47519.0458
[276,] -42712.6641 3717.8309
[277,] 69994.6309 -42712.6641
[278,] 11014.9295 69994.6309
[279,] -52463.7184 11014.9295
[280,] -38389.0381 -52463.7184
[281,] -2902.3382 -38389.0381
[282,] -25904.5990 -2902.3382
[283,] -5399.9761 -25904.5990
[284,] 13225.4431 -5399.9761
[285,] -24549.5861 13225.4431
[286,] 95615.7907 -24549.5861
[287,] -87376.6075 95615.7907
[288,] -6356.2017 -87376.6075
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 13677.8715 2355.0288
2 -3243.7456 13677.8715
3 5564.8420 -3243.7456
4 -7629.8150 5564.8420
5 -2128.9579 -7629.8150
6 48650.4142 -2128.9579
7 -13636.1166 48650.4142
8 10687.7774 -13636.1166
9 -16898.4159 10687.7774
10 -9555.8765 -16898.4159
11 -702.2978 -9555.8765
12 -9679.1943 -702.2978
13 -6352.8345 -9679.1943
14 -6808.0907 -6352.8345
15 1903.1345 -6808.0907
16 2892.1168 1903.1345
17 20223.9669 2892.1168
18 -2639.8582 20223.9669
19 10796.3416 -2639.8582
20 -2312.8652 10796.3416
21 -14300.8541 -2312.8652
22 2364.1832 -14300.8541
23 -2460.8157 2364.1832
24 -103855.4234 -2460.8157
25 25983.1382 -103855.4234
26 68753.9708 25983.1382
27 -15148.1371 68753.9708
28 61253.7002 -15148.1371
29 -18105.3690 61253.7002
30 22242.4642 -18105.3690
31 -59059.8719 22242.4642
32 -18224.6106 -59059.8719
33 -25081.8417 -18224.6106
34 -33282.7875 -25081.8417
35 -2886.1817 -33282.7875
36 -7289.9013 -2886.1817
37 2786.0526 -7289.9013
38 61240.5056 2786.0526
39 -6235.2594 61240.5056
40 17548.3892 -6235.2594
41 -26166.6653 17548.3892
42 12672.5805 -26166.6653
43 8567.1686 12672.5805
44 -145645.5446 8567.1686
45 2256.5480 -145645.5446
46 8627.2496 2256.5480
47 -21471.9062 8627.2496
48 10167.3725 -21471.9062
49 -10096.2483 10167.3725
50 15235.0193 -10096.2483
51 -42715.8559 15235.0193
52 27150.2432 -42715.8559
53 153079.7188 27150.2432
54 56521.1877 153079.7188
55 31132.7744 56521.1877
56 30192.0237 31132.7744
57 24380.7181 30192.0237
58 -28031.2173 24380.7181
59 -36779.8506 -28031.2173
60 -16690.9177 -36779.8506
61 -10850.6935 -16690.9177
62 1276.1932 -10850.6935
63 -1437.2367 1276.1932
64 -14260.9528 -1437.2367
65 22094.5891 -14260.9528
66 -19653.4238 22094.5891
67 -146.7942 -19653.4238
68 -20089.1833 -146.7942
69 -9872.8997 -20089.1833
70 -3571.6530 -9872.8997
71 -3995.0035 -3571.6530
72 -6707.1345 -3995.0035
73 -5091.8876 -6707.1345
74 -7293.3526 -5091.8876
75 17153.6177 -7293.3526
76 -16782.2326 17153.6177
77 -33707.5443 -16782.2326
78 -12343.9437 -33707.5443
79 -10347.2033 -12343.9437
80 4767.5248 -10347.2033
81 -26331.4982 4767.5248
82 -2650.8228 -26331.4982
83 79196.3640 -2650.8228
84 -8035.4208 79196.3640
85 43478.0319 -8035.4208
86 10364.9721 43478.0319
87 38365.6553 10364.9721
88 13803.5330 38365.6553
89 6647.7187 13803.5330
90 30215.8906 6647.7187
91 -27328.2547 30215.8906
92 -13122.1008 -27328.2547
93 -22811.4801 -13122.1008
94 27355.4477 -22811.4801
95 -23305.3541 27355.4477
96 -3514.0341 -23305.3541
97 -34075.1426 -3514.0341
98 20696.4727 -34075.1426
99 -26982.8310 20696.4727
100 -6698.8771 -26982.8310
101 -22499.1919 -6698.8771
102 -16044.9402 -22499.1919
103 -36563.4909 -16044.9402
104 3270.6423 -36563.4909
105 16589.0480 3270.6423
106 61267.4197 16589.0480
107 -3194.4860 61267.4197
108 456.7992 -3194.4860
109 -26429.4573 456.7992
110 -56504.2413 -26429.4573
111 11880.1891 -56504.2413
112 -4799.0026 11880.1891
113 19514.6426 -4799.0026
114 179.3858 19514.6426
115 310.1202 179.3858
116 4563.4845 310.1202
117 35727.5684 4563.4845
118 -27968.5402 35727.5684
119 -9447.9087 -27968.5402
120 44690.9095 -9447.9087
121 32679.6206 44690.9095
122 -6940.9532 32679.6206
123 -21134.7774 -6940.9532
124 -21203.0905 -21134.7774
125 -39922.4663 -21203.0905
126 3444.6048 -39922.4663
127 18898.1939 3444.6048
128 2247.9926 18898.1939
129 11226.3705 2247.9926
130 -18545.5010 11226.3705
131 -13386.6968 -18545.5010
132 -26575.1964 -13386.6968
133 -32564.8722 -26575.1964
134 -32750.6187 -32564.8722
135 12801.1116 -32750.6187
136 8994.6943 12801.1116
137 5213.5884 8994.6943
138 13540.4420 5213.5884
139 377.7720 13540.4420
140 135896.7115 377.7720
141 -9140.6020 135896.7115
142 71982.7638 -9140.6020
143 129591.5312 71982.7638
144 -36246.3361 129591.5312
145 6078.4838 -36246.3361
146 19378.4327 6078.4838
147 -10156.9736 19378.4327
148 -40725.8804 -10156.9736
149 -41683.9317 -40725.8804
150 440.9151 -41683.9317
151 18398.2658 440.9151
152 24264.1481 18398.2658
153 9539.9235 24264.1481
154 -6445.0481 9539.9235
155 -10322.6732 -6445.0481
156 92083.5781 -10322.6732
157 -53605.2041 92083.5781
158 30039.0146 -53605.2041
159 11782.0894 30039.0146
160 1845.1877 11782.0894
161 -8226.5762 1845.1877
162 1805.5966 -8226.5762
163 -32754.1827 1805.5966
164 33104.7581 -32754.1827
165 5020.7409 33104.7581
166 -58955.7357 5020.7409
167 23790.6840 -58955.7357
168 12831.3486 23790.6840
169 -10605.3856 12831.3486
170 2720.6640 -10605.3856
171 6735.5023 2720.6640
172 31018.5467 6735.5023
173 -4792.4641 31018.5467
174 -54797.0392 -4792.4641
175 -15699.6796 -54797.0392
176 -5585.2214 -15699.6796
177 -39572.8579 -5585.2214
178 -16958.1269 -39572.8579
179 -57466.7917 -16958.1269
180 -7015.0344 -57466.7917
181 28989.0636 -7015.0344
182 -84124.0614 28989.0636
183 12702.6719 -84124.0614
184 -19063.2697 12702.6719
185 -49365.1105 -19063.2697
186 -18544.3531 -49365.1105
187 5425.2565 -18544.3531
188 -12451.3836 5425.2565
189 -14137.9616 -12451.3836
190 -53176.4233 -14137.9616
191 1081.5176 -53176.4233
192 40090.7323 1081.5176
193 -5297.9592 40090.7323
194 101087.5393 -5297.9592
195 -12582.5485 101087.5393
196 -2553.3206 -12582.5485
197 -8792.0328 -2553.3206
198 -19073.3023 -8792.0328
199 -133002.0205 -19073.3023
200 -21799.1355 -133002.0205
201 -22468.5157 -21799.1355
202 -5486.8261 -22468.5157
203 664.0490 -5486.8261
204 7480.6969 664.0490
205 31862.6186 7480.6969
206 19255.1157 31862.6186
207 26104.4063 19255.1157
208 -631.9995 26104.4063
209 -45564.4296 -631.9995
210 16699.5018 -45564.4296
211 -35053.1265 16699.5018
212 74080.9117 -35053.1265
213 41502.6981 74080.9117
214 3467.0784 41502.6981
215 -55992.7728 3467.0784
216 813.8630 -55992.7728
217 10835.9131 813.8630
218 80832.6136 10835.9131
219 11534.6135 80832.6136
220 16501.7526 11534.6135
221 19508.0542 16501.7526
222 12386.9522 19508.0542
223 -16354.7068 12386.9522
224 -78860.0683 -16354.7068
225 -28698.3091 -78860.0683
226 65244.1790 -28698.3091
227 644.4465 65244.1790
228 21213.4329 644.4465
229 37975.5976 21213.4329
230 11855.1450 37975.5976
231 19805.1940 11855.1450
232 69558.3294 19805.1940
233 -12980.1623 69558.3294
234 -25998.7983 -12980.1623
235 -2724.6155 -25998.7983
236 -271.7253 -2724.6155
237 -32756.7095 -271.7253
238 -68605.3980 -32756.7095
239 -16205.3857 -68605.3980
240 68817.7955 -16205.3857
241 -29924.8068 68817.7955
242 56401.2851 -29924.8068
243 -2640.3245 56401.2851
244 -60378.4673 -2640.3245
245 1399.0658 -60378.4673
246 -14573.3137 1399.0658
247 45002.9792 -14573.3137
248 17139.0493 45002.9792
249 7133.0257 17139.0493
250 57170.0093 7133.0257
251 -11160.4286 57170.0093
252 -18814.9604 -11160.4286
253 -6313.2639 -18814.9604
254 28075.4344 -6313.2639
255 44823.8293 28075.4344
256 4764.6329 44823.8293
257 57888.9659 4764.6329
258 65717.1710 57888.9659
259 -9216.5379 65717.1710
260 -40122.8092 -9216.5379
261 -8940.9487 -40122.8092
262 -1661.4986 -8940.9487
263 -48361.0254 -1661.4986
264 -15635.1489 -48361.0254
265 -27220.1859 -15635.1489
266 -21386.0629 -27220.1859
267 44125.0537 -21386.0629
268 -30200.0084 44125.0537
269 33414.1774 -30200.0084
270 -33637.5549 33414.1774
271 39286.6916 -33637.5549
272 2387.1506 39286.6916
273 32868.1253 2387.1506
274 -47519.0458 32868.1253
275 3717.8309 -47519.0458
276 -42712.6641 3717.8309
277 69994.6309 -42712.6641
278 11014.9295 69994.6309
279 -52463.7184 11014.9295
280 -38389.0381 -52463.7184
281 -2902.3382 -38389.0381
282 -25904.5990 -2902.3382
283 -5399.9761 -25904.5990
284 13225.4431 -5399.9761
285 -24549.5861 13225.4431
286 95615.7907 -24549.5861
287 -87376.6075 95615.7907
288 -6356.2017 -87376.6075
> 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/7uabt1323954524.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/8gwuy1323954524.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/9bgm81323954524.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/10ernv1323954524.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/11x04l1323954524.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/12kwx61323954524.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/13v1mn1323954524.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/14qgol1323954524.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/15r0f81323954524.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/16w9681323954524.tab")
+ }
>
> try(system("convert tmp/1hfmo1323954524.ps tmp/1hfmo1323954524.png",intern=TRUE))
character(0)
> try(system("convert tmp/2xx4e1323954524.ps tmp/2xx4e1323954524.png",intern=TRUE))
character(0)
> try(system("convert tmp/3h5d91323954524.ps tmp/3h5d91323954524.png",intern=TRUE))
character(0)
> try(system("convert tmp/4okrc1323954524.ps tmp/4okrc1323954524.png",intern=TRUE))
character(0)
> try(system("convert tmp/5lz3f1323954524.ps tmp/5lz3f1323954524.png",intern=TRUE))
character(0)
> try(system("convert tmp/6spt01323954524.ps tmp/6spt01323954524.png",intern=TRUE))
character(0)
> try(system("convert tmp/7uabt1323954524.ps tmp/7uabt1323954524.png",intern=TRUE))
character(0)
> try(system("convert tmp/8gwuy1323954524.ps tmp/8gwuy1323954524.png",intern=TRUE))
character(0)
> try(system("convert tmp/9bgm81323954524.ps tmp/9bgm81323954524.png",intern=TRUE))
character(0)
> try(system("convert tmp/10ernv1323954524.ps tmp/10ernv1323954524.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
8.347 0.766 10.441