R version 3.0.2 (2013-09-25) -- "Frisbee Sailing"
Copyright (C) 2013 The R Foundation for Statistical Computing
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 <- c(221.102
+ ,220.892
+ ,225.537
+ ,219.334
+ ,216.126
+ ,224.573
+ ,213.991
+ ,221.865
+ ,215.382
+ ,213.962
+ ,217.009
+ ,318.110
+ ,219.662
+ ,212.650
+ ,209.307
+ ,210.541
+ ,210.609
+ ,208.910
+ ,207.541
+ ,207.699
+ ,205.005
+ ,205.747
+ ,203.644
+ ,229.937
+ ,214.446
+ ,210.194
+ ,206.535
+ ,216.524
+ ,198.243
+ ,208.274
+ ,207.493
+ ,215.525
+ ,207.562
+ ,213.355
+ ,209.048
+ ,220.497
+ ,214.563
+ ,211.571
+ ,216.385
+ ,211.496
+ ,209.683
+ ,206.304
+ ,213.925
+ ,204.829
+ ,205.729
+ ,200.296
+ ,207.960
+ ,207.729
+ ,208.327
+ ,207.794
+ ,213.700
+ ,213.136
+ ,210.354
+ ,202.205
+ ,220.425
+ ,210.038
+ ,205.201
+ ,197.441
+ ,200.021
+ ,203.165
+ ,201.699
+ ,201.003
+ ,232.308
+ ,211.412
+ ,209.661
+ ,213.587
+ ,193.611
+ ,192.543
+ ,196.040
+ ,191.407
+ ,189.726
+ ,191.272
+ ,187.924
+ ,198.213
+ ,199.352
+ ,199.289
+ ,195.475
+ ,198.045
+ ,197.615
+ ,189.015
+ ,189.668
+ ,189.120
+ ,194.168
+ ,192.304
+ ,185.913
+ ,197.599
+ ,186.085
+ ,190.566
+ ,187.054
+ ,193.222
+ ,189.856
+ ,190.608
+ ,190.588
+ ,186.773
+ ,183.510
+ ,180.106
+ ,180.150
+ ,178.412
+ ,182.353
+ ,203.805
+ ,186.054
+ ,184.290
+ ,187.483
+ ,187.111
+ ,189.561
+ ,184.439
+ ,182.985
+ ,183.828
+ ,184.036
+ ,183.214
+ ,183.464
+ ,173.718
+ ,180.210
+ ,171.252
+ ,172.705
+ ,174.006
+ ,172.043
+ ,169.445
+ ,169.449
+ ,177.073
+ ,170.799
+ ,171.648
+ ,172.220
+ ,165.795
+ ,167.466
+ ,165.528
+ ,162.851
+ ,165.864
+ ,162.094
+ ,162.385
+ ,164.293
+ ,165.983
+ ,159.680
+ ,161.739
+ ,159.302
+ ,167.795
+ ,164.242
+ ,159.743
+ ,160.887
+ ,163.844
+ ,161.172
+ ,159.330
+ ,155.570
+ ,156.749
+ ,155.012
+ ,163.419
+ ,153.630
+ ,154.535
+ ,151.543
+ ,152.955
+ ,150.166
+ ,151.416
+ ,150.332
+ ,152.196
+ ,153.422
+ ,147.435)
> par2 = '12'
> par1 = '750'
> par2 <- '12'
> par1 <- '750'
> #'GNU S' R Code compiled by R2WASP v. 1.2.291 ()
> #Author: root
> #To cite this work: Wessa P., (2012), Blocked Bootstrap Plot for Central Tendency (v1.0.4) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_bootstrapplot.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #
> 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 is 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* 191.7280   0.01923042    5.554739
t2* 192.8825   0.56048133    8.414909
t3* 232.7725 -13.20773867   20.945466
> postscript(file="/var/fisher/rcomp/tmp/16wi01385120623.ps",horizontal=F,onefile=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/fisher/rcomp/tmp/2td3o1385120623.ps",horizontal=F,onefile=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/fisher/rcomp/tmp/3be481385120623.ps",horizontal=F,onefile=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/fisher/rcomp/tmp/4n5nv1385120623.ps",horizontal=F,onefile=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/fisher/rcomp/tmp/54itu1385120623.ps",horizontal=F,onefile=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/fisher/rcomp/tmp/6gz551385120623.ps",horizontal=F,onefile=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]))
> colnames(z) <- list('mean','median','midrange')
> postscript(file="/var/fisher/rcomp/tmp/7qem21385120623.ps",horizontal=F,onefile=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(176.769044871795, 187.964673076923, 191.787740384615,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> 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,'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/fisher/rcomp/tmp/8fxk81385120623.tab") 
> 
> try(system("convert tmp/16wi01385120623.ps tmp/16wi01385120623.png",intern=TRUE))
character(0)
> try(system("convert tmp/2td3o1385120623.ps tmp/2td3o1385120623.png",intern=TRUE))
character(0)
> try(system("convert tmp/3be481385120623.ps tmp/3be481385120623.png",intern=TRUE))
character(0)
> try(system("convert tmp/4n5nv1385120623.ps tmp/4n5nv1385120623.png",intern=TRUE))
character(0)
> try(system("convert tmp/54itu1385120623.ps tmp/54itu1385120623.png",intern=TRUE))
character(0)
> try(system("convert tmp/6gz551385120623.ps tmp/6gz551385120623.png",intern=TRUE))
character(0)
> try(system("convert tmp/7qem21385120623.ps tmp/7qem21385120623.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  5.314   0.889   6.193