R version 2.9.0 (2009-04-17)
Copyright (C) 2009 The R Foundation for Statistical Computing
ISBN 3-900051-07-0

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 <- c(10926126
+ ,204514194
+ ,20655116
+ ,7365385
+ ,3918657
+ ,38494822
+ ,90838469
+ ,8652220
+ ,47482511
+ ,5438191
+ ,5986958
+ ,59447386
+ ,37113055
+ ,23669000
+ ,5572851
+ ,7784176
+ ,5270786
+ ,43583746
+ ,19903814
+ ,7799064
+ ,35989363
+ ,29552113
+ ,21506498
+ ,7859684
+ ,13006528
+ ,89426565
+ ,40100910
+ ,38650820
+ ,17818636
+ ,51985948
+ ,24571878
+ ,29266469
+ ,45758250
+ ,22712439
+ ,86459749
+ ,49703065
+ ,65187751
+ ,39091325
+ ,34955804
+ ,33411147
+ ,15575837
+ ,36116313
+ ,27405375
+ ,55382004
+ ,7478767
+ ,72881801
+ ,37479848
+ ,27030251
+ ,109851111
+ ,28356456
+ ,48560722
+ ,23550558
+ ,52319234
+ ,105318425
+ ,17345510
+ ,71519979
+ ,59398481
+ ,36287185
+ ,44559709
+ ,114256076
+ ,23174724
+ ,38795625
+ ,94992548
+ ,48389051
+ ,46785833
+ ,56063795
+ ,27622103
+ ,11848002
+ ,55125816
+ ,42484133
+ ,49639784
+ ,30006747
+ ,19776143
+ ,41176416
+ ,34933018
+ ,39696259
+ ,38610419
+ ,44404281
+ ,44628459
+ ,20428394
+ ,33494104
+ ,15246622
+ ,10078039
+ ,36293305
+ ,11594321
+ ,34818494
+ ,36718379
+ ,19690322
+ ,44432615
+ ,36047082
+ ,40379563
+ ,23458900
+ ,24477773
+ ,21472076
+ ,27623825
+ ,42834376
+ ,21692408
+ ,2490957
+ ,5804447
+ ,4149313
+ ,22772412
+ ,8743945
+ ,6818743
+ ,30796414
+ ,17697315
+ ,62514975
+ ,17618542
+ ,37171201
+ ,29747133
+ ,39734547
+ ,18636718
+ ,30751381
+ ,58198919
+ ,35172258
+ ,30568494
+ ,32732200
+ ,33931175
+ ,38495426
+ ,39639276
+ ,18386312
+ ,49064295
+ ,34517203
+ ,54053990
+ ,56632146
+ ,46131933
+ ,31138610
+ ,37573144
+ ,32922160
+ ,71766977
+ ,21574297
+ ,46733768
+ ,36375026
+ ,50499581
+ ,44157910
+ ,26780541
+ ,71675809
+ ,17411180
+ ,138402202
+ ,89520475
+ ,53764838
+ ,79645460
+ ,45232765
+ ,60335913
+ ,48862499
+ ,56443228
+ ,149401818
+ ,55893936
+ ,39185258
+ ,54999796
+ ,62938759
+ ,130610415
+ ,3578335
+ ,119330610
+ ,43612479
+ ,67573324
+ ,94235304
+ ,52529046
+ ,38148241
+ ,41788652
+ ,21748127
+ ,16815679
+ ,29140007
+ ,80020386
+ ,10008346
+ ,13150180
+ ,33070836
+ ,44343752
+ ,38761463
+ ,63242438
+ ,10776650
+ ,13079087
+ ,42254065
+ ,45334695
+ ,75653353
+ ,37723883
+ ,35383003
+ ,42169175
+ ,37836247
+ ,16157265
+ ,25483137
+ ,24758200
+ ,25931415
+ ,35341058
+ ,59793935
+ ,46240034
+ ,25935511
+ ,29348057
+ ,16621841
+ ,34421545
+ ,30201523
+ ,68504357
+ ,21760376
+ ,68420735
+ ,34576094
+ ,80009790
+ ,23896896
+ ,43959624
+ ,31004430
+ ,96339703
+ ,78115427
+ ,13773684
+ ,79714473
+ ,30300608
+ ,47400000
+ ,29685182
+ ,72574056
+ ,26667187
+ ,56654962
+ ,34032684
+ ,37819725
+ ,25253403
+ ,55674674
+ ,40214078
+ ,32066001
+ ,26916456
+ ,31925699
+ ,36677072
+ ,22865662
+ ,33057141
+ ,87224667
+ ,30046064
+ ,33705054
+ ,94643577
+ ,74500253
+ ,13705457
+ ,81920371
+ ,75534971
+ ,17984250
+ ,26564680
+ ,75039636
+ ,31345912
+ ,156181130
+ ,37274013
+ ,60344095
+ ,19480849
+ ,15213215
+ ,36646854
+ ,37800225
+ ,35192619
+ ,51894705
+ ,43116866
+ ,14343902
+ ,21908691
+ ,34532932
+ ,38781483
+ ,32523543
+ ,46755676
+ ,68059999
+ ,34481666
+ ,23934252
+ ,56464201
+ ,30381482
+ ,41957576
+ ,11196762
+ ,26125759
+ ,45626240
+ ,33816982
+ ,150181287
+ ,33424206
+ ,67497380
+ ,54983068
+ ,83797137
+ ,44798290
+ ,39918703
+ ,22068324
+ ,78296142
+ ,41187726
+ ,87854574
+ ,56492262
+ ,36294788
+ ,102238796
+ ,42832378
+ ,22490585
+ ,53384345
+ ,22461067
+ ,24979034
+ ,106897458
+ ,51351879
+ ,39339354
+ ,32067212
+ ,32690247
+ ,53227096
+ ,59342371
+ ,64960155
+ ,54764596
+ ,102374548
+ ,76913992
+ ,42777761
+ ,40847167
+ ,86951363
+ ,43659732
+ ,76698629
+ ,78268805
+ ,91261974
+ ,65606487
+ ,28139063
+ ,75895280
+ ,51122832
+ ,53066011
+ ,1354614
+ ,30021100
+ ,33899609
+ ,58995770
+ ,57876331
+ ,87561067
+ ,26724300
+ ,76278583
+ ,50626962)
> par2 = '12'
> par1 = '500'
> #'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!)
> par1 <- as.numeric(par1)
> par2 <- as.numeric(par2)
> if (par1 < 10) par1 = 10
> if (par1 > 5000) par1 = 5000
> if (par2 < 3) par2 = 3
> if (par2 > length(x)) par2 = length(x)
> library(lattice)
> library(boot)

