R version 2.13.0 (2011-04-13)
Copyright (C) 2011 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i486-pc-linux-gnu (32-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> x <- array(list(1418
+ ,210907
+ ,56
+ ,3
+ ,79
+ ,30
+ ,112285
+ ,869
+ ,120982
+ ,56
+ ,4
+ ,58
+ ,28
+ ,84786
+ ,1530
+ ,176508
+ ,54
+ ,12
+ ,60
+ ,38
+ ,83123
+ ,2172
+ ,179321
+ ,89
+ ,2
+ ,108
+ ,30
+ ,101193
+ ,901
+ ,123185
+ ,40
+ ,1
+ ,49
+ ,22
+ ,38361
+ ,463
+ ,52746
+ ,25
+ ,3
+ ,0
+ ,26
+ ,68504
+ ,3201
+ ,385534
+ ,92
+ ,0
+ ,121
+ ,25
+ ,119182
+ ,371
+ ,33170
+ ,18
+ ,0
+ ,1
+ ,18
+ ,22807
+ ,1192
+ ,101645
+ ,63
+ ,0
+ ,20
+ ,11
+ ,17140
+ ,1583
+ ,149061
+ ,44
+ ,5
+ ,43
+ ,26
+ ,116174
+ ,1439
+ ,165446
+ ,33
+ ,0
+ ,69
+ ,25
+ ,57635
+ ,1764
+ ,237213
+ ,84
+ ,0
+ ,78
+ ,38
+ ,66198
+ ,1495
+ ,173326
+ ,88
+ ,7
+ ,86
+ ,44
+ ,71701
+ ,1373
+ ,133131
+ ,55
+ ,7
+ ,44
+ ,30
+ ,57793
+ ,2187
+ ,258873
+ ,60
+ ,3
+ ,104
+ ,40
+ ,80444
+ ,1491
+ ,180083
+ ,66
+ ,9
+ ,63
+ ,34
+ ,53855
+ ,4041
+ ,324799
+ ,154
+ ,0
+ ,158
+ ,47
+ ,97668
+ ,1706
+ ,230964
+ ,53
+ ,4
+ ,102
+ ,30
+ ,133824
+ ,2152
+ ,236785
+ ,119
+ ,3
+ ,77
+ ,31
+ ,101481
+ ,1036
+ ,135473
+ ,41
+ ,0
+ ,82
+ ,23
+ ,99645
+ ,1882
+ ,202925
+ ,61
+ ,7
+ ,115
+ ,36
+ ,114789
+ ,1929
+ ,215147
+ ,58
+ ,0
+ ,101
+ ,36
+ ,99052
+ ,2242
+ ,344297
+ ,75
+ ,1
+ ,80
+ ,30
+ ,67654
+ ,1220
+ ,153935
+ ,33
+ ,5
+ ,50
+ ,25
+ ,65553
+ ,1289
+ ,132943
+ ,40
+ ,7
+ ,83
+ ,39
+ ,97500
+ ,2515
+ ,174724
+ ,92
+ ,0
+ ,123
+ ,34
+ ,69112
+ ,2147
+ ,174415
+ ,100
+ ,0
+ ,73
+ ,31
+ ,82753
+ ,2352
+ ,225548
+ ,112
+ ,5
+ ,81
+ ,31
+ ,85323
+ ,1638
+ ,223632
+ ,73
+ ,0
+ ,105
+ ,33
+ ,72654
+ ,1222
+ ,124817
+ ,40
+ ,0
+ ,47
+ ,25
+ ,30727
+ ,1812
+ ,221698
+ ,45
+ ,0
+ ,105
+ ,33
+ ,77873
+ ,1677
+ ,210767
+ ,60
+ ,3
+ ,94
+ ,35
+ ,117478
+ ,1579
+ ,170266
+ ,62
+ ,4
+ ,44
+ ,42
+ ,74007
+ ,1731
+ ,260561
+ ,75
+ ,1
+ ,114
+ ,43
+ ,90183
+ ,807
+ ,84853
+ ,31
+ ,4
+ ,38
+ ,30
+ ,61542
+ ,2452
+ ,294424
+ ,77
+ ,2
+ ,107
+ ,33
+ ,101494
+ ,829
+ ,101011
+ ,34
+ ,0
+ ,30
+ ,13
+ ,27570
+ ,1940
+ ,215641
+ ,46
+ ,0
+ ,71
+ ,32
+ ,55813
+ ,2662
+ ,325107
+ ,99
+ ,0
+ ,84
+ ,36
+ ,79215
+ ,186
+ ,7176
+ ,17
+ ,0
+ ,0
+ ,0
+ ,1423
+ ,1499
+ ,167542
+ ,66
+ ,2
+ ,59
+ ,28
+ ,55461
+ ,865
+ ,106408
+ ,30
+ ,1
+ ,33
+ ,14
+ ,31081
+ ,1793
+ ,96560
+ ,76
+ ,0
+ ,42
+ ,17
+ ,22996
+ ,2527
+ ,265769
+ ,146
+ ,2
+ ,96
+ ,32
+ ,83122
+ ,2747
+ ,269651
+ ,67
+ ,10
+ ,106
+ ,30
+ ,70106
+ ,1324
+ ,149112
+ ,56
+ ,6
+ ,56
+ ,35
+ ,60578
+ ,2702
+ ,175824
+ ,107
+ ,0
+ ,57
+ ,20
+ ,39992
+ ,1383
+ ,152871
+ ,58
+ ,5
+ ,59
+ ,28
+ ,79892
+ ,1179
+ ,111665
+ ,34
+ ,4
+ ,39
+ ,28
+ ,49810
+ ,2099
+ ,116408
+ ,61
+ ,1
+ ,34
+ ,39
+ ,71570
+ ,4308
+ ,362301
+ ,119
+ ,2
+ ,76
+ ,34
+ ,100708
+ ,918
+ ,78800
+ ,42
+ ,2
+ ,20
+ ,26
+ ,33032
+ ,1831
+ ,183167
+ ,66
+ ,0
+ ,91
+ ,39
+ ,82875
+ ,3373
+ ,277965
+ ,89
+ ,8
+ ,115
+ ,39
+ ,139077
+ ,1713
+ ,150629
+ ,44
+ ,3
+ ,85
+ ,33
+ ,71595
+ ,1438
+ ,168809
+ ,66
+ ,0
+ ,76
+ ,28
+ ,72260
+ ,496
+ ,24188
+ ,24
+ ,0
+ ,8
+ ,4
+ ,5950
+ ,2253
+ ,329267
+ ,259
+ ,8
+ ,79
+ ,39
+ ,115762
+ ,744
+ ,65029
+ ,17
+ ,5
+ ,21
+ ,18
+ ,32551
+ ,1161
+ ,101097
+ ,64
+ ,3
+ ,30
+ ,14
+ ,31701
+ ,2352
+ ,218946
+ ,41
+ ,1
+ ,76
+ ,29
+ ,80670
+ ,2144
+ ,244052
+ ,68
+ ,5
+ ,101
+ ,44
+ ,143558
+ ,4691
+ ,341570
+ ,168
+ ,1
+ ,94
+ ,21
+ ,117105
+ ,1112
+ ,103597
+ ,43
+ ,1
+ ,27
+ ,16
+ ,23789
+ ,2694
+ ,233328
+ ,132
+ ,5
+ ,92
+ ,28
+ ,120733
+ ,1973
+ ,256462
+ ,105
+ ,0
+ ,123
+ ,35
+ ,105195
+ ,1769
+ ,206161
+ ,71
+ ,12
+ ,75
+ ,28
+ ,73107
+ ,3148
+ ,311473
+ ,112
+ ,8
+ ,128
+ ,38
+ ,132068
+ ,2474
+ ,235800
+ ,94
+ ,8
+ ,105
+ ,23
+ ,149193
+ ,2084
+ ,177939
+ ,82
+ ,8
+ ,55
+ ,36
+ ,46821
+ ,1954
+ ,207176
+ ,70
+ ,8
+ ,56
+ ,32
+ ,87011
+ ,1226
+ ,196553
+ ,57
+ ,2
+ ,41
+ ,29
+ ,95260
+ ,1389
+ ,174184
+ ,53
+ ,0
+ ,72
+ ,25
+ ,55183
+ ,1496
+ ,143246
+ ,103
+ ,5
+ ,67
+ ,27
+ ,106671
+ ,2269
+ ,187559
+ ,121
+ ,8
+ ,75
+ ,36
+ ,73511
+ ,1833
+ ,187681
+ ,62
+ ,2
+ ,114
+ ,28
+ ,92945
+ ,1268
+ ,119016
+ ,52
+ ,5
+ ,118
+ ,23
+ ,78664
+ ,1943
+ ,182192
+ ,52
+ ,12
+ ,77
+ ,40
+ ,70054
+ ,893
+ ,73566
+ ,32
+ ,6
+ ,22
+ ,23
+ ,22618
+ ,1762
+ ,194979
+ ,62
+ ,7
+ ,66
+ ,40
+ ,74011
+ ,1403
+ ,167488
+ ,45
+ ,2
+ ,69
+ ,28
+ ,83737
+ ,1425
+ ,143756
+ ,46
+ ,0
+ ,105
+ ,34
+ ,69094
+ ,1857
+ ,275541
+ ,63
+ ,4
+ ,116
+ ,33
+ ,93133
+ ,1840
+ ,243199
+ ,75
+ ,3
+ ,88
+ ,28
+ ,95536
+ ,1502
+ ,182999
+ ,88
+ ,6
+ ,73
+ ,34
+ ,225920
+ ,1441
+ ,135649
+ ,46
+ ,2
+ ,99
+ ,30
+ ,62133
+ ,1420
+ ,152299
+ ,53
+ ,0
+ ,62
+ ,33
+ ,61370
+ ,1416
+ ,120221
+ ,37
+ ,1
+ ,53
+ ,22
+ ,43836
+ ,2970
+ ,346485
+ ,90
+ ,0
+ ,118
+ ,38
+ ,106117
+ ,1317
+ ,145790
+ ,63
+ ,5
+ ,30
+ ,26
+ ,38692
+ ,1644
+ ,193339
+ ,78
+ ,2
+ ,100
+ ,35
+ ,84651
+ ,870
+ ,80953
+ ,25
+ ,0
+ ,49
+ ,8
+ ,56622
+ ,1654
+ ,122774
+ ,45
+ ,0
+ ,24
+ ,24
+ ,15986
+ ,1054
+ ,130585
+ ,46
+ ,5
+ ,67
+ ,29
+ ,95364
+ ,937
+ ,112611
+ ,41
+ ,0
+ ,46
+ ,20
+ ,26706
+ ,3004
+ ,286468
+ ,144
+ ,1
+ ,57
+ ,29
+ ,89691
+ ,2008
+ ,241066
+ ,82
+ ,0
+ ,75
+ ,45
+ ,67267
+ ,2547
+ ,148446
+ ,91
+ ,1
+ ,135
+ ,37
+ ,126846
+ ,1885
+ ,204713
+ ,71
+ ,1
+ ,68
+ ,33
+ ,41140
+ ,1626
+ ,182079
+ ,63
+ ,2
+ ,124
+ ,33
+ ,102860
+ ,1468
+ ,140344
+ ,53
+ ,6
+ ,33
+ ,25
+ ,51715
+ ,2445
+ ,220516
+ ,62
+ ,1
+ ,98
+ ,32
+ ,55801
+ ,1964
+ ,243060
+ ,63
+ ,4
+ ,58
+ ,29
+ ,111813
+ ,1381
+ ,162765
+ ,32
+ ,2
+ ,68
+ ,28
+ ,120293
+ ,1369
+ ,182613
+ ,39
+ ,3
+ ,81
+ ,28
+ ,138599
+ ,1659
+ ,232138
+ ,62
+ ,0
+ ,131
+ ,31
+ ,161647
+ ,2888
+ ,265318
+ ,117
+ ,10
+ ,110
+ ,52
+ ,115929
+ ,1290
+ ,85574
+ ,34
+ ,0
+ ,37
+ ,21
+ ,24266
+ ,2845
+ ,310839
+ ,92
+ ,9
+ ,130
+ ,24
+ ,162901
+ ,1982
+ ,225060
+ ,93
+ ,7
+ ,93
+ ,41
+ ,109825
+ ,1904
+ ,232317
+ ,54
+ ,0
+ ,118
+ ,33
+ ,129838
+ ,1391
+ ,144966
+ ,144
+ ,0
+ ,39
+ ,32
+ ,37510
+ ,602
+ ,43287
+ ,14
+ ,4
+ ,13
+ ,19
+ ,43750
+ ,1743
+ ,155754
+ ,61
+ ,4
+ ,74
+ ,20
+ ,40652
+ ,1559
+ ,164709
+ ,109
+ ,0
+ ,81
+ ,31
+ ,87771
+ ,2014
+ ,201940
+ ,38
+ ,0
+ ,109
+ ,31
+ ,85872
+ ,2143
+ ,235454
+ ,73
+ ,0
+ ,151
+ ,32
+ ,89275
+ ,2146
+ ,220801
+ ,75
+ ,1
+ ,51
+ ,18
+ ,44418
+ ,874
+ ,99466
+ ,50
+ ,0
+ ,28
+ ,23
+ ,192565
+ ,1590
+ ,92661
+ ,61
+ ,1
+ ,40
+ ,17
+ ,35232
+ ,1590
+ ,133328
+ ,55
+ ,0
+ ,56
+ ,20
+ ,40909
+ ,1210
+ ,61361
+ ,77
+ ,0
+ ,27
+ ,12
+ ,13294
+ ,2072
+ ,125930
+ ,75
+ ,4
+ ,37
+ ,17
+ ,32387
+ ,1281
+ ,100750
+ ,72
+ ,0
+ ,83
+ ,30
+ ,140867
+ ,1401
+ ,224549
+ ,50
+ ,4
+ ,54
+ ,31
+ ,120662
+ ,834
+ ,82316
+ ,32
+ ,4
+ ,27
+ ,10
+ ,21233
+ ,1105
+ ,102010
+ ,53
+ ,3
+ ,28
+ ,13
+ ,44332
+ ,1272
+ ,101523
+ ,42
+ ,0
+ ,59
+ ,22
+ ,61056
+ ,1944
+ ,243511
+ ,71
+ ,0
+ ,133
+ ,42
+ ,101338
+ ,391
+ ,22938
+ ,10
+ ,0
+ ,12
+ ,1
+ ,1168
+ ,761
+ ,41566
+ ,35
+ ,5
+ ,0
+ ,9
+ ,13497
+ ,1605
+ ,152474
+ ,65
+ ,0
+ ,106
+ ,32
+ ,65567
+ ,530
+ ,61857
+ ,25
+ ,4
+ ,23
+ ,11
+ ,25162
+ ,1988
+ ,99923
+ ,66
+ ,0
+ ,44
+ ,25
+ ,32334
+ ,1386
+ ,132487
+ ,41
+ ,0
+ ,71
+ ,36
+ ,40735
+ ,2395
+ ,317394
+ ,86
+ ,1
+ ,116
+ ,31
+ ,91413
+ ,387
+ ,21054
+ ,16
+ ,0
+ ,4
+ ,0
+ ,855
+ ,1742
+ ,209641
+ ,42
+ ,5
+ ,62
+ ,24
+ ,97068
+ ,620
+ ,22648
+ ,19
+ ,0
+ ,12
+ ,13
+ ,44339
+ ,449
+ ,31414
+ ,19
+ ,0
+ ,18
+ ,8
+ ,14116
+ ,800
+ ,46698
+ ,45
+ ,0
+ ,14
+ ,13
+ ,10288
+ ,1684
+ ,131698
+ ,65
+ ,0
+ ,60
+ ,19
+ ,65622
+ ,1050
+ ,91735
+ ,35
+ ,0
+ ,7
+ ,18
+ ,16563
+ ,2699
+ ,244749
+ ,95
+ ,2
+ ,98
+ ,33
+ ,76643
+ ,1606
+ ,184510
+ ,49
+ ,7
+ ,64
+ ,40
+ ,110681
+ ,1502
+ ,79863
+ ,37
+ ,1
+ ,29
+ ,22
+ ,29011
+ ,1204
+ ,128423
+ ,64
+ ,8
+ ,32
+ ,38
+ ,92696
+ ,1138
+ ,97839
+ ,38
+ ,2
+ ,25
+ ,24
+ ,94785
+ ,568
+ ,38214
+ ,34
+ ,0
+ ,16
+ ,8
+ ,8773
+ ,1459
+ ,151101
+ ,32
+ ,2
+ ,48
+ ,35
+ ,83209
+ ,2158
+ ,272458
+ ,65
+ ,0
+ ,100
+ ,43
+ ,93815
+ ,1111
+ ,172494
+ ,52
+ ,0
+ ,46
+ ,43
+ ,86687
+ ,1421
+ ,108043
+ ,62
+ ,1
+ ,45
+ ,14
+ ,34553
+ ,2833
+ ,328107
+ ,65
+ ,3
+ ,129
+ ,41
+ ,105547
+ ,1955
+ ,250579
+ ,83
+ ,0
+ ,130
+ ,38
+ ,103487
+ ,2922
+ ,351067
+ ,95
+ ,3
+ ,136
+ ,45
+ ,213688
+ ,1002
+ ,158015
+ ,29
+ ,0
+ ,59
+ ,31
+ ,71220
+ ,1060
+ ,98866
+ ,18
+ ,0
+ ,25
+ ,13
+ ,23517
+ ,956
+ ,85439
+ ,33
+ ,0
+ ,32
+ ,28
+ ,56926
+ ,2186
+ ,229242
+ ,247
+ ,4
+ ,63
+ ,31
+ ,91721
+ ,3604
+ ,351619
+ ,139
+ ,4
+ ,95
+ ,40
+ ,115168
+ ,1035
+ ,84207
+ ,29
+ ,11
+ ,14
+ ,30
+ ,111194
+ ,1417
+ ,120445
+ ,118
+ ,0
+ ,36
+ ,16
+ ,51009
+ ,3261
+ ,324598
+ ,110
+ ,0
+ ,113
+ ,37
+ ,135777
+ ,1587
+ ,131069
+ ,67
+ ,4
+ ,47
+ ,30
+ ,51513
+ ,1424
+ ,204271
+ ,42
+ ,0
+ ,92
+ ,35
+ ,74163
+ ,1701
+ ,165543
+ ,65
+ ,1
+ ,70
+ ,32
+ ,51633
+ ,1249
+ ,141722
+ ,94
+ ,0
+ ,19
+ ,27
+ ,75345
+ ,946
+ ,116048
+ ,64
+ ,0
+ ,50
+ ,20
+ ,33416
+ ,1926
+ ,250047
+ ,81
+ ,0
+ ,41
+ ,18
+ ,83305
+ ,3352
+ ,299775
+ ,95
+ ,9
+ ,91
+ ,31
+ ,98952
+ ,1641
+ ,195838
+ ,67
+ ,1
+ ,111
+ ,31
+ ,102372
+ ,2035
+ ,173260
+ ,63
+ ,3
+ ,41
+ ,21
+ ,37238
+ ,2312
+ ,254488
+ ,83
+ ,10
+ ,120
+ ,39
+ ,103772
+ ,1369
+ ,104389
+ ,45
+ ,5
+ ,135
+ ,41
+ ,123969
+ ,1577
+ ,136084
+ ,30
+ ,0
+ ,27
+ ,13
+ ,27142
+ ,2201
+ ,199476
+ ,70
+ ,2
+ ,87
+ ,32
+ ,135400
+ ,961
+ ,92499
+ ,32
+ ,0
+ ,25
+ ,18
+ ,21399
+ ,1900
+ ,224330
+ ,83
+ ,1
+ ,131
+ ,39
+ ,130115
+ ,1254
+ ,135781
+ ,31
+ ,2
+ ,45
+ ,14
+ ,24874
+ ,1335
+ ,74408
+ ,67
+ ,4
+ ,29
+ ,7
+ ,34988
+ ,1597
+ ,81240
+ ,66
+ ,0
+ ,58
+ ,17
+ ,45549
+ ,207
+ ,14688
+ ,10
+ ,0
+ ,4
+ ,0
+ ,6023
+ ,1645
+ ,181633
+ ,70
+ ,2
+ ,47
+ ,30
+ ,64466
+ ,2429
+ ,271856
+ ,103
+ ,1
+ ,109
+ ,37
+ ,54990
+ ,151
+ ,7199
+ ,5
+ ,0
+ ,7
+ ,0
+ ,1644
+ ,474
+ ,46660
+ ,20
+ ,0
+ ,12
+ ,5
+ ,6179
+ ,141
+ ,17547
+ ,5
+ ,0
+ ,0
+ ,1
+ ,3926
+ ,1639
+ ,133368
+ ,36
+ ,1
+ ,37
+ ,16
+ ,32755
+ ,872
+ ,95227
+ ,34
+ ,0
+ ,37
+ ,32
+ ,34777
+ ,1318
+ ,152601
+ ,48
+ ,2
+ ,46
+ ,24
+ ,73224
+ ,1018
+ ,98146
+ ,40
+ ,0
+ ,15
+ ,17
+ ,27114
+ ,1383
+ ,79619
+ ,43
+ ,3
+ ,42
+ ,11
+ ,20760
+ ,1314
+ ,59194
+ ,31
+ ,6
+ ,7
+ ,24
+ ,37636
+ ,1335
+ ,139942
+ ,42
+ ,0
+ ,54
+ ,22
+ ,65461
+ ,1403
+ ,118612
+ ,46
+ ,2
+ ,54
+ ,12
+ ,30080
+ ,910
+ ,72880
+ ,33
+ ,0
+ ,14
+ ,19
+ ,24094
+ ,616
+ ,65475
+ ,18
+ ,2
+ ,16
+ ,13
+ ,69008
+ ,1407
+ ,99643
+ ,55
+ ,1
+ ,33
+ ,17
+ ,54968
+ ,771
+ ,71965
+ ,35
+ ,1
+ ,32
+ ,15
+ ,46090
+ ,766
+ ,77272
+ ,59
+ ,2
+ ,21
+ ,16
+ ,27507
+ ,473
+ ,49289
+ ,19
+ ,1
+ ,15
+ ,24
+ ,10672
+ ,1376
+ ,135131
+ ,66
+ ,0
+ ,38
+ ,15
+ ,34029
+ ,1232
+ ,108446
+ ,60
+ ,1
+ ,22
+ ,17
+ ,46300
+ ,1521
+ ,89746
+ ,36
+ ,3
+ ,28
+ ,18
+ ,24760
+ ,572
+ ,44296
+ ,25
+ ,0
+ ,10
+ ,20
+ ,18779
+ ,1059
+ ,77648
+ ,47
+ ,0
+ ,31
+ ,16
+ ,21280
+ ,1544
+ ,181528
+ ,54
+ ,0
+ ,32
+ ,16
+ ,40662
+ ,1230
+ ,134019
+ ,53
+ ,0
+ ,32
+ ,18
+ ,28987
+ ,1206
+ ,124064
+ ,40
+ ,1
+ ,43
+ ,22
+ ,22827
+ ,1205
+ ,92630
+ ,40
+ ,4
+ ,27
+ ,8
+ ,18513
+ ,1255
+ ,121848
+ ,39
+ ,0
+ ,37
+ ,17
+ ,30594
+ ,613
+ ,52915
+ ,14
+ ,0
+ ,20
+ ,18
+ ,24006
+ ,721
+ ,81872
+ ,45
+ ,0
+ ,32
+ ,16
+ ,27913
+ ,1109
+ ,58981
+ ,36
+ ,7
+ ,0
+ ,23
+ ,42744
+ ,740
+ ,53515
+ ,28
+ ,2
+ ,5
+ ,22
+ ,12934
+ ,1126
+ ,60812
+ ,44
+ ,0
+ ,26
+ ,13
+ ,22574
+ ,728
+ ,56375
+ ,30
+ ,7
+ ,10
+ ,13
+ ,41385
+ ,689
+ ,65490
+ ,22
+ ,3
+ ,27
+ ,16
+ ,18653
+ ,592
+ ,80949
+ ,17
+ ,0
+ ,11
+ ,16
+ ,18472
+ ,995
+ ,76302
+ ,31
+ ,0
+ ,29
+ ,20
+ ,30976
+ ,1613
+ ,104011
+ ,55
+ ,6
+ ,25
+ ,22
+ ,63339
+ ,2048
+ ,98104
+ ,54
+ ,2
+ ,55
+ ,17
+ ,25568
+ ,705
+ ,67989
+ ,21
+ ,0
+ ,23
+ ,18
+ ,33747
+ ,301
+ ,30989
+ ,14
+ ,0
+ ,5
+ ,17
+ ,4154
+ ,1803
+ ,135458
+ ,81
+ ,3
+ ,43
+ ,12
+ ,19474
+ ,799
+ ,73504
+ ,35
+ ,0
+ ,23
+ ,7
+ ,35130
+ ,861
+ ,63123
+ ,43
+ ,1
+ ,34
+ ,17
+ ,39067
+ ,1186
+ ,61254
+ ,46
+ ,1
+ ,36
+ ,14
+ ,13310
+ ,1451
+ ,74914
+ ,30
+ ,0
+ ,35
+ ,23
+ ,65892
+ ,628
+ ,31774
+ ,23
+ ,1
+ ,0
+ ,17
+ ,4143
+ ,1161
+ ,81437
+ ,38
+ ,0
+ ,37
+ ,14
+ ,28579
+ ,1463
+ ,87186
+ ,54
+ ,0
+ ,28
+ ,15
+ ,51776
+ ,742
+ ,50090
+ ,20
+ ,0
+ ,16
+ ,17
+ ,21152
+ ,979
+ ,65745
+ ,53
+ ,0
+ ,26
+ ,21
+ ,38084
+ ,675
+ ,56653
+ ,45
+ ,0
+ ,38
+ ,18
+ ,27717
+ ,1241
+ ,158399
+ ,39
+ ,0
+ ,23
+ ,18
+ ,32928
+ ,676
+ ,46455
+ ,20
+ ,0
+ ,22
+ ,17
+ ,11342
+ ,1049
+ ,73624
+ ,24
+ ,0
+ ,30
+ ,17
+ ,19499
+ ,620
+ ,38395
+ ,31
+ ,0
+ ,16
+ ,16
+ ,16380
+ ,1081
+ ,91899
+ ,35
+ ,0
+ ,18
+ ,15
+ ,36874
+ ,1688
+ ,139526
+ ,151
+ ,0
+ ,28
+ ,21
+ ,48259
+ ,736
+ ,52164
+ ,52
+ ,0
+ ,32
+ ,16
+ ,16734
+ ,617
+ ,51567
+ ,30
+ ,2
+ ,21
+ ,14
+ ,28207
+ ,812
+ ,70551
+ ,31
+ ,0
+ ,23
+ ,15
+ ,30143
+ ,1051
+ ,84856
+ ,29
+ ,1
+ ,29
+ ,17
+ ,41369
+ ,1656
+ ,102538
+ ,57
+ ,1
+ ,50
+ ,15
+ ,45833
+ ,705
+ ,86678
+ ,40
+ ,0
+ ,12
+ ,15
+ ,29156
+ ,945
+ ,85709
+ ,44
+ ,0
+ ,21
+ ,10
+ ,35944
+ ,554
+ ,34662
+ ,25
+ ,0
+ ,18
+ ,6
+ ,36278
+ ,1597
+ ,150580
+ ,77
+ ,0
+ ,27
+ ,22
+ ,45588
+ ,982
+ ,99611
+ ,35
+ ,0
+ ,41
+ ,21
+ ,45097
+ ,222
+ ,19349
+ ,11
+ ,0
+ ,13
+ ,1
+ ,3895
+ ,1212
+ ,99373
+ ,63
+ ,1
+ ,12
+ ,18
+ ,28394
+ ,1143
+ ,86230
+ ,44
+ ,0
+ ,21
+ ,17
+ ,18632
+ ,435
+ ,30837
+ ,19
+ ,0
+ ,8
+ ,4
+ ,2325
+ ,532
+ ,31706
+ ,13
+ ,0
+ ,26
+ ,10
+ ,25139
+ ,882
+ ,89806
+ ,42
+ ,0
+ ,27
+ ,16
+ ,27975
+ ,608
+ ,62088
+ ,38
+ ,1
+ ,13
+ ,16
+ ,14483
+ ,459
+ ,40151
+ ,29
+ ,0
+ ,16
+ ,9
+ ,13127
+ ,578
+ ,27634
+ ,20
+ ,0
+ ,2
+ ,16
+ ,5839
+ ,826
+ ,76990
+ ,27
+ ,0
+ ,42
+ ,17
+ ,24069
+ ,509
+ ,37460
+ ,20
+ ,0
+ ,5
+ ,7
+ ,3738
+ ,717
+ ,54157
+ ,19
+ ,0
+ ,37
+ ,15
+ ,18625
+ ,637
+ ,49862
+ ,37
+ ,0
+ ,17
+ ,14
+ ,36341
+ ,857
+ ,84337
+ ,26
+ ,0
+ ,38
+ ,14
+ ,24548
+ ,830
+ ,64175
+ ,42
+ ,0
+ ,37
+ ,18
+ ,21792
+ ,652
+ ,59382
+ ,49
+ ,0
+ ,29
+ ,12
+ ,26263
+ ,707
+ ,119308
+ ,30
+ ,0
+ ,32
+ ,16
+ ,23686
+ ,954
+ ,76702
+ ,49
+ ,0
+ ,35
+ ,21
+ ,49303
+ ,1461
+ ,103425
+ ,67
+ ,1
+ ,17
+ ,19
+ ,25659
+ ,672
+ ,70344
+ ,28
+ ,0
+ ,20
+ ,16
+ ,28904
+ ,778
+ ,43410
+ ,19
+ ,0
+ ,7
+ ,1
+ ,2781
+ ,1141
+ ,104838
+ ,49
+ ,1
+ ,46
+ ,16
+ ,29236
+ ,680
+ ,62215
+ ,27
+ ,0
+ ,24
+ ,10
+ ,19546
+ ,1090
+ ,69304
+ ,30
+ ,6
+ ,40
+ ,19
+ ,22818
+ ,616
+ ,53117
+ ,22
+ ,3
+ ,3
+ ,12
+ ,32689
+ ,285
+ ,19764
+ ,12
+ ,1
+ ,10
+ ,2
+ ,5752
+ ,1145
+ ,86680
+ ,31
+ ,2
+ ,37
+ ,14
+ ,22197
+ ,733
+ ,84105
+ ,20
+ ,0
+ ,17
+ ,17
+ ,20055
+ ,888
+ ,77945
+ ,20
+ ,0
+ ,28
+ ,19
+ ,25272
+ ,849
+ ,89113
+ ,39
+ ,0
+ ,19
+ ,14
+ ,82206
+ ,1182
+ ,91005
+ ,29
+ ,3
+ ,29
+ ,11
+ ,32073
+ ,528
+ ,40248
+ ,16
+ ,1
+ ,8
+ ,4
+ ,5444
+ ,642
+ ,64187
+ ,27
+ ,0
+ ,10
+ ,16
+ ,20154
+ ,947
+ ,50857
+ ,21
+ ,0
+ ,15
+ ,20
+ ,36944
+ ,819
+ ,56613
+ ,19
+ ,1
+ ,15
+ ,12
+ ,8019
+ ,757
+ ,62792
+ ,35
+ ,0
+ ,28
+ ,15
+ ,30884
+ ,894
+ ,72535
+ ,14
+ ,0
+ ,17
+ ,16
+ ,19540)
+ ,dim=c(7
+ ,289)
+ ,dimnames=list(c('pageviews'
+ ,'time'
+ ,'logins'
+ ,'shared'
+ ,'blogged'
+ ,'reviewed'
+ ,'totsize')
+ ,1:289))
> y <- array(NA,dim=c(7,289),dimnames=list(c('pageviews','time','logins','shared','blogged','reviewed','totsize'),1:289))
> for (i in 1:dim(x)[1])
+ {
+ for (j in 1:dim(x)[2])
+ {
+ y[i,j] <- as.numeric(x[i,j])
+ }
+ }
> par3 = 'No Linear Trend'
> par2 = 'Do not include Seasonal Dummies'
> par1 = '7'
> library(lattice)
> library(lmtest)
Loading required package: zoo
> n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test
> par1 <- as.numeric(par1)
> x <- t(y)
> k <- length(x[1,])
> n <- length(x[,1])
> x1 <- cbind(x[,par1], x[,1:k!=par1])
> mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1])
> colnames(x1) <- mycolnames #colnames(x)[par1]
> x <- x1
> if (par3 == 'First Differences'){
+ x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep='')))
+ for (i in 1:n-1) {
+ for (j in 1:k) {
+ x2[i,j] <- x[i+1,j] - x[i,j]
+ }
+ }
+ x <- x2
+ }
> if (par2 == 'Include Monthly Dummies'){
+ x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep ='')))
+ for (i in 1:11){
+ x2[seq(i,n,12),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> if (par2 == 'Include Quarterly Dummies'){
+ x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep ='')))
+ for (i in 1:3){
+ x2[seq(i,n,4),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> k <- length(x[1,])
> if (par3 == 'Linear Trend'){
+ x <- cbind(x, c(1:n))
+ colnames(x)[k+1] <- 't'
+ }
> x
totsize pageviews time logins shared blogged reviewed
1 112285 1418 210907 56 3 79 30
2 84786 869 120982 56 4 58 28
3 83123 1530 176508 54 12 60 38
4 101193 2172 179321 89 2 108 30
5 38361 901 123185 40 1 49 22
6 68504 463 52746 25 3 0 26
7 119182 3201 385534 92 0 121 25
8 22807 371 33170 18 0 1 18
9 17140 1192 101645 63 0 20 11
10 116174 1583 149061 44 5 43 26
11 57635 1439 165446 33 0 69 25
12 66198 1764 237213 84 0 78 38
13 71701 1495 173326 88 7 86 44
14 57793 1373 133131 55 7 44 30
15 80444 2187 258873 60 3 104 40
16 53855 1491 180083 66 9 63 34
17 97668 4041 324799 154 0 158 47
18 133824 1706 230964 53 4 102 30
19 101481 2152 236785 119 3 77 31
20 99645 1036 135473 41 0 82 23
21 114789 1882 202925 61 7 115 36
22 99052 1929 215147 58 0 101 36
23 67654 2242 344297 75 1 80 30
24 65553 1220 153935 33 5 50 25
25 97500 1289 132943 40 7 83 39
26 69112 2515 174724 92 0 123 34
27 82753 2147 174415 100 0 73 31
28 85323 2352 225548 112 5 81 31
29 72654 1638 223632 73 0 105 33
30 30727 1222 124817 40 0 47 25
31 77873 1812 221698 45 0 105 33
32 117478 1677 210767 60 3 94 35
33 74007 1579 170266 62 4 44 42
34 90183 1731 260561 75 1 114 43
35 61542 807 84853 31 4 38 30
36 101494 2452 294424 77 2 107 33
37 27570 829 101011 34 0 30 13
38 55813 1940 215641 46 0 71 32
39 79215 2662 325107 99 0 84 36
40 1423 186 7176 17 0 0 0
41 55461 1499 167542 66 2 59 28
42 31081 865 106408 30 1 33 14
43 22996 1793 96560 76 0 42 17
44 83122 2527 265769 146 2 96 32
45 70106 2747 269651 67 10 106 30
46 60578 1324 149112 56 6 56 35
47 39992 2702 175824 107 0 57 20
48 79892 1383 152871 58 5 59 28
49 49810 1179 111665 34 4 39 28
50 71570 2099 116408 61 1 34 39
51 100708 4308 362301 119 2 76 34
52 33032 918 78800 42 2 20 26
53 82875 1831 183167 66 0 91 39
54 139077 3373 277965 89 8 115 39
55 71595 1713 150629 44 3 85 33
56 72260 1438 168809 66 0 76 28
57 5950 496 24188 24 0 8 4
58 115762 2253 329267 259 8 79 39
59 32551 744 65029 17 5 21 18
60 31701 1161 101097 64 3 30 14
61 80670 2352 218946 41 1 76 29
62 143558 2144 244052 68 5 101 44
63 117105 4691 341570 168 1 94 21
64 23789 1112 103597 43 1 27 16
65 120733 2694 233328 132 5 92 28
66 105195 1973 256462 105 0 123 35
67 73107 1769 206161 71 12 75 28
68 132068 3148 311473 112 8 128 38
69 149193 2474 235800 94 8 105 23
70 46821 2084 177939 82 8 55 36
71 87011 1954 207176 70 8 56 32
72 95260 1226 196553 57 2 41 29
73 55183 1389 174184 53 0 72 25
74 106671 1496 143246 103 5 67 27
75 73511 2269 187559 121 8 75 36
76 92945 1833 187681 62 2 114 28
77 78664 1268 119016 52 5 118 23
78 70054 1943 182192 52 12 77 40
79 22618 893 73566 32 6 22 23
80 74011 1762 194979 62 7 66 40
81 83737 1403 167488 45 2 69 28
82 69094 1425 143756 46 0 105 34
83 93133 1857 275541 63 4 116 33
84 95536 1840 243199 75 3 88 28
85 225920 1502 182999 88 6 73 34
86 62133 1441 135649 46 2 99 30
87 61370 1420 152299 53 0 62 33
88 43836 1416 120221 37 1 53 22
89 106117 2970 346485 90 0 118 38
90 38692 1317 145790 63 5 30 26
91 84651 1644 193339 78 2 100 35
92 56622 870 80953 25 0 49 8
93 15986 1654 122774 45 0 24 24
94 95364 1054 130585 46 5 67 29
95 26706 937 112611 41 0 46 20
96 89691 3004 286468 144 1 57 29
97 67267 2008 241066 82 0 75 45
98 126846 2547 148446 91 1 135 37
99 41140 1885 204713 71 1 68 33
100 102860 1626 182079 63 2 124 33
101 51715 1468 140344 53 6 33 25
102 55801 2445 220516 62 1 98 32
103 111813 1964 243060 63 4 58 29
104 120293 1381 162765 32 2 68 28
105 138599 1369 182613 39 3 81 28
106 161647 1659 232138 62 0 131 31
107 115929 2888 265318 117 10 110 52
108 24266 1290 85574 34 0 37 21
109 162901 2845 310839 92 9 130 24
110 109825 1982 225060 93 7 93 41
111 129838 1904 232317 54 0 118 33
112 37510 1391 144966 144 0 39 32
113 43750 602 43287 14 4 13 19
114 40652 1743 155754 61 4 74 20
115 87771 1559 164709 109 0 81 31
116 85872 2014 201940 38 0 109 31
117 89275 2143 235454 73 0 151 32
118 44418 2146 220801 75 1 51 18
119 192565 874 99466 50 0 28 23
120 35232 1590 92661 61 1 40 17
121 40909 1590 133328 55 0 56 20
122 13294 1210 61361 77 0 27 12
123 32387 2072 125930 75 4 37 17
124 140867 1281 100750 72 0 83 30
125 120662 1401 224549 50 4 54 31
126 21233 834 82316 32 4 27 10
127 44332 1105 102010 53 3 28 13
128 61056 1272 101523 42 0 59 22
129 101338 1944 243511 71 0 133 42
130 1168 391 22938 10 0 12 1
131 13497 761 41566 35 5 0 9
132 65567 1605 152474 65 0 106 32
133 25162 530 61857 25 4 23 11
134 32334 1988 99923 66 0 44 25
135 40735 1386 132487 41 0 71 36
136 91413 2395 317394 86 1 116 31
137 855 387 21054 16 0 4 0
138 97068 1742 209641 42 5 62 24
139 44339 620 22648 19 0 12 13
140 14116 449 31414 19 0 18 8
141 10288 800 46698 45 0 14 13
142 65622 1684 131698 65 0 60 19
143 16563 1050 91735 35 0 7 18
144 76643 2699 244749 95 2 98 33
145 110681 1606 184510 49 7 64 40
146 29011 1502 79863 37 1 29 22
147 92696 1204 128423 64 8 32 38
148 94785 1138 97839 38 2 25 24
149 8773 568 38214 34 0 16 8
150 83209 1459 151101 32 2 48 35
151 93815 2158 272458 65 0 100 43
152 86687 1111 172494 52 0 46 43
153 34553 1421 108043 62 1 45 14
154 105547 2833 328107 65 3 129 41
155 103487 1955 250579 83 0 130 38
156 213688 2922 351067 95 3 136 45
157 71220 1002 158015 29 0 59 31
158 23517 1060 98866 18 0 25 13
159 56926 956 85439 33 0 32 28
160 91721 2186 229242 247 4 63 31
161 115168 3604 351619 139 4 95 40
162 111194 1035 84207 29 11 14 30
163 51009 1417 120445 118 0 36 16
164 135777 3261 324598 110 0 113 37
165 51513 1587 131069 67 4 47 30
166 74163 1424 204271 42 0 92 35
167 51633 1701 165543 65 1 70 32
168 75345 1249 141722 94 0 19 27
169 33416 946 116048 64 0 50 20
170 83305 1926 250047 81 0 41 18
171 98952 3352 299775 95 9 91 31
172 102372 1641 195838 67 1 111 31
173 37238 2035 173260 63 3 41 21
174 103772 2312 254488 83 10 120 39
175 123969 1369 104389 45 5 135 41
176 27142 1577 136084 30 0 27 13
177 135400 2201 199476 70 2 87 32
178 21399 961 92499 32 0 25 18
179 130115 1900 224330 83 1 131 39
180 24874 1254 135781 31 2 45 14
181 34988 1335 74408 67 4 29 7
182 45549 1597 81240 66 0 58 17
183 6023 207 14688 10 0 4 0
184 64466 1645 181633 70 2 47 30
185 54990 2429 271856 103 1 109 37
186 1644 151 7199 5 0 7 0
187 6179 474 46660 20 0 12 5
188 3926 141 17547 5 0 0 1
189 32755 1639 133368 36 1 37 16
190 34777 872 95227 34 0 37 32
191 73224 1318 152601 48 2 46 24
192 27114 1018 98146 40 0 15 17
193 20760 1383 79619 43 3 42 11
194 37636 1314 59194 31 6 7 24
195 65461 1335 139942 42 0 54 22
196 30080 1403 118612 46 2 54 12
197 24094 910 72880 33 0 14 19
198 69008 616 65475 18 2 16 13
199 54968 1407 99643 55 1 33 17
200 46090 771 71965 35 1 32 15
201 27507 766 77272 59 2 21 16
202 10672 473 49289 19 1 15 24
203 34029 1376 135131 66 0 38 15
204 46300 1232 108446 60 1 22 17
205 24760 1521 89746 36 3 28 18
206 18779 572 44296 25 0 10 20
207 21280 1059 77648 47 0 31 16
208 40662 1544 181528 54 0 32 16
209 28987 1230 134019 53 0 32 18
210 22827 1206 124064 40 1 43 22
211 18513 1205 92630 40 4 27 8
212 30594 1255 121848 39 0 37 17
213 24006 613 52915 14 0 20 18
214 27913 721 81872 45 0 32 16
215 42744 1109 58981 36 7 0 23
216 12934 740 53515 28 2 5 22
217 22574 1126 60812 44 0 26 13
218 41385 728 56375 30 7 10 13
219 18653 689 65490 22 3 27 16
220 18472 592 80949 17 0 11 16
221 30976 995 76302 31 0 29 20
222 63339 1613 104011 55 6 25 22
223 25568 2048 98104 54 2 55 17
224 33747 705 67989 21 0 23 18
225 4154 301 30989 14 0 5 17
226 19474 1803 135458 81 3 43 12
227 35130 799 73504 35 0 23 7
228 39067 861 63123 43 1 34 17
229 13310 1186 61254 46 1 36 14
230 65892 1451 74914 30 0 35 23
231 4143 628 31774 23 1 0 17
232 28579 1161 81437 38 0 37 14
233 51776 1463 87186 54 0 28 15
234 21152 742 50090 20 0 16 17
235 38084 979 65745 53 0 26 21
236 27717 675 56653 45 0 38 18
237 32928 1241 158399 39 0 23 18
238 11342 676 46455 20 0 22 17
239 19499 1049 73624 24 0 30 17
240 16380 620 38395 31 0 16 16
241 36874 1081 91899 35 0 18 15
242 48259 1688 139526 151 0 28 21
243 16734 736 52164 52 0 32 16
244 28207 617 51567 30 2 21 14
245 30143 812 70551 31 0 23 15
246 41369 1051 84856 29 1 29 17
247 45833 1656 102538 57 1 50 15
248 29156 705 86678 40 0 12 15
249 35944 945 85709 44 0 21 10
250 36278 554 34662 25 0 18 6
251 45588 1597 150580 77 0 27 22
252 45097 982 99611 35 0 41 21
253 3895 222 19349 11 0 13 1
254 28394 1212 99373 63 1 12 18
255 18632 1143 86230 44 0 21 17
256 2325 435 30837 19 0 8 4
257 25139 532 31706 13 0 26 10
258 27975 882 89806 42 0 27 16
259 14483 608 62088 38 1 13 16
260 13127 459 40151 29 0 16 9
261 5839 578 27634 20 0 2 16
262 24069 826 76990 27 0 42 17
263 3738 509 37460 20 0 5 7
264 18625 717 54157 19 0 37 15
265 36341 637 49862 37 0 17 14
266 24548 857 84337 26 0 38 14
267 21792 830 64175 42 0 37 18
268 26263 652 59382 49 0 29 12
269 23686 707 119308 30 0 32 16
270 49303 954 76702 49 0 35 21
271 25659 1461 103425 67 1 17 19
272 28904 672 70344 28 0 20 16
273 2781 778 43410 19 0 7 1
274 29236 1141 104838 49 1 46 16
275 19546 680 62215 27 0 24 10
276 22818 1090 69304 30 6 40 19
277 32689 616 53117 22 3 3 12
278 5752 285 19764 12 1 10 2
279 22197 1145 86680 31 2 37 14
280 20055 733 84105 20 0 17 17
281 25272 888 77945 20 0 28 19
282 82206 849 89113 39 0 19 14
283 32073 1182 91005 29 3 29 11
284 5444 528 40248 16 1 8 4
285 20154 642 64187 27 0 10 16
286 36944 947 50857 21 0 15 20
287 8019 819 56613 19 1 15 12
288 30884 757 62792 35 0 28 15
289 19540 894 72535 14 0 17 16
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) pageviews time logins shared blogged
-278.3956 -10.6874 0.1734 71.4210 2235.7235 419.8750
reviewed
810.2645
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-51260 -13652 -4339 9173 150970
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -278.39564 3712.36512 -0.075 0.940275
pageviews -10.68741 4.53671 -2.356 0.019170 *
time 0.17344 0.04568 3.797 0.000180 ***
logins 71.42100 63.22201 1.130 0.259569
shared 2235.72347 551.72201 4.052 6.56e-05 ***
blogged 419.87499 73.25188 5.732 2.55e-08 ***
reviewed 810.26447 224.28649 3.613 0.000359 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 23250 on 282 degrees of freedom
Multiple R-squared: 0.6762, Adjusted R-squared: 0.6693
F-statistic: 98.13 on 6 and 282 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.84948248 3.010350e-01 1.505175e-01
[2,] 0.79181726 4.163655e-01 2.081827e-01
[3,] 0.89432472 2.113506e-01 1.056753e-01
[4,] 0.86795538 2.640892e-01 1.320446e-01
[5,] 0.83899971 3.220006e-01 1.610003e-01
[6,] 0.84399278 3.120144e-01 1.560072e-01
[7,] 0.86701425 2.659715e-01 1.329858e-01
[8,] 0.81961114 3.607777e-01 1.803889e-01
[9,] 0.83346049 3.330790e-01 1.665395e-01
[10,] 0.85834944 2.833011e-01 1.416506e-01
[11,] 0.83230394 3.353921e-01 1.676961e-01
[12,] 0.77948923 4.410215e-01 2.205108e-01
[13,] 0.72672417 5.465517e-01 2.732758e-01
[14,] 0.71916891 5.616622e-01 2.808311e-01
[15,] 0.67198452 6.560310e-01 3.280155e-01
[16,] 0.60758936 7.848213e-01 3.924106e-01
[17,] 0.61201960 7.759608e-01 3.879804e-01
[18,] 0.60109296 7.978141e-01 3.989070e-01
[19,] 0.53855561 9.228888e-01 4.614444e-01
[20,] 0.52089156 9.582169e-01 4.791084e-01
[21,] 0.53681511 9.263698e-01 4.631849e-01
[22,] 0.49424070 9.884814e-01 5.057593e-01
[23,] 0.51223133 9.755373e-01 4.877687e-01
[24,] 0.48260904 9.652181e-01 5.173910e-01
[25,] 0.43600508 8.720102e-01 5.639949e-01
[26,] 0.37987031 7.597406e-01 6.201297e-01
[27,] 0.32744243 6.548849e-01 6.725576e-01
[28,] 0.32407875 6.481575e-01 6.759212e-01
[29,] 0.29764725 5.952945e-01 7.023527e-01
[30,] 0.25695451 5.139090e-01 7.430455e-01
[31,] 0.26098488 5.219698e-01 7.390151e-01
[32,] 0.22759494 4.551899e-01 7.724051e-01
[33,] 0.20704134 4.140827e-01 7.929587e-01
[34,] 0.18422358 3.684472e-01 8.157764e-01
[35,] 0.15439957 3.087991e-01 8.456004e-01
[36,] 0.23113538 4.622708e-01 7.688646e-01
[37,] 0.20973037 4.194607e-01 7.902696e-01
[38,] 0.17623723 3.524745e-01 8.237628e-01
[39,] 0.15071291 3.014258e-01 8.492871e-01
[40,] 0.12518674 2.503735e-01 8.748133e-01
[41,] 0.13849498 2.769900e-01 8.615050e-01
[42,] 0.13960607 2.792121e-01 8.603939e-01
[43,] 0.11924985 2.384997e-01 8.807501e-01
[44,] 0.09732787 1.946557e-01 9.026721e-01
[45,] 0.10700610 2.140122e-01 8.929939e-01
[46,] 0.09079193 1.815839e-01 9.092081e-01
[47,] 0.07399957 1.479991e-01 9.260004e-01
[48,] 0.06141552 1.228310e-01 9.385845e-01
[49,] 0.05967224 1.193445e-01 9.403278e-01
[50,] 0.04948595 9.897190e-02 9.505141e-01
[51,] 0.04044943 8.089886e-02 9.595506e-01
[52,] 0.03252109 6.504219e-02 9.674789e-01
[53,] 0.04947715 9.895430e-02 9.505228e-01
[54,] 0.07010536 1.402107e-01 9.298946e-01
[55,] 0.06099032 1.219806e-01 9.390097e-01
[56,] 0.07100836 1.420167e-01 9.289916e-01
[57,] 0.05897843 1.179569e-01 9.410216e-01
[58,] 0.06509637 1.301927e-01 9.349036e-01
[59,] 0.05341489 1.068298e-01 9.465851e-01
[60,] 0.10435347 2.087069e-01 8.956465e-01
[61,] 0.13804482 2.760896e-01 8.619552e-01
[62,] 0.11728542 2.345708e-01 8.827146e-01
[63,] 0.13991306 2.798261e-01 8.600869e-01
[64,] 0.12490919 2.498184e-01 8.750908e-01
[65,] 0.14537605 2.907521e-01 8.546240e-01
[66,] 0.14111564 2.822313e-01 8.588844e-01
[67,] 0.11972997 2.394599e-01 8.802700e-01
[68,] 0.10686930 2.137386e-01 8.931307e-01
[69,] 0.13010442 2.602088e-01 8.698956e-01
[70,] 0.13257275 2.651455e-01 8.674272e-01
[71,] 0.12481556 2.496311e-01 8.751844e-01
[72,] 0.11168141 2.233628e-01 8.883186e-01
[73,] 0.09992297 1.998459e-01 9.000770e-01
[74,] 0.09800610 1.960122e-01 9.019939e-01
[75,] 0.08341926 1.668385e-01 9.165807e-01
[76,] 0.95843511 8.312979e-02 4.156489e-02
[77,] 0.95536721 8.926557e-02 4.463279e-02
[78,] 0.94650285 1.069943e-01 5.349715e-02
[79,] 0.93670278 1.265944e-01 6.329722e-02
[80,] 0.92555464 1.488907e-01 7.444536e-02
[81,] 0.92413065 1.517387e-01 7.586935e-02
[82,] 0.91363333 1.727333e-01 8.636667e-02
[83,] 0.90996218 1.800756e-01 9.003782e-02
[84,] 0.90533914 1.893217e-01 9.466086e-02
[85,] 0.89937809 2.012438e-01 1.006219e-01
[86,] 0.89670310 2.065938e-01 1.032969e-01
[87,] 0.88477010 2.304598e-01 1.152299e-01
[88,] 0.88552437 2.289513e-01 1.144756e-01
[89,] 0.90203574 1.959285e-01 9.796426e-02
[90,] 0.92036593 1.592681e-01 7.963407e-02
[91,] 0.90638448 1.872310e-01 9.361552e-02
[92,] 0.89349749 2.130050e-01 1.065025e-01
[93,] 0.90177960 1.964408e-01 9.822040e-02
[94,] 0.91261825 1.747635e-01 8.738175e-02
[95,] 0.95152761 9.694479e-02 4.847239e-02
[96,] 0.97961653 4.076694e-02 2.038347e-02
[97,] 0.99207990 1.584021e-02 7.920103e-03
[98,] 0.99158138 1.683725e-02 8.418624e-03
[99,] 0.98997620 2.004759e-02 1.002380e-02
[100,] 0.99275045 1.449909e-02 7.249546e-03
[101,] 0.99098865 1.802270e-02 9.011349e-03
[102,] 0.99233850 1.532300e-02 7.661502e-03
[103,] 0.99299677 1.400646e-02 7.003232e-03
[104,] 0.99155081 1.689838e-02 8.449189e-03
[105,] 0.99285121 1.429758e-02 7.148792e-03
[106,] 0.99126113 1.747774e-02 8.738869e-03
[107,] 0.98907450 2.185100e-02 1.092550e-02
[108,] 0.98916709 2.166582e-02 1.083291e-02
[109,] 0.98766105 2.467790e-02 1.233895e-02
[110,] 1.00000000 6.775322e-09 3.387661e-09
[111,] 0.99999999 1.110723e-08 5.553614e-09
[112,] 0.99999999 1.678053e-08 8.390265e-09
[113,] 0.99999999 2.349884e-08 1.174942e-08
[114,] 0.99999998 3.285877e-08 1.642939e-08
[115,] 1.00000000 1.574214e-10 7.871069e-11
[116,] 1.00000000 7.003448e-11 3.501724e-11
[117,] 1.00000000 9.145779e-11 4.572889e-11
[118,] 1.00000000 1.558635e-10 7.793176e-11
[119,] 1.00000000 2.223486e-10 1.111743e-10
[120,] 1.00000000 3.301074e-10 1.650537e-10
[121,] 1.00000000 5.219081e-10 2.609540e-10
[122,] 1.00000000 8.226319e-10 4.113159e-10
[123,] 1.00000000 1.051608e-09 5.258038e-10
[124,] 1.00000000 1.607918e-09 8.039588e-10
[125,] 1.00000000 2.533686e-09 1.266843e-09
[126,] 1.00000000 1.562930e-09 7.814648e-10
[127,] 1.00000000 1.797782e-09 8.988910e-10
[128,] 1.00000000 2.900227e-09 1.450114e-09
[129,] 1.00000000 3.300296e-09 1.650148e-09
[130,] 1.00000000 1.994605e-09 9.973025e-10
[131,] 1.00000000 3.279412e-09 1.639706e-09
[132,] 1.00000000 5.051400e-09 2.525700e-09
[133,] 1.00000000 6.178056e-09 3.089028e-09
[134,] 1.00000000 9.269383e-09 4.634692e-09
[135,] 0.99999999 1.119697e-08 5.598484e-09
[136,] 0.99999999 1.479934e-08 7.399671e-09
[137,] 0.99999999 2.351270e-08 1.175635e-08
[138,] 0.99999998 3.231024e-08 1.615512e-08
[139,] 1.00000000 2.984762e-09 1.492381e-09
[140,] 1.00000000 4.774844e-09 2.387422e-09
[141,] 1.00000000 6.248456e-09 3.124228e-09
[142,] 1.00000000 8.376713e-09 4.188357e-09
[143,] 0.99999999 1.306021e-08 6.530106e-09
[144,] 0.99999999 2.104242e-08 1.052121e-08
[145,] 0.99999999 1.965883e-08 9.829416e-09
[146,] 0.99999999 2.840397e-08 1.420199e-08
[147,] 1.00000000 3.144051e-11 1.572025e-11
[148,] 1.00000000 5.483835e-11 2.741918e-11
[149,] 1.00000000 9.741911e-11 4.870956e-11
[150,] 1.00000000 1.387374e-10 6.936869e-11
[151,] 1.00000000 2.292367e-10 1.146184e-10
[152,] 1.00000000 4.032128e-10 2.016064e-10
[153,] 1.00000000 1.058287e-11 5.291437e-12
[154,] 1.00000000 1.898176e-11 9.490881e-12
[155,] 1.00000000 1.018756e-11 5.093778e-12
[156,] 1.00000000 1.610513e-11 8.052564e-12
[157,] 1.00000000 2.547983e-11 1.273992e-11
[158,] 1.00000000 2.542445e-11 1.271222e-11
[159,] 1.00000000 1.480859e-11 7.404295e-12
[160,] 1.00000000 1.840949e-11 9.204747e-12
[161,] 1.00000000 8.553456e-12 4.276728e-12
[162,] 1.00000000 1.589520e-11 7.947600e-12
[163,] 1.00000000 2.293582e-11 1.146791e-11
[164,] 1.00000000 3.159204e-11 1.579602e-11
[165,] 1.00000000 3.673766e-11 1.836883e-11
[166,] 1.00000000 3.646198e-11 1.823099e-11
[167,] 1.00000000 6.665253e-11 3.332626e-11
[168,] 1.00000000 2.568627e-13 1.284314e-13
[169,] 1.00000000 4.296936e-13 2.148468e-13
[170,] 1.00000000 1.332711e-14 6.663556e-15
[171,] 1.00000000 1.874474e-14 9.372370e-15
[172,] 1.00000000 3.777222e-14 1.888611e-14
[173,] 1.00000000 6.122139e-14 3.061070e-14
[174,] 1.00000000 1.285326e-13 6.426632e-14
[175,] 1.00000000 2.368695e-13 1.184348e-13
[176,] 1.00000000 6.106095e-14 3.053048e-14
[177,] 1.00000000 1.284749e-13 6.423744e-14
[178,] 1.00000000 2.418825e-13 1.209413e-13
[179,] 1.00000000 5.051198e-13 2.525599e-13
[180,] 1.00000000 9.606477e-13 4.803238e-13
[181,] 1.00000000 1.576248e-12 7.881238e-13
[182,] 1.00000000 1.104312e-12 5.521559e-13
[183,] 1.00000000 2.224826e-12 1.112413e-12
[184,] 1.00000000 3.514609e-12 1.757305e-12
[185,] 1.00000000 7.152932e-12 3.576466e-12
[186,] 1.00000000 5.306189e-12 2.653094e-12
[187,] 1.00000000 9.469291e-12 4.734646e-12
[188,] 1.00000000 1.873358e-11 9.366788e-12
[189,] 1.00000000 2.753961e-13 1.376980e-13
[190,] 1.00000000 2.324873e-13 1.162437e-13
[191,] 1.00000000 1.941714e-13 9.708569e-14
[192,] 1.00000000 4.141448e-13 2.070724e-13
[193,] 1.00000000 4.305446e-13 2.152723e-13
[194,] 1.00000000 9.082222e-13 4.541111e-13
[195,] 1.00000000 1.421728e-12 7.108642e-13
[196,] 1.00000000 2.357224e-12 1.178612e-12
[197,] 1.00000000 4.728244e-12 2.364122e-12
[198,] 1.00000000 8.416445e-12 4.208222e-12
[199,] 1.00000000 1.789892e-11 8.949462e-12
[200,] 1.00000000 3.146734e-11 1.573367e-11
[201,] 1.00000000 2.796900e-11 1.398450e-11
[202,] 1.00000000 4.590102e-11 2.295051e-11
[203,] 1.00000000 8.800381e-11 4.400190e-11
[204,] 1.00000000 1.842830e-10 9.214149e-11
[205,] 1.00000000 3.765850e-10 1.882925e-10
[206,] 1.00000000 6.611194e-10 3.305597e-10
[207,] 1.00000000 7.926798e-10 3.963399e-10
[208,] 1.00000000 1.597178e-09 7.985892e-10
[209,] 1.00000000 1.497823e-09 7.489113e-10
[210,] 1.00000000 2.681205e-09 1.340603e-09
[211,] 1.00000000 4.900663e-09 2.450331e-09
[212,] 1.00000000 9.821430e-09 4.910715e-09
[213,] 1.00000000 3.146493e-09 1.573246e-09
[214,] 1.00000000 3.490311e-09 1.745155e-09
[215,] 1.00000000 6.522027e-09 3.261014e-09
[216,] 1.00000000 9.228755e-09 4.614378e-09
[217,] 1.00000000 6.150428e-09 3.075214e-09
[218,] 1.00000000 8.651221e-09 4.325610e-09
[219,] 0.99999999 1.381731e-08 6.908656e-09
[220,] 0.99999999 1.238066e-08 6.190328e-09
[221,] 1.00000000 3.467680e-09 1.733840e-09
[222,] 1.00000000 5.128204e-09 2.564102e-09
[223,] 0.99999999 1.082968e-08 5.414839e-09
[224,] 1.00000000 8.952023e-09 4.476011e-09
[225,] 0.99999999 1.948569e-08 9.742844e-09
[226,] 0.99999998 3.681562e-08 1.840781e-08
[227,] 0.99999996 7.741960e-08 3.870980e-08
[228,] 0.99999994 1.274619e-07 6.373095e-08
[229,] 0.99999990 2.008534e-07 1.004267e-07
[230,] 0.99999983 3.401285e-07 1.700643e-07
[231,] 0.99999966 6.714194e-07 3.357097e-07
[232,] 0.99999940 1.205540e-06 6.027698e-07
[233,] 0.99999886 2.270129e-06 1.135065e-06
[234,] 0.99999854 2.921175e-06 1.460587e-06
[235,] 0.99999727 5.460007e-06 2.730004e-06
[236,] 0.99999473 1.054970e-05 5.274851e-06
[237,] 0.99999277 1.446243e-05 7.231217e-06
[238,] 0.99998839 2.322884e-05 1.161442e-05
[239,] 0.99997746 4.507605e-05 2.253803e-05
[240,] 0.99996491 7.017438e-05 3.508719e-05
[241,] 0.99997921 4.157462e-05 2.078731e-05
[242,] 0.99995945 8.110965e-05 4.055483e-05
[243,] 0.99993301 1.339778e-04 6.698888e-05
[244,] 0.99987261 2.547728e-04 1.273864e-04
[245,] 0.99978374 4.325275e-04 2.162638e-04
[246,] 0.99969866 6.026816e-04 3.013408e-04
[247,] 0.99949019 1.019611e-03 5.098056e-04
[248,] 0.99934986 1.300272e-03 6.501359e-04
[249,] 0.99887857 2.242854e-03 1.121427e-03
[250,] 0.99875561 2.488778e-03 1.244389e-03
[251,] 0.99787777 4.244456e-03 2.122228e-03
[252,] 0.99750046 4.999072e-03 2.499536e-03
[253,] 0.99561731 8.765389e-03 4.382695e-03
[254,] 0.99399648 1.200705e-02 6.003525e-03
[255,] 0.98973322 2.053355e-02 1.026678e-02
[256,] 0.98402073 3.195853e-02 1.597927e-02
[257,] 0.97394700 5.210599e-02 2.605300e-02
[258,] 0.96311699 7.376601e-02 3.688301e-02
[259,] 0.94548518 1.090296e-01 5.451482e-02
[260,] 0.94036852 1.192630e-01 5.963148e-02
[261,] 0.91795309 1.640938e-01 8.204691e-02
[262,] 0.97277644 5.444713e-02 2.722356e-02
[263,] 0.95237096 9.525808e-02 4.762904e-02
[264,] 0.93540463 1.291907e-01 6.459537e-02
[265,] 0.94463998 1.107200e-01 5.536002e-02
[266,] 0.90443029 1.911394e-01 9.556971e-02
[267,] 0.85090632 2.981874e-01 1.490937e-01
[268,] 0.85247846 2.950431e-01 1.475215e-01
[269,] 0.94906484 1.018703e-01 5.093516e-02
[270,] 0.90042306 1.991539e-01 9.957694e-02
> postscript(file="/var/wessaorg/rcomp/tmp/102hx1324509044.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/2m6561324509044.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/3c8fo1324509044.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/4aulj1324509045.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/53s1o1324509045.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> qqnorm(mysum$resid, main='Residual Normal Q-Q Plot')
> qqline(mysum$resid)
> grid()
> dev.off()
null device
1
> (myerror <- as.ts(mysum$resid))
Time Series:
Start = 1
End = 289
Frequency = 1
1 2 3 4 5 6
22954.54564 13386.53852 -17527.61333 13101.44254 -16588.17515 35023.05749
7 8 9 10 11 12
9173.23445 5007.35006 -9281.01268 54075.42922 -6986.53663 -25351.86787
13 14 15 16 17 18
-35799.83338 -12705.10718 -27872.05612 -40000.90310 -30619.02753 32414.33491
19 20 21 22 23 24
11036.89677 31505.61108 2524.77953 6912.93449 -33310.03682 -2613.65352
25 26 27 28 29 30
3540.56691 -19798.49271 12816.30860 -6685.72606 -24386.72310 -20429.88502
31 32 33 34 35 36
-14972.90224 20304.73341 -4245.98368 -26528.68346 4308.47059 -4721.91061
37 38 39 40 41 42
-6368.65247 -19599.82181 -19951.56435 1230.52205 -13943.19673 -7428.86684
43 44 45 46 47 48
-11147.27157 -16821.71696 -42981.51488 -20140.98439 -9126.59037 5656.73582
49 50 51 52 53 54
-7111.57958 21623.49724 11761.80199 -7480.79802 -3568.37161 23066.78688
55 56 57 58 59 60
-8221.30897 -682.29921 -979.87388 -18141.77078 -6292.42931 -8364.41419
61 62 63 64 65 66
7539.98436 30328.62695 37559.35498 -11623.30222 27413.78600 -5423.08436
67 68 69 70 71 72
-29541.84697 1550.84673 47693.59900 -37493.96743 -85.83245 25296.87518
73 74 75 76 77 78
-14176.43221 29550.05019 -21678.16673 810.13634 -11221.42876 -35784.08481
79 80 81 82 83 84
-23891.86158 -20896.19418 10617.26895 -15251.69442 -23417.35971 1599.70270
85 86 87 88 89 90
132613.07870 -19346.79251 -6145.85002 -6560.33531 -8719.03978 -21580.73672
91 92 93 94 95 96
-11421.41937 23316.78936 -20089.31286 18165.55335 -20980.06396 12439.71474
97 98 99 100 101 102
-26612.77280 33201.44974 -36537.21252 1162.93661 -7970.28216 -29775.26493
103 104 105 106 107 108
29633.36210 49105.63041 55646.98610 54844.87927 -17976.46231 -11489.67373
109 110 111 112 113 114
38952.17023 -2309.02191 30032.46240 -25075.90284 12158.51828 -28030.36180
115 116 117 118 119 120
9231.69167 -947.32947 -22922.92434 -14253.87419 150969.60813 -729.24788
121 122 123 124 125 126
-8589.89797 -10697.24222 -10640.38763 73062.45319 36663.26374 -14519.47004
127 128 129 130 131 132
5945.38980 11722.94613 -14786.14169 -4916.06948 -6271.24646 -18523.31821
133 134 135 136 137 138
-8921.93940 -6916.14734 -29060.73164 -19961.32032 -1204.32800 19947.74299
139 140 141 142 143 144
30386.68001 -1652.13047 -8608.40959 15826.97273 -7870.56579 -15824.41021
145 146 147 148 149 150
17690.49105 -3389.70406 16886.12790 53128.22833 -7134.25644 17603.84951
151 152 153 154 155 156
-11568.34180 11052.98556 -5622.20856 -19536.87823 -10101.79230 77250.40634
157 158 159 160 161 162
2839.71446 -4338.73162 14123.09141 -2550.86772 1811.42967 51079.01578
163 164 165 166 167 168
9034.58690 29328.19999 -11749.88476 -15755.04400 -20818.18137 27823.98268
169 170 171 172 173 174
-18092.15839 23215.54592 -7170.41349 7477.89025 -16221.31189 -25648.09778
175 176 177 178 179 180
16477.15369 -3340.07200 52676.62708 -11461.73643 17025.04907 -21918.43501
181 182 183 184 185 186
5052.75736 5964.24089 3572.55621 -2689.42998 -51259.65616 -1008.59910
187 188 189 190 191 192
-7087.52246 1500.67549 -5887.19446 -16033.09227 14461.69518 -1679.29891
193 194 195 196 197 198
-14315.62554 3677.43402 12237.40452 -15371.95638 -2172.22128 41505.67723
199 200 201 202 203 204
19207.68038 11801.66566 -7896.66118 -21880.12447 -7246.26593 11404.15060
205 206 207 208 209 210
-9890.80066 -4701.47702 -9927.71427 -4298.63690 -12638.83720 -26495.71659
211 212 213 214 215 216
-14014.09676 -8942.98367 -2323.72498 -7916.68224 7788.01869 -14556.75986
217 218 219 220 221 222
-253.26471 7141.51138 -17642.56944 -7759.11381 -1940.83582 17151.77925
223 224 225 226 227 228
-14476.32772 4026.57583 -14599.05671 -24741.52426 13370.70969 4242.40298
229 230 231 232 233 234
-16340.25874 33210.72721 -12030.55970 -2451.67853 24801.71091 -1247.85566
235 236 237 238 239 240
5705.22411 -8370.20885 -8029.86732 -13652.03584 -9865.38996 -5270.75600
241 242 243 244 245 246
10555.46888 2822.33330 -14282.89276 -639.18412 2838.36956 7905.00251
247 248 249 250 251 252
6571.56030 1886.65688 11394.45670 22260.73114 2156.38632 1864.16876
253 254 255 256 257 258
-3864.07649 -1967.79371 -9563.64834 -6052.87404 5656.27809 -5196.41030
259 260 261 262 263 264
-12881.26047 -4734.28788 -7730.40554 -13515.22842 -6240.25962 -11872.82824
265 266 267 268 269 270
13655.27162 -9797.43601 -13309.10582 -2188.54726 -17714.72682 11263.51400
271 272 273 274 275 276
-5939.71436 802.65750 -1261.03174 -14487.75180 -3806.47724 -25021.11006
277 278 279 280 281 282
11077.22647 -3263.52972 -13885.49992 -8760.32989 -7057.57248 53995.88803
283 284 285 286 287 288
-667.26316 -5593.60846 -2929.95068 14519.69785 -12382.35114 1952.19034
289
-4309.21710
> postscript(file="/var/wessaorg/rcomp/tmp/668od1324509045.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 22954.54564 NA
1 13386.53852 22954.54564
2 -17527.61333 13386.53852
3 13101.44254 -17527.61333
4 -16588.17515 13101.44254
5 35023.05749 -16588.17515
6 9173.23445 35023.05749
7 5007.35006 9173.23445
8 -9281.01268 5007.35006
9 54075.42922 -9281.01268
10 -6986.53663 54075.42922
11 -25351.86787 -6986.53663
12 -35799.83338 -25351.86787
13 -12705.10718 -35799.83338
14 -27872.05612 -12705.10718
15 -40000.90310 -27872.05612
16 -30619.02753 -40000.90310
17 32414.33491 -30619.02753
18 11036.89677 32414.33491
19 31505.61108 11036.89677
20 2524.77953 31505.61108
21 6912.93449 2524.77953
22 -33310.03682 6912.93449
23 -2613.65352 -33310.03682
24 3540.56691 -2613.65352
25 -19798.49271 3540.56691
26 12816.30860 -19798.49271
27 -6685.72606 12816.30860
28 -24386.72310 -6685.72606
29 -20429.88502 -24386.72310
30 -14972.90224 -20429.88502
31 20304.73341 -14972.90224
32 -4245.98368 20304.73341
33 -26528.68346 -4245.98368
34 4308.47059 -26528.68346
35 -4721.91061 4308.47059
36 -6368.65247 -4721.91061
37 -19599.82181 -6368.65247
38 -19951.56435 -19599.82181
39 1230.52205 -19951.56435
40 -13943.19673 1230.52205
41 -7428.86684 -13943.19673
42 -11147.27157 -7428.86684
43 -16821.71696 -11147.27157
44 -42981.51488 -16821.71696
45 -20140.98439 -42981.51488
46 -9126.59037 -20140.98439
47 5656.73582 -9126.59037
48 -7111.57958 5656.73582
49 21623.49724 -7111.57958
50 11761.80199 21623.49724
51 -7480.79802 11761.80199
52 -3568.37161 -7480.79802
53 23066.78688 -3568.37161
54 -8221.30897 23066.78688
55 -682.29921 -8221.30897
56 -979.87388 -682.29921
57 -18141.77078 -979.87388
58 -6292.42931 -18141.77078
59 -8364.41419 -6292.42931
60 7539.98436 -8364.41419
61 30328.62695 7539.98436
62 37559.35498 30328.62695
63 -11623.30222 37559.35498
64 27413.78600 -11623.30222
65 -5423.08436 27413.78600
66 -29541.84697 -5423.08436
67 1550.84673 -29541.84697
68 47693.59900 1550.84673
69 -37493.96743 47693.59900
70 -85.83245 -37493.96743
71 25296.87518 -85.83245
72 -14176.43221 25296.87518
73 29550.05019 -14176.43221
74 -21678.16673 29550.05019
75 810.13634 -21678.16673
76 -11221.42876 810.13634
77 -35784.08481 -11221.42876
78 -23891.86158 -35784.08481
79 -20896.19418 -23891.86158
80 10617.26895 -20896.19418
81 -15251.69442 10617.26895
82 -23417.35971 -15251.69442
83 1599.70270 -23417.35971
84 132613.07870 1599.70270
85 -19346.79251 132613.07870
86 -6145.85002 -19346.79251
87 -6560.33531 -6145.85002
88 -8719.03978 -6560.33531
89 -21580.73672 -8719.03978
90 -11421.41937 -21580.73672
91 23316.78936 -11421.41937
92 -20089.31286 23316.78936
93 18165.55335 -20089.31286
94 -20980.06396 18165.55335
95 12439.71474 -20980.06396
96 -26612.77280 12439.71474
97 33201.44974 -26612.77280
98 -36537.21252 33201.44974
99 1162.93661 -36537.21252
100 -7970.28216 1162.93661
101 -29775.26493 -7970.28216
102 29633.36210 -29775.26493
103 49105.63041 29633.36210
104 55646.98610 49105.63041
105 54844.87927 55646.98610
106 -17976.46231 54844.87927
107 -11489.67373 -17976.46231
108 38952.17023 -11489.67373
109 -2309.02191 38952.17023
110 30032.46240 -2309.02191
111 -25075.90284 30032.46240
112 12158.51828 -25075.90284
113 -28030.36180 12158.51828
114 9231.69167 -28030.36180
115 -947.32947 9231.69167
116 -22922.92434 -947.32947
117 -14253.87419 -22922.92434
118 150969.60813 -14253.87419
119 -729.24788 150969.60813
120 -8589.89797 -729.24788
121 -10697.24222 -8589.89797
122 -10640.38763 -10697.24222
123 73062.45319 -10640.38763
124 36663.26374 73062.45319
125 -14519.47004 36663.26374
126 5945.38980 -14519.47004
127 11722.94613 5945.38980
128 -14786.14169 11722.94613
129 -4916.06948 -14786.14169
130 -6271.24646 -4916.06948
131 -18523.31821 -6271.24646
132 -8921.93940 -18523.31821
133 -6916.14734 -8921.93940
134 -29060.73164 -6916.14734
135 -19961.32032 -29060.73164
136 -1204.32800 -19961.32032
137 19947.74299 -1204.32800
138 30386.68001 19947.74299
139 -1652.13047 30386.68001
140 -8608.40959 -1652.13047
141 15826.97273 -8608.40959
142 -7870.56579 15826.97273
143 -15824.41021 -7870.56579
144 17690.49105 -15824.41021
145 -3389.70406 17690.49105
146 16886.12790 -3389.70406
147 53128.22833 16886.12790
148 -7134.25644 53128.22833
149 17603.84951 -7134.25644
150 -11568.34180 17603.84951
151 11052.98556 -11568.34180
152 -5622.20856 11052.98556
153 -19536.87823 -5622.20856
154 -10101.79230 -19536.87823
155 77250.40634 -10101.79230
156 2839.71446 77250.40634
157 -4338.73162 2839.71446
158 14123.09141 -4338.73162
159 -2550.86772 14123.09141
160 1811.42967 -2550.86772
161 51079.01578 1811.42967
162 9034.58690 51079.01578
163 29328.19999 9034.58690
164 -11749.88476 29328.19999
165 -15755.04400 -11749.88476
166 -20818.18137 -15755.04400
167 27823.98268 -20818.18137
168 -18092.15839 27823.98268
169 23215.54592 -18092.15839
170 -7170.41349 23215.54592
171 7477.89025 -7170.41349
172 -16221.31189 7477.89025
173 -25648.09778 -16221.31189
174 16477.15369 -25648.09778
175 -3340.07200 16477.15369
176 52676.62708 -3340.07200
177 -11461.73643 52676.62708
178 17025.04907 -11461.73643
179 -21918.43501 17025.04907
180 5052.75736 -21918.43501
181 5964.24089 5052.75736
182 3572.55621 5964.24089
183 -2689.42998 3572.55621
184 -51259.65616 -2689.42998
185 -1008.59910 -51259.65616
186 -7087.52246 -1008.59910
187 1500.67549 -7087.52246
188 -5887.19446 1500.67549
189 -16033.09227 -5887.19446
190 14461.69518 -16033.09227
191 -1679.29891 14461.69518
192 -14315.62554 -1679.29891
193 3677.43402 -14315.62554
194 12237.40452 3677.43402
195 -15371.95638 12237.40452
196 -2172.22128 -15371.95638
197 41505.67723 -2172.22128
198 19207.68038 41505.67723
199 11801.66566 19207.68038
200 -7896.66118 11801.66566
201 -21880.12447 -7896.66118
202 -7246.26593 -21880.12447
203 11404.15060 -7246.26593
204 -9890.80066 11404.15060
205 -4701.47702 -9890.80066
206 -9927.71427 -4701.47702
207 -4298.63690 -9927.71427
208 -12638.83720 -4298.63690
209 -26495.71659 -12638.83720
210 -14014.09676 -26495.71659
211 -8942.98367 -14014.09676
212 -2323.72498 -8942.98367
213 -7916.68224 -2323.72498
214 7788.01869 -7916.68224
215 -14556.75986 7788.01869
216 -253.26471 -14556.75986
217 7141.51138 -253.26471
218 -17642.56944 7141.51138
219 -7759.11381 -17642.56944
220 -1940.83582 -7759.11381
221 17151.77925 -1940.83582
222 -14476.32772 17151.77925
223 4026.57583 -14476.32772
224 -14599.05671 4026.57583
225 -24741.52426 -14599.05671
226 13370.70969 -24741.52426
227 4242.40298 13370.70969
228 -16340.25874 4242.40298
229 33210.72721 -16340.25874
230 -12030.55970 33210.72721
231 -2451.67853 -12030.55970
232 24801.71091 -2451.67853
233 -1247.85566 24801.71091
234 5705.22411 -1247.85566
235 -8370.20885 5705.22411
236 -8029.86732 -8370.20885
237 -13652.03584 -8029.86732
238 -9865.38996 -13652.03584
239 -5270.75600 -9865.38996
240 10555.46888 -5270.75600
241 2822.33330 10555.46888
242 -14282.89276 2822.33330
243 -639.18412 -14282.89276
244 2838.36956 -639.18412
245 7905.00251 2838.36956
246 6571.56030 7905.00251
247 1886.65688 6571.56030
248 11394.45670 1886.65688
249 22260.73114 11394.45670
250 2156.38632 22260.73114
251 1864.16876 2156.38632
252 -3864.07649 1864.16876
253 -1967.79371 -3864.07649
254 -9563.64834 -1967.79371
255 -6052.87404 -9563.64834
256 5656.27809 -6052.87404
257 -5196.41030 5656.27809
258 -12881.26047 -5196.41030
259 -4734.28788 -12881.26047
260 -7730.40554 -4734.28788
261 -13515.22842 -7730.40554
262 -6240.25962 -13515.22842
263 -11872.82824 -6240.25962
264 13655.27162 -11872.82824
265 -9797.43601 13655.27162
266 -13309.10582 -9797.43601
267 -2188.54726 -13309.10582
268 -17714.72682 -2188.54726
269 11263.51400 -17714.72682
270 -5939.71436 11263.51400
271 802.65750 -5939.71436
272 -1261.03174 802.65750
273 -14487.75180 -1261.03174
274 -3806.47724 -14487.75180
275 -25021.11006 -3806.47724
276 11077.22647 -25021.11006
277 -3263.52972 11077.22647
278 -13885.49992 -3263.52972
279 -8760.32989 -13885.49992
280 -7057.57248 -8760.32989
281 53995.88803 -7057.57248
282 -667.26316 53995.88803
283 -5593.60846 -667.26316
284 -2929.95068 -5593.60846
285 14519.69785 -2929.95068
286 -12382.35114 14519.69785
287 1952.19034 -12382.35114
288 -4309.21710 1952.19034
289 NA -4309.21710
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] 13386.53852 22954.54564
[2,] -17527.61333 13386.53852
[3,] 13101.44254 -17527.61333
[4,] -16588.17515 13101.44254
[5,] 35023.05749 -16588.17515
[6,] 9173.23445 35023.05749
[7,] 5007.35006 9173.23445
[8,] -9281.01268 5007.35006
[9,] 54075.42922 -9281.01268
[10,] -6986.53663 54075.42922
[11,] -25351.86787 -6986.53663
[12,] -35799.83338 -25351.86787
[13,] -12705.10718 -35799.83338
[14,] -27872.05612 -12705.10718
[15,] -40000.90310 -27872.05612
[16,] -30619.02753 -40000.90310
[17,] 32414.33491 -30619.02753
[18,] 11036.89677 32414.33491
[19,] 31505.61108 11036.89677
[20,] 2524.77953 31505.61108
[21,] 6912.93449 2524.77953
[22,] -33310.03682 6912.93449
[23,] -2613.65352 -33310.03682
[24,] 3540.56691 -2613.65352
[25,] -19798.49271 3540.56691
[26,] 12816.30860 -19798.49271
[27,] -6685.72606 12816.30860
[28,] -24386.72310 -6685.72606
[29,] -20429.88502 -24386.72310
[30,] -14972.90224 -20429.88502
[31,] 20304.73341 -14972.90224
[32,] -4245.98368 20304.73341
[33,] -26528.68346 -4245.98368
[34,] 4308.47059 -26528.68346
[35,] -4721.91061 4308.47059
[36,] -6368.65247 -4721.91061
[37,] -19599.82181 -6368.65247
[38,] -19951.56435 -19599.82181
[39,] 1230.52205 -19951.56435
[40,] -13943.19673 1230.52205
[41,] -7428.86684 -13943.19673
[42,] -11147.27157 -7428.86684
[43,] -16821.71696 -11147.27157
[44,] -42981.51488 -16821.71696
[45,] -20140.98439 -42981.51488
[46,] -9126.59037 -20140.98439
[47,] 5656.73582 -9126.59037
[48,] -7111.57958 5656.73582
[49,] 21623.49724 -7111.57958
[50,] 11761.80199 21623.49724
[51,] -7480.79802 11761.80199
[52,] -3568.37161 -7480.79802
[53,] 23066.78688 -3568.37161
[54,] -8221.30897 23066.78688
[55,] -682.29921 -8221.30897
[56,] -979.87388 -682.29921
[57,] -18141.77078 -979.87388
[58,] -6292.42931 -18141.77078
[59,] -8364.41419 -6292.42931
[60,] 7539.98436 -8364.41419
[61,] 30328.62695 7539.98436
[62,] 37559.35498 30328.62695
[63,] -11623.30222 37559.35498
[64,] 27413.78600 -11623.30222
[65,] -5423.08436 27413.78600
[66,] -29541.84697 -5423.08436
[67,] 1550.84673 -29541.84697
[68,] 47693.59900 1550.84673
[69,] -37493.96743 47693.59900
[70,] -85.83245 -37493.96743
[71,] 25296.87518 -85.83245
[72,] -14176.43221 25296.87518
[73,] 29550.05019 -14176.43221
[74,] -21678.16673 29550.05019
[75,] 810.13634 -21678.16673
[76,] -11221.42876 810.13634
[77,] -35784.08481 -11221.42876
[78,] -23891.86158 -35784.08481
[79,] -20896.19418 -23891.86158
[80,] 10617.26895 -20896.19418
[81,] -15251.69442 10617.26895
[82,] -23417.35971 -15251.69442
[83,] 1599.70270 -23417.35971
[84,] 132613.07870 1599.70270
[85,] -19346.79251 132613.07870
[86,] -6145.85002 -19346.79251
[87,] -6560.33531 -6145.85002
[88,] -8719.03978 -6560.33531
[89,] -21580.73672 -8719.03978
[90,] -11421.41937 -21580.73672
[91,] 23316.78936 -11421.41937
[92,] -20089.31286 23316.78936
[93,] 18165.55335 -20089.31286
[94,] -20980.06396 18165.55335
[95,] 12439.71474 -20980.06396
[96,] -26612.77280 12439.71474
[97,] 33201.44974 -26612.77280
[98,] -36537.21252 33201.44974
[99,] 1162.93661 -36537.21252
[100,] -7970.28216 1162.93661
[101,] -29775.26493 -7970.28216
[102,] 29633.36210 -29775.26493
[103,] 49105.63041 29633.36210
[104,] 55646.98610 49105.63041
[105,] 54844.87927 55646.98610
[106,] -17976.46231 54844.87927
[107,] -11489.67373 -17976.46231
[108,] 38952.17023 -11489.67373
[109,] -2309.02191 38952.17023
[110,] 30032.46240 -2309.02191
[111,] -25075.90284 30032.46240
[112,] 12158.51828 -25075.90284
[113,] -28030.36180 12158.51828
[114,] 9231.69167 -28030.36180
[115,] -947.32947 9231.69167
[116,] -22922.92434 -947.32947
[117,] -14253.87419 -22922.92434
[118,] 150969.60813 -14253.87419
[119,] -729.24788 150969.60813
[120,] -8589.89797 -729.24788
[121,] -10697.24222 -8589.89797
[122,] -10640.38763 -10697.24222
[123,] 73062.45319 -10640.38763
[124,] 36663.26374 73062.45319
[125,] -14519.47004 36663.26374
[126,] 5945.38980 -14519.47004
[127,] 11722.94613 5945.38980
[128,] -14786.14169 11722.94613
[129,] -4916.06948 -14786.14169
[130,] -6271.24646 -4916.06948
[131,] -18523.31821 -6271.24646
[132,] -8921.93940 -18523.31821
[133,] -6916.14734 -8921.93940
[134,] -29060.73164 -6916.14734
[135,] -19961.32032 -29060.73164
[136,] -1204.32800 -19961.32032
[137,] 19947.74299 -1204.32800
[138,] 30386.68001 19947.74299
[139,] -1652.13047 30386.68001
[140,] -8608.40959 -1652.13047
[141,] 15826.97273 -8608.40959
[142,] -7870.56579 15826.97273
[143,] -15824.41021 -7870.56579
[144,] 17690.49105 -15824.41021
[145,] -3389.70406 17690.49105
[146,] 16886.12790 -3389.70406
[147,] 53128.22833 16886.12790
[148,] -7134.25644 53128.22833
[149,] 17603.84951 -7134.25644
[150,] -11568.34180 17603.84951
[151,] 11052.98556 -11568.34180
[152,] -5622.20856 11052.98556
[153,] -19536.87823 -5622.20856
[154,] -10101.79230 -19536.87823
[155,] 77250.40634 -10101.79230
[156,] 2839.71446 77250.40634
[157,] -4338.73162 2839.71446
[158,] 14123.09141 -4338.73162
[159,] -2550.86772 14123.09141
[160,] 1811.42967 -2550.86772
[161,] 51079.01578 1811.42967
[162,] 9034.58690 51079.01578
[163,] 29328.19999 9034.58690
[164,] -11749.88476 29328.19999
[165,] -15755.04400 -11749.88476
[166,] -20818.18137 -15755.04400
[167,] 27823.98268 -20818.18137
[168,] -18092.15839 27823.98268
[169,] 23215.54592 -18092.15839
[170,] -7170.41349 23215.54592
[171,] 7477.89025 -7170.41349
[172,] -16221.31189 7477.89025
[173,] -25648.09778 -16221.31189
[174,] 16477.15369 -25648.09778
[175,] -3340.07200 16477.15369
[176,] 52676.62708 -3340.07200
[177,] -11461.73643 52676.62708
[178,] 17025.04907 -11461.73643
[179,] -21918.43501 17025.04907
[180,] 5052.75736 -21918.43501
[181,] 5964.24089 5052.75736
[182,] 3572.55621 5964.24089
[183,] -2689.42998 3572.55621
[184,] -51259.65616 -2689.42998
[185,] -1008.59910 -51259.65616
[186,] -7087.52246 -1008.59910
[187,] 1500.67549 -7087.52246
[188,] -5887.19446 1500.67549
[189,] -16033.09227 -5887.19446
[190,] 14461.69518 -16033.09227
[191,] -1679.29891 14461.69518
[192,] -14315.62554 -1679.29891
[193,] 3677.43402 -14315.62554
[194,] 12237.40452 3677.43402
[195,] -15371.95638 12237.40452
[196,] -2172.22128 -15371.95638
[197,] 41505.67723 -2172.22128
[198,] 19207.68038 41505.67723
[199,] 11801.66566 19207.68038
[200,] -7896.66118 11801.66566
[201,] -21880.12447 -7896.66118
[202,] -7246.26593 -21880.12447
[203,] 11404.15060 -7246.26593
[204,] -9890.80066 11404.15060
[205,] -4701.47702 -9890.80066
[206,] -9927.71427 -4701.47702
[207,] -4298.63690 -9927.71427
[208,] -12638.83720 -4298.63690
[209,] -26495.71659 -12638.83720
[210,] -14014.09676 -26495.71659
[211,] -8942.98367 -14014.09676
[212,] -2323.72498 -8942.98367
[213,] -7916.68224 -2323.72498
[214,] 7788.01869 -7916.68224
[215,] -14556.75986 7788.01869
[216,] -253.26471 -14556.75986
[217,] 7141.51138 -253.26471
[218,] -17642.56944 7141.51138
[219,] -7759.11381 -17642.56944
[220,] -1940.83582 -7759.11381
[221,] 17151.77925 -1940.83582
[222,] -14476.32772 17151.77925
[223,] 4026.57583 -14476.32772
[224,] -14599.05671 4026.57583
[225,] -24741.52426 -14599.05671
[226,] 13370.70969 -24741.52426
[227,] 4242.40298 13370.70969
[228,] -16340.25874 4242.40298
[229,] 33210.72721 -16340.25874
[230,] -12030.55970 33210.72721
[231,] -2451.67853 -12030.55970
[232,] 24801.71091 -2451.67853
[233,] -1247.85566 24801.71091
[234,] 5705.22411 -1247.85566
[235,] -8370.20885 5705.22411
[236,] -8029.86732 -8370.20885
[237,] -13652.03584 -8029.86732
[238,] -9865.38996 -13652.03584
[239,] -5270.75600 -9865.38996
[240,] 10555.46888 -5270.75600
[241,] 2822.33330 10555.46888
[242,] -14282.89276 2822.33330
[243,] -639.18412 -14282.89276
[244,] 2838.36956 -639.18412
[245,] 7905.00251 2838.36956
[246,] 6571.56030 7905.00251
[247,] 1886.65688 6571.56030
[248,] 11394.45670 1886.65688
[249,] 22260.73114 11394.45670
[250,] 2156.38632 22260.73114
[251,] 1864.16876 2156.38632
[252,] -3864.07649 1864.16876
[253,] -1967.79371 -3864.07649
[254,] -9563.64834 -1967.79371
[255,] -6052.87404 -9563.64834
[256,] 5656.27809 -6052.87404
[257,] -5196.41030 5656.27809
[258,] -12881.26047 -5196.41030
[259,] -4734.28788 -12881.26047
[260,] -7730.40554 -4734.28788
[261,] -13515.22842 -7730.40554
[262,] -6240.25962 -13515.22842
[263,] -11872.82824 -6240.25962
[264,] 13655.27162 -11872.82824
[265,] -9797.43601 13655.27162
[266,] -13309.10582 -9797.43601
[267,] -2188.54726 -13309.10582
[268,] -17714.72682 -2188.54726
[269,] 11263.51400 -17714.72682
[270,] -5939.71436 11263.51400
[271,] 802.65750 -5939.71436
[272,] -1261.03174 802.65750
[273,] -14487.75180 -1261.03174
[274,] -3806.47724 -14487.75180
[275,] -25021.11006 -3806.47724
[276,] 11077.22647 -25021.11006
[277,] -3263.52972 11077.22647
[278,] -13885.49992 -3263.52972
[279,] -8760.32989 -13885.49992
[280,] -7057.57248 -8760.32989
[281,] 53995.88803 -7057.57248
[282,] -667.26316 53995.88803
[283,] -5593.60846 -667.26316
[284,] -2929.95068 -5593.60846
[285,] 14519.69785 -2929.95068
[286,] -12382.35114 14519.69785
[287,] 1952.19034 -12382.35114
[288,] -4309.21710 1952.19034
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 13386.53852 22954.54564
2 -17527.61333 13386.53852
3 13101.44254 -17527.61333
4 -16588.17515 13101.44254
5 35023.05749 -16588.17515
6 9173.23445 35023.05749
7 5007.35006 9173.23445
8 -9281.01268 5007.35006
9 54075.42922 -9281.01268
10 -6986.53663 54075.42922
11 -25351.86787 -6986.53663
12 -35799.83338 -25351.86787
13 -12705.10718 -35799.83338
14 -27872.05612 -12705.10718
15 -40000.90310 -27872.05612
16 -30619.02753 -40000.90310
17 32414.33491 -30619.02753
18 11036.89677 32414.33491
19 31505.61108 11036.89677
20 2524.77953 31505.61108
21 6912.93449 2524.77953
22 -33310.03682 6912.93449
23 -2613.65352 -33310.03682
24 3540.56691 -2613.65352
25 -19798.49271 3540.56691
26 12816.30860 -19798.49271
27 -6685.72606 12816.30860
28 -24386.72310 -6685.72606
29 -20429.88502 -24386.72310
30 -14972.90224 -20429.88502
31 20304.73341 -14972.90224
32 -4245.98368 20304.73341
33 -26528.68346 -4245.98368
34 4308.47059 -26528.68346
35 -4721.91061 4308.47059
36 -6368.65247 -4721.91061
37 -19599.82181 -6368.65247
38 -19951.56435 -19599.82181
39 1230.52205 -19951.56435
40 -13943.19673 1230.52205
41 -7428.86684 -13943.19673
42 -11147.27157 -7428.86684
43 -16821.71696 -11147.27157
44 -42981.51488 -16821.71696
45 -20140.98439 -42981.51488
46 -9126.59037 -20140.98439
47 5656.73582 -9126.59037
48 -7111.57958 5656.73582
49 21623.49724 -7111.57958
50 11761.80199 21623.49724
51 -7480.79802 11761.80199
52 -3568.37161 -7480.79802
53 23066.78688 -3568.37161
54 -8221.30897 23066.78688
55 -682.29921 -8221.30897
56 -979.87388 -682.29921
57 -18141.77078 -979.87388
58 -6292.42931 -18141.77078
59 -8364.41419 -6292.42931
60 7539.98436 -8364.41419
61 30328.62695 7539.98436
62 37559.35498 30328.62695
63 -11623.30222 37559.35498
64 27413.78600 -11623.30222
65 -5423.08436 27413.78600
66 -29541.84697 -5423.08436
67 1550.84673 -29541.84697
68 47693.59900 1550.84673
69 -37493.96743 47693.59900
70 -85.83245 -37493.96743
71 25296.87518 -85.83245
72 -14176.43221 25296.87518
73 29550.05019 -14176.43221
74 -21678.16673 29550.05019
75 810.13634 -21678.16673
76 -11221.42876 810.13634
77 -35784.08481 -11221.42876
78 -23891.86158 -35784.08481
79 -20896.19418 -23891.86158
80 10617.26895 -20896.19418
81 -15251.69442 10617.26895
82 -23417.35971 -15251.69442
83 1599.70270 -23417.35971
84 132613.07870 1599.70270
85 -19346.79251 132613.07870
86 -6145.85002 -19346.79251
87 -6560.33531 -6145.85002
88 -8719.03978 -6560.33531
89 -21580.73672 -8719.03978
90 -11421.41937 -21580.73672
91 23316.78936 -11421.41937
92 -20089.31286 23316.78936
93 18165.55335 -20089.31286
94 -20980.06396 18165.55335
95 12439.71474 -20980.06396
96 -26612.77280 12439.71474
97 33201.44974 -26612.77280
98 -36537.21252 33201.44974
99 1162.93661 -36537.21252
100 -7970.28216 1162.93661
101 -29775.26493 -7970.28216
102 29633.36210 -29775.26493
103 49105.63041 29633.36210
104 55646.98610 49105.63041
105 54844.87927 55646.98610
106 -17976.46231 54844.87927
107 -11489.67373 -17976.46231
108 38952.17023 -11489.67373
109 -2309.02191 38952.17023
110 30032.46240 -2309.02191
111 -25075.90284 30032.46240
112 12158.51828 -25075.90284
113 -28030.36180 12158.51828
114 9231.69167 -28030.36180
115 -947.32947 9231.69167
116 -22922.92434 -947.32947
117 -14253.87419 -22922.92434
118 150969.60813 -14253.87419
119 -729.24788 150969.60813
120 -8589.89797 -729.24788
121 -10697.24222 -8589.89797
122 -10640.38763 -10697.24222
123 73062.45319 -10640.38763
124 36663.26374 73062.45319
125 -14519.47004 36663.26374
126 5945.38980 -14519.47004
127 11722.94613 5945.38980
128 -14786.14169 11722.94613
129 -4916.06948 -14786.14169
130 -6271.24646 -4916.06948
131 -18523.31821 -6271.24646
132 -8921.93940 -18523.31821
133 -6916.14734 -8921.93940
134 -29060.73164 -6916.14734
135 -19961.32032 -29060.73164
136 -1204.32800 -19961.32032
137 19947.74299 -1204.32800
138 30386.68001 19947.74299
139 -1652.13047 30386.68001
140 -8608.40959 -1652.13047
141 15826.97273 -8608.40959
142 -7870.56579 15826.97273
143 -15824.41021 -7870.56579
144 17690.49105 -15824.41021
145 -3389.70406 17690.49105
146 16886.12790 -3389.70406
147 53128.22833 16886.12790
148 -7134.25644 53128.22833
149 17603.84951 -7134.25644
150 -11568.34180 17603.84951
151 11052.98556 -11568.34180
152 -5622.20856 11052.98556
153 -19536.87823 -5622.20856
154 -10101.79230 -19536.87823
155 77250.40634 -10101.79230
156 2839.71446 77250.40634
157 -4338.73162 2839.71446
158 14123.09141 -4338.73162
159 -2550.86772 14123.09141
160 1811.42967 -2550.86772
161 51079.01578 1811.42967
162 9034.58690 51079.01578
163 29328.19999 9034.58690
164 -11749.88476 29328.19999
165 -15755.04400 -11749.88476
166 -20818.18137 -15755.04400
167 27823.98268 -20818.18137
168 -18092.15839 27823.98268
169 23215.54592 -18092.15839
170 -7170.41349 23215.54592
171 7477.89025 -7170.41349
172 -16221.31189 7477.89025
173 -25648.09778 -16221.31189
174 16477.15369 -25648.09778
175 -3340.07200 16477.15369
176 52676.62708 -3340.07200
177 -11461.73643 52676.62708
178 17025.04907 -11461.73643
179 -21918.43501 17025.04907
180 5052.75736 -21918.43501
181 5964.24089 5052.75736
182 3572.55621 5964.24089
183 -2689.42998 3572.55621
184 -51259.65616 -2689.42998
185 -1008.59910 -51259.65616
186 -7087.52246 -1008.59910
187 1500.67549 -7087.52246
188 -5887.19446 1500.67549
189 -16033.09227 -5887.19446
190 14461.69518 -16033.09227
191 -1679.29891 14461.69518
192 -14315.62554 -1679.29891
193 3677.43402 -14315.62554
194 12237.40452 3677.43402
195 -15371.95638 12237.40452
196 -2172.22128 -15371.95638
197 41505.67723 -2172.22128
198 19207.68038 41505.67723
199 11801.66566 19207.68038
200 -7896.66118 11801.66566
201 -21880.12447 -7896.66118
202 -7246.26593 -21880.12447
203 11404.15060 -7246.26593
204 -9890.80066 11404.15060
205 -4701.47702 -9890.80066
206 -9927.71427 -4701.47702
207 -4298.63690 -9927.71427
208 -12638.83720 -4298.63690
209 -26495.71659 -12638.83720
210 -14014.09676 -26495.71659
211 -8942.98367 -14014.09676
212 -2323.72498 -8942.98367
213 -7916.68224 -2323.72498
214 7788.01869 -7916.68224
215 -14556.75986 7788.01869
216 -253.26471 -14556.75986
217 7141.51138 -253.26471
218 -17642.56944 7141.51138
219 -7759.11381 -17642.56944
220 -1940.83582 -7759.11381
221 17151.77925 -1940.83582
222 -14476.32772 17151.77925
223 4026.57583 -14476.32772
224 -14599.05671 4026.57583
225 -24741.52426 -14599.05671
226 13370.70969 -24741.52426
227 4242.40298 13370.70969
228 -16340.25874 4242.40298
229 33210.72721 -16340.25874
230 -12030.55970 33210.72721
231 -2451.67853 -12030.55970
232 24801.71091 -2451.67853
233 -1247.85566 24801.71091
234 5705.22411 -1247.85566
235 -8370.20885 5705.22411
236 -8029.86732 -8370.20885
237 -13652.03584 -8029.86732
238 -9865.38996 -13652.03584
239 -5270.75600 -9865.38996
240 10555.46888 -5270.75600
241 2822.33330 10555.46888
242 -14282.89276 2822.33330
243 -639.18412 -14282.89276
244 2838.36956 -639.18412
245 7905.00251 2838.36956
246 6571.56030 7905.00251
247 1886.65688 6571.56030
248 11394.45670 1886.65688
249 22260.73114 11394.45670
250 2156.38632 22260.73114
251 1864.16876 2156.38632
252 -3864.07649 1864.16876
253 -1967.79371 -3864.07649
254 -9563.64834 -1967.79371
255 -6052.87404 -9563.64834
256 5656.27809 -6052.87404
257 -5196.41030 5656.27809
258 -12881.26047 -5196.41030
259 -4734.28788 -12881.26047
260 -7730.40554 -4734.28788
261 -13515.22842 -7730.40554
262 -6240.25962 -13515.22842
263 -11872.82824 -6240.25962
264 13655.27162 -11872.82824
265 -9797.43601 13655.27162
266 -13309.10582 -9797.43601
267 -2188.54726 -13309.10582
268 -17714.72682 -2188.54726
269 11263.51400 -17714.72682
270 -5939.71436 11263.51400
271 802.65750 -5939.71436
272 -1261.03174 802.65750
273 -14487.75180 -1261.03174
274 -3806.47724 -14487.75180
275 -25021.11006 -3806.47724
276 11077.22647 -25021.11006
277 -3263.52972 11077.22647
278 -13885.49992 -3263.52972
279 -8760.32989 -13885.49992
280 -7057.57248 -8760.32989
281 53995.88803 -7057.57248
282 -667.26316 53995.88803
283 -5593.60846 -667.26316
284 -2929.95068 -5593.60846
285 14519.69785 -2929.95068
286 -12382.35114 14519.69785
287 1952.19034 -12382.35114
288 -4309.21710 1952.19034
> 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/7vx2n1324509045.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/8njvn1324509045.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/9mcl81324509045.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/10iqfv1324509045.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/11zvk91324509045.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/12q5st1324509045.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/13xuh31324509045.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/1436hy1324509045.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/153jlq1324509045.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/162hep1324509045.tab")
+ }
>
> try(system("convert tmp/102hx1324509044.ps tmp/102hx1324509044.png",intern=TRUE))
character(0)
> try(system("convert tmp/2m6561324509044.ps tmp/2m6561324509044.png",intern=TRUE))
character(0)
> try(system("convert tmp/3c8fo1324509044.ps tmp/3c8fo1324509044.png",intern=TRUE))
character(0)
> try(system("convert tmp/4aulj1324509045.ps tmp/4aulj1324509045.png",intern=TRUE))
character(0)
> try(system("convert tmp/53s1o1324509045.ps tmp/53s1o1324509045.png",intern=TRUE))
character(0)
> try(system("convert tmp/668od1324509045.ps tmp/668od1324509045.png",intern=TRUE))
character(0)
> try(system("convert tmp/7vx2n1324509045.ps tmp/7vx2n1324509045.png",intern=TRUE))
character(0)
> try(system("convert tmp/8njvn1324509045.ps tmp/8njvn1324509045.png",intern=TRUE))
character(0)
> try(system("convert tmp/9mcl81324509045.ps tmp/9mcl81324509045.png",intern=TRUE))
character(0)
> try(system("convert tmp/10iqfv1324509045.ps tmp/10iqfv1324509045.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
8.686 0.734 9.458