R version 2.15.2 (2012-10-26) -- "Trick or Treat"
Copyright (C) 2012 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i686-pc-linux-gnu (32-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> x <- array(list(1418
+ ,210907
+ ,56
+ ,396
+ ,81
+ ,3
+ ,869
+ ,120982
+ ,56
+ ,297
+ ,55
+ ,4
+ ,1530
+ ,176508
+ ,54
+ ,559
+ ,50
+ ,12
+ ,2172
+ ,179321
+ ,89
+ ,967
+ ,125
+ ,2
+ ,901
+ ,123185
+ ,40
+ ,270
+ ,40
+ ,1
+ ,463
+ ,52746
+ ,25
+ ,143
+ ,37
+ ,3
+ ,3201
+ ,385534
+ ,92
+ ,1562
+ ,63
+ ,0
+ ,371
+ ,33170
+ ,18
+ ,109
+ ,44
+ ,0
+ ,1192
+ ,101645
+ ,63
+ ,371
+ ,88
+ ,0
+ ,1583
+ ,149061
+ ,44
+ ,656
+ ,66
+ ,5
+ ,1439
+ ,165446
+ ,33
+ ,511
+ ,57
+ ,0
+ ,1764
+ ,237213
+ ,84
+ ,655
+ ,74
+ ,0
+ ,1495
+ ,173326
+ ,88
+ ,465
+ ,49
+ ,7
+ ,1373
+ ,133131
+ ,55
+ ,525
+ ,52
+ ,7
+ ,2187
+ ,258873
+ ,60
+ ,885
+ ,88
+ ,3
+ ,1491
+ ,180083
+ ,66
+ ,497
+ ,36
+ ,9
+ ,4041
+ ,324799
+ ,154
+ ,1436
+ ,108
+ ,0
+ ,1706
+ ,230964
+ ,53
+ ,612
+ ,43
+ ,4
+ ,2152
+ ,236785
+ ,119
+ ,865
+ ,75
+ ,3
+ ,1036
+ ,135473
+ ,41
+ ,385
+ ,32
+ ,0
+ ,1882
+ ,202925
+ ,61
+ ,567
+ ,44
+ ,7
+ ,1929
+ ,215147
+ ,58
+ ,639
+ ,85
+ ,0
+ ,2242
+ ,344297
+ ,75
+ ,963
+ ,86
+ ,1
+ ,1220
+ ,153935
+ ,33
+ ,398
+ ,56
+ ,5
+ ,1289
+ ,132943
+ ,40
+ ,410
+ ,50
+ ,7
+ ,2515
+ ,174724
+ ,92
+ ,966
+ ,135
+ ,0
+ ,2147
+ ,174415
+ ,100
+ ,801
+ ,63
+ ,0
+ ,2352
+ ,225548
+ ,112
+ ,892
+ ,81
+ ,5
+ ,1638
+ ,223632
+ ,73
+ ,513
+ ,52
+ ,0
+ ,1222
+ ,124817
+ ,40
+ ,469
+ ,44
+ ,0
+ ,1812
+ ,221698
+ ,45
+ ,683
+ ,113
+ ,0
+ ,1677
+ ,210767
+ ,60
+ ,643
+ ,39
+ ,3
+ ,1579
+ ,170266
+ ,62
+ ,535
+ ,73
+ ,4
+ ,1731
+ ,260561
+ ,75
+ ,625
+ ,48
+ ,1
+ ,807
+ ,84853
+ ,31
+ ,264
+ ,33
+ ,4
+ ,2452
+ ,294424
+ ,77
+ ,992
+ ,59
+ ,2
+ ,829
+ ,101011
+ ,34
+ ,238
+ ,41
+ ,0
+ ,1940
+ ,215641
+ ,46
+ ,818
+ ,69
+ ,0
+ ,2662
+ ,325107
+ ,99
+ ,937
+ ,64
+ ,0
+ ,186
+ ,7176
+ ,17
+ ,70
+ ,1
+ ,0
+ ,1499
+ ,167542
+ ,66
+ ,507
+ ,59
+ ,2
+ ,865
+ ,106408
+ ,30
+ ,260
+ ,32
+ ,1
+ ,1793
+ ,96560
+ ,76
+ ,503
+ ,129
+ ,0
+ ,2527
+ ,265769
+ ,146
+ ,927
+ ,37
+ ,2
+ ,2747
+ ,269651
+ ,67
+ ,1269
+ ,31
+ ,10
+ ,1324
+ ,149112
+ ,56
+ ,537
+ ,65
+ ,6
+ ,2702
+ ,175824
+ ,107
+ ,910
+ ,107
+ ,0
+ ,1383
+ ,152871
+ ,58
+ ,532
+ ,74
+ ,5
+ ,1179
+ ,111665
+ ,34
+ ,345
+ ,54
+ ,4
+ ,2099
+ ,116408
+ ,61
+ ,918
+ ,76
+ ,1
+ ,4308
+ ,362301
+ ,119
+ ,1635
+ ,715
+ ,2
+ ,918
+ ,78800
+ ,42
+ ,330
+ ,57
+ ,2
+ ,1831
+ ,183167
+ ,66
+ ,557
+ ,66
+ ,0
+ ,3373
+ ,277965
+ ,89
+ ,1178
+ ,106
+ ,8
+ ,1713
+ ,150629
+ ,44
+ ,740
+ ,54
+ ,3
+ ,1438
+ ,168809
+ ,66
+ ,452
+ ,32
+ ,0
+ ,496
+ ,24188
+ ,24
+ ,218
+ ,20
+ ,0
+ ,2253
+ ,329267
+ ,259
+ ,764
+ ,71
+ ,8
+ ,744
+ ,65029
+ ,17
+ ,255
+ ,21
+ ,5
+ ,1161
+ ,101097
+ ,64
+ ,454
+ ,70
+ ,3
+ ,2352
+ ,218946
+ ,41
+ ,866
+ ,112
+ ,1
+ ,2144
+ ,244052
+ ,68
+ ,574
+ ,66
+ ,5
+ ,4691
+ ,341570
+ ,168
+ ,1276
+ ,190
+ ,1
+ ,1112
+ ,103597
+ ,43
+ ,379
+ ,66
+ ,1
+ ,2694
+ ,233328
+ ,132
+ ,825
+ ,165
+ ,5
+ ,1973
+ ,256462
+ ,105
+ ,798
+ ,56
+ ,0
+ ,1769
+ ,206161
+ ,71
+ ,663
+ ,61
+ ,12
+ ,3148
+ ,311473
+ ,112
+ ,1069
+ ,53
+ ,8
+ ,2474
+ ,235800
+ ,94
+ ,921
+ ,127
+ ,8
+ ,2084
+ ,177939
+ ,82
+ ,858
+ ,63
+ ,8
+ ,1954
+ ,207176
+ ,70
+ ,711
+ ,38
+ ,8
+ ,1226
+ ,196553
+ ,57
+ ,503
+ ,50
+ ,2
+ ,1389
+ ,174184
+ ,53
+ ,382
+ ,52
+ ,0
+ ,1496
+ ,143246
+ ,103
+ ,464
+ ,42
+ ,5
+ ,2269
+ ,187559
+ ,121
+ ,717
+ ,76
+ ,8
+ ,1833
+ ,187681
+ ,62
+ ,690
+ ,67
+ ,2
+ ,1268
+ ,119016
+ ,52
+ ,462
+ ,50
+ ,5
+ ,1943
+ ,182192
+ ,52
+ ,657
+ ,53
+ ,12
+ ,893
+ ,73566
+ ,32
+ ,385
+ ,39
+ ,6
+ ,1762
+ ,194979
+ ,62
+ ,577
+ ,50
+ ,7
+ ,1403
+ ,167488
+ ,45
+ ,619
+ ,77
+ ,2
+ ,1425
+ ,143756
+ ,46
+ ,479
+ ,57
+ ,0
+ ,1857
+ ,275541
+ ,63
+ ,817
+ ,73
+ ,4
+ ,1840
+ ,243199
+ ,75
+ ,752
+ ,34
+ ,3
+ ,1502
+ ,182999
+ ,88
+ ,430
+ ,39
+ ,6
+ ,1441
+ ,135649
+ ,46
+ ,451
+ ,46
+ ,2
+ ,1420
+ ,152299
+ ,53
+ ,537
+ ,63
+ ,0
+ ,1416
+ ,120221
+ ,37
+ ,519
+ ,35
+ ,1
+ ,2970
+ ,346485
+ ,90
+ ,1000
+ ,106
+ ,0
+ ,1317
+ ,145790
+ ,63
+ ,637
+ ,43
+ ,5
+ ,1644
+ ,193339
+ ,78
+ ,465
+ ,47
+ ,2
+ ,870
+ ,80953
+ ,25
+ ,437
+ ,31
+ ,0
+ ,1654
+ ,122774
+ ,45
+ ,711
+ ,162
+ ,0
+ ,1054
+ ,130585
+ ,46
+ ,299
+ ,57
+ ,5
+ ,937
+ ,112611
+ ,41
+ ,248
+ ,36
+ ,0
+ ,3004
+ ,286468
+ ,144
+ ,1162
+ ,263
+ ,1
+ ,2008
+ ,241066
+ ,82
+ ,714
+ ,78
+ ,0
+ ,2547
+ ,148446
+ ,91
+ ,905
+ ,63
+ ,1
+ ,1885
+ ,204713
+ ,71
+ ,649
+ ,54
+ ,1
+ ,1626
+ ,182079
+ ,63
+ ,512
+ ,63
+ ,2
+ ,1468
+ ,140344
+ ,53
+ ,472
+ ,77
+ ,6
+ ,2445
+ ,220516
+ ,62
+ ,905
+ ,79
+ ,1
+ ,1964
+ ,243060
+ ,63
+ ,786
+ ,110
+ ,4
+ ,1381
+ ,162765
+ ,32
+ ,489
+ ,56
+ ,2
+ ,1369
+ ,182613
+ ,39
+ ,479
+ ,56
+ ,3
+ ,1659
+ ,232138
+ ,62
+ ,617
+ ,43
+ ,0
+ ,2888
+ ,265318
+ ,117
+ ,925
+ ,111
+ ,10
+ ,1290
+ ,85574
+ ,34
+ ,351
+ ,71
+ ,0
+ ,2845
+ ,310839
+ ,92
+ ,1144
+ ,62
+ ,9
+ ,1982
+ ,225060
+ ,93
+ ,669
+ ,56
+ ,7
+ ,1904
+ ,232317
+ ,54
+ ,707
+ ,74
+ ,0
+ ,1391
+ ,144966
+ ,144
+ ,458
+ ,60
+ ,0
+ ,602
+ ,43287
+ ,14
+ ,214
+ ,43
+ ,4
+ ,1743
+ ,155754
+ ,61
+ ,599
+ ,68
+ ,4
+ ,1559
+ ,164709
+ ,109
+ ,572
+ ,53
+ ,0
+ ,2014
+ ,201940
+ ,38
+ ,897
+ ,87
+ ,0
+ ,2143
+ ,235454
+ ,73
+ ,819
+ ,46
+ ,0
+ ,2146
+ ,220801
+ ,75
+ ,720
+ ,105
+ ,1
+ ,874
+ ,99466
+ ,50
+ ,273
+ ,32
+ ,0
+ ,1590
+ ,92661
+ ,61
+ ,508
+ ,133
+ ,1
+ ,1590
+ ,133328
+ ,55
+ ,506
+ ,79
+ ,0
+ ,1210
+ ,61361
+ ,77
+ ,451
+ ,51
+ ,0
+ ,2072
+ ,125930
+ ,75
+ ,699
+ ,207
+ ,4
+ ,1281
+ ,100750
+ ,72
+ ,407
+ ,67
+ ,0
+ ,1401
+ ,224549
+ ,50
+ ,465
+ ,47
+ ,4
+ ,834
+ ,82316
+ ,32
+ ,245
+ ,34
+ ,4
+ ,1105
+ ,102010
+ ,53
+ ,370
+ ,66
+ ,3
+ ,1272
+ ,101523
+ ,42
+ ,316
+ ,76
+ ,0
+ ,1944
+ ,243511
+ ,71
+ ,603
+ ,65
+ ,0
+ ,391
+ ,22938
+ ,10
+ ,154
+ ,9
+ ,0
+ ,761
+ ,41566
+ ,35
+ ,229
+ ,42
+ ,5
+ ,1605
+ ,152474
+ ,65
+ ,577
+ ,45
+ ,0
+ ,530
+ ,61857
+ ,25
+ ,192
+ ,25
+ ,4
+ ,1988
+ ,99923
+ ,66
+ ,617
+ ,115
+ ,0
+ ,1386
+ ,132487
+ ,41
+ ,411
+ ,97
+ ,0
+ ,2395
+ ,317394
+ ,86
+ ,975
+ ,53
+ ,1
+ ,387
+ ,21054
+ ,16
+ ,146
+ ,2
+ ,0
+ ,1742
+ ,209641
+ ,42
+ ,705
+ ,52
+ ,5
+ ,620
+ ,22648
+ ,19
+ ,184
+ ,44
+ ,0
+ ,449
+ ,31414
+ ,19
+ ,200
+ ,22
+ ,0
+ ,800
+ ,46698
+ ,45
+ ,274
+ ,35
+ ,0
+ ,1684
+ ,131698
+ ,65
+ ,502
+ ,74
+ ,0
+ ,1050
+ ,91735
+ ,35
+ ,382
+ ,103
+ ,0
+ ,2699
+ ,244749
+ ,95
+ ,964
+ ,144
+ ,2
+ ,1606
+ ,184510
+ ,49
+ ,537
+ ,60
+ ,7
+ ,1502
+ ,79863
+ ,37
+ ,438
+ ,134
+ ,1
+ ,1204
+ ,128423
+ ,64
+ ,369
+ ,89
+ ,8
+ ,1138
+ ,97839
+ ,38
+ ,417
+ ,42
+ ,2
+ ,568
+ ,38214
+ ,34
+ ,276
+ ,52
+ ,0
+ ,1459
+ ,151101
+ ,32
+ ,514
+ ,98
+ ,2
+ ,2158
+ ,272458
+ ,65
+ ,822
+ ,99
+ ,0
+ ,1111
+ ,172494
+ ,52
+ ,389
+ ,52
+ ,0
+ ,1421
+ ,108043
+ ,62
+ ,466
+ ,29
+ ,1
+ ,2833
+ ,328107
+ ,65
+ ,1255
+ ,125
+ ,3
+ ,1955
+ ,250579
+ ,83
+ ,694
+ ,106
+ ,0
+ ,2922
+ ,351067
+ ,95
+ ,1024
+ ,95
+ ,3
+ ,1002
+ ,158015
+ ,29
+ ,400
+ ,40
+ ,0
+ ,1060
+ ,98866
+ ,18
+ ,397
+ ,140
+ ,0
+ ,956
+ ,85439
+ ,33
+ ,350
+ ,43
+ ,0
+ ,2186
+ ,229242
+ ,247
+ ,719
+ ,128
+ ,4
+ ,3604
+ ,351619
+ ,139
+ ,1277
+ ,142
+ ,4
+ ,1035
+ ,84207
+ ,29
+ ,356
+ ,73
+ ,11
+ ,1417
+ ,120445
+ ,118
+ ,457
+ ,72
+ ,0
+ ,3261
+ ,324598
+ ,110
+ ,1402
+ ,128
+ ,0
+ ,1587
+ ,131069
+ ,67
+ ,600
+ ,61
+ ,4
+ ,1424
+ ,204271
+ ,42
+ ,480
+ ,73
+ ,0
+ ,1701
+ ,165543
+ ,65
+ ,595
+ ,148
+ ,1
+ ,1249
+ ,141722
+ ,94
+ ,436
+ ,64
+ ,0
+ ,946
+ ,116048
+ ,64
+ ,230
+ ,45
+ ,0
+ ,1926
+ ,250047
+ ,81
+ ,651
+ ,58
+ ,0
+ ,3352
+ ,299775
+ ,95
+ ,1367
+ ,97
+ ,9
+ ,1641
+ ,195838
+ ,67
+ ,564
+ ,50
+ ,1
+ ,2035
+ ,173260
+ ,63
+ ,716
+ ,37
+ ,3
+ ,2312
+ ,254488
+ ,83
+ ,747
+ ,50
+ ,10
+ ,1369
+ ,104389
+ ,45
+ ,467
+ ,105
+ ,5
+ ,1577
+ ,136084
+ ,30
+ ,671
+ ,69
+ ,0
+ ,2201
+ ,199476
+ ,70
+ ,861
+ ,46
+ ,2
+ ,961
+ ,92499
+ ,32
+ ,319
+ ,57
+ ,0
+ ,1900
+ ,224330
+ ,83
+ ,612
+ ,52
+ ,1
+ ,1254
+ ,135781
+ ,31
+ ,433
+ ,98
+ ,2
+ ,1335
+ ,74408
+ ,67
+ ,434
+ ,61
+ ,4
+ ,1597
+ ,81240
+ ,66
+ ,503
+ ,89
+ ,0
+ ,207
+ ,14688
+ ,10
+ ,85
+ ,0
+ ,0
+ ,1645
+ ,181633
+ ,70
+ ,564
+ ,48
+ ,2
+ ,2429
+ ,271856
+ ,103
+ ,824
+ ,91
+ ,1
+ ,151
+ ,7199
+ ,5
+ ,74
+ ,0
+ ,0
+ ,474
+ ,46660
+ ,20
+ ,259
+ ,7
+ ,0
+ ,141
+ ,17547
+ ,5
+ ,69
+ ,3
+ ,0
+ ,1639
+ ,133368
+ ,36
+ ,535
+ ,54
+ ,1
+ ,872
+ ,95227
+ ,34
+ ,239
+ ,70
+ ,0
+ ,1318
+ ,152601
+ ,48
+ ,438
+ ,36
+ ,2
+ ,1018
+ ,98146
+ ,40
+ ,459
+ ,37
+ ,0
+ ,1383
+ ,79619
+ ,43
+ ,426
+ ,123
+ ,3
+ ,1314
+ ,59194
+ ,31
+ ,288
+ ,247
+ ,6
+ ,1335
+ ,139942
+ ,42
+ ,498
+ ,46
+ ,0
+ ,1403
+ ,118612
+ ,46
+ ,454
+ ,72
+ ,2
+ ,910
+ ,72880
+ ,33
+ ,376
+ ,41
+ ,0
+ ,616
+ ,65475
+ ,18
+ ,225
+ ,24
+ ,2
+ ,1407
+ ,99643
+ ,55
+ ,555
+ ,45
+ ,1
+ ,771
+ ,71965
+ ,35
+ ,252
+ ,33
+ ,1
+ ,766
+ ,77272
+ ,59
+ ,208
+ ,27
+ ,2
+ ,473
+ ,49289
+ ,19
+ ,130
+ ,36
+ ,1
+ ,1376
+ ,135131
+ ,66
+ ,481
+ ,87
+ ,0
+ ,1232
+ ,108446
+ ,60
+ ,389
+ ,90
+ ,1
+ ,1521
+ ,89746
+ ,36
+ ,565
+ ,114
+ ,3
+ ,572
+ ,44296
+ ,25
+ ,173
+ ,31
+ ,0
+ ,1059
+ ,77648
+ ,47
+ ,278
+ ,45
+ ,0
+ ,1544
+ ,181528
+ ,54
+ ,609
+ ,69
+ ,0
+ ,1230
+ ,134019
+ ,53
+ ,422
+ ,51
+ ,0
+ ,1206
+ ,124064
+ ,40
+ ,445
+ ,34
+ ,1
+ ,1205
+ ,92630
+ ,40
+ ,387
+ ,60
+ ,4
+ ,1255
+ ,121848
+ ,39
+ ,339
+ ,45
+ ,0
+ ,613
+ ,52915
+ ,14
+ ,181
+ ,54
+ ,0
+ ,721
+ ,81872
+ ,45
+ ,245
+ ,25
+ ,0
+ ,1109
+ ,58981
+ ,36
+ ,384
+ ,38
+ ,7
+ ,740
+ ,53515
+ ,28
+ ,212
+ ,52
+ ,2
+ ,1126
+ ,60812
+ ,44
+ ,399
+ ,67
+ ,0
+ ,728
+ ,56375
+ ,30
+ ,229
+ ,74
+ ,7
+ ,689
+ ,65490
+ ,22
+ ,224
+ ,38
+ ,3
+ ,592
+ ,80949
+ ,17
+ ,203
+ ,30
+ ,0
+ ,995
+ ,76302
+ ,31
+ ,333
+ ,26
+ ,0
+ ,1613
+ ,104011
+ ,55
+ ,384
+ ,67
+ ,6
+ ,2048
+ ,98104
+ ,54
+ ,636
+ ,132
+ ,2
+ ,705
+ ,67989
+ ,21
+ ,185
+ ,42
+ ,0
+ ,301
+ ,30989
+ ,14
+ ,93
+ ,35
+ ,0
+ ,1803
+ ,135458
+ ,81
+ ,581
+ ,118
+ ,3
+ ,799
+ ,73504
+ ,35
+ ,248
+ ,68
+ ,0
+ ,861
+ ,63123
+ ,43
+ ,304
+ ,43
+ ,1
+ ,1186
+ ,61254
+ ,46
+ ,344
+ ,76
+ ,1
+ ,1451
+ ,74914
+ ,30
+ ,407
+ ,64
+ ,0
+ ,628
+ ,31774
+ ,23
+ ,170
+ ,48
+ ,1
+ ,1161
+ ,81437
+ ,38
+ ,312
+ ,64
+ ,0
+ ,1463
+ ,87186
+ ,54
+ ,507
+ ,56
+ ,0
+ ,742
+ ,50090
+ ,20
+ ,224
+ ,71
+ ,0
+ ,979
+ ,65745
+ ,53
+ ,340
+ ,75
+ ,0
+ ,675
+ ,56653
+ ,45
+ ,168
+ ,39
+ ,0
+ ,1241
+ ,158399
+ ,39
+ ,443
+ ,42
+ ,0
+ ,676
+ ,46455
+ ,20
+ ,204
+ ,39
+ ,0
+ ,1049
+ ,73624
+ ,24
+ ,367
+ ,93
+ ,0
+ ,620
+ ,38395
+ ,31
+ ,210
+ ,38
+ ,0
+ ,1081
+ ,91899
+ ,35
+ ,335
+ ,60
+ ,0
+ ,1688
+ ,139526
+ ,151
+ ,364
+ ,71
+ ,0
+ ,736
+ ,52164
+ ,52
+ ,178
+ ,52
+ ,0
+ ,617
+ ,51567
+ ,30
+ ,206
+ ,27
+ ,2
+ ,812
+ ,70551
+ ,31
+ ,279
+ ,59
+ ,0
+ ,1051
+ ,84856
+ ,29
+ ,387
+ ,40
+ ,1
+ ,1656
+ ,102538
+ ,57
+ ,490
+ ,79
+ ,1
+ ,705
+ ,86678
+ ,40
+ ,238
+ ,44
+ ,0
+ ,945
+ ,85709
+ ,44
+ ,343
+ ,65
+ ,0
+ ,554
+ ,34662
+ ,25
+ ,232
+ ,10
+ ,0
+ ,1597
+ ,150580
+ ,77
+ ,530
+ ,124
+ ,0
+ ,982
+ ,99611
+ ,35
+ ,291
+ ,81
+ ,0
+ ,222
+ ,19349
+ ,11
+ ,67
+ ,15
+ ,0
+ ,1212
+ ,99373
+ ,63
+ ,397
+ ,92
+ ,1
+ ,1143
+ ,86230
+ ,44
+ ,467
+ ,42
+ ,0
+ ,435
+ ,30837
+ ,19
+ ,178
+ ,10
+ ,0
+ ,532
+ ,31706
+ ,13
+ ,175
+ ,24
+ ,0
+ ,882
+ ,89806
+ ,42
+ ,299
+ ,64
+ ,0
+ ,608
+ ,62088
+ ,38
+ ,154
+ ,45
+ ,1
+ ,459
+ ,40151
+ ,29
+ ,106
+ ,22
+ ,0
+ ,578
+ ,27634
+ ,20
+ ,189
+ ,56
+ ,0
+ ,826
+ ,76990
+ ,27
+ ,194
+ ,94
+ ,0
+ ,509
+ ,37460
+ ,20
+ ,135
+ ,19
+ ,0
+ ,717
+ ,54157
+ ,19
+ ,201
+ ,35
+ ,0
+ ,637
+ ,49862
+ ,37
+ ,207
+ ,32
+ ,0
+ ,857
+ ,84337
+ ,26
+ ,280
+ ,35
+ ,0
+ ,830
+ ,64175
+ ,42
+ ,260
+ ,48
+ ,0
+ ,652
+ ,59382
+ ,49
+ ,227
+ ,49
+ ,0
+ ,707
+ ,119308
+ ,30
+ ,239
+ ,48
+ ,0
+ ,954
+ ,76702
+ ,49
+ ,333
+ ,62
+ ,0
+ ,1461
+ ,103425
+ ,67
+ ,428
+ ,96
+ ,1
+ ,672
+ ,70344
+ ,28
+ ,230
+ ,45
+ ,0
+ ,778
+ ,43410
+ ,19
+ ,292
+ ,63
+ ,0
+ ,1141
+ ,104838
+ ,49
+ ,350
+ ,71
+ ,1
+ ,680
+ ,62215
+ ,27
+ ,186
+ ,26
+ ,0
+ ,1090
+ ,69304
+ ,30
+ ,326
+ ,48
+ ,6
+ ,616
+ ,53117
+ ,22
+ ,155
+ ,29
+ ,3
+ ,285
+ ,19764
+ ,12
+ ,75
+ ,19
+ ,1
+ ,1145
+ ,86680
+ ,31
+ ,361
+ ,45
+ ,2
+ ,733
+ ,84105
+ ,20
+ ,261
+ ,45
+ ,0
+ ,888
+ ,77945
+ ,20
+ ,299
+ ,67
+ ,0
+ ,849
+ ,89113
+ ,39
+ ,300
+ ,30
+ ,0
+ ,1182
+ ,91005
+ ,29
+ ,450
+ ,36
+ ,3
+ ,528
+ ,40248
+ ,16
+ ,183
+ ,34
+ ,1
+ ,642
+ ,64187
+ ,27
+ ,238
+ ,36
+ ,0
+ ,947
+ ,50857
+ ,21
+ ,165
+ ,34
+ ,0
+ ,819
+ ,56613
+ ,19
+ ,234
+ ,37
+ ,1
+ ,757
+ ,62792
+ ,35
+ ,176
+ ,46
+ ,0
+ ,894
+ ,72535
+ ,14
+ ,329
+ ,44
+ ,0)
+ ,dim=c(6
+ ,289)
+ ,dimnames=list(c('
pageviews'
+ ,'time_in_rfc'
+ ,'logins'
+ ,'compendium_views_info'
+ ,'compendium_views_pr'
+ ,'shared_compendiums
')
+ ,1:289))
> y <- array(NA,dim=c(6,289),dimnames=list(c('
pageviews','time_in_rfc','logins','compendium_views_info','compendium_views_pr','shared_compendiums
'),1:289))
> for (i in 1:dim(x)[1])
+ {
+ for (j in 1:dim(x)[2])
+ {
+ y[i,j] <- as.numeric(x[i,j])
+ }
+ }
> par3 = 'No Linear Trend'
> par2 = 'Do not include Seasonal Dummies'
> par1 = '1'
> par3 <- 'No Linear Trend'
> par2 <- 'Do not include Seasonal Dummies'
> par1 <- '1'
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Prof. Dr. P. Wessa
> #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #Technical description: Write here your technical program description (don't use hard returns!)
> library(lattice)
> library(lmtest)
Loading required package: zoo
Attaching package: 'zoo'
The following object(s) are masked from 'package:base':
as.Date, as.Date.numeric
> n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test
> par1 <- as.numeric(par1)
> x <- t(y)
> k <- length(x[1,])
> n <- length(x[,1])
> x1 <- cbind(x[,par1], x[,1:k!=par1])
> mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1])
> colnames(x1) <- mycolnames #colnames(x)[par1]
> x <- x1
> if (par3 == 'First Differences'){
+ x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep='')))
+ for (i in 1:n-1) {
+ for (j in 1:k) {
+ x2[i,j] <- x[i+1,j] - x[i,j]
+ }
+ }
+ x <- x2
+ }
> if (par2 == 'Include Monthly Dummies'){
+ x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep ='')))
+ for (i in 1:11){
+ x2[seq(i,n,12),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> if (par2 == 'Include Quarterly Dummies'){
+ x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep ='')))
+ for (i in 1:3){
+ x2[seq(i,n,4),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> k <- length(x[1,])
> if (par3 == 'Linear Trend'){
+ x <- cbind(x, c(1:n))
+ colnames(x)[k+1] <- 't'
+ }
> x
\rpageviews time_in_rfc logins compendium_views_info compendium_views_pr
1 1418 210907 56 396 81
2 869 120982 56 297 55
3 1530 176508 54 559 50
4 2172 179321 89 967 125
5 901 123185 40 270 40
6 463 52746 25 143 37
7 3201 385534 92 1562 63
8 371 33170 18 109 44
9 1192 101645 63 371 88
10 1583 149061 44 656 66
11 1439 165446 33 511 57
12 1764 237213 84 655 74
13 1495 173326 88 465 49
14 1373 133131 55 525 52
15 2187 258873 60 885 88
16 1491 180083 66 497 36
17 4041 324799 154 1436 108
18 1706 230964 53 612 43
19 2152 236785 119 865 75
20 1036 135473 41 385 32
21 1882 202925 61 567 44
22 1929 215147 58 639 85
23 2242 344297 75 963 86
24 1220 153935 33 398 56
25 1289 132943 40 410 50
26 2515 174724 92 966 135
27 2147 174415 100 801 63
28 2352 225548 112 892 81
29 1638 223632 73 513 52
30 1222 124817 40 469 44
31 1812 221698 45 683 113
32 1677 210767 60 643 39
33 1579 170266 62 535 73
34 1731 260561 75 625 48
35 807 84853 31 264 33
36 2452 294424 77 992 59
37 829 101011 34 238 41
38 1940 215641 46 818 69
39 2662 325107 99 937 64
40 186 7176 17 70 1
41 1499 167542 66 507 59
42 865 106408 30 260 32
43 1793 96560 76 503 129
44 2527 265769 146 927 37
45 2747 269651 67 1269 31
46 1324 149112 56 537 65
47 2702 175824 107 910 107
48 1383 152871 58 532 74
49 1179 111665 34 345 54
50 2099 116408 61 918 76
51 4308 362301 119 1635 715
52 918 78800 42 330 57
53 1831 183167 66 557 66
54 3373 277965 89 1178 106
55 1713 150629 44 740 54
56 1438 168809 66 452 32
57 496 24188 24 218 20
58 2253 329267 259 764 71
59 744 65029 17 255 21
60 1161 101097 64 454 70
61 2352 218946 41 866 112
62 2144 244052 68 574 66
63 4691 341570 168 1276 190
64 1112 103597 43 379 66
65 2694 233328 132 825 165
66 1973 256462 105 798 56
67 1769 206161 71 663 61
68 3148 311473 112 1069 53
69 2474 235800 94 921 127
70 2084 177939 82 858 63
71 1954 207176 70 711 38
72 1226 196553 57 503 50
73 1389 174184 53 382 52
74 1496 143246 103 464 42
75 2269 187559 121 717 76
76 1833 187681 62 690 67
77 1268 119016 52 462 50
78 1943 182192 52 657 53
79 893 73566 32 385 39
80 1762 194979 62 577 50
81 1403 167488 45 619 77
82 1425 143756 46 479 57
83 1857 275541 63 817 73
84 1840 243199 75 752 34
85 1502 182999 88 430 39
86 1441 135649 46 451 46
87 1420 152299 53 537 63
88 1416 120221 37 519 35
89 2970 346485 90 1000 106
90 1317 145790 63 637 43
91 1644 193339 78 465 47
92 870 80953 25 437 31
93 1654 122774 45 711 162
94 1054 130585 46 299 57
95 937 112611 41 248 36
96 3004 286468 144 1162 263
97 2008 241066 82 714 78
98 2547 148446 91 905 63
99 1885 204713 71 649 54
100 1626 182079 63 512 63
101 1468 140344 53 472 77
102 2445 220516 62 905 79
103 1964 243060 63 786 110
104 1381 162765 32 489 56
105 1369 182613 39 479 56
106 1659 232138 62 617 43
107 2888 265318 117 925 111
108 1290 85574 34 351 71
109 2845 310839 92 1144 62
110 1982 225060 93 669 56
111 1904 232317 54 707 74
112 1391 144966 144 458 60
113 602 43287 14 214 43
114 1743 155754 61 599 68
115 1559 164709 109 572 53
116 2014 201940 38 897 87
117 2143 235454 73 819 46
118 2146 220801 75 720 105
119 874 99466 50 273 32
120 1590 92661 61 508 133
121 1590 133328 55 506 79
122 1210 61361 77 451 51
123 2072 125930 75 699 207
124 1281 100750 72 407 67
125 1401 224549 50 465 47
126 834 82316 32 245 34
127 1105 102010 53 370 66
128 1272 101523 42 316 76
129 1944 243511 71 603 65
130 391 22938 10 154 9
131 761 41566 35 229 42
132 1605 152474 65 577 45
133 530 61857 25 192 25
134 1988 99923 66 617 115
135 1386 132487 41 411 97
136 2395 317394 86 975 53
137 387 21054 16 146 2
138 1742 209641 42 705 52
139 620 22648 19 184 44
140 449 31414 19 200 22
141 800 46698 45 274 35
142 1684 131698 65 502 74
143 1050 91735 35 382 103
144 2699 244749 95 964 144
145 1606 184510 49 537 60
146 1502 79863 37 438 134
147 1204 128423 64 369 89
148 1138 97839 38 417 42
149 568 38214 34 276 52
150 1459 151101 32 514 98
151 2158 272458 65 822 99
152 1111 172494 52 389 52
153 1421 108043 62 466 29
154 2833 328107 65 1255 125
155 1955 250579 83 694 106
156 2922 351067 95 1024 95
157 1002 158015 29 400 40
158 1060 98866 18 397 140
159 956 85439 33 350 43
160 2186 229242 247 719 128
161 3604 351619 139 1277 142
162 1035 84207 29 356 73
163 1417 120445 118 457 72
164 3261 324598 110 1402 128
165 1587 131069 67 600 61
166 1424 204271 42 480 73
167 1701 165543 65 595 148
168 1249 141722 94 436 64
169 946 116048 64 230 45
170 1926 250047 81 651 58
171 3352 299775 95 1367 97
172 1641 195838 67 564 50
173 2035 173260 63 716 37
174 2312 254488 83 747 50
175 1369 104389 45 467 105
176 1577 136084 30 671 69
177 2201 199476 70 861 46
178 961 92499 32 319 57
179 1900 224330 83 612 52
180 1254 135781 31 433 98
181 1335 74408 67 434 61
182 1597 81240 66 503 89
183 207 14688 10 85 0
184 1645 181633 70 564 48
185 2429 271856 103 824 91
186 151 7199 5 74 0
187 474 46660 20 259 7
188 141 17547 5 69 3
189 1639 133368 36 535 54
190 872 95227 34 239 70
191 1318 152601 48 438 36
192 1018 98146 40 459 37
193 1383 79619 43 426 123
194 1314 59194 31 288 247
195 1335 139942 42 498 46
196 1403 118612 46 454 72
197 910 72880 33 376 41
198 616 65475 18 225 24
199 1407 99643 55 555 45
200 771 71965 35 252 33
201 766 77272 59 208 27
202 473 49289 19 130 36
203 1376 135131 66 481 87
204 1232 108446 60 389 90
205 1521 89746 36 565 114
206 572 44296 25 173 31
207 1059 77648 47 278 45
208 1544 181528 54 609 69
209 1230 134019 53 422 51
210 1206 124064 40 445 34
211 1205 92630 40 387 60
212 1255 121848 39 339 45
213 613 52915 14 181 54
214 721 81872 45 245 25
215 1109 58981 36 384 38
216 740 53515 28 212 52
217 1126 60812 44 399 67
218 728 56375 30 229 74
219 689 65490 22 224 38
220 592 80949 17 203 30
221 995 76302 31 333 26
222 1613 104011 55 384 67
223 2048 98104 54 636 132
224 705 67989 21 185 42
225 301 30989 14 93 35
226 1803 135458 81 581 118
227 799 73504 35 248 68
228 861 63123 43 304 43
229 1186 61254 46 344 76
230 1451 74914 30 407 64
231 628 31774 23 170 48
232 1161 81437 38 312 64
233 1463 87186 54 507 56
234 742 50090 20 224 71
235 979 65745 53 340 75
236 675 56653 45 168 39
237 1241 158399 39 443 42
238 676 46455 20 204 39
239 1049 73624 24 367 93
240 620 38395 31 210 38
241 1081 91899 35 335 60
242 1688 139526 151 364 71
243 736 52164 52 178 52
244 617 51567 30 206 27
245 812 70551 31 279 59
246 1051 84856 29 387 40
247 1656 102538 57 490 79
248 705 86678 40 238 44
249 945 85709 44 343 65
250 554 34662 25 232 10
251 1597 150580 77 530 124
252 982 99611 35 291 81
253 222 19349 11 67 15
254 1212 99373 63 397 92
255 1143 86230 44 467 42
256 435 30837 19 178 10
257 532 31706 13 175 24
258 882 89806 42 299 64
259 608 62088 38 154 45
260 459 40151 29 106 22
261 578 27634 20 189 56
262 826 76990 27 194 94
263 509 37460 20 135 19
264 717 54157 19 201 35
265 637 49862 37 207 32
266 857 84337 26 280 35
267 830 64175 42 260 48
268 652 59382 49 227 49
269 707 119308 30 239 48
270 954 76702 49 333 62
271 1461 103425 67 428 96
272 672 70344 28 230 45
273 778 43410 19 292 63
274 1141 104838 49 350 71
275 680 62215 27 186 26
276 1090 69304 30 326 48
277 616 53117 22 155 29
278 285 19764 12 75 19
279 1145 86680 31 361 45
280 733 84105 20 261 45
281 888 77945 20 299 67
282 849 89113 39 300 30
283 1182 91005 29 450 36
284 528 40248 16 183 34
285 642 64187 27 238 36
286 947 50857 21 165 34
287 819 56613 19 234 37
288 757 62792 35 176 46
289 894 72535 14 329 44
shared_compendiums\r
1 3
2 4
3 12
4 2
5 1
6 3
7 0
8 0
9 0
10 5
11 0
12 0
13 7
14 7
15 3
16 9
17 0
18 4
19 3
20 0
21 7
22 0
23 1
24 5
25 7
26 0
27 0
28 5
29 0
30 0
31 0
32 3
33 4
34 1
35 4
36 2
37 0
38 0
39 0
40 0
41 2
42 1
43 0
44 2
45 10
46 6
47 0
48 5
49 4
50 1
51 2
52 2
53 0
54 8
55 3
56 0
57 0
58 8
59 5
60 3
61 1
62 5
63 1
64 1
65 5
66 0
67 12
68 8
69 8
70 8
71 8
72 2
73 0
74 5
75 8
76 2
77 5
78 12
79 6
80 7
81 2
82 0
83 4
84 3
85 6
86 2
87 0
88 1
89 0
90 5
91 2
92 0
93 0
94 5
95 0
96 1
97 0
98 1
99 1
100 2
101 6
102 1
103 4
104 2
105 3
106 0
107 10
108 0
109 9
110 7
111 0
112 0
113 4
114 4
115 0
116 0
117 0
118 1
119 0
120 1
121 0
122 0
123 4
124 0
125 4
126 4
127 3
128 0
129 0
130 0
131 5
132 0
133 4
134 0
135 0
136 1
137 0
138 5
139 0
140 0
141 0
142 0
143 0
144 2
145 7
146 1
147 8
148 2
149 0
150 2
151 0
152 0
153 1
154 3
155 0
156 3
157 0
158 0
159 0
160 4
161 4
162 11
163 0
164 0
165 4
166 0
167 1
168 0
169 0
170 0
171 9
172 1
173 3
174 10
175 5
176 0
177 2
178 0
179 1
180 2
181 4
182 0
183 0
184 2
185 1
186 0
187 0
188 0
189 1
190 0
191 2
192 0
193 3
194 6
195 0
196 2
197 0
198 2
199 1
200 1
201 2
202 1
203 0
204 1
205 3
206 0
207 0
208 0
209 0
210 1
211 4
212 0
213 0
214 0
215 7
216 2
217 0
218 7
219 3
220 0
221 0
222 6
223 2
224 0
225 0
226 3
227 0
228 1
229 1
230 0
231 1
232 0
233 0
234 0
235 0
236 0
237 0
238 0
239 0
240 0
241 0
242 0
243 0
244 2
245 0
246 1
247 1
248 0
249 0
250 0
251 0
252 0
253 0
254 1
255 0
256 0
257 0
258 0
259 1
260 0
261 0
262 0
263 0
264 0
265 0
266 0
267 0
268 0
269 0
270 0
271 1
272 0
273 0
274 1
275 0
276 6
277 3
278 1
279 2
280 0
281 0
282 0
283 3
284 1
285 0
286 0
287 1
288 0
289 0
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) time_in_rfc logins
1.072e+02 1.117e-03 3.188e+00
compendium_views_info compendium_views_pr `shared_compendiums\\r`
1.779e+00 1.353e+00 5.944e+00
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-550.08 -84.09 -10.05 76.26 1134.16
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.072e+02 2.069e+01 5.183 4.16e-07 ***
time_in_rfc 1.117e-03 2.805e-04 3.981 8.73e-05 ***
logins 3.188e+00 4.191e-01 7.607 4.19e-13 ***
compendium_views_info 1.779e+00 8.278e-02 21.485 < 2e-16 ***
compendium_views_pr 1.353e+00 2.214e-01 6.112 3.25e-09 ***
`shared_compendiums\\r` 5.944e+00 3.735e+00 1.592 0.113
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 164.4 on 283 degrees of freedom
Multiple R-squared: 0.951, Adjusted R-squared: 0.9501
F-statistic: 1098 on 5 and 283 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.19640988 3.928198e-01 8.035901e-01
[2,] 0.10678854 2.135771e-01 8.932115e-01
[3,] 0.07380016 1.476003e-01 9.261998e-01
[4,] 0.03515871 7.031743e-02 9.648413e-01
[5,] 0.06569283 1.313857e-01 9.343072e-01
[6,] 0.03681130 7.362259e-02 9.631887e-01
[7,] 0.01835064 3.670128e-02 9.816494e-01
[8,] 0.01107092 2.214184e-02 9.889291e-01
[9,] 0.71688927 5.662215e-01 2.831107e-01
[10,] 0.66878019 6.624396e-01 3.312198e-01
[11,] 0.80314617 3.937077e-01 1.968538e-01
[12,] 0.74617492 5.076502e-01 2.538251e-01
[13,] 0.84554728 3.089054e-01 1.544527e-01
[14,] 0.84865842 3.026832e-01 1.513416e-01
[15,] 0.89914340 2.017132e-01 1.008566e-01
[16,] 0.87715177 2.456965e-01 1.228482e-01
[17,] 0.85940735 2.811853e-01 1.405926e-01
[18,] 0.83196832 3.360634e-01 1.680317e-01
[19,] 0.79322853 4.135429e-01 2.067715e-01
[20,] 0.75835952 4.832810e-01 2.416405e-01
[21,] 0.71593917 5.681217e-01 2.840608e-01
[22,] 0.66597396 6.680521e-01 3.340260e-01
[23,] 0.61250498 7.749900e-01 3.874950e-01
[24,] 0.55662321 8.867536e-01 4.433768e-01
[25,] 0.50078663 9.984267e-01 4.992134e-01
[26,] 0.45330838 9.066168e-01 5.466916e-01
[27,] 0.39808950 7.961790e-01 6.019105e-01
[28,] 0.35153479 7.030696e-01 6.484652e-01
[29,] 0.30485330 6.097066e-01 6.951467e-01
[30,] 0.26433963 5.286793e-01 7.356604e-01
[31,] 0.26325760 5.265152e-01 7.367424e-01
[32,] 0.22695549 4.539110e-01 7.730445e-01
[33,] 0.18893645 3.778729e-01 8.110636e-01
[34,] 0.16186960 3.237392e-01 8.381304e-01
[35,] 0.18870495 3.774099e-01 8.112951e-01
[36,] 0.16261897 3.252379e-01 8.373810e-01
[37,] 0.14957614 2.991523e-01 8.504239e-01
[38,] 0.15772494 3.154499e-01 8.422751e-01
[39,] 0.22954378 4.590876e-01 7.704562e-01
[40,] 0.22175886 4.435177e-01 7.782411e-01
[41,] 0.22597885 4.519577e-01 7.740211e-01
[42,] 0.19497222 3.899444e-01 8.050278e-01
[43,] 0.35056593 7.011319e-01 6.494341e-01
[44,] 0.32089132 6.417826e-01 6.791087e-01
[45,] 0.36877700 7.375540e-01 6.312230e-01
[46,] 0.74303302 5.139340e-01 2.569670e-01
[47,] 0.71090699 5.781860e-01 2.890930e-01
[48,] 0.67922845 6.415431e-01 3.207715e-01
[49,] 0.66108973 6.778205e-01 3.389103e-01
[50,] 0.94326880 1.134624e-01 5.673120e-02
[51,] 0.93021948 1.395610e-01 6.978052e-02
[52,] 0.92982712 1.403458e-01 7.017288e-02
[53,] 0.93316094 1.336781e-01 6.683906e-02
[54,] 0.98257794 3.484411e-02 1.742206e-02
[55,] 1.00000000 3.783459e-10 1.891729e-10
[56,] 1.00000000 6.916433e-10 3.458217e-10
[57,] 1.00000000 6.285202e-10 3.142601e-10
[58,] 1.00000000 3.085211e-10 1.542606e-10
[59,] 1.00000000 4.231350e-10 2.115675e-10
[60,] 1.00000000 9.478987e-11 4.739493e-11
[61,] 1.00000000 1.622335e-10 8.111674e-11
[62,] 1.00000000 2.130455e-10 1.065228e-10
[63,] 1.00000000 3.703640e-10 1.851820e-10
[64,] 1.00000000 1.874185e-10 9.370923e-11
[65,] 1.00000000 1.921913e-10 9.609565e-11
[66,] 1.00000000 3.483013e-10 1.741506e-10
[67,] 1.00000000 4.089784e-10 2.044892e-10
[68,] 1.00000000 7.304069e-10 3.652034e-10
[69,] 1.00000000 1.226367e-09 6.131833e-10
[70,] 1.00000000 1.215587e-09 6.077937e-10
[71,] 1.00000000 1.286094e-09 6.430471e-10
[72,] 1.00000000 1.797921e-09 8.989607e-10
[73,] 1.00000000 9.613434e-10 4.806717e-10
[74,] 1.00000000 1.444835e-09 7.224175e-10
[75,] 1.00000000 2.641091e-10 1.320545e-10
[76,] 1.00000000 2.492819e-10 1.246410e-10
[77,] 1.00000000 4.155377e-10 2.077688e-10
[78,] 1.00000000 4.241892e-10 2.120946e-10
[79,] 1.00000000 6.902207e-10 3.451103e-10
[80,] 1.00000000 1.030762e-09 5.153810e-10
[81,] 1.00000000 4.700343e-10 2.350171e-10
[82,] 1.00000000 4.591814e-11 2.295907e-11
[83,] 1.00000000 4.400960e-11 2.200480e-11
[84,] 1.00000000 2.872121e-11 1.436061e-11
[85,] 1.00000000 1.620047e-11 8.100236e-12
[86,] 1.00000000 2.879277e-11 1.439639e-11
[87,] 1.00000000 4.236556e-11 2.118278e-11
[88,] 1.00000000 4.223005e-12 2.111502e-12
[89,] 1.00000000 7.765132e-12 3.882566e-12
[90,] 1.00000000 1.948259e-12 9.741294e-13
[91,] 1.00000000 2.894553e-12 1.447277e-12
[92,] 1.00000000 4.033854e-12 2.016927e-12
[93,] 1.00000000 6.742738e-12 3.371369e-12
[94,] 1.00000000 6.248816e-12 3.124408e-12
[95,] 1.00000000 3.881947e-12 1.940973e-12
[96,] 1.00000000 6.910307e-12 3.455154e-12
[97,] 1.00000000 1.239609e-11 6.198046e-12
[98,] 1.00000000 2.081572e-11 1.040786e-11
[99,] 1.00000000 9.798392e-12 4.899196e-12
[100,] 1.00000000 3.591792e-12 1.795896e-12
[101,] 1.00000000 5.633019e-12 2.816509e-12
[102,] 1.00000000 1.019383e-11 5.096915e-12
[103,] 1.00000000 1.845892e-11 9.229462e-12
[104,] 1.00000000 1.102201e-11 5.511007e-12
[105,] 1.00000000 1.767523e-11 8.837617e-12
[106,] 1.00000000 2.694801e-11 1.347401e-11
[107,] 1.00000000 2.805668e-11 1.402834e-11
[108,] 1.00000000 2.561506e-11 1.280753e-11
[109,] 1.00000000 4.468773e-11 2.234387e-11
[110,] 1.00000000 5.934375e-11 2.967188e-11
[111,] 1.00000000 1.044122e-10 5.220612e-11
[112,] 1.00000000 1.548752e-10 7.743760e-11
[113,] 1.00000000 1.627393e-10 8.136963e-11
[114,] 1.00000000 2.565317e-10 1.282659e-10
[115,] 1.00000000 3.718779e-10 1.859390e-10
[116,] 1.00000000 6.355692e-10 3.177846e-10
[117,] 1.00000000 1.069683e-09 5.348413e-10
[118,] 1.00000000 1.779995e-09 8.899973e-10
[119,] 1.00000000 2.801295e-09 1.400648e-09
[120,] 1.00000000 1.235570e-09 6.177848e-10
[121,] 1.00000000 9.447896e-10 4.723948e-10
[122,] 1.00000000 1.543967e-09 7.719837e-10
[123,] 1.00000000 2.578738e-09 1.289369e-09
[124,] 1.00000000 4.131539e-09 2.065769e-09
[125,] 1.00000000 4.948648e-09 2.474324e-09
[126,] 1.00000000 1.291440e-09 6.457201e-10
[127,] 1.00000000 1.522640e-09 7.613199e-10
[128,] 1.00000000 1.673822e-09 8.369108e-10
[129,] 1.00000000 2.728825e-09 1.364412e-09
[130,] 1.00000000 3.519981e-09 1.759991e-09
[131,] 1.00000000 5.696610e-09 2.848305e-09
[132,] 1.00000000 6.547102e-09 3.273551e-09
[133,] 0.99999999 1.076292e-08 5.381460e-09
[134,] 1.00000000 5.338060e-09 2.669030e-09
[135,] 1.00000000 6.999299e-09 3.499650e-09
[136,] 0.99999999 1.029918e-08 5.149589e-09
[137,] 0.99999999 1.635270e-08 8.176349e-09
[138,] 0.99999999 1.150216e-08 5.751080e-09
[139,] 0.99999999 1.304125e-08 6.520625e-09
[140,] 0.99999999 2.150068e-08 1.075034e-08
[141,] 1.00000000 9.189194e-09 4.594597e-09
[142,] 0.99999999 1.511742e-08 7.558708e-09
[143,] 0.99999999 2.322532e-08 1.161266e-08
[144,] 0.99999999 2.901782e-08 1.450891e-08
[145,] 0.99999998 3.170271e-08 1.585136e-08
[146,] 1.00000000 4.860000e-09 2.430000e-09
[147,] 1.00000000 6.596456e-09 3.298228e-09
[148,] 1.00000000 8.399867e-09 4.199934e-09
[149,] 1.00000000 8.102231e-09 4.051116e-09
[150,] 1.00000000 6.533729e-09 3.266864e-09
[151,] 0.99999999 1.080588e-08 5.402940e-09
[152,] 1.00000000 9.481924e-11 4.740962e-11
[153,] 1.00000000 1.242493e-10 6.212465e-11
[154,] 1.00000000 8.666475e-11 4.333238e-11
[155,] 1.00000000 8.531343e-11 4.265672e-11
[156,] 1.00000000 1.919163e-11 9.595815e-12
[157,] 1.00000000 2.383370e-11 1.191685e-11
[158,] 1.00000000 4.415635e-11 2.207817e-11
[159,] 1.00000000 4.163371e-11 2.081685e-11
[160,] 1.00000000 1.516086e-11 7.580428e-12
[161,] 1.00000000 2.807702e-11 1.403851e-11
[162,] 1.00000000 5.188821e-11 2.594410e-11
[163,] 1.00000000 2.681122e-11 1.340561e-11
[164,] 1.00000000 4.936604e-11 2.468302e-11
[165,] 1.00000000 4.900258e-11 2.450129e-11
[166,] 1.00000000 6.872521e-11 3.436260e-11
[167,] 1.00000000 8.715897e-11 4.357949e-11
[168,] 1.00000000 1.226684e-10 6.133419e-11
[169,] 1.00000000 2.145981e-10 1.072990e-10
[170,] 1.00000000 3.933919e-10 1.966959e-10
[171,] 1.00000000 6.113371e-10 3.056685e-10
[172,] 1.00000000 9.264423e-10 4.632211e-10
[173,] 1.00000000 1.550289e-09 7.751445e-10
[174,] 1.00000000 1.741015e-09 8.705074e-10
[175,] 1.00000000 2.894911e-09 1.447455e-09
[176,] 1.00000000 5.039946e-09 2.519973e-09
[177,] 1.00000000 8.488887e-09 4.244443e-09
[178,] 0.99999999 1.353115e-08 6.765576e-09
[179,] 1.00000000 8.330039e-09 4.165019e-09
[180,] 0.99999999 1.215278e-08 6.076390e-09
[181,] 1.00000000 5.339959e-09 2.669980e-09
[182,] 1.00000000 9.293674e-09 4.646837e-09
[183,] 0.99999999 1.584134e-08 7.920671e-09
[184,] 1.00000000 8.325138e-09 4.162569e-09
[185,] 0.99999999 1.365999e-08 6.829994e-09
[186,] 0.99999999 1.863781e-08 9.318905e-09
[187,] 0.99999998 3.282192e-08 1.641096e-08
[188,] 0.99999997 5.199536e-08 2.599768e-08
[189,] 0.99999997 6.194636e-08 3.097318e-08
[190,] 0.99999995 9.643306e-08 4.821653e-08
[191,] 0.99999994 1.297571e-07 6.487854e-08
[192,] 0.99999989 2.184106e-07 1.092053e-07
[193,] 0.99999983 3.486239e-07 1.743120e-07
[194,] 0.99999971 5.848536e-07 2.924268e-07
[195,] 0.99999965 6.903086e-07 3.451543e-07
[196,] 0.99999949 1.025393e-06 5.126963e-07
[197,] 0.99999937 1.255059e-06 6.275293e-07
[198,] 0.99999895 2.100986e-06 1.050493e-06
[199,] 0.99999907 1.861994e-06 9.309970e-07
[200,] 0.99999930 1.394864e-06 6.974318e-07
[201,] 0.99999891 2.180390e-06 1.090195e-06
[202,] 0.99999833 3.331637e-06 1.665818e-06
[203,] 0.99999732 5.351035e-06 2.675518e-06
[204,] 0.99999871 2.579843e-06 1.289921e-06
[205,] 0.99999786 4.283654e-06 2.141827e-06
[206,] 0.99999696 6.088015e-06 3.044008e-06
[207,] 0.99999595 8.097954e-06 4.048977e-06
[208,] 0.99999334 1.332170e-05 6.660850e-06
[209,] 0.99998985 2.029757e-05 1.014878e-05
[210,] 0.99999656 6.888281e-06 3.444141e-06
[211,] 0.99999556 8.877977e-06 4.438989e-06
[212,] 0.99999273 1.453656e-05 7.268281e-06
[213,] 0.99999020 1.959451e-05 9.797255e-06
[214,] 0.99999697 6.053563e-06 3.026782e-06
[215,] 0.99999836 3.285616e-06 1.642808e-06
[216,] 0.99999773 4.540920e-06 2.270460e-06
[217,] 0.99999668 6.639531e-06 3.319766e-06
[218,] 0.99999506 9.883380e-06 4.941690e-06
[219,] 0.99999207 1.586755e-05 7.933776e-06
[220,] 0.99998933 2.133261e-05 1.066630e-05
[221,] 0.99998545 2.909822e-05 1.454911e-05
[222,] 0.99999945 1.099163e-06 5.495817e-07
[223,] 0.99999900 1.992410e-06 9.962049e-07
[224,] 0.99999958 8.428548e-07 4.214274e-07
[225,] 0.99999952 9.552604e-07 4.776302e-07
[226,] 0.99999912 1.765897e-06 8.829484e-07
[227,] 0.99999877 2.458967e-06 1.229484e-06
[228,] 0.99999770 4.607349e-06 2.303674e-06
[229,] 0.99999574 8.523011e-06 4.261506e-06
[230,] 0.99999313 1.374696e-05 6.873479e-06
[231,] 0.99998752 2.496009e-05 1.248004e-05
[232,] 0.99997879 4.241662e-05 2.120831e-05
[233,] 0.99997162 5.676027e-05 2.838013e-05
[234,] 0.99998556 2.887962e-05 1.443981e-05
[235,] 0.99997637 4.725027e-05 2.362513e-05
[236,] 0.99996495 7.010507e-05 3.505253e-05
[237,] 0.99994258 1.148451e-04 5.742254e-05
[238,] 0.99989805 2.039031e-04 1.019515e-04
[239,] 0.99998826 2.347013e-05 1.173506e-05
[240,] 0.99998300 3.399704e-05 1.699852e-05
[241,] 0.99997391 5.217448e-05 2.608724e-05
[242,] 0.99995463 9.073810e-05 4.536905e-05
[243,] 0.99991631 1.673889e-04 8.369445e-05
[244,] 0.99984547 3.090657e-04 1.545328e-04
[245,] 0.99977925 4.414984e-04 2.207492e-04
[246,] 0.99963289 7.342260e-04 3.671130e-04
[247,] 0.99937116 1.257689e-03 6.288443e-04
[248,] 0.99908984 1.820315e-03 9.101577e-04
[249,] 0.99841641 3.167178e-03 1.583589e-03
[250,] 0.99771418 4.571639e-03 2.285820e-03
[251,] 0.99644869 7.102612e-03 3.551306e-03
[252,] 0.99408504 1.182992e-02 5.914961e-03
[253,] 0.99148420 1.703159e-02 8.515797e-03
[254,] 0.98636642 2.726717e-02 1.363358e-02
[255,] 0.97901518 4.196963e-02 2.098482e-02
[256,] 0.97183640 5.632720e-02 2.816360e-02
[257,] 0.95845124 8.309752e-02 4.154876e-02
[258,] 0.94049749 1.190050e-01 5.950251e-02
[259,] 0.91304689 1.739062e-01 8.695311e-02
[260,] 0.92522752 1.495450e-01 7.477248e-02
[261,] 0.91756832 1.648634e-01 8.243168e-02
[262,] 0.89800592 2.039882e-01 1.019941e-01
[263,] 0.85805136 2.838973e-01 1.419486e-01
[264,] 0.83002061 3.399588e-01 1.699794e-01
[265,] 0.77535388 4.492922e-01 2.246461e-01
[266,] 0.70551574 5.889685e-01 2.944843e-01
[267,] 0.61255711 7.748858e-01 3.874429e-01
[268,] 0.50451916 9.909617e-01 4.954808e-01
[269,] 0.39301790 7.860358e-01 6.069821e-01
[270,] 0.37153493 7.430699e-01 6.284651e-01
[271,] 0.28567248 5.713450e-01 7.143275e-01
[272,] 0.47876601 9.575320e-01 5.212340e-01
> postscript(file="/var/fisher/rcomp/tmp/1m0d71356095843.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/fisher/rcomp/tmp/2b72l1356095843.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/fisher/rcomp/tmp/3w9sp1356095843.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/fisher/rcomp/tmp/4ukgk1356095843.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/fisher/rcomp/tmp/5ba7g1356095843.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
64.9089301 -178.3322975 -79.7499290 -320.1926827 -11.6325642 -105.0961244
7 8 9 10 11 12
-493.5189133 -84.0883739 -8.5529664 -116.7879067 55.7764574 -141.0722811
13 14 15 16 17 18
-21.3125969 -104.0050870 -111.6181327 -13.9392965 379.8735609 1.3838980
19 20 21 22 23 24
-256.8531226 -81.3096248 244.0475657 145.0196243 -323.9718634 22.2445700
25 26 27 28 29 30
67.2646840 18.5110126 16.2745388 -90.0215954 65.4966649 -45.8647010
31 32 33 34 35 36
-54.0133998 -71.1507690 9.8316237 -88.8551552 -31.8204636 -85.6263155
37 38 39 40 41 42
21.7639792 -102.9924746 122.9215875 -109.3041567 0.7617717 31.5949697
43 44 45 46 47 48
266.4191689 -53.1853346 -233.4102915 -207.0350376 293.9669334 -155.9525791
49 50 51 52 53 54
128.1825094 -74.2581739 -470.8442112 -87.1064717 228.7959786 385.4046771
55 56 57 58 59 60
-109.8100112 84.6026619 -129.5654591 -550.0811443 -1.7484040 -183.2221743
61 62 63 64 65 66
171.7373193 407.4636879 1134.1596358 -17.3763995 185.0084505 -250.4866083
67 68 69 70 71 72
-127.9195692 315.2678047 -53.7672350 -142.2189977 28.6702902 -256.6541621
73 74 75 76 77 78
168.4705701 -11.3908047 140.9049810 -11.2900641 -57.0307575 154.9146905
79 80 81 82 83 84
-171.6219510 103.8274044 -251.8027063 81.4762530 -334.4942157 -179.2858188
85 86 87 88 89 90
56.6143975 159.3306128 -66.6577549 80.1386125 266.8364124 -374.7762772
91 92 93 94 95 96
169.6418322 -226.5506392 -217.6533044 15.6105210 83.4767119 -310.8426564
97 98 99 100 101 102
-5.3511189 283.0484166 89.4613695 106.7822830 55.6868783 171.3436525
103 104 105 106 107 108
-186.1612396 32.5519719 -12.0870607 -60.7312204 256.6028097 258.4187242
109 110 111 112 113 114
-74.7916040 19.6633338 7.5382736 -232.9743161 -60.8133774 86.1643790
115 116 117 118 119 120
-168.7415624 -153.0612903 21.1650668 124.4346731 -32.5820807 95.3369798
121 122 123 124 125 126
151.6348978 -82.3973031 37.8505468 17.1519145 -30.8495364 27.2650708
127 128 129 130 131 132
-50.3619942 252.5883781 178.0015408 -59.8239963 1.8983247 33.1143144
133 134 135 136 137 138
-125.1251541 305.7290963 137.8045225 -152.6673669 -57.1425187 -87.2727207
139 140 141 142 143 144
40.0809054 -139.3898349 -37.5496508 229.4607564 -90.0980400 94.2283225
145 146 147 148 149 150
58.5676852 221.2817162 -74.9992256 -10.0544618 -251.5769505 22.2670670
151 152 153 154 155 156
-56.7323747 -116.9044778 121.4347430 -267.0381471 -74.4893129 152.1497302
157 158 159 160 161 162
-139.7405150 -110.6305435 -32.5663419 -440.4525877 173.7341386 -56.1052750
163 164 165 166 167 168
-111.1729599 -226.2434010 -53.6964598 2.2068471 -62.8509767 -178.2540290
169 170 171 172 173 174
35.1555111 44.9274835 -8.9963925 24.7182393 192.0452671 200.2308820
175 176 177 178 179 180
-0.7165750 -64.6847990 42.3197937 3.9238359 112.8155514 -18.3693367
181 182 183 184 185 186
52.8310485 173.5454576 -99.7057223 31.7825277 95.1358161 -111.8389846
187 188 189 190 191 192
-219.2392894 -128.5635503 237.4692619 30.1937548 47.6793673 -192.8171157
193 194 195 196 197 198
107.7710338 159.6129298 -10.4183393 99.8314699 -108.0763893 -66.3035143
199 200 201 202 203 204
-40.8210049 -27.0007674 -33.9966859 -35.7452439 -65.8060126 -7.2517602
205 206 207 208 209 210
21.7314371 -14.0622807 159.8595901 -114.6690975 -15.4618486 -10.7469805
211 212 213 214 215 216
73.4937775 223.5023713 7.0166050 -90.7191432 45.1041867 24.4012006
217 218 219 220 221 222
10.2389619 -86.9035219 -29.1856665 -61.5008583 76.2590205 404.9407732
223 224 225 226 227 228
337.3207383 68.9954799 -98.2630168 75.3630850 -35.0331624 -58.6486297
229 230 231 232 233 234
143.0666939 353.9480888 38.6799210 200.1280956 108.7075129 20.5568711
235 236 237 238 239 240
-76.8502145 9.4490530 -12.2322994 37.5003406 4.4069572 -53.8800935
241 242 243 244 245 246
82.5130325 200.0870020 17.7676366 -58.2895008 -48.9394573 8.1434659
247 248 249 250 251 252
268.1600706 -109.4147438 -96.2587683 -97.8161848 -34.3545617 24.7346102
253 254 255 256 257 258
-81.3867402 -43.6173598 -88.2556819 -97.3744156 4.1677343 -77.8475933
259 260 261 262 263 264
-30.4725899 -3.8363881 -35.8102220 74.4277657 30.3391288 83.8358677
265 266 267 268 269 270
-55.3559764 27.3073820 -10.2008591 -147.8213090 -119.1732965 -71.2912152
271 272 273 274 275 276
127.5562775 -73.0453410 -42.9113251 35.9234368 51.1955986 129.2687862
277 278 279 280 281 282
46.5373293 -47.6250547 127.2625328 -57.0493595 7.4669030 -56.2703767
283 284 285 286 287 288
13.7497825 -52.6433927 -95.0285918 376.5289732 115.7479121 92.7660084
289
16.4075216
> postscript(file="/var/fisher/rcomp/tmp/62fgh1356095843.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 64.9089301 NA
1 -178.3322975 64.9089301
2 -79.7499290 -178.3322975
3 -320.1926827 -79.7499290
4 -11.6325642 -320.1926827
5 -105.0961244 -11.6325642
6 -493.5189133 -105.0961244
7 -84.0883739 -493.5189133
8 -8.5529664 -84.0883739
9 -116.7879067 -8.5529664
10 55.7764574 -116.7879067
11 -141.0722811 55.7764574
12 -21.3125969 -141.0722811
13 -104.0050870 -21.3125969
14 -111.6181327 -104.0050870
15 -13.9392965 -111.6181327
16 379.8735609 -13.9392965
17 1.3838980 379.8735609
18 -256.8531226 1.3838980
19 -81.3096248 -256.8531226
20 244.0475657 -81.3096248
21 145.0196243 244.0475657
22 -323.9718634 145.0196243
23 22.2445700 -323.9718634
24 67.2646840 22.2445700
25 18.5110126 67.2646840
26 16.2745388 18.5110126
27 -90.0215954 16.2745388
28 65.4966649 -90.0215954
29 -45.8647010 65.4966649
30 -54.0133998 -45.8647010
31 -71.1507690 -54.0133998
32 9.8316237 -71.1507690
33 -88.8551552 9.8316237
34 -31.8204636 -88.8551552
35 -85.6263155 -31.8204636
36 21.7639792 -85.6263155
37 -102.9924746 21.7639792
38 122.9215875 -102.9924746
39 -109.3041567 122.9215875
40 0.7617717 -109.3041567
41 31.5949697 0.7617717
42 266.4191689 31.5949697
43 -53.1853346 266.4191689
44 -233.4102915 -53.1853346
45 -207.0350376 -233.4102915
46 293.9669334 -207.0350376
47 -155.9525791 293.9669334
48 128.1825094 -155.9525791
49 -74.2581739 128.1825094
50 -470.8442112 -74.2581739
51 -87.1064717 -470.8442112
52 228.7959786 -87.1064717
53 385.4046771 228.7959786
54 -109.8100112 385.4046771
55 84.6026619 -109.8100112
56 -129.5654591 84.6026619
57 -550.0811443 -129.5654591
58 -1.7484040 -550.0811443
59 -183.2221743 -1.7484040
60 171.7373193 -183.2221743
61 407.4636879 171.7373193
62 1134.1596358 407.4636879
63 -17.3763995 1134.1596358
64 185.0084505 -17.3763995
65 -250.4866083 185.0084505
66 -127.9195692 -250.4866083
67 315.2678047 -127.9195692
68 -53.7672350 315.2678047
69 -142.2189977 -53.7672350
70 28.6702902 -142.2189977
71 -256.6541621 28.6702902
72 168.4705701 -256.6541621
73 -11.3908047 168.4705701
74 140.9049810 -11.3908047
75 -11.2900641 140.9049810
76 -57.0307575 -11.2900641
77 154.9146905 -57.0307575
78 -171.6219510 154.9146905
79 103.8274044 -171.6219510
80 -251.8027063 103.8274044
81 81.4762530 -251.8027063
82 -334.4942157 81.4762530
83 -179.2858188 -334.4942157
84 56.6143975 -179.2858188
85 159.3306128 56.6143975
86 -66.6577549 159.3306128
87 80.1386125 -66.6577549
88 266.8364124 80.1386125
89 -374.7762772 266.8364124
90 169.6418322 -374.7762772
91 -226.5506392 169.6418322
92 -217.6533044 -226.5506392
93 15.6105210 -217.6533044
94 83.4767119 15.6105210
95 -310.8426564 83.4767119
96 -5.3511189 -310.8426564
97 283.0484166 -5.3511189
98 89.4613695 283.0484166
99 106.7822830 89.4613695
100 55.6868783 106.7822830
101 171.3436525 55.6868783
102 -186.1612396 171.3436525
103 32.5519719 -186.1612396
104 -12.0870607 32.5519719
105 -60.7312204 -12.0870607
106 256.6028097 -60.7312204
107 258.4187242 256.6028097
108 -74.7916040 258.4187242
109 19.6633338 -74.7916040
110 7.5382736 19.6633338
111 -232.9743161 7.5382736
112 -60.8133774 -232.9743161
113 86.1643790 -60.8133774
114 -168.7415624 86.1643790
115 -153.0612903 -168.7415624
116 21.1650668 -153.0612903
117 124.4346731 21.1650668
118 -32.5820807 124.4346731
119 95.3369798 -32.5820807
120 151.6348978 95.3369798
121 -82.3973031 151.6348978
122 37.8505468 -82.3973031
123 17.1519145 37.8505468
124 -30.8495364 17.1519145
125 27.2650708 -30.8495364
126 -50.3619942 27.2650708
127 252.5883781 -50.3619942
128 178.0015408 252.5883781
129 -59.8239963 178.0015408
130 1.8983247 -59.8239963
131 33.1143144 1.8983247
132 -125.1251541 33.1143144
133 305.7290963 -125.1251541
134 137.8045225 305.7290963
135 -152.6673669 137.8045225
136 -57.1425187 -152.6673669
137 -87.2727207 -57.1425187
138 40.0809054 -87.2727207
139 -139.3898349 40.0809054
140 -37.5496508 -139.3898349
141 229.4607564 -37.5496508
142 -90.0980400 229.4607564
143 94.2283225 -90.0980400
144 58.5676852 94.2283225
145 221.2817162 58.5676852
146 -74.9992256 221.2817162
147 -10.0544618 -74.9992256
148 -251.5769505 -10.0544618
149 22.2670670 -251.5769505
150 -56.7323747 22.2670670
151 -116.9044778 -56.7323747
152 121.4347430 -116.9044778
153 -267.0381471 121.4347430
154 -74.4893129 -267.0381471
155 152.1497302 -74.4893129
156 -139.7405150 152.1497302
157 -110.6305435 -139.7405150
158 -32.5663419 -110.6305435
159 -440.4525877 -32.5663419
160 173.7341386 -440.4525877
161 -56.1052750 173.7341386
162 -111.1729599 -56.1052750
163 -226.2434010 -111.1729599
164 -53.6964598 -226.2434010
165 2.2068471 -53.6964598
166 -62.8509767 2.2068471
167 -178.2540290 -62.8509767
168 35.1555111 -178.2540290
169 44.9274835 35.1555111
170 -8.9963925 44.9274835
171 24.7182393 -8.9963925
172 192.0452671 24.7182393
173 200.2308820 192.0452671
174 -0.7165750 200.2308820
175 -64.6847990 -0.7165750
176 42.3197937 -64.6847990
177 3.9238359 42.3197937
178 112.8155514 3.9238359
179 -18.3693367 112.8155514
180 52.8310485 -18.3693367
181 173.5454576 52.8310485
182 -99.7057223 173.5454576
183 31.7825277 -99.7057223
184 95.1358161 31.7825277
185 -111.8389846 95.1358161
186 -219.2392894 -111.8389846
187 -128.5635503 -219.2392894
188 237.4692619 -128.5635503
189 30.1937548 237.4692619
190 47.6793673 30.1937548
191 -192.8171157 47.6793673
192 107.7710338 -192.8171157
193 159.6129298 107.7710338
194 -10.4183393 159.6129298
195 99.8314699 -10.4183393
196 -108.0763893 99.8314699
197 -66.3035143 -108.0763893
198 -40.8210049 -66.3035143
199 -27.0007674 -40.8210049
200 -33.9966859 -27.0007674
201 -35.7452439 -33.9966859
202 -65.8060126 -35.7452439
203 -7.2517602 -65.8060126
204 21.7314371 -7.2517602
205 -14.0622807 21.7314371
206 159.8595901 -14.0622807
207 -114.6690975 159.8595901
208 -15.4618486 -114.6690975
209 -10.7469805 -15.4618486
210 73.4937775 -10.7469805
211 223.5023713 73.4937775
212 7.0166050 223.5023713
213 -90.7191432 7.0166050
214 45.1041867 -90.7191432
215 24.4012006 45.1041867
216 10.2389619 24.4012006
217 -86.9035219 10.2389619
218 -29.1856665 -86.9035219
219 -61.5008583 -29.1856665
220 76.2590205 -61.5008583
221 404.9407732 76.2590205
222 337.3207383 404.9407732
223 68.9954799 337.3207383
224 -98.2630168 68.9954799
225 75.3630850 -98.2630168
226 -35.0331624 75.3630850
227 -58.6486297 -35.0331624
228 143.0666939 -58.6486297
229 353.9480888 143.0666939
230 38.6799210 353.9480888
231 200.1280956 38.6799210
232 108.7075129 200.1280956
233 20.5568711 108.7075129
234 -76.8502145 20.5568711
235 9.4490530 -76.8502145
236 -12.2322994 9.4490530
237 37.5003406 -12.2322994
238 4.4069572 37.5003406
239 -53.8800935 4.4069572
240 82.5130325 -53.8800935
241 200.0870020 82.5130325
242 17.7676366 200.0870020
243 -58.2895008 17.7676366
244 -48.9394573 -58.2895008
245 8.1434659 -48.9394573
246 268.1600706 8.1434659
247 -109.4147438 268.1600706
248 -96.2587683 -109.4147438
249 -97.8161848 -96.2587683
250 -34.3545617 -97.8161848
251 24.7346102 -34.3545617
252 -81.3867402 24.7346102
253 -43.6173598 -81.3867402
254 -88.2556819 -43.6173598
255 -97.3744156 -88.2556819
256 4.1677343 -97.3744156
257 -77.8475933 4.1677343
258 -30.4725899 -77.8475933
259 -3.8363881 -30.4725899
260 -35.8102220 -3.8363881
261 74.4277657 -35.8102220
262 30.3391288 74.4277657
263 83.8358677 30.3391288
264 -55.3559764 83.8358677
265 27.3073820 -55.3559764
266 -10.2008591 27.3073820
267 -147.8213090 -10.2008591
268 -119.1732965 -147.8213090
269 -71.2912152 -119.1732965
270 127.5562775 -71.2912152
271 -73.0453410 127.5562775
272 -42.9113251 -73.0453410
273 35.9234368 -42.9113251
274 51.1955986 35.9234368
275 129.2687862 51.1955986
276 46.5373293 129.2687862
277 -47.6250547 46.5373293
278 127.2625328 -47.6250547
279 -57.0493595 127.2625328
280 7.4669030 -57.0493595
281 -56.2703767 7.4669030
282 13.7497825 -56.2703767
283 -52.6433927 13.7497825
284 -95.0285918 -52.6433927
285 376.5289732 -95.0285918
286 115.7479121 376.5289732
287 92.7660084 115.7479121
288 16.4075216 92.7660084
289 NA 16.4075216
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] -178.3322975 64.9089301
[2,] -79.7499290 -178.3322975
[3,] -320.1926827 -79.7499290
[4,] -11.6325642 -320.1926827
[5,] -105.0961244 -11.6325642
[6,] -493.5189133 -105.0961244
[7,] -84.0883739 -493.5189133
[8,] -8.5529664 -84.0883739
[9,] -116.7879067 -8.5529664
[10,] 55.7764574 -116.7879067
[11,] -141.0722811 55.7764574
[12,] -21.3125969 -141.0722811
[13,] -104.0050870 -21.3125969
[14,] -111.6181327 -104.0050870
[15,] -13.9392965 -111.6181327
[16,] 379.8735609 -13.9392965
[17,] 1.3838980 379.8735609
[18,] -256.8531226 1.3838980
[19,] -81.3096248 -256.8531226
[20,] 244.0475657 -81.3096248
[21,] 145.0196243 244.0475657
[22,] -323.9718634 145.0196243
[23,] 22.2445700 -323.9718634
[24,] 67.2646840 22.2445700
[25,] 18.5110126 67.2646840
[26,] 16.2745388 18.5110126
[27,] -90.0215954 16.2745388
[28,] 65.4966649 -90.0215954
[29,] -45.8647010 65.4966649
[30,] -54.0133998 -45.8647010
[31,] -71.1507690 -54.0133998
[32,] 9.8316237 -71.1507690
[33,] -88.8551552 9.8316237
[34,] -31.8204636 -88.8551552
[35,] -85.6263155 -31.8204636
[36,] 21.7639792 -85.6263155
[37,] -102.9924746 21.7639792
[38,] 122.9215875 -102.9924746
[39,] -109.3041567 122.9215875
[40,] 0.7617717 -109.3041567
[41,] 31.5949697 0.7617717
[42,] 266.4191689 31.5949697
[43,] -53.1853346 266.4191689
[44,] -233.4102915 -53.1853346
[45,] -207.0350376 -233.4102915
[46,] 293.9669334 -207.0350376
[47,] -155.9525791 293.9669334
[48,] 128.1825094 -155.9525791
[49,] -74.2581739 128.1825094
[50,] -470.8442112 -74.2581739
[51,] -87.1064717 -470.8442112
[52,] 228.7959786 -87.1064717
[53,] 385.4046771 228.7959786
[54,] -109.8100112 385.4046771
[55,] 84.6026619 -109.8100112
[56,] -129.5654591 84.6026619
[57,] -550.0811443 -129.5654591
[58,] -1.7484040 -550.0811443
[59,] -183.2221743 -1.7484040
[60,] 171.7373193 -183.2221743
[61,] 407.4636879 171.7373193
[62,] 1134.1596358 407.4636879
[63,] -17.3763995 1134.1596358
[64,] 185.0084505 -17.3763995
[65,] -250.4866083 185.0084505
[66,] -127.9195692 -250.4866083
[67,] 315.2678047 -127.9195692
[68,] -53.7672350 315.2678047
[69,] -142.2189977 -53.7672350
[70,] 28.6702902 -142.2189977
[71,] -256.6541621 28.6702902
[72,] 168.4705701 -256.6541621
[73,] -11.3908047 168.4705701
[74,] 140.9049810 -11.3908047
[75,] -11.2900641 140.9049810
[76,] -57.0307575 -11.2900641
[77,] 154.9146905 -57.0307575
[78,] -171.6219510 154.9146905
[79,] 103.8274044 -171.6219510
[80,] -251.8027063 103.8274044
[81,] 81.4762530 -251.8027063
[82,] -334.4942157 81.4762530
[83,] -179.2858188 -334.4942157
[84,] 56.6143975 -179.2858188
[85,] 159.3306128 56.6143975
[86,] -66.6577549 159.3306128
[87,] 80.1386125 -66.6577549
[88,] 266.8364124 80.1386125
[89,] -374.7762772 266.8364124
[90,] 169.6418322 -374.7762772
[91,] -226.5506392 169.6418322
[92,] -217.6533044 -226.5506392
[93,] 15.6105210 -217.6533044
[94,] 83.4767119 15.6105210
[95,] -310.8426564 83.4767119
[96,] -5.3511189 -310.8426564
[97,] 283.0484166 -5.3511189
[98,] 89.4613695 283.0484166
[99,] 106.7822830 89.4613695
[100,] 55.6868783 106.7822830
[101,] 171.3436525 55.6868783
[102,] -186.1612396 171.3436525
[103,] 32.5519719 -186.1612396
[104,] -12.0870607 32.5519719
[105,] -60.7312204 -12.0870607
[106,] 256.6028097 -60.7312204
[107,] 258.4187242 256.6028097
[108,] -74.7916040 258.4187242
[109,] 19.6633338 -74.7916040
[110,] 7.5382736 19.6633338
[111,] -232.9743161 7.5382736
[112,] -60.8133774 -232.9743161
[113,] 86.1643790 -60.8133774
[114,] -168.7415624 86.1643790
[115,] -153.0612903 -168.7415624
[116,] 21.1650668 -153.0612903
[117,] 124.4346731 21.1650668
[118,] -32.5820807 124.4346731
[119,] 95.3369798 -32.5820807
[120,] 151.6348978 95.3369798
[121,] -82.3973031 151.6348978
[122,] 37.8505468 -82.3973031
[123,] 17.1519145 37.8505468
[124,] -30.8495364 17.1519145
[125,] 27.2650708 -30.8495364
[126,] -50.3619942 27.2650708
[127,] 252.5883781 -50.3619942
[128,] 178.0015408 252.5883781
[129,] -59.8239963 178.0015408
[130,] 1.8983247 -59.8239963
[131,] 33.1143144 1.8983247
[132,] -125.1251541 33.1143144
[133,] 305.7290963 -125.1251541
[134,] 137.8045225 305.7290963
[135,] -152.6673669 137.8045225
[136,] -57.1425187 -152.6673669
[137,] -87.2727207 -57.1425187
[138,] 40.0809054 -87.2727207
[139,] -139.3898349 40.0809054
[140,] -37.5496508 -139.3898349
[141,] 229.4607564 -37.5496508
[142,] -90.0980400 229.4607564
[143,] 94.2283225 -90.0980400
[144,] 58.5676852 94.2283225
[145,] 221.2817162 58.5676852
[146,] -74.9992256 221.2817162
[147,] -10.0544618 -74.9992256
[148,] -251.5769505 -10.0544618
[149,] 22.2670670 -251.5769505
[150,] -56.7323747 22.2670670
[151,] -116.9044778 -56.7323747
[152,] 121.4347430 -116.9044778
[153,] -267.0381471 121.4347430
[154,] -74.4893129 -267.0381471
[155,] 152.1497302 -74.4893129
[156,] -139.7405150 152.1497302
[157,] -110.6305435 -139.7405150
[158,] -32.5663419 -110.6305435
[159,] -440.4525877 -32.5663419
[160,] 173.7341386 -440.4525877
[161,] -56.1052750 173.7341386
[162,] -111.1729599 -56.1052750
[163,] -226.2434010 -111.1729599
[164,] -53.6964598 -226.2434010
[165,] 2.2068471 -53.6964598
[166,] -62.8509767 2.2068471
[167,] -178.2540290 -62.8509767
[168,] 35.1555111 -178.2540290
[169,] 44.9274835 35.1555111
[170,] -8.9963925 44.9274835
[171,] 24.7182393 -8.9963925
[172,] 192.0452671 24.7182393
[173,] 200.2308820 192.0452671
[174,] -0.7165750 200.2308820
[175,] -64.6847990 -0.7165750
[176,] 42.3197937 -64.6847990
[177,] 3.9238359 42.3197937
[178,] 112.8155514 3.9238359
[179,] -18.3693367 112.8155514
[180,] 52.8310485 -18.3693367
[181,] 173.5454576 52.8310485
[182,] -99.7057223 173.5454576
[183,] 31.7825277 -99.7057223
[184,] 95.1358161 31.7825277
[185,] -111.8389846 95.1358161
[186,] -219.2392894 -111.8389846
[187,] -128.5635503 -219.2392894
[188,] 237.4692619 -128.5635503
[189,] 30.1937548 237.4692619
[190,] 47.6793673 30.1937548
[191,] -192.8171157 47.6793673
[192,] 107.7710338 -192.8171157
[193,] 159.6129298 107.7710338
[194,] -10.4183393 159.6129298
[195,] 99.8314699 -10.4183393
[196,] -108.0763893 99.8314699
[197,] -66.3035143 -108.0763893
[198,] -40.8210049 -66.3035143
[199,] -27.0007674 -40.8210049
[200,] -33.9966859 -27.0007674
[201,] -35.7452439 -33.9966859
[202,] -65.8060126 -35.7452439
[203,] -7.2517602 -65.8060126
[204,] 21.7314371 -7.2517602
[205,] -14.0622807 21.7314371
[206,] 159.8595901 -14.0622807
[207,] -114.6690975 159.8595901
[208,] -15.4618486 -114.6690975
[209,] -10.7469805 -15.4618486
[210,] 73.4937775 -10.7469805
[211,] 223.5023713 73.4937775
[212,] 7.0166050 223.5023713
[213,] -90.7191432 7.0166050
[214,] 45.1041867 -90.7191432
[215,] 24.4012006 45.1041867
[216,] 10.2389619 24.4012006
[217,] -86.9035219 10.2389619
[218,] -29.1856665 -86.9035219
[219,] -61.5008583 -29.1856665
[220,] 76.2590205 -61.5008583
[221,] 404.9407732 76.2590205
[222,] 337.3207383 404.9407732
[223,] 68.9954799 337.3207383
[224,] -98.2630168 68.9954799
[225,] 75.3630850 -98.2630168
[226,] -35.0331624 75.3630850
[227,] -58.6486297 -35.0331624
[228,] 143.0666939 -58.6486297
[229,] 353.9480888 143.0666939
[230,] 38.6799210 353.9480888
[231,] 200.1280956 38.6799210
[232,] 108.7075129 200.1280956
[233,] 20.5568711 108.7075129
[234,] -76.8502145 20.5568711
[235,] 9.4490530 -76.8502145
[236,] -12.2322994 9.4490530
[237,] 37.5003406 -12.2322994
[238,] 4.4069572 37.5003406
[239,] -53.8800935 4.4069572
[240,] 82.5130325 -53.8800935
[241,] 200.0870020 82.5130325
[242,] 17.7676366 200.0870020
[243,] -58.2895008 17.7676366
[244,] -48.9394573 -58.2895008
[245,] 8.1434659 -48.9394573
[246,] 268.1600706 8.1434659
[247,] -109.4147438 268.1600706
[248,] -96.2587683 -109.4147438
[249,] -97.8161848 -96.2587683
[250,] -34.3545617 -97.8161848
[251,] 24.7346102 -34.3545617
[252,] -81.3867402 24.7346102
[253,] -43.6173598 -81.3867402
[254,] -88.2556819 -43.6173598
[255,] -97.3744156 -88.2556819
[256,] 4.1677343 -97.3744156
[257,] -77.8475933 4.1677343
[258,] -30.4725899 -77.8475933
[259,] -3.8363881 -30.4725899
[260,] -35.8102220 -3.8363881
[261,] 74.4277657 -35.8102220
[262,] 30.3391288 74.4277657
[263,] 83.8358677 30.3391288
[264,] -55.3559764 83.8358677
[265,] 27.3073820 -55.3559764
[266,] -10.2008591 27.3073820
[267,] -147.8213090 -10.2008591
[268,] -119.1732965 -147.8213090
[269,] -71.2912152 -119.1732965
[270,] 127.5562775 -71.2912152
[271,] -73.0453410 127.5562775
[272,] -42.9113251 -73.0453410
[273,] 35.9234368 -42.9113251
[274,] 51.1955986 35.9234368
[275,] 129.2687862 51.1955986
[276,] 46.5373293 129.2687862
[277,] -47.6250547 46.5373293
[278,] 127.2625328 -47.6250547
[279,] -57.0493595 127.2625328
[280,] 7.4669030 -57.0493595
[281,] -56.2703767 7.4669030
[282,] 13.7497825 -56.2703767
[283,] -52.6433927 13.7497825
[284,] -95.0285918 -52.6433927
[285,] 376.5289732 -95.0285918
[286,] 115.7479121 376.5289732
[287,] 92.7660084 115.7479121
[288,] 16.4075216 92.7660084
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 -178.3322975 64.9089301
2 -79.7499290 -178.3322975
3 -320.1926827 -79.7499290
4 -11.6325642 -320.1926827
5 -105.0961244 -11.6325642
6 -493.5189133 -105.0961244
7 -84.0883739 -493.5189133
8 -8.5529664 -84.0883739
9 -116.7879067 -8.5529664
10 55.7764574 -116.7879067
11 -141.0722811 55.7764574
12 -21.3125969 -141.0722811
13 -104.0050870 -21.3125969
14 -111.6181327 -104.0050870
15 -13.9392965 -111.6181327
16 379.8735609 -13.9392965
17 1.3838980 379.8735609
18 -256.8531226 1.3838980
19 -81.3096248 -256.8531226
20 244.0475657 -81.3096248
21 145.0196243 244.0475657
22 -323.9718634 145.0196243
23 22.2445700 -323.9718634
24 67.2646840 22.2445700
25 18.5110126 67.2646840
26 16.2745388 18.5110126
27 -90.0215954 16.2745388
28 65.4966649 -90.0215954
29 -45.8647010 65.4966649
30 -54.0133998 -45.8647010
31 -71.1507690 -54.0133998
32 9.8316237 -71.1507690
33 -88.8551552 9.8316237
34 -31.8204636 -88.8551552
35 -85.6263155 -31.8204636
36 21.7639792 -85.6263155
37 -102.9924746 21.7639792
38 122.9215875 -102.9924746
39 -109.3041567 122.9215875
40 0.7617717 -109.3041567
41 31.5949697 0.7617717
42 266.4191689 31.5949697
43 -53.1853346 266.4191689
44 -233.4102915 -53.1853346
45 -207.0350376 -233.4102915
46 293.9669334 -207.0350376
47 -155.9525791 293.9669334
48 128.1825094 -155.9525791
49 -74.2581739 128.1825094
50 -470.8442112 -74.2581739
51 -87.1064717 -470.8442112
52 228.7959786 -87.1064717
53 385.4046771 228.7959786
54 -109.8100112 385.4046771
55 84.6026619 -109.8100112
56 -129.5654591 84.6026619
57 -550.0811443 -129.5654591
58 -1.7484040 -550.0811443
59 -183.2221743 -1.7484040
60 171.7373193 -183.2221743
61 407.4636879 171.7373193
62 1134.1596358 407.4636879
63 -17.3763995 1134.1596358
64 185.0084505 -17.3763995
65 -250.4866083 185.0084505
66 -127.9195692 -250.4866083
67 315.2678047 -127.9195692
68 -53.7672350 315.2678047
69 -142.2189977 -53.7672350
70 28.6702902 -142.2189977
71 -256.6541621 28.6702902
72 168.4705701 -256.6541621
73 -11.3908047 168.4705701
74 140.9049810 -11.3908047
75 -11.2900641 140.9049810
76 -57.0307575 -11.2900641
77 154.9146905 -57.0307575
78 -171.6219510 154.9146905
79 103.8274044 -171.6219510
80 -251.8027063 103.8274044
81 81.4762530 -251.8027063
82 -334.4942157 81.4762530
83 -179.2858188 -334.4942157
84 56.6143975 -179.2858188
85 159.3306128 56.6143975
86 -66.6577549 159.3306128
87 80.1386125 -66.6577549
88 266.8364124 80.1386125
89 -374.7762772 266.8364124
90 169.6418322 -374.7762772
91 -226.5506392 169.6418322
92 -217.6533044 -226.5506392
93 15.6105210 -217.6533044
94 83.4767119 15.6105210
95 -310.8426564 83.4767119
96 -5.3511189 -310.8426564
97 283.0484166 -5.3511189
98 89.4613695 283.0484166
99 106.7822830 89.4613695
100 55.6868783 106.7822830
101 171.3436525 55.6868783
102 -186.1612396 171.3436525
103 32.5519719 -186.1612396
104 -12.0870607 32.5519719
105 -60.7312204 -12.0870607
106 256.6028097 -60.7312204
107 258.4187242 256.6028097
108 -74.7916040 258.4187242
109 19.6633338 -74.7916040
110 7.5382736 19.6633338
111 -232.9743161 7.5382736
112 -60.8133774 -232.9743161
113 86.1643790 -60.8133774
114 -168.7415624 86.1643790
115 -153.0612903 -168.7415624
116 21.1650668 -153.0612903
117 124.4346731 21.1650668
118 -32.5820807 124.4346731
119 95.3369798 -32.5820807
120 151.6348978 95.3369798
121 -82.3973031 151.6348978
122 37.8505468 -82.3973031
123 17.1519145 37.8505468
124 -30.8495364 17.1519145
125 27.2650708 -30.8495364
126 -50.3619942 27.2650708
127 252.5883781 -50.3619942
128 178.0015408 252.5883781
129 -59.8239963 178.0015408
130 1.8983247 -59.8239963
131 33.1143144 1.8983247
132 -125.1251541 33.1143144
133 305.7290963 -125.1251541
134 137.8045225 305.7290963
135 -152.6673669 137.8045225
136 -57.1425187 -152.6673669
137 -87.2727207 -57.1425187
138 40.0809054 -87.2727207
139 -139.3898349 40.0809054
140 -37.5496508 -139.3898349
141 229.4607564 -37.5496508
142 -90.0980400 229.4607564
143 94.2283225 -90.0980400
144 58.5676852 94.2283225
145 221.2817162 58.5676852
146 -74.9992256 221.2817162
147 -10.0544618 -74.9992256
148 -251.5769505 -10.0544618
149 22.2670670 -251.5769505
150 -56.7323747 22.2670670
151 -116.9044778 -56.7323747
152 121.4347430 -116.9044778
153 -267.0381471 121.4347430
154 -74.4893129 -267.0381471
155 152.1497302 -74.4893129
156 -139.7405150 152.1497302
157 -110.6305435 -139.7405150
158 -32.5663419 -110.6305435
159 -440.4525877 -32.5663419
160 173.7341386 -440.4525877
161 -56.1052750 173.7341386
162 -111.1729599 -56.1052750
163 -226.2434010 -111.1729599
164 -53.6964598 -226.2434010
165 2.2068471 -53.6964598
166 -62.8509767 2.2068471
167 -178.2540290 -62.8509767
168 35.1555111 -178.2540290
169 44.9274835 35.1555111
170 -8.9963925 44.9274835
171 24.7182393 -8.9963925
172 192.0452671 24.7182393
173 200.2308820 192.0452671
174 -0.7165750 200.2308820
175 -64.6847990 -0.7165750
176 42.3197937 -64.6847990
177 3.9238359 42.3197937
178 112.8155514 3.9238359
179 -18.3693367 112.8155514
180 52.8310485 -18.3693367
181 173.5454576 52.8310485
182 -99.7057223 173.5454576
183 31.7825277 -99.7057223
184 95.1358161 31.7825277
185 -111.8389846 95.1358161
186 -219.2392894 -111.8389846
187 -128.5635503 -219.2392894
188 237.4692619 -128.5635503
189 30.1937548 237.4692619
190 47.6793673 30.1937548
191 -192.8171157 47.6793673
192 107.7710338 -192.8171157
193 159.6129298 107.7710338
194 -10.4183393 159.6129298
195 99.8314699 -10.4183393
196 -108.0763893 99.8314699
197 -66.3035143 -108.0763893
198 -40.8210049 -66.3035143
199 -27.0007674 -40.8210049
200 -33.9966859 -27.0007674
201 -35.7452439 -33.9966859
202 -65.8060126 -35.7452439
203 -7.2517602 -65.8060126
204 21.7314371 -7.2517602
205 -14.0622807 21.7314371
206 159.8595901 -14.0622807
207 -114.6690975 159.8595901
208 -15.4618486 -114.6690975
209 -10.7469805 -15.4618486
210 73.4937775 -10.7469805
211 223.5023713 73.4937775
212 7.0166050 223.5023713
213 -90.7191432 7.0166050
214 45.1041867 -90.7191432
215 24.4012006 45.1041867
216 10.2389619 24.4012006
217 -86.9035219 10.2389619
218 -29.1856665 -86.9035219
219 -61.5008583 -29.1856665
220 76.2590205 -61.5008583
221 404.9407732 76.2590205
222 337.3207383 404.9407732
223 68.9954799 337.3207383
224 -98.2630168 68.9954799
225 75.3630850 -98.2630168
226 -35.0331624 75.3630850
227 -58.6486297 -35.0331624
228 143.0666939 -58.6486297
229 353.9480888 143.0666939
230 38.6799210 353.9480888
231 200.1280956 38.6799210
232 108.7075129 200.1280956
233 20.5568711 108.7075129
234 -76.8502145 20.5568711
235 9.4490530 -76.8502145
236 -12.2322994 9.4490530
237 37.5003406 -12.2322994
238 4.4069572 37.5003406
239 -53.8800935 4.4069572
240 82.5130325 -53.8800935
241 200.0870020 82.5130325
242 17.7676366 200.0870020
243 -58.2895008 17.7676366
244 -48.9394573 -58.2895008
245 8.1434659 -48.9394573
246 268.1600706 8.1434659
247 -109.4147438 268.1600706
248 -96.2587683 -109.4147438
249 -97.8161848 -96.2587683
250 -34.3545617 -97.8161848
251 24.7346102 -34.3545617
252 -81.3867402 24.7346102
253 -43.6173598 -81.3867402
254 -88.2556819 -43.6173598
255 -97.3744156 -88.2556819
256 4.1677343 -97.3744156
257 -77.8475933 4.1677343
258 -30.4725899 -77.8475933
259 -3.8363881 -30.4725899
260 -35.8102220 -3.8363881
261 74.4277657 -35.8102220
262 30.3391288 74.4277657
263 83.8358677 30.3391288
264 -55.3559764 83.8358677
265 27.3073820 -55.3559764
266 -10.2008591 27.3073820
267 -147.8213090 -10.2008591
268 -119.1732965 -147.8213090
269 -71.2912152 -119.1732965
270 127.5562775 -71.2912152
271 -73.0453410 127.5562775
272 -42.9113251 -73.0453410
273 35.9234368 -42.9113251
274 51.1955986 35.9234368
275 129.2687862 51.1955986
276 46.5373293 129.2687862
277 -47.6250547 46.5373293
278 127.2625328 -47.6250547
279 -57.0493595 127.2625328
280 7.4669030 -57.0493595
281 -56.2703767 7.4669030
282 13.7497825 -56.2703767
283 -52.6433927 13.7497825
284 -95.0285918 -52.6433927
285 376.5289732 -95.0285918
286 115.7479121 376.5289732
287 92.7660084 115.7479121
288 16.4075216 92.7660084
> 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/fisher/rcomp/tmp/7lrlu1356095843.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/fisher/rcomp/tmp/8lnot1356095843.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/fisher/rcomp/tmp/9qri11356095843.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/fisher/rcomp/tmp/100glj1356095843.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/fisher/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/fisher/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/fisher/rcomp/tmp/11p7p21356095843.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/fisher/rcomp/tmp/126mza1356095844.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/fisher/rcomp/tmp/13uxss1356095844.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/fisher/rcomp/tmp/14ypl51356095844.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/fisher/rcomp/tmp/15rkkb1356095844.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/fisher/rcomp/tmp/1627x91356095844.tab")
+ }
>
> try(system("convert tmp/1m0d71356095843.ps tmp/1m0d71356095843.png",intern=TRUE))
character(0)
> try(system("convert tmp/2b72l1356095843.ps tmp/2b72l1356095843.png",intern=TRUE))
character(0)
> try(system("convert tmp/3w9sp1356095843.ps tmp/3w9sp1356095843.png",intern=TRUE))
character(0)
> try(system("convert tmp/4ukgk1356095843.ps tmp/4ukgk1356095843.png",intern=TRUE))
character(0)
> try(system("convert tmp/5ba7g1356095843.ps tmp/5ba7g1356095843.png",intern=TRUE))
character(0)
> try(system("convert tmp/62fgh1356095843.ps tmp/62fgh1356095843.png",intern=TRUE))
character(0)
> try(system("convert tmp/7lrlu1356095843.ps tmp/7lrlu1356095843.png",intern=TRUE))
character(0)
> try(system("convert tmp/8lnot1356095843.ps tmp/8lnot1356095843.png",intern=TRUE))
character(0)
> try(system("convert tmp/9qri11356095843.ps tmp/9qri11356095843.png",intern=TRUE))
character(0)
> try(system("convert tmp/100glj1356095843.ps tmp/100glj1356095843.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
11.709 1.802 13.550