Attaching package: 'boot'


	The following object(s) are masked from package:lattice :

	 melanoma 

> boot.stat <- function(s)
+ {
+ s.mean <- mean(s)
+ s.median <- median(s)
+ s.midrange <- (max(s) + min(s)) / 2
+ c(s.mean, s.median, s.midrange)
+ }
> (r <- tsboot(x, boot.stat, R=par1, l=12, sim='fixed'))

BLOCK BOOTSTRAP FOR TIME SERIES

Fixed Block Length of 12 

Call:
tsboot(tseries = x, statistic = boot.stat, R = par1, l = 12, 
    sim = "fixed")


Bootstrap Statistics :
     original      bias    std. error
t1*  43903412    29167.66     2310872
t2*  37809975   146567.22     1597169
t3* 102934404 -8226763.20    12146808
> postscript(file="/var/www/html/rcomp/tmp/1ob681290781095.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> plot(r$t[,1],type='p',ylab='simulated values',main='Simulation of Mean')
> grid()
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/2ob681290781095.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> plot(r$t[,2],type='p',ylab='simulated values',main='Simulation of Median')
> grid()
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/3g2nt1290781095.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> plot(r$t[,3],type='p',ylab='simulated values',main='Simulation of Midrange')
> grid()
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/4g2nt1290781095.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> densityplot(~r$t[,1],col='black',main='Density Plot',xlab='mean')
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5g2nt1290781095.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> densityplot(~r$t[,2],col='black',main='Density Plot',xlab='median')
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6ru4e1290781095.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> densityplot(~r$t[,3],col='black',main='Density Plot',xlab='midrange')
> dev.off()
null device 
          1 
> z <- data.frame(cbind(r$t[,1],r$t[,2],r$t[,3]))
Warning message:
In data.row.names(row.names, rowsi, i) :
  some row.names duplicated: 2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,80,81,82,83,84,85,86,87,88,89,90,91,92,93,94,95,96,97,98,99,100,101,102,103,104,105,106,107,108,109,110,111,112,113,114,115,116,117,118,119,120,121,122,123,124,125,126,127,128,129,130,131,132,133,134,135,136,137,138,139,140,141,142,143,144,145,146,147,148,149,150,151,152,153,154,155,156,157,158,159,160,161,162,163,164,165,166,167,168,169,170,171,172,173,174,175,176,177,178,179,180,181,182,183,184,185,186,187,188,189,190,191,192,193,194,195,196,197,198,199,200,201,202,203,204,205,206,207,208,209,210,211,212,213,214,215,216,217,218,219,220,221,222,223,224,225,226,227,228,229,230,231,232,233,234,235,236,237,238,239,240,241,242,243,244,245,246,247,248,249,250,251,252,253,254,255,256,257,258,259,260,261,262,263,264,265,266,267,268,269,270,271 [... truncated]
> colnames(z) <- list('mean','median','midrange')
> postscript(file="/var/www/html/rcomp/tmp/7ru4e1290781095.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> boxplot(z,notch=TRUE,ylab='simulated values',main='Bootstrap Simulation - Central Tendency')
Warning message:
In bxp(list(stats = c(37675147.1948052, 42349247.4220779, 43975979.1996753,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> grid()
> dev.off()
null device 
          1 
> 
> #Note: the /var/www/html/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/www/html/rcomp/createtable")
> 
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,'Estimation Results of Blocked Bootstrap',6,TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'statistic',header=TRUE)
> a<-table.element(a,'Q1',header=TRUE)
> a<-table.element(a,'Estimate',header=TRUE)
> a<-table.element(a,'Q3',header=TRUE)
> a<-table.element(a,'S.D.',header=TRUE)
> a<-table.element(a,'IQR',header=TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'mean',header=TRUE)
> q1 <- quantile(r$t[,1],0.25)[[1]]
> q3 <- quantile(r$t[,1],0.75)[[1]]
> a<-table.element(a,q1)
> a<-table.element(a,r$t0[1])
> a<-table.element(a,q3)
> a<-table.element(a,sqrt(var(r$t[,1])))
> a<-table.element(a,q3-q1)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'median',header=TRUE)
> q1 <- quantile(r$t[,2],0.25)[[1]]
> q3 <- quantile(r$t[,2],0.75)[[1]]
> a<-table.element(a,q1)
> a<-table.element(a,r$t0[2])
> a<-table.element(a,q3)
> a<-table.element(a,sqrt(var(r$t[,2])))
> a<-table.element(a,q3-q1)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'midrange',header=TRUE)
> q1 <- quantile(r$t[,3],0.25)[[1]]
> q3 <- quantile(r$t[,3],0.75)[[1]]
> a<-table.element(a,q1)
> a<-table.element(a,r$t0[3])
> a<-table.element(a,q3)
> a<-table.element(a,sqrt(var(r$t[,3])))
> a<-table.element(a,q3-q1)
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/8gvjp1290781095.tab") 
> 
> try(system("convert tmp/1ob681290781095.ps tmp/1ob681290781095.png",intern=TRUE))
character(0)
> try(system("convert tmp/2ob681290781095.ps tmp/2ob681290781095.png",intern=TRUE))
character(0)
> try(system("convert tmp/3g2nt1290781095.ps tmp/3g2nt1290781095.png",intern=TRUE))
character(0)
> try(system("convert tmp/4g2nt1290781095.ps tmp/4g2nt1290781095.png",intern=TRUE))
character(0)
> try(system("convert tmp/5g2nt1290781095.ps tmp/5g2nt1290781095.png",intern=TRUE))
character(0)
> try(system("convert tmp/6ru4e1290781095.ps tmp/6ru4e1290781095.png",intern=TRUE))
character(0)
> try(system("convert tmp/7ru4e1290781095.ps tmp/7ru4e1290781095.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  2.416   1.253   5.313