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 <- c(414.89
+ ,444.50
+ ,481.29
+ ,491.09
+ ,419.70
+ ,432.88
+ ,438.01
+ ,412.84
+ ,402.91
+ ,416.20
+ ,411.80
+ ,393.60
+ ,381.70
+ ,388.34
+ ,370.89
+ ,385.96
+ ,394.26
+ ,381.37
+ ,377.40
+ ,377.70
+ ,347.71
+ ,347.70
+ ,340.20
+ ,341.00
+ ,342.00
+ ,319.54
+ ,302.79
+ ,299.10
+ ,313.50
+ ,326.80
+ ,316.00
+ ,316.50
+ ,317.20
+ ,330.40
+ ,323.35
+ ,325.85
+ ,321.50
+ ,321.90
+ ,347.48
+ ,338.89
+ ,345.70
+ ,340.44
+ ,342.40
+ ,342.70
+ ,348.34
+ ,376.66
+ ,417.73
+ ,423.51
+ ,397.56
+ ,390.92
+ ,408.26
+ ,401.12
+ ,408.91
+ ,438.35
+ ,460.23
+ ,449.59
+ ,450.52
+ ,461.15
+ ,460.38
+ ,465.35
+ ,467.57
+ ,486.24
+ ,476.58
+ ,442.07
+ ,443.61
+ ,451.55
+ ,451.07
+ ,451.33
+ ,437.63
+ ,431.28
+ ,413.46
+ ,406.78
+ ,420.17
+ ,419.05
+ ,404.01
+ ,387.51
+ ,390.15
+ ,384.40
+ ,371.05
+ ,367.60
+ ,375.04
+ ,365.14
+ ,361.75
+ ,366.88
+ ,394.26
+ ,409.39
+ ,410.11
+ ,416.81
+ ,393.06
+ ,374.24
+ ,369.05
+ ,352.33
+ ,362.53
+ ,394.73
+ ,389.32
+ ,380.74
+ ,381.73
+ ,376.95
+ ,383.64
+ ,363.83
+ ,363.34
+ ,358.38
+ ,356.95
+ ,366.72
+ ,367.69
+ ,356.31
+ ,348.74
+ ,358.69
+ ,360.17
+ ,361.73
+ ,354.45
+ ,353.91
+ ,344.34
+ ,338.62
+ ,337.24
+ ,340.81
+ ,352.72
+ ,343.06
+ ,345.43
+ ,344.38
+ ,335.02
+ ,334.82
+ ,329.01
+ ,329.31
+ ,330.08
+ ,342.15
+ ,367.18
+ ,371.89
+ ,392.19
+ ,378.84
+ ,355.28
+ ,364.18
+ ,373.83
+ ,383.30
+ ,386.88
+ ,381.91
+ ,384.13
+ ,377.27
+ ,381.43
+ ,385.64
+ ,385.49
+ ,380.36
+ ,391.58
+ ,389.77
+ ,384.39
+ ,379.29
+ ,378.55
+ ,376.64
+ ,382.12
+ ,391.03
+ ,385.22
+ ,387.56
+ ,386.23
+ ,383.67
+ ,383.06
+ ,383.14
+ ,385.31
+ ,387.44
+ ,399.45
+ ,404.76
+ ,396.21
+ ,392.85
+ ,391.93
+ ,385.27
+ ,383.47
+ ,387.35
+ ,383.14
+ ,381.07
+ ,377.85
+ ,369.00
+ ,355.11
+ ,346.58
+ ,351.81
+ ,344.47
+ ,343.84
+ ,340.76
+ ,324.10
+ ,324.01
+ ,322.82
+ ,324.87
+ ,306.04
+ ,288.74
+ ,289.10
+ ,297.49
+ ,295.94
+ ,308.29
+ ,299.10
+ ,292.32
+ ,292.87
+ ,284.11
+ ,288.98
+ ,295.93
+ ,294.12
+ ,291.68
+ ,287.08
+ ,287.33
+ ,285.96
+ ,282.62
+ ,276.44
+ ,261.31
+ ,256.08
+ ,256.69
+ ,264.74
+ ,310.72
+ ,293.18
+ ,283.07
+ ,284.32
+ ,299.86
+ ,286.39
+ ,279.69
+ ,275.19
+ ,285.73
+ ,281.59
+ ,274.47
+ ,273.68
+ ,270.00
+ ,266.01
+ ,271.45
+ ,265.49
+ ,261.87
+ ,263.03
+ ,260.48
+ ,272.36
+ ,270.23
+ ,267.53
+ ,272.39
+ ,283.42
+ ,283.06
+ ,276.16
+ ,275.85
+ ,281.51
+ ,295.50
+ ,294.06
+ ,302.68
+ ,314.49
+ ,321.18
+ ,313.29
+ ,310.26
+ ,319.14
+ ,316.56
+ ,319.07
+ ,331.92
+ ,356.86
+ ,358.97
+ ,340.55
+ ,328.18
+ ,355.68
+ ,356.35
+ ,351.02
+ ,359.77
+ ,378.95
+ ,378.92
+ ,389.91
+ ,406.95
+ ,413.79
+ ,404.88
+ ,406.67
+ ,403.26
+ ,383.78
+ ,392.37
+ ,398.09
+ ,400.51
+ ,405.28
+ ,420.46
+ ,439.38
+ ,442.08
+ ,424.03
+ ,423.35
+ ,433.85
+ ,429.23
+ ,421.87
+ ,430.66
+ ,424.48
+ ,437.93
+ ,456.05
+ ,469.90
+ ,476.67
+ ,510.10
+ ,549.86
+ ,555.00
+ ,557.09
+ ,610.65
+ ,675.39
+ ,596.15
+ ,633.71
+ ,632.59
+ ,598.19
+ ,585.78
+ ,627.83
+ ,629.79
+ ,631.17
+ ,664.75
+ ,654.90
+ ,679.37
+ ,667.31
+ ,655.66
+ ,665.38
+ ,665.41
+ ,712.65
+ ,754.60
+ ,806.25
+ ,803.20
+ ,889.60
+ ,922.30
+ ,968.43
+ ,909.71
+ ,888.66
+ ,889.49
+ ,939.77
+ ,839.03
+ ,829.93
+ ,806.62
+ ,760.86
+ ,816.09
+ ,858.69
+ ,943.00
+ ,924.27
+ ,890.20
+ ,928.65
+ ,945.67
+ ,934.23
+ ,949.38
+ ,996.59
+ ,1043.16
+ ,1127.04
+ ,1134.72
+ ,1117.96
+ ,1095.41
+ ,1113.34
+ ,1148.69
+ ,1205.43
+ ,1232.92
+ ,1192.97
+ ,1215.81
+ ,1270.98
+ ,1342.02
+ ,1369.89
+ ,1390.55
+ ,1356.40
+ ,1372.73
+ ,1424.00
+ ,1479.76
+ ,1512.60
+ ,1528.66
+ ,1572.21
+ ,1757.21
+ ,1770.95
+ ,1665.21
+ ,1738.11
+ ,1641.84
+ ,1652.21
+ ,1742.14
+ ,1673.77
+ ,1649.69
+ ,1591.19
+ ,1598.76
+ ,1589.90
+ ,1630.31
+ ,1744.81
+ ,1746.58
+ ,1721.64)
> par2 = '12'
> par1 = '500'
> par2 <- '12'
> par1 <- '500'
> #'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(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*  527.1657   5.237252    58.69898
t2*  384.4000   2.696540    13.07892
t3* 1013.5150 -14.677280    54.79128
> postscript(file="/var/wessaorg/rcomp/tmp/1rrhc1355657277.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/wessaorg/rcomp/tmp/2nho51355657277.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/wessaorg/rcomp/tmp/3ib5n1355657277.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/wessaorg/rcomp/tmp/4v4ey1355657277.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/wessaorg/rcomp/tmp/5jbfn1355657277.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/wessaorg/rcomp/tmp/6dmt71355657277.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/wessaorg/rcomp/tmp/735nf1355657277.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(392.798088642659, 484.22012465374, 532.065027700831,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> grid()
> dev.off()
null device 
          1 
> 
> #Note: the /var/wessaorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/wessaorg/rcomp/createtable")
> 
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,'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/wessaorg/rcomp/tmp/8ly491355657277.tab") 
> 
> try(system("convert tmp/1rrhc1355657277.ps tmp/1rrhc1355657277.png",intern=TRUE))
character(0)
> try(system("convert tmp/2nho51355657277.ps tmp/2nho51355657277.png",intern=TRUE))
character(0)
> try(system("convert tmp/3ib5n1355657277.ps tmp/3ib5n1355657277.png",intern=TRUE))
character(0)
> try(system("convert tmp/4v4ey1355657277.ps tmp/4v4ey1355657277.png",intern=TRUE))
character(0)
> try(system("convert tmp/5jbfn1355657277.ps tmp/5jbfn1355657277.png",intern=TRUE))
character(0)
> try(system("convert tmp/6dmt71355657277.ps tmp/6dmt71355657277.png",intern=TRUE))
character(0)
> try(system("convert tmp/735nf1355657277.ps tmp/735nf1355657277.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  5.460   0.872   6.322