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(1
+ ,1
+ ,1
+ ,1418
+ ,1418
+ ,56
+ ,56
+ ,30
+ ,30
+ ,145
+ ,145
+ ,3
+ ,3
+ ,210907
+ ,1
+ ,2
+ ,2
+ ,869
+ ,869
+ ,56
+ ,56
+ ,28
+ ,28
+ ,101
+ ,101
+ ,4
+ ,4
+ ,120982
+ ,1
+ ,3
+ ,3
+ ,1530
+ ,1530
+ ,54
+ ,54
+ ,38
+ ,38
+ ,98
+ ,98
+ ,12
+ ,12
+ ,176508
+ ,1
+ ,4
+ ,4
+ ,2172
+ ,2172
+ ,89
+ ,89
+ ,30
+ ,30
+ ,132
+ ,132
+ ,2
+ ,2
+ ,179321
+ ,1
+ ,5
+ ,5
+ ,901
+ ,901
+ ,40
+ ,40
+ ,22
+ ,22
+ ,60
+ ,60
+ ,1
+ ,1
+ ,123185
+ ,1
+ ,6
+ ,6
+ ,463
+ ,463
+ ,25
+ ,25
+ ,26
+ ,26
+ ,38
+ ,38
+ ,3
+ ,3
+ ,52746
+ ,1
+ ,7
+ ,7
+ ,3201
+ ,3201
+ ,92
+ ,92
+ ,25
+ ,25
+ ,144
+ ,144
+ ,0
+ ,0
+ ,385534
+ ,1
+ ,8
+ ,8
+ ,371
+ ,371
+ ,18
+ ,18
+ ,18
+ ,18
+ ,5
+ ,5
+ ,0
+ ,0
+ ,33170
+ ,1
+ ,9
+ ,9
+ ,1583
+ ,1583
+ ,44
+ ,44
+ ,26
+ ,26
+ ,84
+ ,84
+ ,5
+ ,5
+ ,149061
+ ,1
+ ,10
+ ,10
+ ,1439
+ ,1439
+ ,33
+ ,33
+ ,25
+ ,25
+ ,79
+ ,79
+ ,0
+ ,0
+ ,165446
+ ,1
+ ,11
+ ,11
+ ,1764
+ ,1764
+ ,84
+ ,84
+ ,38
+ ,38
+ ,127
+ ,127
+ ,0
+ ,0
+ ,237213
+ ,1
+ ,12
+ ,12
+ ,1495
+ ,1495
+ ,88
+ ,88
+ ,44
+ ,44
+ ,78
+ ,78
+ ,7
+ ,7
+ ,173326
+ ,1
+ ,13
+ ,13
+ ,1373
+ ,1373
+ ,55
+ ,55
+ ,30
+ ,30
+ ,60
+ ,60
+ ,7
+ ,7
+ ,133131
+ ,1
+ ,14
+ ,14
+ ,2187
+ ,2187
+ ,60
+ ,60
+ ,40
+ ,40
+ ,131
+ ,131
+ ,3
+ ,3
+ ,258873
+ ,1
+ ,15
+ ,15
+ ,1491
+ ,1491
+ ,66
+ ,66
+ ,34
+ ,34
+ ,84
+ ,84
+ ,9
+ ,9
+ ,180083
+ ,1
+ ,16
+ ,16
+ ,4041
+ ,4041
+ ,154
+ ,154
+ ,47
+ ,47
+ ,133
+ ,133
+ ,0
+ ,0
+ ,324799
+ ,1
+ ,17
+ ,17
+ ,1706
+ ,1706
+ ,53
+ ,53
+ ,30
+ ,30
+ ,150
+ ,150
+ ,4
+ ,4
+ ,230964
+ ,1
+ ,18
+ ,18
+ ,2152
+ ,2152
+ ,119
+ ,119
+ ,31
+ ,31
+ ,91
+ ,91
+ ,3
+ ,3
+ ,236785
+ ,1
+ ,19
+ ,19
+ ,1036
+ ,1036
+ ,41
+ ,41
+ ,23
+ ,23
+ ,132
+ ,132
+ ,0
+ ,0
+ ,135473
+ ,1
+ ,20
+ ,20
+ ,1882
+ ,1882
+ ,61
+ ,61
+ ,36
+ ,36
+ ,136
+ ,136
+ ,7
+ ,7
+ ,202925
+ ,1
+ ,21
+ ,21
+ ,1929
+ ,1929
+ ,58
+ ,58
+ ,36
+ ,36
+ ,124
+ ,124
+ ,0
+ ,0
+ ,215147
+ ,1
+ ,22
+ ,22
+ ,2242
+ ,2242
+ ,75
+ ,75
+ ,30
+ ,30
+ ,118
+ ,118
+ ,1
+ ,1
+ ,344297
+ ,1
+ ,23
+ ,23
+ ,1220
+ ,1220
+ ,33
+ ,33
+ ,25
+ ,25
+ ,70
+ ,70
+ ,5
+ ,5
+ ,153935
+ ,1
+ ,24
+ ,24
+ ,1289
+ ,1289
+ ,40
+ ,40
+ ,39
+ ,39
+ ,107
+ ,107
+ ,7
+ ,7
+ ,132943
+ ,1
+ ,25
+ ,25
+ ,2515
+ ,2515
+ ,92
+ ,92
+ ,34
+ ,34
+ ,119
+ ,119
+ ,0
+ ,0
+ ,174724
+ ,1
+ ,26
+ ,26
+ ,2147
+ ,2147
+ ,100
+ ,100
+ ,31
+ ,31
+ ,89
+ ,89
+ ,0
+ ,0
+ ,174415
+ ,1
+ ,27
+ ,27
+ ,2352
+ ,2352
+ ,112
+ ,112
+ ,31
+ ,31
+ ,112
+ ,112
+ ,5
+ ,5
+ ,225548
+ ,1
+ ,28
+ ,28
+ ,1638
+ ,1638
+ ,73
+ ,73
+ ,33
+ ,33
+ ,108
+ ,108
+ ,0
+ ,0
+ ,223632
+ ,1
+ ,29
+ ,29
+ ,1222
+ ,1222
+ ,40
+ ,40
+ ,25
+ ,25
+ ,52
+ ,52
+ ,0
+ ,0
+ ,124817
+ ,1
+ ,30
+ ,30
+ ,1812
+ ,1812
+ ,45
+ ,45
+ ,33
+ ,33
+ ,112
+ ,112
+ ,0
+ ,0
+ ,221698
+ ,1
+ ,31
+ ,31
+ ,1677
+ ,1677
+ ,60
+ ,60
+ ,35
+ ,35
+ ,116
+ ,116
+ ,3
+ ,3
+ ,210767
+ ,1
+ ,32
+ ,32
+ ,1579
+ ,1579
+ ,62
+ ,62
+ ,42
+ ,42
+ ,123
+ ,123
+ ,4
+ ,4
+ ,170266
+ ,1
+ ,33
+ ,33
+ ,1731
+ ,1731
+ ,75
+ ,75
+ ,43
+ ,43
+ ,125
+ ,125
+ ,1
+ ,1
+ ,260561
+ ,1
+ ,34
+ ,34
+ ,807
+ ,807
+ ,31
+ ,31
+ ,30
+ ,30
+ ,27
+ ,27
+ ,4
+ ,4
+ ,84853
+ ,1
+ ,35
+ ,35
+ ,2452
+ ,2452
+ ,77
+ ,77
+ ,33
+ ,33
+ ,162
+ ,162
+ ,2
+ ,2
+ ,294424
+ ,1
+ ,36
+ ,36
+ ,1940
+ ,1940
+ ,46
+ ,46
+ ,32
+ ,32
+ ,64
+ ,64
+ ,0
+ ,0
+ ,215641
+ ,1
+ ,37
+ ,37
+ ,2662
+ ,2662
+ ,99
+ ,99
+ ,36
+ ,36
+ ,92
+ ,92
+ ,0
+ ,0
+ ,325107
+ ,1
+ ,38
+ ,38
+ ,1499
+ ,1499
+ ,66
+ ,66
+ ,28
+ ,28
+ ,83
+ ,83
+ ,2
+ ,2
+ ,167542
+ ,1
+ ,39
+ ,39
+ ,865
+ ,865
+ ,30
+ ,30
+ ,14
+ ,14
+ ,41
+ ,41
+ ,1
+ ,1
+ ,106408
+ ,1
+ ,40
+ ,40
+ ,2527
+ ,2527
+ ,146
+ ,146
+ ,32
+ ,32
+ ,120
+ ,120
+ ,2
+ ,2
+ ,265769
+ ,1
+ ,41
+ ,41
+ ,2747
+ ,2747
+ ,67
+ ,67
+ ,30
+ ,30
+ ,105
+ ,105
+ ,10
+ ,10
+ ,269651
+ ,1
+ ,42
+ ,42
+ ,1324
+ ,1324
+ ,56
+ ,56
+ ,35
+ ,35
+ ,79
+ ,79
+ ,6
+ ,6
+ ,149112
+ ,1
+ ,43
+ ,43
+ ,1383
+ ,1383
+ ,58
+ ,58
+ ,28
+ ,28
+ ,70
+ ,70
+ ,5
+ ,5
+ ,152871
+ ,1
+ ,44
+ ,44
+ ,1179
+ ,1179
+ ,34
+ ,34
+ ,28
+ ,28
+ ,55
+ ,55
+ ,4
+ ,4
+ ,111665
+ ,1
+ ,45
+ ,45
+ ,2099
+ ,2099
+ ,61
+ ,61
+ ,39
+ ,39
+ ,39
+ ,39
+ ,1
+ ,1
+ ,116408
+ ,1
+ ,46
+ ,46
+ ,4308
+ ,4308
+ ,119
+ ,119
+ ,34
+ ,34
+ ,67
+ ,67
+ ,2
+ ,2
+ ,362301
+ ,1
+ ,47
+ ,47
+ ,918
+ ,918
+ ,42
+ ,42
+ ,26
+ ,26
+ ,21
+ ,21
+ ,2
+ ,2
+ ,78800
+ ,1
+ ,48
+ ,48
+ ,1831
+ ,1831
+ ,66
+ ,66
+ ,39
+ ,39
+ ,127
+ ,127
+ ,0
+ ,0
+ ,183167
+ ,1
+ ,49
+ ,49
+ ,3373
+ ,3373
+ ,89
+ ,89
+ ,39
+ ,39
+ ,152
+ ,152
+ ,8
+ ,8
+ ,277965
+ ,1
+ ,50
+ ,50
+ ,1713
+ ,1713
+ ,44
+ ,44
+ ,33
+ ,33
+ ,113
+ ,113
+ ,3
+ ,3
+ ,150629
+ ,1
+ ,51
+ ,51
+ ,1438
+ ,1438
+ ,66
+ ,66
+ ,28
+ ,28
+ ,99
+ ,99
+ ,0
+ ,0
+ ,168809
+ ,1
+ ,52
+ ,52
+ ,496
+ ,496
+ ,24
+ ,24
+ ,4
+ ,4
+ ,7
+ ,7
+ ,0
+ ,0
+ ,24188
+ ,1
+ ,53
+ ,53
+ ,2253
+ ,2253
+ ,259
+ ,259
+ ,39
+ ,39
+ ,141
+ ,141
+ ,8
+ ,8
+ ,329267
+ ,1
+ ,54
+ ,54
+ ,744
+ ,744
+ ,17
+ ,17
+ ,18
+ ,18
+ ,21
+ ,21
+ ,5
+ ,5
+ ,65029
+ ,1
+ ,55
+ ,55
+ ,1161
+ ,1161
+ ,64
+ ,64
+ ,14
+ ,14
+ ,35
+ ,35
+ ,3
+ ,3
+ ,101097
+ ,1
+ ,56
+ ,56
+ ,2352
+ ,2352
+ ,41
+ ,41
+ ,29
+ ,29
+ ,109
+ ,109
+ ,1
+ ,1
+ ,218946
+ ,1
+ ,57
+ ,57
+ ,2144
+ ,2144
+ ,68
+ ,68
+ ,44
+ ,44
+ ,133
+ ,133
+ ,5
+ ,5
+ ,244052
+ ,1
+ ,58
+ ,58
+ ,2694
+ ,2694
+ ,132
+ ,132
+ ,28
+ ,28
+ ,230
+ ,230
+ ,5
+ ,5
+ ,233328
+ ,1
+ ,59
+ ,59
+ ,1973
+ ,1973
+ ,105
+ ,105
+ ,35
+ ,35
+ ,166
+ ,166
+ ,0
+ ,0
+ ,256462
+ ,1
+ ,60
+ ,60
+ ,1769
+ ,1769
+ ,71
+ ,71
+ ,28
+ ,28
+ ,68
+ ,68
+ ,12
+ ,12
+ ,206161
+ ,1
+ ,61
+ ,61
+ ,3148
+ ,3148
+ ,112
+ ,112
+ ,38
+ ,38
+ ,147
+ ,147
+ ,8
+ ,8
+ ,311473
+ ,1
+ ,62
+ ,62
+ ,2474
+ ,2474
+ ,94
+ ,94
+ ,23
+ ,23
+ ,179
+ ,179
+ ,8
+ ,8
+ ,235800
+ ,1
+ ,63
+ ,63
+ ,2084
+ ,2084
+ ,82
+ ,82
+ ,36
+ ,36
+ ,61
+ ,61
+ ,8
+ ,8
+ ,177939
+ ,1
+ ,64
+ ,64
+ ,1954
+ ,1954
+ ,70
+ ,70
+ ,32
+ ,32
+ ,101
+ ,101
+ ,8
+ ,8
+ ,207176
+ ,1
+ ,65
+ ,65
+ ,1226
+ ,1226
+ ,57
+ ,57
+ ,29
+ ,29
+ ,108
+ ,108
+ ,2
+ ,2
+ ,196553
+ ,1
+ ,66
+ ,66
+ ,1389
+ ,1389
+ ,53
+ ,53
+ ,25
+ ,25
+ ,90
+ ,90
+ ,0
+ ,0
+ ,174184
+ ,1
+ ,67
+ ,67
+ ,1496
+ ,1496
+ ,103
+ ,103
+ ,27
+ ,27
+ ,114
+ ,114
+ ,5
+ ,5
+ ,143246
+ ,1
+ ,68
+ ,68
+ ,2269
+ ,2269
+ ,121
+ ,121
+ ,36
+ ,36
+ ,103
+ ,103
+ ,8
+ ,8
+ ,187559
+ ,1
+ ,69
+ ,69
+ ,1833
+ ,1833
+ ,62
+ ,62
+ ,28
+ ,28
+ ,142
+ ,142
+ ,2
+ ,2
+ ,187681
+ ,1
+ ,70
+ ,70
+ ,1268
+ ,1268
+ ,52
+ ,52
+ ,23
+ ,23
+ ,79
+ ,79
+ ,5
+ ,5
+ ,119016
+ ,1
+ ,71
+ ,71
+ ,1943
+ ,1943
+ ,52
+ ,52
+ ,40
+ ,40
+ ,88
+ ,88
+ ,12
+ ,12
+ ,182192
+ ,1
+ ,72
+ ,72
+ ,893
+ ,893
+ ,32
+ ,32
+ ,23
+ ,23
+ ,25
+ ,25
+ ,6
+ ,6
+ ,73566
+ ,1
+ ,73
+ ,73
+ ,1762
+ ,1762
+ ,62
+ ,62
+ ,40
+ ,40
+ ,83
+ ,83
+ ,7
+ ,7
+ ,194979
+ ,1
+ ,74
+ ,74
+ ,1403
+ ,1403
+ ,45
+ ,45
+ ,28
+ ,28
+ ,113
+ ,113
+ ,2
+ ,2
+ ,167488
+ ,1
+ ,75
+ ,75
+ ,1425
+ ,1425
+ ,46
+ ,46
+ ,34
+ ,34
+ ,118
+ ,118
+ ,0
+ ,0
+ ,143756
+ ,1
+ ,76
+ ,76
+ ,1857
+ ,1857
+ ,63
+ ,63
+ ,33
+ ,33
+ ,110
+ ,110
+ ,4
+ ,4
+ ,275541
+ ,1
+ ,77
+ ,77
+ ,1840
+ ,1840
+ ,75
+ ,75
+ ,28
+ ,28
+ ,129
+ ,129
+ ,3
+ ,3
+ ,243199
+ ,1
+ ,78
+ ,78
+ ,1502
+ ,1502
+ ,88
+ ,88
+ ,34
+ ,34
+ ,51
+ ,51
+ ,6
+ ,6
+ ,182999
+ ,1
+ ,79
+ ,79
+ ,1441
+ ,1441
+ ,46
+ ,46
+ ,30
+ ,30
+ ,93
+ ,93
+ ,2
+ ,2
+ ,135649
+ ,1
+ ,80
+ ,80
+ ,1420
+ ,1420
+ ,53
+ ,53
+ ,33
+ ,33
+ ,76
+ ,76
+ ,0
+ ,0
+ ,152299
+ ,1
+ ,81
+ ,81
+ ,1416
+ ,1416
+ ,37
+ ,37
+ ,22
+ ,22
+ ,49
+ ,49
+ ,1
+ ,1
+ ,120221
+ ,1
+ ,82
+ ,82
+ ,2970
+ ,2970
+ ,90
+ ,90
+ ,38
+ ,38
+ ,118
+ ,118
+ ,0
+ ,0
+ ,346485
+ ,1
+ ,83
+ ,83
+ ,1317
+ ,1317
+ ,63
+ ,63
+ ,26
+ ,26
+ ,38
+ ,38
+ ,5
+ ,5
+ ,145790
+ ,1
+ ,84
+ ,84
+ ,1644
+ ,1644
+ ,78
+ ,78
+ ,35
+ ,35
+ ,141
+ ,141
+ ,2
+ ,2
+ ,193339
+ ,1
+ ,85
+ ,85
+ ,870
+ ,870
+ ,25
+ ,25
+ ,8
+ ,8
+ ,58
+ ,58
+ ,0
+ ,0
+ ,80953
+ ,1
+ ,86
+ ,86
+ ,1654
+ ,1654
+ ,45
+ ,45
+ ,24
+ ,24
+ ,27
+ ,27
+ ,0
+ ,0
+ ,122774
+ ,1
+ ,87
+ ,87
+ ,1054
+ ,1054
+ ,46
+ ,46
+ ,29
+ ,29
+ ,91
+ ,91
+ ,5
+ ,5
+ ,130585
+ ,1
+ ,88
+ ,88
+ ,3004
+ ,3004
+ ,144
+ ,144
+ ,29
+ ,29
+ ,63
+ ,63
+ ,1
+ ,1
+ ,286468
+ ,1
+ ,89
+ ,89
+ ,2008
+ ,2008
+ ,82
+ ,82
+ ,45
+ ,45
+ ,56
+ ,56
+ ,0
+ ,0
+ ,241066
+ ,1
+ ,90
+ ,90
+ ,2547
+ ,2547
+ ,91
+ ,91
+ ,37
+ ,37
+ ,144
+ ,144
+ ,1
+ ,1
+ ,148446
+ ,1
+ ,91
+ ,91
+ ,1885
+ ,1885
+ ,71
+ ,71
+ ,33
+ ,33
+ ,73
+ ,73
+ ,1
+ ,1
+ ,204713
+ ,1
+ ,92
+ ,92
+ ,1626
+ ,1626
+ ,63
+ ,63
+ ,33
+ ,33
+ ,168
+ ,168
+ ,2
+ ,2
+ ,182079
+ ,1
+ ,93
+ ,93
+ ,1468
+ ,1468
+ ,53
+ ,53
+ ,25
+ ,25
+ ,64
+ ,64
+ ,6
+ ,6
+ ,140344
+ ,1
+ ,94
+ ,94
+ ,2445
+ ,2445
+ ,62
+ ,62
+ ,32
+ ,32
+ ,97
+ ,97
+ ,1
+ ,1
+ ,220516
+ ,1
+ ,95
+ ,95
+ ,1964
+ ,1964
+ ,63
+ ,63
+ ,29
+ ,29
+ ,117
+ ,117
+ ,4
+ ,4
+ ,243060
+ ,1
+ ,96
+ ,96
+ ,1381
+ ,1381
+ ,32
+ ,32
+ ,28
+ ,28
+ ,100
+ ,100
+ ,2
+ ,2
+ ,162765
+ ,1
+ ,97
+ ,97
+ ,1369
+ ,1369
+ ,39
+ ,39
+ ,28
+ ,28
+ ,149
+ ,149
+ ,3
+ ,3
+ ,182613
+ ,1
+ ,98
+ ,98
+ ,1659
+ ,1659
+ ,62
+ ,62
+ ,31
+ ,31
+ ,187
+ ,187
+ ,0
+ ,0
+ ,232138
+ ,1
+ ,99
+ ,99
+ ,2888
+ ,2888
+ ,117
+ ,117
+ ,52
+ ,52
+ ,127
+ ,127
+ ,10
+ ,10
+ ,265318
+ ,1
+ ,100
+ ,100
+ ,2845
+ ,2845
+ ,92
+ ,92
+ ,24
+ ,24
+ ,245
+ ,245
+ ,9
+ ,9
+ ,310839
+ ,1
+ ,101
+ ,101
+ ,1982
+ ,1982
+ ,93
+ ,93
+ ,41
+ ,41
+ ,87
+ ,87
+ ,7
+ ,7
+ ,225060
+ ,1
+ ,102
+ ,102
+ ,1904
+ ,1904
+ ,54
+ ,54
+ ,33
+ ,33
+ ,177
+ ,177
+ ,0
+ ,0
+ ,232317
+ ,1
+ ,103
+ ,103
+ ,1391
+ ,1391
+ ,144
+ ,144
+ ,32
+ ,32
+ ,49
+ ,49
+ ,0
+ ,0
+ ,144966
+ ,1
+ ,104
+ ,104
+ ,602
+ ,602
+ ,14
+ ,14
+ ,19
+ ,19
+ ,49
+ ,49
+ ,4
+ ,4
+ ,43287
+ ,1
+ ,105
+ ,105
+ ,1743
+ ,1743
+ ,61
+ ,61
+ ,20
+ ,20
+ ,73
+ ,73
+ ,4
+ ,4
+ ,155754
+ ,1
+ ,106
+ ,106
+ ,1559
+ ,1559
+ ,109
+ ,109
+ ,31
+ ,31
+ ,177
+ ,177
+ ,0
+ ,0
+ ,164709
+ ,1
+ ,107
+ ,107
+ ,2014
+ ,2014
+ ,38
+ ,38
+ ,31
+ ,31
+ ,94
+ ,94
+ ,0
+ ,0
+ ,201940
+ ,1
+ ,108
+ ,108
+ ,2143
+ ,2143
+ ,73
+ ,73
+ ,32
+ ,32
+ ,117
+ ,117
+ ,0
+ ,0
+ ,235454
+ ,1
+ ,109
+ ,109
+ ,874
+ ,874
+ ,50
+ ,50
+ ,23
+ ,23
+ ,55
+ ,55
+ ,0
+ ,0
+ ,99466
+ ,1
+ ,110
+ ,110
+ ,1281
+ ,1281
+ ,72
+ ,72
+ ,30
+ ,30
+ ,58
+ ,58
+ ,0
+ ,0
+ ,100750
+ ,1
+ ,111
+ ,111
+ ,1401
+ ,1401
+ ,50
+ ,50
+ ,31
+ ,31
+ ,95
+ ,95
+ ,4
+ ,4
+ ,224549
+ ,1
+ ,112
+ ,112
+ ,1944
+ ,1944
+ ,71
+ ,71
+ ,42
+ ,42
+ ,129
+ ,129
+ ,0
+ ,0
+ ,243511
+ ,1
+ ,113
+ ,113
+ ,391
+ ,391
+ ,10
+ ,10
+ ,1
+ ,1
+ ,11
+ ,11
+ ,0
+ ,0
+ ,22938
+ ,1
+ ,114
+ ,114
+ ,1605
+ ,1605
+ ,65
+ ,65
+ ,32
+ ,32
+ ,101
+ ,101
+ ,0
+ ,0
+ ,152474
+ ,1
+ ,115
+ ,115
+ ,530
+ ,530
+ ,25
+ ,25
+ ,11
+ ,11
+ ,28
+ ,28
+ ,4
+ ,4
+ ,61857
+ ,1
+ ,116
+ ,116
+ ,1386
+ ,1386
+ ,41
+ ,41
+ ,36
+ ,36
+ ,89
+ ,89
+ ,0
+ ,0
+ ,132487
+ ,1
+ ,117
+ ,117
+ ,2395
+ ,2395
+ ,86
+ ,86
+ ,31
+ ,31
+ ,193
+ ,193
+ ,1
+ ,1
+ ,317394
+ ,1
+ ,118
+ ,118
+ ,387
+ ,387
+ ,16
+ ,16
+ ,0
+ ,0
+ ,4
+ ,4
+ ,0
+ ,0
+ ,21054
+ ,1
+ ,119
+ ,119
+ ,1742
+ ,1742
+ ,42
+ ,42
+ ,24
+ ,24
+ ,84
+ ,84
+ ,5
+ ,5
+ ,209641
+ ,1
+ ,120
+ ,120
+ ,449
+ ,449
+ ,19
+ ,19
+ ,8
+ ,8
+ ,39
+ ,39
+ ,0
+ ,0
+ ,31414
+ ,1
+ ,121
+ ,121
+ ,2699
+ ,2699
+ ,95
+ ,95
+ ,33
+ ,33
+ ,101
+ ,101
+ ,2
+ ,2
+ ,244749
+ ,1
+ ,122
+ ,122
+ ,1606
+ ,1606
+ ,49
+ ,49
+ ,40
+ ,40
+ ,82
+ ,82
+ ,7
+ ,7
+ ,184510
+ ,1
+ ,123
+ ,123
+ ,1204
+ ,1204
+ ,64
+ ,64
+ ,38
+ ,38
+ ,36
+ ,36
+ ,8
+ ,8
+ ,128423
+ ,1
+ ,124
+ ,124
+ ,1138
+ ,1138
+ ,38
+ ,38
+ ,24
+ ,24
+ ,75
+ ,75
+ ,2
+ ,2
+ ,97839
+ ,1
+ ,125
+ ,125
+ ,568
+ ,568
+ ,34
+ ,34
+ ,8
+ ,8
+ ,16
+ ,16
+ ,0
+ ,0
+ ,38214
+ ,1
+ ,126
+ ,126
+ ,1459
+ ,1459
+ ,32
+ ,32
+ ,35
+ ,35
+ ,55
+ ,55
+ ,2
+ ,2
+ ,151101
+ ,1
+ ,127
+ ,127
+ ,2158
+ ,2158
+ ,65
+ ,65
+ ,43
+ ,43
+ ,131
+ ,131
+ ,0
+ ,0
+ ,272458
+ ,1
+ ,128
+ ,128
+ ,1111
+ ,1111
+ ,52
+ ,52
+ ,43
+ ,43
+ ,131
+ ,131
+ ,0
+ ,0
+ ,172494
+ ,1
+ ,129
+ ,129
+ ,2833
+ ,2833
+ ,65
+ ,65
+ ,41
+ ,41
+ ,144
+ ,144
+ ,3
+ ,3
+ ,328107
+ ,1
+ ,130
+ ,130
+ ,1955
+ ,1955
+ ,83
+ ,83
+ ,38
+ ,38
+ ,139
+ ,139
+ ,0
+ ,0
+ ,250579
+ ,1
+ ,131
+ ,131
+ ,2922
+ ,2922
+ ,95
+ ,95
+ ,45
+ ,45
+ ,211
+ ,211
+ ,3
+ ,3
+ ,351067
+ ,1
+ ,132
+ ,132
+ ,1002
+ ,1002
+ ,29
+ ,29
+ ,31
+ ,31
+ ,78
+ ,78
+ ,0
+ ,0
+ ,158015
+ ,1
+ ,133
+ ,133
+ ,956
+ ,956
+ ,33
+ ,33
+ ,28
+ ,28
+ ,39
+ ,39
+ ,0
+ ,0
+ ,85439
+ ,1
+ ,134
+ ,134
+ ,2186
+ ,2186
+ ,247
+ ,247
+ ,31
+ ,31
+ ,90
+ ,90
+ ,4
+ ,4
+ ,229242
+ ,1
+ ,135
+ ,135
+ ,3604
+ ,3604
+ ,139
+ ,139
+ ,40
+ ,40
+ ,166
+ ,166
+ ,4
+ ,4
+ ,351619
+ ,1
+ ,136
+ ,136
+ ,1035
+ ,1035
+ ,29
+ ,29
+ ,30
+ ,30
+ ,12
+ ,12
+ ,11
+ ,11
+ ,84207
+ ,1
+ ,137
+ ,137
+ ,3261
+ ,3261
+ ,110
+ ,110
+ ,37
+ ,37
+ ,133
+ ,133
+ ,0
+ ,0
+ ,324598
+ ,1
+ ,138
+ ,138
+ ,1587
+ ,1587
+ ,67
+ ,67
+ ,30
+ ,30
+ ,69
+ ,69
+ ,4
+ ,4
+ ,131069
+ ,1
+ ,139
+ ,139
+ ,1424
+ ,1424
+ ,42
+ ,42
+ ,35
+ ,35
+ ,119
+ ,119
+ ,0
+ ,0
+ ,204271
+ ,1
+ ,140
+ ,140
+ ,1701
+ ,1701
+ ,65
+ ,65
+ ,32
+ ,32
+ ,119
+ ,119
+ ,1
+ ,1
+ ,165543
+ ,1
+ ,141
+ ,141
+ ,1249
+ ,1249
+ ,94
+ ,94
+ ,27
+ ,27
+ ,65
+ ,65
+ ,0
+ ,0
+ ,141722
+ ,1
+ ,142
+ ,142
+ ,3352
+ ,3352
+ ,95
+ ,95
+ ,31
+ ,31
+ ,101
+ ,101
+ ,9
+ ,9
+ ,299775
+ ,1
+ ,143
+ ,143
+ ,1641
+ ,1641
+ ,67
+ ,67
+ ,31
+ ,31
+ ,196
+ ,196
+ ,1
+ ,1
+ ,195838
+ ,1
+ ,144
+ ,144
+ ,2035
+ ,2035
+ ,63
+ ,63
+ ,21
+ ,21
+ ,15
+ ,15
+ ,3
+ ,3
+ ,173260
+ ,1
+ ,145
+ ,145
+ ,2312
+ ,2312
+ ,83
+ ,83
+ ,39
+ ,39
+ ,136
+ ,136
+ ,10
+ ,10
+ ,254488
+ ,1
+ ,146
+ ,146
+ ,1369
+ ,1369
+ ,45
+ ,45
+ ,41
+ ,41
+ ,89
+ ,89
+ ,5
+ ,5
+ ,104389
+ ,1
+ ,147
+ ,147
+ ,2201
+ ,2201
+ ,70
+ ,70
+ ,32
+ ,32
+ ,123
+ ,123
+ ,2
+ ,2
+ ,199476
+ ,1
+ ,148
+ ,148
+ ,1900
+ ,1900
+ ,83
+ ,83
+ ,39
+ ,39
+ ,163
+ ,163
+ ,1
+ ,1
+ ,224330
+ ,1
+ ,149
+ ,149
+ ,207
+ ,207
+ ,10
+ ,10
+ ,0
+ ,0
+ ,5
+ ,5
+ ,0
+ ,0
+ ,14688
+ ,1
+ ,150
+ ,150
+ ,1645
+ ,1645
+ ,70
+ ,70
+ ,30
+ ,30
+ ,96
+ ,96
+ ,2
+ ,2
+ ,181633
+ ,1
+ ,151
+ ,151
+ ,2429
+ ,2429
+ ,103
+ ,103
+ ,37
+ ,37
+ ,151
+ ,151
+ ,1
+ ,1
+ ,271856
+ ,1
+ ,152
+ ,152
+ ,151
+ ,151
+ ,5
+ ,5
+ ,0
+ ,0
+ ,6
+ ,6
+ ,0
+ ,0
+ ,7199
+ ,1
+ ,153
+ ,153
+ ,474
+ ,474
+ ,20
+ ,20
+ ,5
+ ,5
+ ,13
+ ,13
+ ,0
+ ,0
+ ,46660
+ ,1
+ ,154
+ ,154
+ ,141
+ ,141
+ ,5
+ ,5
+ ,1
+ ,1
+ ,3
+ ,3
+ ,0
+ ,0
+ ,17547
+ ,1
+ ,155
+ ,155
+ ,872
+ ,872
+ ,34
+ ,34
+ ,32
+ ,32
+ ,23
+ ,23
+ ,0
+ ,0
+ ,95227
+ ,1
+ ,156
+ ,156
+ ,1318
+ ,1318
+ ,48
+ ,48
+ ,24
+ ,24
+ ,57
+ ,57
+ ,2
+ ,2
+ ,152601
+ ,0
+ ,157
+ ,0
+ ,1018
+ ,0
+ ,40
+ ,0
+ ,17
+ ,0
+ ,14
+ ,0
+ ,0
+ ,0
+ ,98146
+ ,0
+ ,158
+ ,0
+ ,1383
+ ,0
+ ,43
+ ,0
+ ,11
+ ,0
+ ,43
+ ,0
+ ,3
+ ,0
+ ,79619
+ ,0
+ ,159
+ ,0
+ ,1314
+ ,0
+ ,31
+ ,0
+ ,24
+ ,0
+ ,20
+ ,0
+ ,6
+ ,0
+ ,59194
+ ,0
+ ,160
+ ,0
+ ,1335
+ ,0
+ ,42
+ ,0
+ ,22
+ ,0
+ ,72
+ ,0
+ ,0
+ ,0
+ ,139942
+ ,0
+ ,161
+ ,0
+ ,1403
+ ,0
+ ,46
+ ,0
+ ,12
+ ,0
+ ,87
+ ,0
+ ,2
+ ,0
+ ,118612
+ ,0
+ ,162
+ ,0
+ ,910
+ ,0
+ ,33
+ ,0
+ ,19
+ ,0
+ ,21
+ ,0
+ ,0
+ ,0
+ ,72880
+ ,0
+ ,163
+ ,0
+ ,616
+ ,0
+ ,18
+ ,0
+ ,13
+ ,0
+ ,56
+ ,0
+ ,2
+ ,0
+ ,65475
+ ,0
+ ,164
+ ,0
+ ,1407
+ ,0
+ ,55
+ ,0
+ ,17
+ ,0
+ ,59
+ ,0
+ ,1
+ ,0
+ ,99643
+ ,0
+ ,165
+ ,0
+ ,771
+ ,0
+ ,35
+ ,0
+ ,15
+ ,0
+ ,82
+ ,0
+ ,1
+ ,0
+ ,71965
+ ,0
+ ,166
+ ,0
+ ,766
+ ,0
+ ,59
+ ,0
+ ,16
+ ,0
+ ,43
+ ,0
+ ,2
+ ,0
+ ,77272
+ ,0
+ ,167
+ ,0
+ ,473
+ ,0
+ ,19
+ ,0
+ ,24
+ ,0
+ ,25
+ ,0
+ ,1
+ ,0
+ ,49289
+ ,0
+ ,168
+ ,0
+ ,1376
+ ,0
+ ,66
+ ,0
+ ,15
+ ,0
+ ,38
+ ,0
+ ,0
+ ,0
+ ,135131
+ ,0
+ ,169
+ ,0
+ ,1232
+ ,0
+ ,60
+ ,0
+ ,17
+ ,0
+ ,25
+ ,0
+ ,1
+ ,0
+ ,108446
+ ,0
+ ,170
+ ,0
+ ,1521
+ ,0
+ ,36
+ ,0
+ ,18
+ ,0
+ ,38
+ ,0
+ ,3
+ ,0
+ ,89746
+ ,0
+ ,171
+ ,0
+ ,572
+ ,0
+ ,25
+ ,0
+ ,20
+ ,0
+ ,12
+ ,0
+ ,0
+ ,0
+ ,44296
+ ,0
+ ,172
+ ,0
+ ,1059
+ ,0
+ ,47
+ ,0
+ ,16
+ ,0
+ ,29
+ ,0
+ ,0
+ ,0
+ ,77648
+ ,0
+ ,173
+ ,0
+ ,1544
+ ,0
+ ,54
+ ,0
+ ,16
+ ,0
+ ,47
+ ,0
+ ,0
+ ,0
+ ,181528
+ ,0
+ ,174
+ ,0
+ ,1230
+ ,0
+ ,53
+ ,0
+ ,18
+ ,0
+ ,45
+ ,0
+ ,0
+ ,0
+ ,134019
+ ,0
+ ,175
+ ,0
+ ,1206
+ ,0
+ ,40
+ ,0
+ ,22
+ ,0
+ ,40
+ ,0
+ ,1
+ ,0
+ ,124064
+ ,0
+ ,176
+ ,0
+ ,1205
+ ,0
+ ,40
+ ,0
+ ,8
+ ,0
+ ,30
+ ,0
+ ,4
+ ,0
+ ,92630
+ ,0
+ ,177
+ ,0
+ ,1255
+ ,0
+ ,39
+ ,0
+ ,17
+ ,0
+ ,41
+ ,0
+ ,0
+ ,0
+ ,121848
+ ,0
+ ,178
+ ,0
+ ,613
+ ,0
+ ,14
+ ,0
+ ,18
+ ,0
+ ,25
+ ,0
+ ,0
+ ,0
+ ,52915
+ ,0
+ ,179
+ ,0
+ ,721
+ ,0
+ ,45
+ ,0
+ ,16
+ ,0
+ ,23
+ ,0
+ ,0
+ ,0
+ ,81872
+ ,0
+ ,180
+ ,0
+ ,1109
+ ,0
+ ,36
+ ,0
+ ,23
+ ,0
+ ,14
+ ,0
+ ,7
+ ,0
+ ,58981
+ ,0
+ ,181
+ ,0
+ ,740
+ ,0
+ ,28
+ ,0
+ ,22
+ ,0
+ ,16
+ ,0
+ ,2
+ ,0
+ ,53515
+ ,0
+ ,182
+ ,0
+ ,1126
+ ,0
+ ,44
+ ,0
+ ,13
+ ,0
+ ,26
+ ,0
+ ,0
+ ,0
+ ,60812
+ ,0
+ ,183
+ ,0
+ ,728
+ ,0
+ ,30
+ ,0
+ ,13
+ ,0
+ ,21
+ ,0
+ ,7
+ ,0
+ ,56375
+ ,0
+ ,184
+ ,0
+ ,689
+ ,0
+ ,22
+ ,0
+ ,16
+ ,0
+ ,27
+ ,0
+ ,3
+ ,0
+ ,65490
+ ,0
+ ,185
+ ,0
+ ,592
+ ,0
+ ,17
+ ,0
+ ,16
+ ,0
+ ,9
+ ,0
+ ,0
+ ,0
+ ,80949
+ ,0
+ ,186
+ ,0
+ ,995
+ ,0
+ ,31
+ ,0
+ ,20
+ ,0
+ ,33
+ ,0
+ ,0
+ ,0
+ ,76302
+ ,0
+ ,187
+ ,0
+ ,1613
+ ,0
+ ,55
+ ,0
+ ,22
+ ,0
+ ,42
+ ,0
+ ,6
+ ,0
+ ,104011
+ ,0
+ ,188
+ ,0
+ ,2048
+ ,0
+ ,54
+ ,0
+ ,17
+ ,0
+ ,68
+ ,0
+ ,2
+ ,0
+ ,98104
+ ,0
+ ,189
+ ,0
+ ,705
+ ,0
+ ,21
+ ,0
+ ,18
+ ,0
+ ,32
+ ,0
+ ,0
+ ,0
+ ,67989
+ ,0
+ ,190
+ ,0
+ ,301
+ ,0
+ ,14
+ ,0
+ ,17
+ ,0
+ ,6
+ ,0
+ ,0
+ ,0
+ ,30989
+ ,0
+ ,191
+ ,0
+ ,1803
+ ,0
+ ,81
+ ,0
+ ,12
+ ,0
+ ,67
+ ,0
+ ,3
+ ,0
+ ,135458
+ ,0
+ ,192
+ ,0
+ ,799
+ ,0
+ ,35
+ ,0
+ ,7
+ ,0
+ ,33
+ ,0
+ ,0
+ ,0
+ ,73504
+ ,0
+ ,193
+ ,0
+ ,861
+ ,0
+ ,43
+ ,0
+ ,17
+ ,0
+ ,77
+ ,0
+ ,1
+ ,0
+ ,63123
+ ,0
+ ,194
+ ,0
+ ,1186
+ ,0
+ ,46
+ ,0
+ ,14
+ ,0
+ ,46
+ ,0
+ ,1
+ ,0
+ ,61254
+ ,0
+ ,195
+ ,0
+ ,1451
+ ,0
+ ,30
+ ,0
+ ,23
+ ,0
+ ,30
+ ,0
+ ,0
+ ,0
+ ,74914
+ ,0
+ ,196
+ ,0
+ ,628
+ ,0
+ ,23
+ ,0
+ ,17
+ ,0
+ ,0
+ ,0
+ ,1
+ ,0
+ ,31774
+ ,0
+ ,197
+ ,0
+ ,1161
+ ,0
+ ,38
+ ,0
+ ,14
+ ,0
+ ,36
+ ,0
+ ,0
+ ,0
+ ,81437
+ ,0
+ ,198
+ ,0
+ ,1463
+ ,0
+ ,54
+ ,0
+ ,15
+ ,0
+ ,46
+ ,0
+ ,0
+ ,0
+ ,87186
+ ,0
+ ,199
+ ,0
+ ,742
+ ,0
+ ,20
+ ,0
+ ,17
+ ,0
+ ,18
+ ,0
+ ,0
+ ,0
+ ,50090
+ ,0
+ ,200
+ ,0
+ ,979
+ ,0
+ ,53
+ ,0
+ ,21
+ ,0
+ ,48
+ ,0
+ ,0
+ ,0
+ ,65745
+ ,0
+ ,201
+ ,0
+ ,675
+ ,0
+ ,45
+ ,0
+ ,18
+ ,0
+ ,29
+ ,0
+ ,0
+ ,0
+ ,56653
+ ,0
+ ,202
+ ,0
+ ,1241
+ ,0
+ ,39
+ ,0
+ ,18
+ ,0
+ ,28
+ ,0
+ ,0
+ ,0
+ ,158399
+ ,0
+ ,203
+ ,0
+ ,676
+ ,0
+ ,20
+ ,0
+ ,17
+ ,0
+ ,34
+ ,0
+ ,0
+ ,0
+ ,46455
+ ,0
+ ,204
+ ,0
+ ,1049
+ ,0
+ ,24
+ ,0
+ ,17
+ ,0
+ ,33
+ ,0
+ ,0
+ ,0
+ ,73624
+ ,0
+ ,205
+ ,0
+ ,620
+ ,0
+ ,31
+ ,0
+ ,16
+ ,0
+ ,34
+ ,0
+ ,0
+ ,0
+ ,38395
+ ,0
+ ,206
+ ,0
+ ,1081
+ ,0
+ ,35
+ ,0
+ ,15
+ ,0
+ ,33
+ ,0
+ ,0
+ ,0
+ ,91899
+ ,0
+ ,207
+ ,0
+ ,1688
+ ,0
+ ,151
+ ,0
+ ,21
+ ,0
+ ,80
+ ,0
+ ,0
+ ,0
+ ,139526
+ ,0
+ ,208
+ ,0
+ ,736
+ ,0
+ ,52
+ ,0
+ ,16
+ ,0
+ ,32
+ ,0
+ ,0
+ ,0
+ ,52164
+ ,0
+ ,209
+ ,0
+ ,617
+ ,0
+ ,30
+ ,0
+ ,14
+ ,0
+ ,30
+ ,0
+ ,2
+ ,0
+ ,51567
+ ,0
+ ,210
+ ,0
+ ,812
+ ,0
+ ,31
+ ,0
+ ,15
+ ,0
+ ,41
+ ,0
+ ,0
+ ,0
+ ,70551
+ ,0
+ ,211
+ ,0
+ ,1051
+ ,0
+ ,29
+ ,0
+ ,17
+ ,0
+ ,41
+ ,0
+ ,1
+ ,0
+ ,84856
+ ,0
+ ,212
+ ,0
+ ,1656
+ ,0
+ ,57
+ ,0
+ ,15
+ ,0
+ ,51
+ ,0
+ ,1
+ ,0
+ ,102538
+ ,0
+ ,213
+ ,0
+ ,705
+ ,0
+ ,40
+ ,0
+ ,15
+ ,0
+ ,18
+ ,0
+ ,0
+ ,0
+ ,86678
+ ,0
+ ,214
+ ,0
+ ,945
+ ,0
+ ,44
+ ,0
+ ,10
+ ,0
+ ,34
+ ,0
+ ,0
+ ,0
+ ,85709
+ ,0
+ ,215
+ ,0
+ ,554
+ ,0
+ ,25
+ ,0
+ ,6
+ ,0
+ ,31
+ ,0
+ ,0
+ ,0
+ ,34662
+ ,0
+ ,216
+ ,0
+ ,1597
+ ,0
+ ,77
+ ,0
+ ,22
+ ,0
+ ,39
+ ,0
+ ,0
+ ,0
+ ,150580
+ ,0
+ ,217
+ ,0
+ ,982
+ ,0
+ ,35
+ ,0
+ ,21
+ ,0
+ ,54
+ ,0
+ ,0
+ ,0
+ ,99611
+ ,0
+ ,218
+ ,0
+ ,222
+ ,0
+ ,11
+ ,0
+ ,1
+ ,0
+ ,14
+ ,0
+ ,0
+ ,0
+ ,19349
+ ,0
+ ,219
+ ,0
+ ,1212
+ ,0
+ ,63
+ ,0
+ ,18
+ ,0
+ ,24
+ ,0
+ ,1
+ ,0
+ ,99373
+ ,0
+ ,220
+ ,0
+ ,1143
+ ,0
+ ,44
+ ,0
+ ,17
+ ,0
+ ,24
+ ,0
+ ,0
+ ,0
+ ,86230
+ ,0
+ ,221
+ ,0
+ ,435
+ ,0
+ ,19
+ ,0
+ ,4
+ ,0
+ ,8
+ ,0
+ ,0
+ ,0
+ ,30837
+ ,0
+ ,222
+ ,0
+ ,532
+ ,0
+ ,13
+ ,0
+ ,10
+ ,0
+ ,26
+ ,0
+ ,0
+ ,0
+ ,31706
+ ,0
+ ,223
+ ,0
+ ,882
+ ,0
+ ,42
+ ,0
+ ,16
+ ,0
+ ,19
+ ,0
+ ,0
+ ,0
+ ,89806
+ ,0
+ ,224
+ ,0
+ ,608
+ ,0
+ ,38
+ ,0
+ ,16
+ ,0
+ ,11
+ ,0
+ ,1
+ ,0
+ ,62088
+ ,0
+ ,225
+ ,0
+ ,459
+ ,0
+ ,29
+ ,0
+ ,9
+ ,0
+ ,14
+ ,0
+ ,0
+ ,0
+ ,40151
+ ,0
+ ,226
+ ,0
+ ,578
+ ,0
+ ,20
+ ,0
+ ,16
+ ,0
+ ,1
+ ,0
+ ,0
+ ,0
+ ,27634
+ ,0
+ ,227
+ ,0
+ ,826
+ ,0
+ ,27
+ ,0
+ ,17
+ ,0
+ ,39
+ ,0
+ ,0
+ ,0
+ ,76990
+ ,0
+ ,228
+ ,0
+ ,509
+ ,0
+ ,20
+ ,0
+ ,7
+ ,0
+ ,5
+ ,0
+ ,0
+ ,0
+ ,37460
+ ,0
+ ,229
+ ,0
+ ,717
+ ,0
+ ,19
+ ,0
+ ,15
+ ,0
+ ,37
+ ,0
+ ,0
+ ,0
+ ,54157
+ ,0
+ ,230
+ ,0
+ ,637
+ ,0
+ ,37
+ ,0
+ ,14
+ ,0
+ ,32
+ ,0
+ ,0
+ ,0
+ ,49862
+ ,0
+ ,231
+ ,0
+ ,857
+ ,0
+ ,26
+ ,0
+ ,14
+ ,0
+ ,38
+ ,0
+ ,0
+ ,0
+ ,84337
+ ,0
+ ,232
+ ,0
+ ,830
+ ,0
+ ,42
+ ,0
+ ,18
+ ,0
+ ,47
+ ,0
+ ,0
+ ,0
+ ,64175
+ ,0
+ ,233
+ ,0
+ ,652
+ ,0
+ ,49
+ ,0
+ ,12
+ ,0
+ ,47
+ ,0
+ ,0
+ ,0
+ ,59382
+ ,0
+ ,234
+ ,0
+ ,707
+ ,0
+ ,30
+ ,0
+ ,16
+ ,0
+ ,37
+ ,0
+ ,0
+ ,0
+ ,119308
+ ,0
+ ,235
+ ,0
+ ,954
+ ,0
+ ,49
+ ,0
+ ,21
+ ,0
+ ,51
+ ,0
+ ,0
+ ,0
+ ,76702
+ ,0
+ ,236
+ ,0
+ ,1461
+ ,0
+ ,67
+ ,0
+ ,19
+ ,0
+ ,45
+ ,0
+ ,1
+ ,0
+ ,103425
+ ,0
+ ,237
+ ,0
+ ,672
+ ,0
+ ,28
+ ,0
+ ,16
+ ,0
+ ,21
+ ,0
+ ,0
+ ,0
+ ,70344
+ ,0
+ ,238
+ ,0
+ ,778
+ ,0
+ ,19
+ ,0
+ ,1
+ ,0
+ ,1
+ ,0
+ ,0
+ ,0
+ ,43410
+ ,0
+ ,239
+ ,0
+ ,1141
+ ,0
+ ,49
+ ,0
+ ,16
+ ,0
+ ,42
+ ,0
+ ,1
+ ,0
+ ,104838
+ ,0
+ ,240
+ ,0
+ ,680
+ ,0
+ ,27
+ ,0
+ ,10
+ ,0
+ ,26
+ ,0
+ ,0
+ ,0
+ ,62215
+ ,0
+ ,241
+ ,0
+ ,1090
+ ,0
+ ,30
+ ,0
+ ,19
+ ,0
+ ,21
+ ,0
+ ,6
+ ,0
+ ,69304
+ ,0
+ ,242
+ ,0
+ ,616
+ ,0
+ ,22
+ ,0
+ ,12
+ ,0
+ ,4
+ ,0
+ ,3
+ ,0
+ ,53117
+ ,0
+ ,243
+ ,0
+ ,285
+ ,0
+ ,12
+ ,0
+ ,2
+ ,0
+ ,10
+ ,0
+ ,1
+ ,0
+ ,19764
+ ,0
+ ,244
+ ,0
+ ,1145
+ ,0
+ ,31
+ ,0
+ ,14
+ ,0
+ ,43
+ ,0
+ ,2
+ ,0
+ ,86680
+ ,0
+ ,245
+ ,0
+ ,733
+ ,0
+ ,20
+ ,0
+ ,17
+ ,0
+ ,34
+ ,0
+ ,0
+ ,0
+ ,84105
+ ,0
+ ,246
+ ,0
+ ,888
+ ,0
+ ,20
+ ,0
+ ,19
+ ,0
+ ,31
+ ,0
+ ,0
+ ,0
+ ,77945
+ ,0
+ ,247
+ ,0
+ ,849
+ ,0
+ ,39
+ ,0
+ ,14
+ ,0
+ ,19
+ ,0
+ ,0
+ ,0
+ ,89113
+ ,0
+ ,248
+ ,0
+ ,1182
+ ,0
+ ,29
+ ,0
+ ,11
+ ,0
+ ,34
+ ,0
+ ,3
+ ,0
+ ,91005
+ ,0
+ ,249
+ ,0
+ ,528
+ ,0
+ ,16
+ ,0
+ ,4
+ ,0
+ ,6
+ ,0
+ ,1
+ ,0
+ ,40248
+ ,0
+ ,250
+ ,0
+ ,642
+ ,0
+ ,27
+ ,0
+ ,16
+ ,0
+ ,11
+ ,0
+ ,0
+ ,0
+ ,64187
+ ,0
+ ,251
+ ,0
+ ,947
+ ,0
+ ,21
+ ,0
+ ,20
+ ,0
+ ,24
+ ,0
+ ,0
+ ,0
+ ,50857
+ ,0
+ ,252
+ ,0
+ ,819
+ ,0
+ ,19
+ ,0
+ ,12
+ ,0
+ ,16
+ ,0
+ ,1
+ ,0
+ ,56613
+ ,0
+ ,253
+ ,0
+ ,757
+ ,0
+ ,35
+ ,0
+ ,15
+ ,0
+ ,72
+ ,0
+ ,0
+ ,0
+ ,62792
+ ,0
+ ,254
+ ,0
+ ,894
+ ,0
+ ,14
+ ,0
+ ,16
+ ,0
+ ,21
+ ,0
+ ,0
+ ,0
+ ,72535
+ ,0
+ ,255
+ ,0
+ ,1192
+ ,0
+ ,63
+ ,0
+ ,11
+ ,0
+ ,28
+ ,0
+ ,0
+ ,0
+ ,101645
+ ,0
+ ,256
+ ,0
+ ,186
+ ,0
+ ,17
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,7176
+ ,0
+ ,257
+ ,0
+ ,1793
+ ,0
+ ,76
+ ,0
+ ,17
+ ,0
+ ,47
+ ,0
+ ,0
+ ,0
+ ,96560
+ ,0
+ ,258
+ ,0
+ ,829
+ ,0
+ ,34
+ ,0
+ ,13
+ ,0
+ ,32
+ ,0
+ ,0
+ ,0
+ ,101011
+ ,0
+ ,259
+ ,0
+ ,2702
+ ,0
+ ,107
+ ,0
+ ,20
+ ,0
+ ,65
+ ,0
+ ,0
+ ,0
+ ,175824
+ ,0
+ ,260
+ ,0
+ ,4691
+ ,0
+ ,168
+ ,0
+ ,21
+ ,0
+ ,123
+ ,0
+ ,1
+ ,0
+ ,341570
+ ,0
+ ,261
+ ,0
+ ,1112
+ ,0
+ ,43
+ ,0
+ ,16
+ ,0
+ ,26
+ ,0
+ ,1
+ ,0
+ ,103597
+ ,0
+ ,262
+ ,0
+ ,937
+ ,0
+ ,41
+ ,0
+ ,20
+ ,0
+ ,48
+ ,0
+ ,0
+ ,0
+ ,112611
+ ,0
+ ,263
+ ,0
+ ,1290
+ ,0
+ ,34
+ ,0
+ ,21
+ ,0
+ ,37
+ ,0
+ ,0
+ ,0
+ ,85574
+ ,0
+ ,264
+ ,0
+ ,2146
+ ,0
+ ,75
+ ,0
+ ,18
+ ,0
+ ,60
+ ,0
+ ,1
+ ,0
+ ,220801
+ ,0
+ ,265
+ ,0
+ ,1590
+ ,0
+ ,61
+ ,0
+ ,17
+ ,0
+ ,39
+ ,0
+ ,1
+ ,0
+ ,92661
+ ,0
+ ,266
+ ,0
+ ,1590
+ ,0
+ ,55
+ ,0
+ ,20
+ ,0
+ ,64
+ ,0
+ ,0
+ ,0
+ ,133328
+ ,0
+ ,267
+ ,0
+ ,1210
+ ,0
+ ,77
+ ,0
+ ,12
+ ,0
+ ,26
+ ,0
+ ,0
+ ,0
+ ,61361
+ ,0
+ ,268
+ ,0
+ ,2072
+ ,0
+ ,75
+ ,0
+ ,17
+ ,0
+ ,64
+ ,0
+ ,4
+ ,0
+ ,125930
+ ,0
+ ,269
+ ,0
+ ,834
+ ,0
+ ,32
+ ,0
+ ,10
+ ,0
+ ,25
+ ,0
+ ,4
+ ,0
+ ,82316
+ ,0
+ ,270
+ ,0
+ ,1105
+ ,0
+ ,53
+ ,0
+ ,13
+ ,0
+ ,26
+ ,0
+ ,3
+ ,0
+ ,102010
+ ,0
+ ,271
+ ,0
+ ,1272
+ ,0
+ ,42
+ ,0
+ ,22
+ ,0
+ ,76
+ ,0
+ ,0
+ ,0
+ ,101523
+ ,0
+ ,272
+ ,0
+ ,761
+ ,0
+ ,35
+ ,0
+ ,9
+ ,0
+ ,2
+ ,0
+ ,5
+ ,0
+ ,41566
+ ,0
+ ,273
+ ,0
+ ,1988
+ ,0
+ ,66
+ ,0
+ ,25
+ ,0
+ ,36
+ ,0
+ ,0
+ ,0
+ ,99923
+ ,0
+ ,274
+ ,0
+ ,620
+ ,0
+ ,19
+ ,0
+ ,13
+ ,0
+ ,23
+ ,0
+ ,0
+ ,0
+ ,22648
+ ,0
+ ,275
+ ,0
+ ,800
+ ,0
+ ,45
+ ,0
+ ,13
+ ,0
+ ,14
+ ,0
+ ,0
+ ,0
+ ,46698
+ ,0
+ ,276
+ ,0
+ ,1684
+ ,0
+ ,65
+ ,0
+ ,19
+ ,0
+ ,78
+ ,0
+ ,0
+ ,0
+ ,131698
+ ,0
+ ,277
+ ,0
+ ,1050
+ ,0
+ ,35
+ ,0
+ ,18
+ ,0
+ ,14
+ ,0
+ ,0
+ ,0
+ ,91735
+ ,0
+ ,278
+ ,0
+ ,1639
+ ,0
+ ,36
+ ,0
+ ,16
+ ,0
+ ,56
+ ,0
+ ,1
+ ,0
+ ,133368
+ ,0
+ ,279
+ ,0
+ ,1254
+ ,0
+ ,31
+ ,0
+ ,14
+ ,0
+ ,29
+ ,0
+ ,2
+ ,0
+ ,135781
+ ,0
+ ,280
+ ,0
+ ,1335
+ ,0
+ ,67
+ ,0
+ ,7
+ ,0
+ ,35
+ ,0
+ ,4
+ ,0
+ ,74408
+ ,0
+ ,281
+ ,0
+ ,1597
+ ,0
+ ,66
+ ,0
+ ,17
+ ,0
+ ,13
+ ,0
+ ,0
+ ,0
+ ,81240
+ ,0
+ ,282
+ ,0
+ ,961
+ ,0
+ ,32
+ ,0
+ ,18
+ ,0
+ ,21
+ ,0
+ ,0
+ ,0
+ ,92499
+ ,0
+ ,283
+ ,0
+ ,1577
+ ,0
+ ,30
+ ,0
+ ,13
+ ,0
+ ,40
+ ,0
+ ,0
+ ,0
+ ,136084
+ ,0
+ ,284
+ ,0
+ ,946
+ ,0
+ ,64
+ ,0
+ ,20
+ ,0
+ ,61
+ ,0
+ ,0
+ ,0
+ ,116048
+ ,0
+ ,285
+ ,0
+ ,1926
+ ,0
+ ,81
+ ,0
+ ,18
+ ,0
+ ,49
+ ,0
+ ,0
+ ,0
+ ,250047
+ ,0
+ ,286
+ ,0
+ ,1417
+ ,0
+ ,118
+ ,0
+ ,16
+ ,0
+ ,57
+ ,0
+ ,0
+ ,0
+ ,120445
+ ,0
+ ,287
+ ,0
+ ,1060
+ ,0
+ ,18
+ ,0
+ ,13
+ ,0
+ ,50
+ ,0
+ ,0
+ ,0
+ ,98866
+ ,0
+ ,288
+ ,0
+ ,1421
+ ,0
+ ,62
+ ,0
+ ,14
+ ,0
+ ,39
+ ,0
+ ,1
+ ,0
+ ,108043
+ ,0
+ ,289
+ ,0
+ ,1502
+ ,0
+ ,37
+ ,0
+ ,22
+ ,0
+ ,24
+ ,0
+ ,1
+ ,0
+ ,79863)
+ ,dim=c(14
+ ,289)
+ ,dimnames=list(c('Pop'
+ ,'t'
+ ,'pop_t'
+ ,'Pageviews'
+ ,'Pageviews_p'
+ ,'Logins'
+ ,'Logins_p'
+ ,'Compendiums_Reviewed'
+ ,'Compendiums_Reviewed_p'
+ ,'Totblogs'
+ ,'Totblogs_p'
+ ,'Shared_Compendiums'
+ ,'Shared_Compendiums_p'
+ ,'Time_in_rfc')
+ ,1:289))
> y <- array(NA,dim=c(14,289),dimnames=list(c('Pop','t','pop_t','Pageviews','Pageviews_p','Logins','Logins_p','Compendiums_Reviewed','Compendiums_Reviewed_p','Totblogs','Totblogs_p','Shared_Compendiums','Shared_Compendiums_p','Time_in_rfc'),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 = '14'
> 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 Pop t pop_t Pageviews Pageviews_p Logins Logins_p
1 210907 1 1 1 1418 1418 56 56
2 120982 1 2 2 869 869 56 56
3 176508 1 3 3 1530 1530 54 54
4 179321 1 4 4 2172 2172 89 89
5 123185 1 5 5 901 901 40 40
6 52746 1 6 6 463 463 25 25
7 385534 1 7 7 3201 3201 92 92
8 33170 1 8 8 371 371 18 18
9 149061 1 9 9 1583 1583 44 44
10 165446 1 10 10 1439 1439 33 33
11 237213 1 11 11 1764 1764 84 84
12 173326 1 12 12 1495 1495 88 88
13 133131 1 13 13 1373 1373 55 55
14 258873 1 14 14 2187 2187 60 60
15 180083 1 15 15 1491 1491 66 66
16 324799 1 16 16 4041 4041 154 154
17 230964 1 17 17 1706 1706 53 53
18 236785 1 18 18 2152 2152 119 119
19 135473 1 19 19 1036 1036 41 41
20 202925 1 20 20 1882 1882 61 61
21 215147 1 21 21 1929 1929 58 58
22 344297 1 22 22 2242 2242 75 75
23 153935 1 23 23 1220 1220 33 33
24 132943 1 24 24 1289 1289 40 40
25 174724 1 25 25 2515 2515 92 92
26 174415 1 26 26 2147 2147 100 100
27 225548 1 27 27 2352 2352 112 112
28 223632 1 28 28 1638 1638 73 73
29 124817 1 29 29 1222 1222 40 40
30 221698 1 30 30 1812 1812 45 45
31 210767 1 31 31 1677 1677 60 60
32 170266 1 32 32 1579 1579 62 62
33 260561 1 33 33 1731 1731 75 75
34 84853 1 34 34 807 807 31 31
35 294424 1 35 35 2452 2452 77 77
36 215641 1 36 36 1940 1940 46 46
37 325107 1 37 37 2662 2662 99 99
38 167542 1 38 38 1499 1499 66 66
39 106408 1 39 39 865 865 30 30
40 265769 1 40 40 2527 2527 146 146
41 269651 1 41 41 2747 2747 67 67
42 149112 1 42 42 1324 1324 56 56
43 152871 1 43 43 1383 1383 58 58
44 111665 1 44 44 1179 1179 34 34
45 116408 1 45 45 2099 2099 61 61
46 362301 1 46 46 4308 4308 119 119
47 78800 1 47 47 918 918 42 42
48 183167 1 48 48 1831 1831 66 66
49 277965 1 49 49 3373 3373 89 89
50 150629 1 50 50 1713 1713 44 44
51 168809 1 51 51 1438 1438 66 66
52 24188 1 52 52 496 496 24 24
53 329267 1 53 53 2253 2253 259 259
54 65029 1 54 54 744 744 17 17
55 101097 1 55 55 1161 1161 64 64
56 218946 1 56 56 2352 2352 41 41
57 244052 1 57 57 2144 2144 68 68
58 233328 1 58 58 2694 2694 132 132
59 256462 1 59 59 1973 1973 105 105
60 206161 1 60 60 1769 1769 71 71
61 311473 1 61 61 3148 3148 112 112
62 235800 1 62 62 2474 2474 94 94
63 177939 1 63 63 2084 2084 82 82
64 207176 1 64 64 1954 1954 70 70
65 196553 1 65 65 1226 1226 57 57
66 174184 1 66 66 1389 1389 53 53
67 143246 1 67 67 1496 1496 103 103
68 187559 1 68 68 2269 2269 121 121
69 187681 1 69 69 1833 1833 62 62
70 119016 1 70 70 1268 1268 52 52
71 182192 1 71 71 1943 1943 52 52
72 73566 1 72 72 893 893 32 32
73 194979 1 73 73 1762 1762 62 62
74 167488 1 74 74 1403 1403 45 45
75 143756 1 75 75 1425 1425 46 46
76 275541 1 76 76 1857 1857 63 63
77 243199 1 77 77 1840 1840 75 75
78 182999 1 78 78 1502 1502 88 88
79 135649 1 79 79 1441 1441 46 46
80 152299 1 80 80 1420 1420 53 53
81 120221 1 81 81 1416 1416 37 37
82 346485 1 82 82 2970 2970 90 90
83 145790 1 83 83 1317 1317 63 63
84 193339 1 84 84 1644 1644 78 78
85 80953 1 85 85 870 870 25 25
86 122774 1 86 86 1654 1654 45 45
87 130585 1 87 87 1054 1054 46 46
88 286468 1 88 88 3004 3004 144 144
89 241066 1 89 89 2008 2008 82 82
90 148446 1 90 90 2547 2547 91 91
91 204713 1 91 91 1885 1885 71 71
92 182079 1 92 92 1626 1626 63 63
93 140344 1 93 93 1468 1468 53 53
94 220516 1 94 94 2445 2445 62 62
95 243060 1 95 95 1964 1964 63 63
96 162765 1 96 96 1381 1381 32 32
97 182613 1 97 97 1369 1369 39 39
98 232138 1 98 98 1659 1659 62 62
99 265318 1 99 99 2888 2888 117 117
100 310839 1 100 100 2845 2845 92 92
101 225060 1 101 101 1982 1982 93 93
102 232317 1 102 102 1904 1904 54 54
103 144966 1 103 103 1391 1391 144 144
104 43287 1 104 104 602 602 14 14
105 155754 1 105 105 1743 1743 61 61
106 164709 1 106 106 1559 1559 109 109
107 201940 1 107 107 2014 2014 38 38
108 235454 1 108 108 2143 2143 73 73
109 99466 1 109 109 874 874 50 50
110 100750 1 110 110 1281 1281 72 72
111 224549 1 111 111 1401 1401 50 50
112 243511 1 112 112 1944 1944 71 71
113 22938 1 113 113 391 391 10 10
114 152474 1 114 114 1605 1605 65 65
115 61857 1 115 115 530 530 25 25
116 132487 1 116 116 1386 1386 41 41
117 317394 1 117 117 2395 2395 86 86
118 21054 1 118 118 387 387 16 16
119 209641 1 119 119 1742 1742 42 42
120 31414 1 120 120 449 449 19 19
121 244749 1 121 121 2699 2699 95 95
122 184510 1 122 122 1606 1606 49 49
123 128423 1 123 123 1204 1204 64 64
124 97839 1 124 124 1138 1138 38 38
125 38214 1 125 125 568 568 34 34
126 151101 1 126 126 1459 1459 32 32
127 272458 1 127 127 2158 2158 65 65
128 172494 1 128 128 1111 1111 52 52
129 328107 1 129 129 2833 2833 65 65
130 250579 1 130 130 1955 1955 83 83
131 351067 1 131 131 2922 2922 95 95
132 158015 1 132 132 1002 1002 29 29
133 85439 1 133 133 956 956 33 33
134 229242 1 134 134 2186 2186 247 247
135 351619 1 135 135 3604 3604 139 139
136 84207 1 136 136 1035 1035 29 29
137 324598 1 137 137 3261 3261 110 110
138 131069 1 138 138 1587 1587 67 67
139 204271 1 139 139 1424 1424 42 42
140 165543 1 140 140 1701 1701 65 65
141 141722 1 141 141 1249 1249 94 94
142 299775 1 142 142 3352 3352 95 95
143 195838 1 143 143 1641 1641 67 67
144 173260 1 144 144 2035 2035 63 63
145 254488 1 145 145 2312 2312 83 83
146 104389 1 146 146 1369 1369 45 45
147 199476 1 147 147 2201 2201 70 70
148 224330 1 148 148 1900 1900 83 83
149 14688 1 149 149 207 207 10 10
150 181633 1 150 150 1645 1645 70 70
151 271856 1 151 151 2429 2429 103 103
152 7199 1 152 152 151 151 5 5
153 46660 1 153 153 474 474 20 20
154 17547 1 154 154 141 141 5 5
155 95227 1 155 155 872 872 34 34
156 152601 1 156 156 1318 1318 48 48
157 98146 0 157 0 1018 0 40 0
158 79619 0 158 0 1383 0 43 0
159 59194 0 159 0 1314 0 31 0
160 139942 0 160 0 1335 0 42 0
161 118612 0 161 0 1403 0 46 0
162 72880 0 162 0 910 0 33 0
163 65475 0 163 0 616 0 18 0
164 99643 0 164 0 1407 0 55 0
165 71965 0 165 0 771 0 35 0
166 77272 0 166 0 766 0 59 0
167 49289 0 167 0 473 0 19 0
168 135131 0 168 0 1376 0 66 0
169 108446 0 169 0 1232 0 60 0
170 89746 0 170 0 1521 0 36 0
171 44296 0 171 0 572 0 25 0
172 77648 0 172 0 1059 0 47 0
173 181528 0 173 0 1544 0 54 0
174 134019 0 174 0 1230 0 53 0
175 124064 0 175 0 1206 0 40 0
176 92630 0 176 0 1205 0 40 0
177 121848 0 177 0 1255 0 39 0
178 52915 0 178 0 613 0 14 0
179 81872 0 179 0 721 0 45 0
180 58981 0 180 0 1109 0 36 0
181 53515 0 181 0 740 0 28 0
182 60812 0 182 0 1126 0 44 0
183 56375 0 183 0 728 0 30 0
184 65490 0 184 0 689 0 22 0
185 80949 0 185 0 592 0 17 0
186 76302 0 186 0 995 0 31 0
187 104011 0 187 0 1613 0 55 0
188 98104 0 188 0 2048 0 54 0
189 67989 0 189 0 705 0 21 0
190 30989 0 190 0 301 0 14 0
191 135458 0 191 0 1803 0 81 0
192 73504 0 192 0 799 0 35 0
193 63123 0 193 0 861 0 43 0
194 61254 0 194 0 1186 0 46 0
195 74914 0 195 0 1451 0 30 0
196 31774 0 196 0 628 0 23 0
197 81437 0 197 0 1161 0 38 0
198 87186 0 198 0 1463 0 54 0
199 50090 0 199 0 742 0 20 0
200 65745 0 200 0 979 0 53 0
201 56653 0 201 0 675 0 45 0
202 158399 0 202 0 1241 0 39 0
203 46455 0 203 0 676 0 20 0
204 73624 0 204 0 1049 0 24 0
205 38395 0 205 0 620 0 31 0
206 91899 0 206 0 1081 0 35 0
207 139526 0 207 0 1688 0 151 0
208 52164 0 208 0 736 0 52 0
209 51567 0 209 0 617 0 30 0
210 70551 0 210 0 812 0 31 0
211 84856 0 211 0 1051 0 29 0
212 102538 0 212 0 1656 0 57 0
213 86678 0 213 0 705 0 40 0
214 85709 0 214 0 945 0 44 0
215 34662 0 215 0 554 0 25 0
216 150580 0 216 0 1597 0 77 0
217 99611 0 217 0 982 0 35 0
218 19349 0 218 0 222 0 11 0
219 99373 0 219 0 1212 0 63 0
220 86230 0 220 0 1143 0 44 0
221 30837 0 221 0 435 0 19 0
222 31706 0 222 0 532 0 13 0
223 89806 0 223 0 882 0 42 0
224 62088 0 224 0 608 0 38 0
225 40151 0 225 0 459 0 29 0
226 27634 0 226 0 578 0 20 0
227 76990 0 227 0 826 0 27 0
228 37460 0 228 0 509 0 20 0
229 54157 0 229 0 717 0 19 0
230 49862 0 230 0 637 0 37 0
231 84337 0 231 0 857 0 26 0
232 64175 0 232 0 830 0 42 0
233 59382 0 233 0 652 0 49 0
234 119308 0 234 0 707 0 30 0
235 76702 0 235 0 954 0 49 0
236 103425 0 236 0 1461 0 67 0
237 70344 0 237 0 672 0 28 0
238 43410 0 238 0 778 0 19 0
239 104838 0 239 0 1141 0 49 0
240 62215 0 240 0 680 0 27 0
241 69304 0 241 0 1090 0 30 0
242 53117 0 242 0 616 0 22 0
243 19764 0 243 0 285 0 12 0
244 86680 0 244 0 1145 0 31 0
245 84105 0 245 0 733 0 20 0
246 77945 0 246 0 888 0 20 0
247 89113 0 247 0 849 0 39 0
248 91005 0 248 0 1182 0 29 0
249 40248 0 249 0 528 0 16 0
250 64187 0 250 0 642 0 27 0
251 50857 0 251 0 947 0 21 0
252 56613 0 252 0 819 0 19 0
253 62792 0 253 0 757 0 35 0
254 72535 0 254 0 894 0 14 0
255 101645 0 255 0 1192 0 63 0
256 7176 0 256 0 186 0 17 0
257 96560 0 257 0 1793 0 76 0
258 101011 0 258 0 829 0 34 0
259 175824 0 259 0 2702 0 107 0
260 341570 0 260 0 4691 0 168 0
261 103597 0 261 0 1112 0 43 0
262 112611 0 262 0 937 0 41 0
263 85574 0 263 0 1290 0 34 0
264 220801 0 264 0 2146 0 75 0
265 92661 0 265 0 1590 0 61 0
266 133328 0 266 0 1590 0 55 0
267 61361 0 267 0 1210 0 77 0
268 125930 0 268 0 2072 0 75 0
269 82316 0 269 0 834 0 32 0
270 102010 0 270 0 1105 0 53 0
271 101523 0 271 0 1272 0 42 0
272 41566 0 272 0 761 0 35 0
273 99923 0 273 0 1988 0 66 0
274 22648 0 274 0 620 0 19 0
275 46698 0 275 0 800 0 45 0
276 131698 0 276 0 1684 0 65 0
277 91735 0 277 0 1050 0 35 0
278 133368 0 278 0 1639 0 36 0
279 135781 0 279 0 1254 0 31 0
280 74408 0 280 0 1335 0 67 0
281 81240 0 281 0 1597 0 66 0
282 92499 0 282 0 961 0 32 0
283 136084 0 283 0 1577 0 30 0
284 116048 0 284 0 946 0 64 0
285 250047 0 285 0 1926 0 81 0
286 120445 0 286 0 1417 0 118 0
287 98866 0 287 0 1060 0 18 0
288 108043 0 288 0 1421 0 62 0
289 79863 0 289 0 1502 0 37 0
Compendiums_Reviewed Compendiums_Reviewed_p Totblogs Totblogs_p
1 30 30 145 145
2 28 28 101 101
3 38 38 98 98
4 30 30 132 132
5 22 22 60 60
6 26 26 38 38
7 25 25 144 144
8 18 18 5 5
9 26 26 84 84
10 25 25 79 79
11 38 38 127 127
12 44 44 78 78
13 30 30 60 60
14 40 40 131 131
15 34 34 84 84
16 47 47 133 133
17 30 30 150 150
18 31 31 91 91
19 23 23 132 132
20 36 36 136 136
21 36 36 124 124
22 30 30 118 118
23 25 25 70 70
24 39 39 107 107
25 34 34 119 119
26 31 31 89 89
27 31 31 112 112
28 33 33 108 108
29 25 25 52 52
30 33 33 112 112
31 35 35 116 116
32 42 42 123 123
33 43 43 125 125
34 30 30 27 27
35 33 33 162 162
36 32 32 64 64
37 36 36 92 92
38 28 28 83 83
39 14 14 41 41
40 32 32 120 120
41 30 30 105 105
42 35 35 79 79
43 28 28 70 70
44 28 28 55 55
45 39 39 39 39
46 34 34 67 67
47 26 26 21 21
48 39 39 127 127
49 39 39 152 152
50 33 33 113 113
51 28 28 99 99
52 4 4 7 7
53 39 39 141 141
54 18 18 21 21
55 14 14 35 35
56 29 29 109 109
57 44 44 133 133
58 28 28 230 230
59 35 35 166 166
60 28 28 68 68
61 38 38 147 147
62 23 23 179 179
63 36 36 61 61
64 32 32 101 101
65 29 29 108 108
66 25 25 90 90
67 27 27 114 114
68 36 36 103 103
69 28 28 142 142
70 23 23 79 79
71 40 40 88 88
72 23 23 25 25
73 40 40 83 83
74 28 28 113 113
75 34 34 118 118
76 33 33 110 110
77 28 28 129 129
78 34 34 51 51
79 30 30 93 93
80 33 33 76 76
81 22 22 49 49
82 38 38 118 118
83 26 26 38 38
84 35 35 141 141
85 8 8 58 58
86 24 24 27 27
87 29 29 91 91
88 29 29 63 63
89 45 45 56 56
90 37 37 144 144
91 33 33 73 73
92 33 33 168 168
93 25 25 64 64
94 32 32 97 97
95 29 29 117 117
96 28 28 100 100
97 28 28 149 149
98 31 31 187 187
99 52 52 127 127
100 24 24 245 245
101 41 41 87 87
102 33 33 177 177
103 32 32 49 49
104 19 19 49 49
105 20 20 73 73
106 31 31 177 177
107 31 31 94 94
108 32 32 117 117
109 23 23 55 55
110 30 30 58 58
111 31 31 95 95
112 42 42 129 129
113 1 1 11 11
114 32 32 101 101
115 11 11 28 28
116 36 36 89 89
117 31 31 193 193
118 0 0 4 4
119 24 24 84 84
120 8 8 39 39
121 33 33 101 101
122 40 40 82 82
123 38 38 36 36
124 24 24 75 75
125 8 8 16 16
126 35 35 55 55
127 43 43 131 131
128 43 43 131 131
129 41 41 144 144
130 38 38 139 139
131 45 45 211 211
132 31 31 78 78
133 28 28 39 39
134 31 31 90 90
135 40 40 166 166
136 30 30 12 12
137 37 37 133 133
138 30 30 69 69
139 35 35 119 119
140 32 32 119 119
141 27 27 65 65
142 31 31 101 101
143 31 31 196 196
144 21 21 15 15
145 39 39 136 136
146 41 41 89 89
147 32 32 123 123
148 39 39 163 163
149 0 0 5 5
150 30 30 96 96
151 37 37 151 151
152 0 0 6 6
153 5 5 13 13
154 1 1 3 3
155 32 32 23 23
156 24 24 57 57
157 17 0 14 0
158 11 0 43 0
159 24 0 20 0
160 22 0 72 0
161 12 0 87 0
162 19 0 21 0
163 13 0 56 0
164 17 0 59 0
165 15 0 82 0
166 16 0 43 0
167 24 0 25 0
168 15 0 38 0
169 17 0 25 0
170 18 0 38 0
171 20 0 12 0
172 16 0 29 0
173 16 0 47 0
174 18 0 45 0
175 22 0 40 0
176 8 0 30 0
177 17 0 41 0
178 18 0 25 0
179 16 0 23 0
180 23 0 14 0
181 22 0 16 0
182 13 0 26 0
183 13 0 21 0
184 16 0 27 0
185 16 0 9 0
186 20 0 33 0
187 22 0 42 0
188 17 0 68 0
189 18 0 32 0
190 17 0 6 0
191 12 0 67 0
192 7 0 33 0
193 17 0 77 0
194 14 0 46 0
195 23 0 30 0
196 17 0 0 0
197 14 0 36 0
198 15 0 46 0
199 17 0 18 0
200 21 0 48 0
201 18 0 29 0
202 18 0 28 0
203 17 0 34 0
204 17 0 33 0
205 16 0 34 0
206 15 0 33 0
207 21 0 80 0
208 16 0 32 0
209 14 0 30 0
210 15 0 41 0
211 17 0 41 0
212 15 0 51 0
213 15 0 18 0
214 10 0 34 0
215 6 0 31 0
216 22 0 39 0
217 21 0 54 0
218 1 0 14 0
219 18 0 24 0
220 17 0 24 0
221 4 0 8 0
222 10 0 26 0
223 16 0 19 0
224 16 0 11 0
225 9 0 14 0
226 16 0 1 0
227 17 0 39 0
228 7 0 5 0
229 15 0 37 0
230 14 0 32 0
231 14 0 38 0
232 18 0 47 0
233 12 0 47 0
234 16 0 37 0
235 21 0 51 0
236 19 0 45 0
237 16 0 21 0
238 1 0 1 0
239 16 0 42 0
240 10 0 26 0
241 19 0 21 0
242 12 0 4 0
243 2 0 10 0
244 14 0 43 0
245 17 0 34 0
246 19 0 31 0
247 14 0 19 0
248 11 0 34 0
249 4 0 6 0
250 16 0 11 0
251 20 0 24 0
252 12 0 16 0
253 15 0 72 0
254 16 0 21 0
255 11 0 28 0
256 0 0 0 0
257 17 0 47 0
258 13 0 32 0
259 20 0 65 0
260 21 0 123 0
261 16 0 26 0
262 20 0 48 0
263 21 0 37 0
264 18 0 60 0
265 17 0 39 0
266 20 0 64 0
267 12 0 26 0
268 17 0 64 0
269 10 0 25 0
270 13 0 26 0
271 22 0 76 0
272 9 0 2 0
273 25 0 36 0
274 13 0 23 0
275 13 0 14 0
276 19 0 78 0
277 18 0 14 0
278 16 0 56 0
279 14 0 29 0
280 7 0 35 0
281 17 0 13 0
282 18 0 21 0
283 13 0 40 0
284 20 0 61 0
285 18 0 49 0
286 16 0 57 0
287 13 0 50 0
288 14 0 39 0
289 22 0 24 0
Shared_Compendiums Shared_Compendiums_p
1 3 3
2 4 4
3 12 12
4 2 2
5 1 1
6 3 3
7 0 0
8 0 0
9 5 5
10 0 0
11 0 0
12 7 7
13 7 7
14 3 3
15 9 9
16 0 0
17 4 4
18 3 3
19 0 0
20 7 7
21 0 0
22 1 1
23 5 5
24 7 7
25 0 0
26 0 0
27 5 5
28 0 0
29 0 0
30 0 0
31 3 3
32 4 4
33 1 1
34 4 4
35 2 2
36 0 0
37 0 0
38 2 2
39 1 1
40 2 2
41 10 10
42 6 6
43 5 5
44 4 4
45 1 1
46 2 2
47 2 2
48 0 0
49 8 8
50 3 3
51 0 0
52 0 0
53 8 8
54 5 5
55 3 3
56 1 1
57 5 5
58 5 5
59 0 0
60 12 12
61 8 8
62 8 8
63 8 8
64 8 8
65 2 2
66 0 0
67 5 5
68 8 8
69 2 2
70 5 5
71 12 12
72 6 6
73 7 7
74 2 2
75 0 0
76 4 4
77 3 3
78 6 6
79 2 2
80 0 0
81 1 1
82 0 0
83 5 5
84 2 2
85 0 0
86 0 0
87 5 5
88 1 1
89 0 0
90 1 1
91 1 1
92 2 2
93 6 6
94 1 1
95 4 4
96 2 2
97 3 3
98 0 0
99 10 10
100 9 9
101 7 7
102 0 0
103 0 0
104 4 4
105 4 4
106 0 0
107 0 0
108 0 0
109 0 0
110 0 0
111 4 4
112 0 0
113 0 0
114 0 0
115 4 4
116 0 0
117 1 1
118 0 0
119 5 5
120 0 0
121 2 2
122 7 7
123 8 8
124 2 2
125 0 0
126 2 2
127 0 0
128 0 0
129 3 3
130 0 0
131 3 3
132 0 0
133 0 0
134 4 4
135 4 4
136 11 11
137 0 0
138 4 4
139 0 0
140 1 1
141 0 0
142 9 9
143 1 1
144 3 3
145 10 10
146 5 5
147 2 2
148 1 1
149 0 0
150 2 2
151 1 1
152 0 0
153 0 0
154 0 0
155 0 0
156 2 2
157 0 0
158 3 0
159 6 0
160 0 0
161 2 0
162 0 0
163 2 0
164 1 0
165 1 0
166 2 0
167 1 0
168 0 0
169 1 0
170 3 0
171 0 0
172 0 0
173 0 0
174 0 0
175 1 0
176 4 0
177 0 0
178 0 0
179 0 0
180 7 0
181 2 0
182 0 0
183 7 0
184 3 0
185 0 0
186 0 0
187 6 0
188 2 0
189 0 0
190 0 0
191 3 0
192 0 0
193 1 0
194 1 0
195 0 0
196 1 0
197 0 0
198 0 0
199 0 0
200 0 0
201 0 0
202 0 0
203 0 0
204 0 0
205 0 0
206 0 0
207 0 0
208 0 0
209 2 0
210 0 0
211 1 0
212 1 0
213 0 0
214 0 0
215 0 0
216 0 0
217 0 0
218 0 0
219 1 0
220 0 0
221 0 0
222 0 0
223 0 0
224 1 0
225 0 0
226 0 0
227 0 0
228 0 0
229 0 0
230 0 0
231 0 0
232 0 0
233 0 0
234 0 0
235 0 0
236 1 0
237 0 0
238 0 0
239 1 0
240 0 0
241 6 0
242 3 0
243 1 0
244 2 0
245 0 0
246 0 0
247 0 0
248 3 0
249 1 0
250 0 0
251 0 0
252 1 0
253 0 0
254 0 0
255 0 0
256 0 0
257 0 0
258 0 0
259 0 0
260 1 0
261 1 0
262 0 0
263 0 0
264 1 0
265 1 0
266 0 0
267 0 0
268 4 0
269 4 0
270 3 0
271 0 0
272 5 0
273 0 0
274 0 0
275 0 0
276 0 0
277 0 0
278 1 0
279 2 0
280 4 0
281 0 0
282 0 0
283 0 0
284 0 0
285 0 0
286 0 0
287 0 0
288 1 0
289 1 0
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) Pop t
2869.07 -4793.20 15.57
pop_t Pageviews Pageviews_p
-57.63 56.86 12.91
Logins Logins_p Compendiums_Reviewed
127.85 -14.15 380.42
Compendiums_Reviewed_p Totblogs Totblogs_p
467.34 283.91 120.38
Shared_Compendiums Shared_Compendiums_p
-3045.40 2074.97
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-122511 -14940 -1226 14356 110025
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 2869.066 18173.224 0.158 0.8747
Pop -4793.199 20243.083 -0.237 0.8130
t 15.574 67.835 0.230 0.8186
pop_t -57.632 84.553 -0.682 0.4961
Pageviews 56.858 8.460 6.721 1.04e-10 ***
Pageviews_p 12.912 9.765 1.322 0.1872
Logins 127.852 169.812 0.753 0.4522
Logins_p -14.148 189.750 -0.075 0.9406
Compendiums_Reviewed 380.421 569.415 0.668 0.5046
Compendiums_Reviewed_p 467.339 651.607 0.717 0.4739
Totblogs 283.911 164.605 1.725 0.0857 .
Totblogs_p 120.377 175.696 0.685 0.4938
Shared_Compendiums -3045.403 1647.321 -1.849 0.0656 .
Shared_Compendiums_p 2074.970 1808.479 1.147 0.2522
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 27950 on 275 degrees of freedom
Multiple R-squared: 0.89, Adjusted R-squared: 0.8848
F-statistic: 171.1 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.9972629530 5.474094e-03 2.737047e-03
[2,] 0.9951400589 9.719882e-03 4.859941e-03
[3,] 0.9964510510 7.097898e-03 3.548949e-03
[4,] 0.9943282402 1.134352e-02 5.671760e-03
[5,] 0.9889001516 2.219970e-02 1.109985e-02
[6,] 0.9998498765 3.002470e-04 1.501235e-04
[7,] 0.9996979395 6.041210e-04 3.020605e-04
[8,] 0.9996460591 7.078818e-04 3.539409e-04
[9,] 0.9999898866 2.022673e-05 1.011337e-05
[10,] 0.9999879736 2.405285e-05 1.202643e-05
[11,] 0.9999759927 4.801456e-05 2.400728e-05
[12,] 0.9999838468 3.230641e-05 1.615321e-05
[13,] 0.9999678238 6.435237e-05 3.217619e-05
[14,] 0.9999466329 1.067342e-04 5.336709e-05
[15,] 0.9999120898 1.758204e-04 8.791019e-05
[16,] 0.9998636949 2.726102e-04 1.363051e-04
[17,] 0.9999558318 8.833646e-05 4.416823e-05
[18,] 0.9999205424 1.589152e-04 7.945762e-05
[19,] 0.9998747568 2.504865e-04 1.252432e-04
[20,] 0.9998336744 3.326512e-04 1.663256e-04
[21,] 0.9999654972 6.900556e-05 3.450278e-05
[22,] 0.9999424050 1.151900e-04 5.759500e-05
[23,] 0.9999113499 1.773003e-04 8.865014e-05
[24,] 0.9998557717 2.884565e-04 1.442283e-04
[25,] 0.9997840691 4.318618e-04 2.159309e-04
[26,] 0.9996593225 6.813550e-04 3.406775e-04
[27,] 0.9994752381 1.049524e-03 5.247619e-04
[28,] 0.9993407967 1.318407e-03 6.592033e-04
[29,] 0.9998811193 2.377614e-04 1.188807e-04
[30,] 0.9998135469 3.729061e-04 1.864531e-04
[31,] 0.9997149815 5.700369e-04 2.850185e-04
[32,] 0.9997382272 5.235455e-04 2.617728e-04
[33,] 0.9998980676 2.038647e-04 1.019324e-04
[34,] 0.9999256515 1.486969e-04 7.434846e-05
[35,] 0.9998827888 2.344223e-04 1.172112e-04
[36,] 0.9998627655 2.744689e-04 1.372345e-04
[37,] 0.9999566661 8.666783e-05 4.333392e-05
[38,] 0.9999316395 1.367210e-04 6.836052e-05
[39,] 0.9999005820 1.988359e-04 9.941797e-05
[40,] 0.9998547225 2.905550e-04 1.452775e-04
[41,] 0.9997989024 4.021951e-04 2.010976e-04
[42,] 0.9999914441 1.711178e-05 8.555888e-06
[43,] 0.9999883913 2.321731e-05 1.160866e-05
[44,] 0.9999932582 1.348370e-05 6.741849e-06
[45,] 0.9999893302 2.133961e-05 1.066980e-05
[46,] 0.9999875312 2.493770e-05 1.246885e-05
[47,] 0.9999820643 3.587138e-05 1.793569e-05
[48,] 0.9999744618 5.107645e-05 2.553823e-05
[49,] 0.9999863146 2.737085e-05 1.368543e-05
[50,] 0.9999827763 3.444730e-05 1.722365e-05
[51,] 0.9999831453 3.370946e-05 1.685473e-05
[52,] 0.9999872156 2.556885e-05 1.278442e-05
[53,] 0.9999825668 3.486631e-05 1.743315e-05
[54,] 0.9999750134 4.997326e-05 2.498663e-05
[55,] 0.9999649962 7.000758e-05 3.500379e-05
[56,] 0.9999489010 1.021980e-04 5.109898e-05
[57,] 0.9999344229 1.311541e-04 6.557707e-05
[58,] 0.9999051788 1.896425e-04 9.482125e-05
[59,] 0.9998960099 2.079801e-04 1.039901e-04
[60,] 0.9999914980 1.700394e-05 8.501971e-06
[61,] 0.9999942101 1.157981e-05 5.789906e-06
[62,] 0.9999947368 1.052639e-05 5.263193e-06
[63,] 0.9999936067 1.278655e-05 6.393274e-06
[64,] 0.9999901708 1.965832e-05 9.829160e-06
[65,] 0.9999860418 2.791643e-05 1.395821e-05
[66,] 0.9999956017 8.796503e-06 4.398251e-06
[67,] 0.9999946622 1.067569e-05 5.337847e-06
[68,] 0.9999920454 1.590913e-05 7.954566e-06
[69,] 0.9999881003 2.379934e-05 1.189967e-05
[70,] 0.9999853595 2.928093e-05 1.464046e-05
[71,] 0.9999782737 4.345266e-05 2.172633e-05
[72,] 0.9999726524 5.469516e-05 2.734758e-05
[73,] 0.9999804717 3.905650e-05 1.952825e-05
[74,] 0.9999999832 3.355196e-08 1.677598e-08
[75,] 0.9999999764 4.724104e-08 2.362052e-08
[76,] 0.9999999707 5.856934e-08 2.928467e-08
[77,] 0.9999999515 9.690634e-08 4.845317e-08
[78,] 0.9999999342 1.316577e-07 6.582887e-08
[79,] 0.9999999572 8.560098e-08 4.280049e-08
[80,] 0.9999999337 1.325564e-07 6.627821e-08
[81,] 0.9999999002 1.995475e-07 9.977375e-08
[82,] 0.9999998629 2.742485e-07 1.371243e-07
[83,] 0.9999998548 2.903526e-07 1.451763e-07
[84,] 0.9999997801 4.397174e-07 2.198587e-07
[85,] 0.9999997266 5.468352e-07 2.734176e-07
[86,] 0.9999995714 8.572548e-07 4.286274e-07
[87,] 0.9999993960 1.207977e-06 6.039887e-07
[88,] 0.9999993193 1.361308e-06 6.806542e-07
[89,] 0.9999990065 1.987055e-06 9.935274e-07
[90,] 0.9999995426 9.148308e-07 4.574154e-07
[91,] 0.9999993172 1.365579e-06 6.827896e-07
[92,] 0.9999990045 1.991071e-06 9.955353e-07
[93,] 0.9999984605 3.079089e-06 1.539544e-06
[94,] 0.9999988876 2.224704e-06 1.112352e-06
[95,] 0.9999998459 3.082911e-07 1.541455e-07
[96,] 0.9999998016 3.967188e-07 1.983594e-07
[97,] 0.9999996911 6.178630e-07 3.089315e-07
[98,] 0.9999997017 5.965840e-07 2.982920e-07
[99,] 0.9999995624 8.752973e-07 4.376486e-07
[100,] 0.9999996049 7.901414e-07 3.950707e-07
[101,] 0.9999997588 4.824017e-07 2.412008e-07
[102,] 0.9999996233 7.534525e-07 3.767263e-07
[103,] 0.9999997680 4.639905e-07 2.319953e-07
[104,] 0.9999996886 6.228048e-07 3.114024e-07
[105,] 0.9999996220 7.559313e-07 3.779656e-07
[106,] 0.9999994764 1.047209e-06 5.236044e-07
[107,] 0.9999992341 1.531724e-06 7.658618e-07
[108,] 0.9999992807 1.438669e-06 7.193343e-07
[109,] 0.9999991582 1.683677e-06 8.418385e-07
[110,] 0.9999987846 2.430775e-06 1.215388e-06
[111,] 0.9999987126 2.574894e-06 1.287447e-06
[112,] 0.9999980559 3.888261e-06 1.944130e-06
[113,] 0.9999984367 3.126595e-06 1.563298e-06
[114,] 0.9999981478 3.704356e-06 1.852178e-06
[115,] 0.9999980464 3.907220e-06 1.953610e-06
[116,] 0.9999988593 2.281395e-06 1.140697e-06
[117,] 0.9999983613 3.277359e-06 1.638680e-06
[118,] 0.9999975040 4.992020e-06 2.496010e-06
[119,] 0.9999962610 7.477977e-06 3.738989e-06
[120,] 0.9999949065 1.018707e-05 5.093537e-06
[121,] 0.9999924474 1.510527e-05 7.552636e-06
[122,] 0.9999915200 1.696002e-05 8.480011e-06
[123,] 0.9999977829 4.434209e-06 2.217104e-06
[124,] 0.9999972462 5.507567e-06 2.753783e-06
[125,] 0.9999960592 7.881612e-06 3.940806e-06
[126,] 0.9999941033 1.179336e-05 5.896679e-06
[127,] 0.9999932641 1.347171e-05 6.735853e-06
[128,] 0.9999904910 1.901802e-05 9.509012e-06
[129,] 0.9999877416 2.451683e-05 1.225841e-05
[130,] 0.9999897506 2.049871e-05 1.024936e-05
[131,] 0.9999866730 2.665402e-05 1.332701e-05
[132,] 0.9999805263 3.894739e-05 1.947370e-05
[133,] 0.9999722616 5.547680e-05 2.773840e-05
[134,] 0.9999598777 8.024457e-05 4.012229e-05
[135,] 0.9999426195 1.147611e-04 5.738053e-05
[136,] 0.9999179736 1.640529e-04 8.202645e-05
[137,] 0.9998857361 2.285279e-04 1.142639e-04
[138,] 0.9998433087 3.133825e-04 1.566913e-04
[139,] 0.9997800196 4.399608e-04 2.199804e-04
[140,] 0.9997134652 5.730696e-04 2.865348e-04
[141,] 0.9996467242 7.065516e-04 3.532758e-04
[142,] 0.9995523478 8.953044e-04 4.476522e-04
[143,] 0.9994740200 1.051960e-03 5.259800e-04
[144,] 0.9993978297 1.204341e-03 6.021703e-04
[145,] 0.9991829945 1.634011e-03 8.170055e-04
[146,] 0.9988920347 2.215931e-03 1.107965e-03
[147,] 0.9985386679 2.922664e-03 1.461332e-03
[148,] 0.9981219991 3.756002e-03 1.878001e-03
[149,] 0.9975936478 4.812704e-03 2.406352e-03
[150,] 0.9969677823 6.064435e-03 3.032218e-03
[151,] 0.9960321924 7.935615e-03 3.967808e-03
[152,] 0.9960793894 7.841221e-03 3.920611e-03
[153,] 0.9952349002 9.530200e-03 4.765100e-03
[154,] 0.9941811759 1.163765e-02 5.818824e-03
[155,] 0.9926564833 1.468703e-02 7.343517e-03
[156,] 0.9907854723 1.842906e-02 9.214528e-03
[157,] 0.9967516994 6.496601e-03 3.248301e-03
[158,] 0.9971208456 5.758309e-03 2.879154e-03
[159,] 0.9972258084 5.548383e-03 2.774192e-03
[160,] 0.9967396921 6.520616e-03 3.260308e-03
[161,] 0.9966812107 6.637579e-03 3.318789e-03
[162,] 0.9957097569 8.580486e-03 4.290243e-03
[163,] 0.9952611173 9.477765e-03 4.738883e-03
[164,] 0.9939093133 1.218137e-02 6.090687e-03
[165,] 0.9921733563 1.565329e-02 7.826644e-03
[166,] 0.9922300546 1.553989e-02 7.769945e-03
[167,] 0.9911958967 1.760821e-02 8.804103e-03
[168,] 0.9896725039 2.065499e-02 1.032750e-02
[169,] 0.9913900483 1.721990e-02 8.609952e-03
[170,] 0.9894844951 2.103101e-02 1.051550e-02
[171,] 0.9867336328 2.653273e-02 1.326637e-02
[172,] 0.9916480471 1.670391e-02 8.351953e-03
[173,] 0.9895157145 2.096857e-02 1.048429e-02
[174,] 0.9868202392 2.635952e-02 1.317976e-02
[175,] 0.9834993988 3.300120e-02 1.650060e-02
[176,] 0.9801130649 3.977387e-02 1.988694e-02
[177,] 0.9776158501 4.476830e-02 2.238415e-02
[178,] 0.9781324482 4.373510e-02 2.186755e-02
[179,] 0.9781805137 4.363897e-02 2.181949e-02
[180,] 0.9737257696 5.254846e-02 2.627423e-02
[181,] 0.9674846969 6.503061e-02 3.251530e-02
[182,] 0.9653618632 6.927627e-02 3.463814e-02
[183,] 0.9580268491 8.394630e-02 4.197315e-02
[184,] 0.9543246521 9.135070e-02 4.567535e-02
[185,] 0.9443817362 1.112365e-01 5.561826e-02
[186,] 0.9826607376 3.467852e-02 1.733926e-02
[187,] 0.9790437785 4.191244e-02 2.095622e-02
[188,] 0.9738950795 5.220984e-02 2.610492e-02
[189,] 0.9708751021 5.824980e-02 2.912490e-02
[190,] 0.9643394355 7.132113e-02 3.566056e-02
[191,] 0.9564218435 8.715631e-02 4.357816e-02
[192,] 0.9492593409 1.014813e-01 5.074066e-02
[193,] 0.9381237537 1.237525e-01 6.187625e-02
[194,] 0.9249622690 1.500755e-01 7.503773e-02
[195,] 0.9103050586 1.793899e-01 8.969494e-02
[196,] 0.9037769070 1.924462e-01 9.622309e-02
[197,] 0.9036541423 1.926917e-01 9.634586e-02
[198,] 0.8863937898 2.272124e-01 1.136062e-01
[199,] 0.8729522429 2.540955e-01 1.270478e-01
[200,] 0.8773623624 2.452753e-01 1.226376e-01
[201,] 0.8596958329 2.806083e-01 1.403042e-01
[202,] 0.8347974547 3.304051e-01 1.652025e-01
[203,] 0.8119658477 3.760683e-01 1.880342e-01
[204,] 0.7813117782 4.373764e-01 2.186882e-01
[205,] 0.7480168543 5.039663e-01 2.519831e-01
[206,] 0.7259791840 5.480416e-01 2.740208e-01
[207,] 0.7143707404 5.712585e-01 2.856293e-01
[208,] 0.6921314854 6.157370e-01 3.078685e-01
[209,] 0.6516448902 6.967102e-01 3.483551e-01
[210,] 0.6164018402 7.671963e-01 3.835982e-01
[211,] 0.5749960827 8.500078e-01 4.250039e-01
[212,] 0.5296918497 9.406163e-01 4.703082e-01
[213,] 0.4920579457 9.841159e-01 5.079421e-01
[214,] 0.4506929888 9.013860e-01 5.493070e-01
[215,] 0.4119813654 8.239627e-01 5.880186e-01
[216,] 0.3778544167 7.557088e-01 6.221456e-01
[217,] 0.3398173019 6.796346e-01 6.601827e-01
[218,] 0.4522498874 9.044998e-01 5.477501e-01
[219,] 0.4092747157 8.185494e-01 5.907253e-01
[220,] 0.3651893705 7.303787e-01 6.348106e-01
[221,] 0.3318233944 6.636468e-01 6.681766e-01
[222,] 0.2913267267 5.826535e-01 7.086733e-01
[223,] 0.2607890117 5.215780e-01 7.392110e-01
[224,] 0.2224175147 4.448350e-01 7.775825e-01
[225,] 0.1875181229 3.750362e-01 8.124819e-01
[226,] 0.1647209549 3.294419e-01 8.352790e-01
[227,] 0.1366232876 2.732466e-01 8.633767e-01
[228,] 0.1106641107 2.213282e-01 8.893359e-01
[229,] 0.0975576580 1.951153e-01 9.024423e-01
[230,] 0.0780629858 1.561260e-01 9.219370e-01
[231,] 0.0754905982 1.509812e-01 9.245094e-01
[232,] 0.0592338898 1.184678e-01 9.407661e-01
[233,] 0.0447443706 8.948874e-02 9.552556e-01
[234,] 0.0397446816 7.948936e-02 9.602553e-01
[235,] 0.0309593319 6.191866e-02 9.690407e-01
[236,] 0.0222172381 4.443448e-02 9.777828e-01
[237,] 0.0192778940 3.855579e-02 9.807221e-01
[238,] 0.0136145969 2.722919e-02 9.863854e-01
[239,] 0.0102385256 2.047705e-02 9.897615e-01
[240,] 0.0071560188 1.431204e-02 9.928440e-01
[241,] 0.0076634418 1.532688e-02 9.923366e-01
[242,] 0.0068696875 1.373937e-02 9.931303e-01
[243,] 0.0048917028 9.783406e-03 9.951083e-01
[244,] 0.0033371168 6.674234e-03 9.966629e-01
[245,] 0.0027023363 5.404673e-03 9.972977e-01
[246,] 0.0031424394 6.284879e-03 9.968576e-01
[247,] 0.0018130628 3.626126e-03 9.981869e-01
[248,] 0.0080814564 1.616291e-02 9.919185e-01
[249,] 0.0048694016 9.738803e-03 9.951306e-01
[250,] 0.0031141070 6.228214e-03 9.968859e-01
[251,] 0.0018993527 3.798705e-03 9.981006e-01
[252,] 0.0012287489 2.457498e-03 9.987713e-01
[253,] 0.0008046153 1.609231e-03 9.991954e-01
[254,] 0.0006834352 1.366870e-03 9.993166e-01
[255,] 0.0002638110 5.276220e-04 9.997362e-01
[256,] 0.0001363975 2.727951e-04 9.998636e-01
> postscript(file="/var/wessaorg/rcomp/tmp/1756t1353439957.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/2sa1l1353439957.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/3e22s1353439957.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/4w5sm1353439957.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/58gtf1353439957.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
26428.27684 -4696.16715 5479.91615 -57105.19577 15970.70964
6 7 8 9 10
-14717.14856 74545.78641 -9781.94403 -15235.39230 10506.35337
11 12 13 14 15
23414.47216 -600.48554 -9343.40404 18015.35866 17055.98749
16 17 18 19 20
-65670.85000 26354.48554 15630.41454 -11612.12797 -11262.74909
21 22 23 24 25
-3878.35131 110025.39208 23312.54944 -28133.91430 -85167.42773
26 27 28 29 30
-45996.62294 -14935.37275 32510.53014 -4063.54144 20087.17737
31 32 33 34 35
16510.28521 -25132.49030 48553.70982 -4089.13216 26458.45180
36 37 38 39 40
25491.84421 63888.39179 3621.97758 18736.01233 2763.30103
41 42 43 44 45
15843.81691 -1728.39900 6331.29350 -12776.74148 -81018.32699
46 47 48 49 50
-1911.33484 -14714.70643 -32551.07734 -50255.62391 -40612.69433
51 52 53 54 55
1282.23261 -15256.82238 64474.80019 -3515.30902 -6053.36940
56 57 58 59 60
-13218.07509 4834.59514 -77149.82608 14488.53542 39528.44019
61 62 63 64 65
-290.60990 -27070.59051 -19629.53912 7303.59050 42884.17448
66 67 68 69 70
18366.86946 -32225.94758 -44121.30039 -21636.68289 -17082.13737
71 72 73 74 75
-12216.46130 -11208.10244 9315.30316 2039.13048 -32348.33020
76 77 78 79 80
75368.83380 38477.42678 29782.91495 -25964.30067 -6214.23589
81 82 83 84 85
-14940.25270 54486.08972 19601.68606 -9510.53209 -7321.42608
86 87 88 89 90
-23463.50915 876.87506 17045.39655 36521.64182 -122510.65607
91 92 93 94 95
14355.88520 -26692.91851 -3515.30740 -16618.09332 36782.54510
96 97 98 99 100
6510.25173 7601.92196 13502.88799 -29117.66968 -2651.09629
101 102 103 104 105
19234.75368 13.41169 -9140.04748 -26044.19559 -9037.71206
106 107 108 109 110
-47913.81304 -757.29441 9672.38077 -2424.24157 -39143.32909
111 112 113 114 115
66902.08143 18680.22956 -4097.42888 -28140.64838 12033.29892
116 117 118 119 120
-28574.54489 44022.92625 -2496.49976 40800.45998 -17651.42793
121 122 123 124 125
-14217.86223 13673.73924 5234.09540 -27467.92748 -11350.71213
126 127 128 129 130
2924.58018 32353.15929 6958.78174 40342.44884 23721.51387
131 132 133 134 135
23287.53322 34468.56626 -17000.22716 -2585.38361 -5176.26558
136 137 138 139 140
-3268.04826 7118.50995 -28993.31106 30130.92328 -26982.83925
141 142 143 144 145
2576.89831 4619.96001 -22885.31622 11138.78397 13422.85690
146 147 148 149 150
-54066.15741 -28855.78003 -7513.33032 5277.87325 4830.94525
151 152 153 154 155
7503.30867 1986.40945 10179.31201 13481.33322 2537.62105
156 157 158 159 160
22221.44904 19394.29135 -17099.74534 -21362.19657 24494.89045
161 162 163 164 165
4407.90596 -1662.15495 7987.78640 -12983.91194 -7727.78430
166 167 168 169 170
8517.00366 1313.56289 26475.59749 14705.11111 -15354.48315
171 172 173 174 175
-7970.67398 -8441.70497 61841.16131 32104.87629 28104.24463
176 177 178 179 180
14012.74159 21771.69084 -3315.50644 16850.50093 -5756.29396
181 182 183 184 185
-4648.81467 -26866.38653 15837.88003 13151.24707 30723.36189
186 187 188 189 190
-6978.51630 -2535.48770 -50724.69543 3473.89417 -1913.91984
191 192 193 194 195
2292.25053 5708.30303 -22487.18182 -33291.70202 -34595.67091
196 197 198 199 200
-16216.77126 -10917.44829 -27620.34008 -12201.56707 -24295.64258
201 202 203 204 205
-8559.95448 62039.78441 -16688.80113 -10970.95796 -22621.84059
206 207 208 209 210
4807.90699 -12550.70419 -17612.17048 -1226.41748 -3067.42889
211 212 213 214 215
173.16567 -22217.70119 24475.95957 6693.55291 -17334.88722
216 217 218 219 220
24258.18266 9732.89461 -5299.18834 5510.52293 -3960.63805
221 222 223 224 225
-6429.30438 -17716.93207 16464.29566 10137.95114 -3426.28300
226 227 228 229 230
-20546.41344 2629.19094 -4540.20448 -11685.94112 -11949.21833
231 232 233 234 235
9704.31931 -15060.59538 -8360.85598 52168.98176 -12802.59537
236 237 238 239 240
-11713.86750 9946.57567 -10494.76038 12141.34411 2306.81170
241 242 243 244 245
1953.02801 12077.25356 -5182.81996 -498.29222 17066.18856
246 247 248 249 250
2168.49389 18418.24224 8658.23019 1254.60436 8259.82371
251 252 253 254 255
-26872.90055 -5238.86963 -17681.55002 1040.23893 6840.93166
256 257 258 259 260
-12429.03136 -41785.86753 28610.93029 -24452.13181 6586.63497
261 262 263 264 265
17516.29541 25907.47054 -17578.46503 61377.18219 -27032.78723
266 267 268 269 270
3101.30768 -36255.92258 -20967.62469 25026.28378 22140.69032
271 272 273 274 275
-13206.37347 -2047.57041 -48401.21391 -33644.92214 -20613.90708
276 277 278 279 280
-8901.94575 9553.82222 9435.93190 45834.84461 -17711.61453
281 282 283 284 285
-35403.94824 13696.50624 19004.81516 21858.69225 102115.43370
286 287 288 289
-4802.30560 9815.29041 -1386.54546 -29776.03775
> postscript(file="/var/wessaorg/rcomp/tmp/6z4eu1353439957.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 26428.27684 NA
1 -4696.16715 26428.27684
2 5479.91615 -4696.16715
3 -57105.19577 5479.91615
4 15970.70964 -57105.19577
5 -14717.14856 15970.70964
6 74545.78641 -14717.14856
7 -9781.94403 74545.78641
8 -15235.39230 -9781.94403
9 10506.35337 -15235.39230
10 23414.47216 10506.35337
11 -600.48554 23414.47216
12 -9343.40404 -600.48554
13 18015.35866 -9343.40404
14 17055.98749 18015.35866
15 -65670.85000 17055.98749
16 26354.48554 -65670.85000
17 15630.41454 26354.48554
18 -11612.12797 15630.41454
19 -11262.74909 -11612.12797
20 -3878.35131 -11262.74909
21 110025.39208 -3878.35131
22 23312.54944 110025.39208
23 -28133.91430 23312.54944
24 -85167.42773 -28133.91430
25 -45996.62294 -85167.42773
26 -14935.37275 -45996.62294
27 32510.53014 -14935.37275
28 -4063.54144 32510.53014
29 20087.17737 -4063.54144
30 16510.28521 20087.17737
31 -25132.49030 16510.28521
32 48553.70982 -25132.49030
33 -4089.13216 48553.70982
34 26458.45180 -4089.13216
35 25491.84421 26458.45180
36 63888.39179 25491.84421
37 3621.97758 63888.39179
38 18736.01233 3621.97758
39 2763.30103 18736.01233
40 15843.81691 2763.30103
41 -1728.39900 15843.81691
42 6331.29350 -1728.39900
43 -12776.74148 6331.29350
44 -81018.32699 -12776.74148
45 -1911.33484 -81018.32699
46 -14714.70643 -1911.33484
47 -32551.07734 -14714.70643
48 -50255.62391 -32551.07734
49 -40612.69433 -50255.62391
50 1282.23261 -40612.69433
51 -15256.82238 1282.23261
52 64474.80019 -15256.82238
53 -3515.30902 64474.80019
54 -6053.36940 -3515.30902
55 -13218.07509 -6053.36940
56 4834.59514 -13218.07509
57 -77149.82608 4834.59514
58 14488.53542 -77149.82608
59 39528.44019 14488.53542
60 -290.60990 39528.44019
61 -27070.59051 -290.60990
62 -19629.53912 -27070.59051
63 7303.59050 -19629.53912
64 42884.17448 7303.59050
65 18366.86946 42884.17448
66 -32225.94758 18366.86946
67 -44121.30039 -32225.94758
68 -21636.68289 -44121.30039
69 -17082.13737 -21636.68289
70 -12216.46130 -17082.13737
71 -11208.10244 -12216.46130
72 9315.30316 -11208.10244
73 2039.13048 9315.30316
74 -32348.33020 2039.13048
75 75368.83380 -32348.33020
76 38477.42678 75368.83380
77 29782.91495 38477.42678
78 -25964.30067 29782.91495
79 -6214.23589 -25964.30067
80 -14940.25270 -6214.23589
81 54486.08972 -14940.25270
82 19601.68606 54486.08972
83 -9510.53209 19601.68606
84 -7321.42608 -9510.53209
85 -23463.50915 -7321.42608
86 876.87506 -23463.50915
87 17045.39655 876.87506
88 36521.64182 17045.39655
89 -122510.65607 36521.64182
90 14355.88520 -122510.65607
91 -26692.91851 14355.88520
92 -3515.30740 -26692.91851
93 -16618.09332 -3515.30740
94 36782.54510 -16618.09332
95 6510.25173 36782.54510
96 7601.92196 6510.25173
97 13502.88799 7601.92196
98 -29117.66968 13502.88799
99 -2651.09629 -29117.66968
100 19234.75368 -2651.09629
101 13.41169 19234.75368
102 -9140.04748 13.41169
103 -26044.19559 -9140.04748
104 -9037.71206 -26044.19559
105 -47913.81304 -9037.71206
106 -757.29441 -47913.81304
107 9672.38077 -757.29441
108 -2424.24157 9672.38077
109 -39143.32909 -2424.24157
110 66902.08143 -39143.32909
111 18680.22956 66902.08143
112 -4097.42888 18680.22956
113 -28140.64838 -4097.42888
114 12033.29892 -28140.64838
115 -28574.54489 12033.29892
116 44022.92625 -28574.54489
117 -2496.49976 44022.92625
118 40800.45998 -2496.49976
119 -17651.42793 40800.45998
120 -14217.86223 -17651.42793
121 13673.73924 -14217.86223
122 5234.09540 13673.73924
123 -27467.92748 5234.09540
124 -11350.71213 -27467.92748
125 2924.58018 -11350.71213
126 32353.15929 2924.58018
127 6958.78174 32353.15929
128 40342.44884 6958.78174
129 23721.51387 40342.44884
130 23287.53322 23721.51387
131 34468.56626 23287.53322
132 -17000.22716 34468.56626
133 -2585.38361 -17000.22716
134 -5176.26558 -2585.38361
135 -3268.04826 -5176.26558
136 7118.50995 -3268.04826
137 -28993.31106 7118.50995
138 30130.92328 -28993.31106
139 -26982.83925 30130.92328
140 2576.89831 -26982.83925
141 4619.96001 2576.89831
142 -22885.31622 4619.96001
143 11138.78397 -22885.31622
144 13422.85690 11138.78397
145 -54066.15741 13422.85690
146 -28855.78003 -54066.15741
147 -7513.33032 -28855.78003
148 5277.87325 -7513.33032
149 4830.94525 5277.87325
150 7503.30867 4830.94525
151 1986.40945 7503.30867
152 10179.31201 1986.40945
153 13481.33322 10179.31201
154 2537.62105 13481.33322
155 22221.44904 2537.62105
156 19394.29135 22221.44904
157 -17099.74534 19394.29135
158 -21362.19657 -17099.74534
159 24494.89045 -21362.19657
160 4407.90596 24494.89045
161 -1662.15495 4407.90596
162 7987.78640 -1662.15495
163 -12983.91194 7987.78640
164 -7727.78430 -12983.91194
165 8517.00366 -7727.78430
166 1313.56289 8517.00366
167 26475.59749 1313.56289
168 14705.11111 26475.59749
169 -15354.48315 14705.11111
170 -7970.67398 -15354.48315
171 -8441.70497 -7970.67398
172 61841.16131 -8441.70497
173 32104.87629 61841.16131
174 28104.24463 32104.87629
175 14012.74159 28104.24463
176 21771.69084 14012.74159
177 -3315.50644 21771.69084
178 16850.50093 -3315.50644
179 -5756.29396 16850.50093
180 -4648.81467 -5756.29396
181 -26866.38653 -4648.81467
182 15837.88003 -26866.38653
183 13151.24707 15837.88003
184 30723.36189 13151.24707
185 -6978.51630 30723.36189
186 -2535.48770 -6978.51630
187 -50724.69543 -2535.48770
188 3473.89417 -50724.69543
189 -1913.91984 3473.89417
190 2292.25053 -1913.91984
191 5708.30303 2292.25053
192 -22487.18182 5708.30303
193 -33291.70202 -22487.18182
194 -34595.67091 -33291.70202
195 -16216.77126 -34595.67091
196 -10917.44829 -16216.77126
197 -27620.34008 -10917.44829
198 -12201.56707 -27620.34008
199 -24295.64258 -12201.56707
200 -8559.95448 -24295.64258
201 62039.78441 -8559.95448
202 -16688.80113 62039.78441
203 -10970.95796 -16688.80113
204 -22621.84059 -10970.95796
205 4807.90699 -22621.84059
206 -12550.70419 4807.90699
207 -17612.17048 -12550.70419
208 -1226.41748 -17612.17048
209 -3067.42889 -1226.41748
210 173.16567 -3067.42889
211 -22217.70119 173.16567
212 24475.95957 -22217.70119
213 6693.55291 24475.95957
214 -17334.88722 6693.55291
215 24258.18266 -17334.88722
216 9732.89461 24258.18266
217 -5299.18834 9732.89461
218 5510.52293 -5299.18834
219 -3960.63805 5510.52293
220 -6429.30438 -3960.63805
221 -17716.93207 -6429.30438
222 16464.29566 -17716.93207
223 10137.95114 16464.29566
224 -3426.28300 10137.95114
225 -20546.41344 -3426.28300
226 2629.19094 -20546.41344
227 -4540.20448 2629.19094
228 -11685.94112 -4540.20448
229 -11949.21833 -11685.94112
230 9704.31931 -11949.21833
231 -15060.59538 9704.31931
232 -8360.85598 -15060.59538
233 52168.98176 -8360.85598
234 -12802.59537 52168.98176
235 -11713.86750 -12802.59537
236 9946.57567 -11713.86750
237 -10494.76038 9946.57567
238 12141.34411 -10494.76038
239 2306.81170 12141.34411
240 1953.02801 2306.81170
241 12077.25356 1953.02801
242 -5182.81996 12077.25356
243 -498.29222 -5182.81996
244 17066.18856 -498.29222
245 2168.49389 17066.18856
246 18418.24224 2168.49389
247 8658.23019 18418.24224
248 1254.60436 8658.23019
249 8259.82371 1254.60436
250 -26872.90055 8259.82371
251 -5238.86963 -26872.90055
252 -17681.55002 -5238.86963
253 1040.23893 -17681.55002
254 6840.93166 1040.23893
255 -12429.03136 6840.93166
256 -41785.86753 -12429.03136
257 28610.93029 -41785.86753
258 -24452.13181 28610.93029
259 6586.63497 -24452.13181
260 17516.29541 6586.63497
261 25907.47054 17516.29541
262 -17578.46503 25907.47054
263 61377.18219 -17578.46503
264 -27032.78723 61377.18219
265 3101.30768 -27032.78723
266 -36255.92258 3101.30768
267 -20967.62469 -36255.92258
268 25026.28378 -20967.62469
269 22140.69032 25026.28378
270 -13206.37347 22140.69032
271 -2047.57041 -13206.37347
272 -48401.21391 -2047.57041
273 -33644.92214 -48401.21391
274 -20613.90708 -33644.92214
275 -8901.94575 -20613.90708
276 9553.82222 -8901.94575
277 9435.93190 9553.82222
278 45834.84461 9435.93190
279 -17711.61453 45834.84461
280 -35403.94824 -17711.61453
281 13696.50624 -35403.94824
282 19004.81516 13696.50624
283 21858.69225 19004.81516
284 102115.43370 21858.69225
285 -4802.30560 102115.43370
286 9815.29041 -4802.30560
287 -1386.54546 9815.29041
288 -29776.03775 -1386.54546
289 NA -29776.03775
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] -4696.16715 26428.27684
[2,] 5479.91615 -4696.16715
[3,] -57105.19577 5479.91615
[4,] 15970.70964 -57105.19577
[5,] -14717.14856 15970.70964
[6,] 74545.78641 -14717.14856
[7,] -9781.94403 74545.78641
[8,] -15235.39230 -9781.94403
[9,] 10506.35337 -15235.39230
[10,] 23414.47216 10506.35337
[11,] -600.48554 23414.47216
[12,] -9343.40404 -600.48554
[13,] 18015.35866 -9343.40404
[14,] 17055.98749 18015.35866
[15,] -65670.85000 17055.98749
[16,] 26354.48554 -65670.85000
[17,] 15630.41454 26354.48554
[18,] -11612.12797 15630.41454
[19,] -11262.74909 -11612.12797
[20,] -3878.35131 -11262.74909
[21,] 110025.39208 -3878.35131
[22,] 23312.54944 110025.39208
[23,] -28133.91430 23312.54944
[24,] -85167.42773 -28133.91430
[25,] -45996.62294 -85167.42773
[26,] -14935.37275 -45996.62294
[27,] 32510.53014 -14935.37275
[28,] -4063.54144 32510.53014
[29,] 20087.17737 -4063.54144
[30,] 16510.28521 20087.17737
[31,] -25132.49030 16510.28521
[32,] 48553.70982 -25132.49030
[33,] -4089.13216 48553.70982
[34,] 26458.45180 -4089.13216
[35,] 25491.84421 26458.45180
[36,] 63888.39179 25491.84421
[37,] 3621.97758 63888.39179
[38,] 18736.01233 3621.97758
[39,] 2763.30103 18736.01233
[40,] 15843.81691 2763.30103
[41,] -1728.39900 15843.81691
[42,] 6331.29350 -1728.39900
[43,] -12776.74148 6331.29350
[44,] -81018.32699 -12776.74148
[45,] -1911.33484 -81018.32699
[46,] -14714.70643 -1911.33484
[47,] -32551.07734 -14714.70643
[48,] -50255.62391 -32551.07734
[49,] -40612.69433 -50255.62391
[50,] 1282.23261 -40612.69433
[51,] -15256.82238 1282.23261
[52,] 64474.80019 -15256.82238
[53,] -3515.30902 64474.80019
[54,] -6053.36940 -3515.30902
[55,] -13218.07509 -6053.36940
[56,] 4834.59514 -13218.07509
[57,] -77149.82608 4834.59514
[58,] 14488.53542 -77149.82608
[59,] 39528.44019 14488.53542
[60,] -290.60990 39528.44019
[61,] -27070.59051 -290.60990
[62,] -19629.53912 -27070.59051
[63,] 7303.59050 -19629.53912
[64,] 42884.17448 7303.59050
[65,] 18366.86946 42884.17448
[66,] -32225.94758 18366.86946
[67,] -44121.30039 -32225.94758
[68,] -21636.68289 -44121.30039
[69,] -17082.13737 -21636.68289
[70,] -12216.46130 -17082.13737
[71,] -11208.10244 -12216.46130
[72,] 9315.30316 -11208.10244
[73,] 2039.13048 9315.30316
[74,] -32348.33020 2039.13048
[75,] 75368.83380 -32348.33020
[76,] 38477.42678 75368.83380
[77,] 29782.91495 38477.42678
[78,] -25964.30067 29782.91495
[79,] -6214.23589 -25964.30067
[80,] -14940.25270 -6214.23589
[81,] 54486.08972 -14940.25270
[82,] 19601.68606 54486.08972
[83,] -9510.53209 19601.68606
[84,] -7321.42608 -9510.53209
[85,] -23463.50915 -7321.42608
[86,] 876.87506 -23463.50915
[87,] 17045.39655 876.87506
[88,] 36521.64182 17045.39655
[89,] -122510.65607 36521.64182
[90,] 14355.88520 -122510.65607
[91,] -26692.91851 14355.88520
[92,] -3515.30740 -26692.91851
[93,] -16618.09332 -3515.30740
[94,] 36782.54510 -16618.09332
[95,] 6510.25173 36782.54510
[96,] 7601.92196 6510.25173
[97,] 13502.88799 7601.92196
[98,] -29117.66968 13502.88799
[99,] -2651.09629 -29117.66968
[100,] 19234.75368 -2651.09629
[101,] 13.41169 19234.75368
[102,] -9140.04748 13.41169
[103,] -26044.19559 -9140.04748
[104,] -9037.71206 -26044.19559
[105,] -47913.81304 -9037.71206
[106,] -757.29441 -47913.81304
[107,] 9672.38077 -757.29441
[108,] -2424.24157 9672.38077
[109,] -39143.32909 -2424.24157
[110,] 66902.08143 -39143.32909
[111,] 18680.22956 66902.08143
[112,] -4097.42888 18680.22956
[113,] -28140.64838 -4097.42888
[114,] 12033.29892 -28140.64838
[115,] -28574.54489 12033.29892
[116,] 44022.92625 -28574.54489
[117,] -2496.49976 44022.92625
[118,] 40800.45998 -2496.49976
[119,] -17651.42793 40800.45998
[120,] -14217.86223 -17651.42793
[121,] 13673.73924 -14217.86223
[122,] 5234.09540 13673.73924
[123,] -27467.92748 5234.09540
[124,] -11350.71213 -27467.92748
[125,] 2924.58018 -11350.71213
[126,] 32353.15929 2924.58018
[127,] 6958.78174 32353.15929
[128,] 40342.44884 6958.78174
[129,] 23721.51387 40342.44884
[130,] 23287.53322 23721.51387
[131,] 34468.56626 23287.53322
[132,] -17000.22716 34468.56626
[133,] -2585.38361 -17000.22716
[134,] -5176.26558 -2585.38361
[135,] -3268.04826 -5176.26558
[136,] 7118.50995 -3268.04826
[137,] -28993.31106 7118.50995
[138,] 30130.92328 -28993.31106
[139,] -26982.83925 30130.92328
[140,] 2576.89831 -26982.83925
[141,] 4619.96001 2576.89831
[142,] -22885.31622 4619.96001
[143,] 11138.78397 -22885.31622
[144,] 13422.85690 11138.78397
[145,] -54066.15741 13422.85690
[146,] -28855.78003 -54066.15741
[147,] -7513.33032 -28855.78003
[148,] 5277.87325 -7513.33032
[149,] 4830.94525 5277.87325
[150,] 7503.30867 4830.94525
[151,] 1986.40945 7503.30867
[152,] 10179.31201 1986.40945
[153,] 13481.33322 10179.31201
[154,] 2537.62105 13481.33322
[155,] 22221.44904 2537.62105
[156,] 19394.29135 22221.44904
[157,] -17099.74534 19394.29135
[158,] -21362.19657 -17099.74534
[159,] 24494.89045 -21362.19657
[160,] 4407.90596 24494.89045
[161,] -1662.15495 4407.90596
[162,] 7987.78640 -1662.15495
[163,] -12983.91194 7987.78640
[164,] -7727.78430 -12983.91194
[165,] 8517.00366 -7727.78430
[166,] 1313.56289 8517.00366
[167,] 26475.59749 1313.56289
[168,] 14705.11111 26475.59749
[169,] -15354.48315 14705.11111
[170,] -7970.67398 -15354.48315
[171,] -8441.70497 -7970.67398
[172,] 61841.16131 -8441.70497
[173,] 32104.87629 61841.16131
[174,] 28104.24463 32104.87629
[175,] 14012.74159 28104.24463
[176,] 21771.69084 14012.74159
[177,] -3315.50644 21771.69084
[178,] 16850.50093 -3315.50644
[179,] -5756.29396 16850.50093
[180,] -4648.81467 -5756.29396
[181,] -26866.38653 -4648.81467
[182,] 15837.88003 -26866.38653
[183,] 13151.24707 15837.88003
[184,] 30723.36189 13151.24707
[185,] -6978.51630 30723.36189
[186,] -2535.48770 -6978.51630
[187,] -50724.69543 -2535.48770
[188,] 3473.89417 -50724.69543
[189,] -1913.91984 3473.89417
[190,] 2292.25053 -1913.91984
[191,] 5708.30303 2292.25053
[192,] -22487.18182 5708.30303
[193,] -33291.70202 -22487.18182
[194,] -34595.67091 -33291.70202
[195,] -16216.77126 -34595.67091
[196,] -10917.44829 -16216.77126
[197,] -27620.34008 -10917.44829
[198,] -12201.56707 -27620.34008
[199,] -24295.64258 -12201.56707
[200,] -8559.95448 -24295.64258
[201,] 62039.78441 -8559.95448
[202,] -16688.80113 62039.78441
[203,] -10970.95796 -16688.80113
[204,] -22621.84059 -10970.95796
[205,] 4807.90699 -22621.84059
[206,] -12550.70419 4807.90699
[207,] -17612.17048 -12550.70419
[208,] -1226.41748 -17612.17048
[209,] -3067.42889 -1226.41748
[210,] 173.16567 -3067.42889
[211,] -22217.70119 173.16567
[212,] 24475.95957 -22217.70119
[213,] 6693.55291 24475.95957
[214,] -17334.88722 6693.55291
[215,] 24258.18266 -17334.88722
[216,] 9732.89461 24258.18266
[217,] -5299.18834 9732.89461
[218,] 5510.52293 -5299.18834
[219,] -3960.63805 5510.52293
[220,] -6429.30438 -3960.63805
[221,] -17716.93207 -6429.30438
[222,] 16464.29566 -17716.93207
[223,] 10137.95114 16464.29566
[224,] -3426.28300 10137.95114
[225,] -20546.41344 -3426.28300
[226,] 2629.19094 -20546.41344
[227,] -4540.20448 2629.19094
[228,] -11685.94112 -4540.20448
[229,] -11949.21833 -11685.94112
[230,] 9704.31931 -11949.21833
[231,] -15060.59538 9704.31931
[232,] -8360.85598 -15060.59538
[233,] 52168.98176 -8360.85598
[234,] -12802.59537 52168.98176
[235,] -11713.86750 -12802.59537
[236,] 9946.57567 -11713.86750
[237,] -10494.76038 9946.57567
[238,] 12141.34411 -10494.76038
[239,] 2306.81170 12141.34411
[240,] 1953.02801 2306.81170
[241,] 12077.25356 1953.02801
[242,] -5182.81996 12077.25356
[243,] -498.29222 -5182.81996
[244,] 17066.18856 -498.29222
[245,] 2168.49389 17066.18856
[246,] 18418.24224 2168.49389
[247,] 8658.23019 18418.24224
[248,] 1254.60436 8658.23019
[249,] 8259.82371 1254.60436
[250,] -26872.90055 8259.82371
[251,] -5238.86963 -26872.90055
[252,] -17681.55002 -5238.86963
[253,] 1040.23893 -17681.55002
[254,] 6840.93166 1040.23893
[255,] -12429.03136 6840.93166
[256,] -41785.86753 -12429.03136
[257,] 28610.93029 -41785.86753
[258,] -24452.13181 28610.93029
[259,] 6586.63497 -24452.13181
[260,] 17516.29541 6586.63497
[261,] 25907.47054 17516.29541
[262,] -17578.46503 25907.47054
[263,] 61377.18219 -17578.46503
[264,] -27032.78723 61377.18219
[265,] 3101.30768 -27032.78723
[266,] -36255.92258 3101.30768
[267,] -20967.62469 -36255.92258
[268,] 25026.28378 -20967.62469
[269,] 22140.69032 25026.28378
[270,] -13206.37347 22140.69032
[271,] -2047.57041 -13206.37347
[272,] -48401.21391 -2047.57041
[273,] -33644.92214 -48401.21391
[274,] -20613.90708 -33644.92214
[275,] -8901.94575 -20613.90708
[276,] 9553.82222 -8901.94575
[277,] 9435.93190 9553.82222
[278,] 45834.84461 9435.93190
[279,] -17711.61453 45834.84461
[280,] -35403.94824 -17711.61453
[281,] 13696.50624 -35403.94824
[282,] 19004.81516 13696.50624
[283,] 21858.69225 19004.81516
[284,] 102115.43370 21858.69225
[285,] -4802.30560 102115.43370
[286,] 9815.29041 -4802.30560
[287,] -1386.54546 9815.29041
[288,] -29776.03775 -1386.54546
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 -4696.16715 26428.27684
2 5479.91615 -4696.16715
3 -57105.19577 5479.91615
4 15970.70964 -57105.19577
5 -14717.14856 15970.70964
6 74545.78641 -14717.14856
7 -9781.94403 74545.78641
8 -15235.39230 -9781.94403
9 10506.35337 -15235.39230
10 23414.47216 10506.35337
11 -600.48554 23414.47216
12 -9343.40404 -600.48554
13 18015.35866 -9343.40404
14 17055.98749 18015.35866
15 -65670.85000 17055.98749
16 26354.48554 -65670.85000
17 15630.41454 26354.48554
18 -11612.12797 15630.41454
19 -11262.74909 -11612.12797
20 -3878.35131 -11262.74909
21 110025.39208 -3878.35131
22 23312.54944 110025.39208
23 -28133.91430 23312.54944
24 -85167.42773 -28133.91430
25 -45996.62294 -85167.42773
26 -14935.37275 -45996.62294
27 32510.53014 -14935.37275
28 -4063.54144 32510.53014
29 20087.17737 -4063.54144
30 16510.28521 20087.17737
31 -25132.49030 16510.28521
32 48553.70982 -25132.49030
33 -4089.13216 48553.70982
34 26458.45180 -4089.13216
35 25491.84421 26458.45180
36 63888.39179 25491.84421
37 3621.97758 63888.39179
38 18736.01233 3621.97758
39 2763.30103 18736.01233
40 15843.81691 2763.30103
41 -1728.39900 15843.81691
42 6331.29350 -1728.39900
43 -12776.74148 6331.29350
44 -81018.32699 -12776.74148
45 -1911.33484 -81018.32699
46 -14714.70643 -1911.33484
47 -32551.07734 -14714.70643
48 -50255.62391 -32551.07734
49 -40612.69433 -50255.62391
50 1282.23261 -40612.69433
51 -15256.82238 1282.23261
52 64474.80019 -15256.82238
53 -3515.30902 64474.80019
54 -6053.36940 -3515.30902
55 -13218.07509 -6053.36940
56 4834.59514 -13218.07509
57 -77149.82608 4834.59514
58 14488.53542 -77149.82608
59 39528.44019 14488.53542
60 -290.60990 39528.44019
61 -27070.59051 -290.60990
62 -19629.53912 -27070.59051
63 7303.59050 -19629.53912
64 42884.17448 7303.59050
65 18366.86946 42884.17448
66 -32225.94758 18366.86946
67 -44121.30039 -32225.94758
68 -21636.68289 -44121.30039
69 -17082.13737 -21636.68289
70 -12216.46130 -17082.13737
71 -11208.10244 -12216.46130
72 9315.30316 -11208.10244
73 2039.13048 9315.30316
74 -32348.33020 2039.13048
75 75368.83380 -32348.33020
76 38477.42678 75368.83380
77 29782.91495 38477.42678
78 -25964.30067 29782.91495
79 -6214.23589 -25964.30067
80 -14940.25270 -6214.23589
81 54486.08972 -14940.25270
82 19601.68606 54486.08972
83 -9510.53209 19601.68606
84 -7321.42608 -9510.53209
85 -23463.50915 -7321.42608
86 876.87506 -23463.50915
87 17045.39655 876.87506
88 36521.64182 17045.39655
89 -122510.65607 36521.64182
90 14355.88520 -122510.65607
91 -26692.91851 14355.88520
92 -3515.30740 -26692.91851
93 -16618.09332 -3515.30740
94 36782.54510 -16618.09332
95 6510.25173 36782.54510
96 7601.92196 6510.25173
97 13502.88799 7601.92196
98 -29117.66968 13502.88799
99 -2651.09629 -29117.66968
100 19234.75368 -2651.09629
101 13.41169 19234.75368
102 -9140.04748 13.41169
103 -26044.19559 -9140.04748
104 -9037.71206 -26044.19559
105 -47913.81304 -9037.71206
106 -757.29441 -47913.81304
107 9672.38077 -757.29441
108 -2424.24157 9672.38077
109 -39143.32909 -2424.24157
110 66902.08143 -39143.32909
111 18680.22956 66902.08143
112 -4097.42888 18680.22956
113 -28140.64838 -4097.42888
114 12033.29892 -28140.64838
115 -28574.54489 12033.29892
116 44022.92625 -28574.54489
117 -2496.49976 44022.92625
118 40800.45998 -2496.49976
119 -17651.42793 40800.45998
120 -14217.86223 -17651.42793
121 13673.73924 -14217.86223
122 5234.09540 13673.73924
123 -27467.92748 5234.09540
124 -11350.71213 -27467.92748
125 2924.58018 -11350.71213
126 32353.15929 2924.58018
127 6958.78174 32353.15929
128 40342.44884 6958.78174
129 23721.51387 40342.44884
130 23287.53322 23721.51387
131 34468.56626 23287.53322
132 -17000.22716 34468.56626
133 -2585.38361 -17000.22716
134 -5176.26558 -2585.38361
135 -3268.04826 -5176.26558
136 7118.50995 -3268.04826
137 -28993.31106 7118.50995
138 30130.92328 -28993.31106
139 -26982.83925 30130.92328
140 2576.89831 -26982.83925
141 4619.96001 2576.89831
142 -22885.31622 4619.96001
143 11138.78397 -22885.31622
144 13422.85690 11138.78397
145 -54066.15741 13422.85690
146 -28855.78003 -54066.15741
147 -7513.33032 -28855.78003
148 5277.87325 -7513.33032
149 4830.94525 5277.87325
150 7503.30867 4830.94525
151 1986.40945 7503.30867
152 10179.31201 1986.40945
153 13481.33322 10179.31201
154 2537.62105 13481.33322
155 22221.44904 2537.62105
156 19394.29135 22221.44904
157 -17099.74534 19394.29135
158 -21362.19657 -17099.74534
159 24494.89045 -21362.19657
160 4407.90596 24494.89045
161 -1662.15495 4407.90596
162 7987.78640 -1662.15495
163 -12983.91194 7987.78640
164 -7727.78430 -12983.91194
165 8517.00366 -7727.78430
166 1313.56289 8517.00366
167 26475.59749 1313.56289
168 14705.11111 26475.59749
169 -15354.48315 14705.11111
170 -7970.67398 -15354.48315
171 -8441.70497 -7970.67398
172 61841.16131 -8441.70497
173 32104.87629 61841.16131
174 28104.24463 32104.87629
175 14012.74159 28104.24463
176 21771.69084 14012.74159
177 -3315.50644 21771.69084
178 16850.50093 -3315.50644
179 -5756.29396 16850.50093
180 -4648.81467 -5756.29396
181 -26866.38653 -4648.81467
182 15837.88003 -26866.38653
183 13151.24707 15837.88003
184 30723.36189 13151.24707
185 -6978.51630 30723.36189
186 -2535.48770 -6978.51630
187 -50724.69543 -2535.48770
188 3473.89417 -50724.69543
189 -1913.91984 3473.89417
190 2292.25053 -1913.91984
191 5708.30303 2292.25053
192 -22487.18182 5708.30303
193 -33291.70202 -22487.18182
194 -34595.67091 -33291.70202
195 -16216.77126 -34595.67091
196 -10917.44829 -16216.77126
197 -27620.34008 -10917.44829
198 -12201.56707 -27620.34008
199 -24295.64258 -12201.56707
200 -8559.95448 -24295.64258
201 62039.78441 -8559.95448
202 -16688.80113 62039.78441
203 -10970.95796 -16688.80113
204 -22621.84059 -10970.95796
205 4807.90699 -22621.84059
206 -12550.70419 4807.90699
207 -17612.17048 -12550.70419
208 -1226.41748 -17612.17048
209 -3067.42889 -1226.41748
210 173.16567 -3067.42889
211 -22217.70119 173.16567
212 24475.95957 -22217.70119
213 6693.55291 24475.95957
214 -17334.88722 6693.55291
215 24258.18266 -17334.88722
216 9732.89461 24258.18266
217 -5299.18834 9732.89461
218 5510.52293 -5299.18834
219 -3960.63805 5510.52293
220 -6429.30438 -3960.63805
221 -17716.93207 -6429.30438
222 16464.29566 -17716.93207
223 10137.95114 16464.29566
224 -3426.28300 10137.95114
225 -20546.41344 -3426.28300
226 2629.19094 -20546.41344
227 -4540.20448 2629.19094
228 -11685.94112 -4540.20448
229 -11949.21833 -11685.94112
230 9704.31931 -11949.21833
231 -15060.59538 9704.31931
232 -8360.85598 -15060.59538
233 52168.98176 -8360.85598
234 -12802.59537 52168.98176
235 -11713.86750 -12802.59537
236 9946.57567 -11713.86750
237 -10494.76038 9946.57567
238 12141.34411 -10494.76038
239 2306.81170 12141.34411
240 1953.02801 2306.81170
241 12077.25356 1953.02801
242 -5182.81996 12077.25356
243 -498.29222 -5182.81996
244 17066.18856 -498.29222
245 2168.49389 17066.18856
246 18418.24224 2168.49389
247 8658.23019 18418.24224
248 1254.60436 8658.23019
249 8259.82371 1254.60436
250 -26872.90055 8259.82371
251 -5238.86963 -26872.90055
252 -17681.55002 -5238.86963
253 1040.23893 -17681.55002
254 6840.93166 1040.23893
255 -12429.03136 6840.93166
256 -41785.86753 -12429.03136
257 28610.93029 -41785.86753
258 -24452.13181 28610.93029
259 6586.63497 -24452.13181
260 17516.29541 6586.63497
261 25907.47054 17516.29541
262 -17578.46503 25907.47054
263 61377.18219 -17578.46503
264 -27032.78723 61377.18219
265 3101.30768 -27032.78723
266 -36255.92258 3101.30768
267 -20967.62469 -36255.92258
268 25026.28378 -20967.62469
269 22140.69032 25026.28378
270 -13206.37347 22140.69032
271 -2047.57041 -13206.37347
272 -48401.21391 -2047.57041
273 -33644.92214 -48401.21391
274 -20613.90708 -33644.92214
275 -8901.94575 -20613.90708
276 9553.82222 -8901.94575
277 9435.93190 9553.82222
278 45834.84461 9435.93190
279 -17711.61453 45834.84461
280 -35403.94824 -17711.61453
281 13696.50624 -35403.94824
282 19004.81516 13696.50624
283 21858.69225 19004.81516
284 102115.43370 21858.69225
285 -4802.30560 102115.43370
286 9815.29041 -4802.30560
287 -1386.54546 9815.29041
288 -29776.03775 -1386.54546
> 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/7b4tz1353439957.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/8jqm81353439957.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/9fg9a1353439957.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/10ogky1353439957.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/111gen1353439957.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/12ause1353439957.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/13kama1353439957.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/14b2h21353439957.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/15yhia1353439957.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/16ft321353439957.tab")
+ }
>
> try(system("convert tmp/1756t1353439957.ps tmp/1756t1353439957.png",intern=TRUE))
character(0)
> try(system("convert tmp/2sa1l1353439957.ps tmp/2sa1l1353439957.png",intern=TRUE))
character(0)
> try(system("convert tmp/3e22s1353439957.ps tmp/3e22s1353439957.png",intern=TRUE))
character(0)
> try(system("convert tmp/4w5sm1353439957.ps tmp/4w5sm1353439957.png",intern=TRUE))
character(0)
> try(system("convert tmp/58gtf1353439957.ps tmp/58gtf1353439957.png",intern=TRUE))
character(0)
> try(system("convert tmp/6z4eu1353439957.ps tmp/6z4eu1353439957.png",intern=TRUE))
character(0)
> try(system("convert tmp/7b4tz1353439957.ps tmp/7b4tz1353439957.png",intern=TRUE))
character(0)
> try(system("convert tmp/8jqm81353439957.ps tmp/8jqm81353439957.png",intern=TRUE))
character(0)
> try(system("convert tmp/9fg9a1353439957.ps tmp/9fg9a1353439957.png",intern=TRUE))
character(0)
> try(system("convert tmp/10ogky1353439957.ps tmp/10ogky1353439957.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
21.204 1.880 23.194