Home » date » 2010 » Nov » 30 »

Blocked Bootstrap Plot 50

*Unverified author*

R Software Module: /rwasp_bootstrapplot.wasp (opens new window with default values)

Title produced by software: Blocked Bootstrap Plot - Central Tendency

Date of computation: Tue, 30 Nov 2010 22:15:56 +0000

Cite this page as follows:

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL http://www.freestatistics.org/blog/date/2010/Nov/30/t1291155347erdc6asw3pey1gt.htm/, Retrieved Wed, 22 May 2013 22:40:07 +0000

Original text written by user:

IsPrivate?

No (this computation is public)

User-defined keywords:

KDGP2W22

System-generated keywords (parent):

(pk = 0)

Estimated Impact

Dataseries X:

» Textfile « » CSV « » Stem and Leaf « » Histogram « » Kernel Density « » Harrell-Davis Quantiles « » Central Tendency « » Variability «

Output produced by software:

Summary of computational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	6 seconds
R Server	'George Udny Yule' @ 72.249.76.132

Estimation Results of Blocked Bootstrap
statistic	Q1	Estimate	Q3	S.D.	IQR
mean	88.8923611111111	92.4305555555556	93.8645833333333	4.64013626011906	4.97222222222221
median	83.5	95.5	97.75	8.82265751994747	14.25
midrange	99	109.5	109.5	9.215913477485	10.5

Charts produced by software:

http://www.freestatistics.org/blog/date/2010/Nov/30/t1291155347erdc6asw3pey1gt/1al5l1291155349.png (opens in new window)

http://www.freestatistics.org/blog/date/2010/Nov/30/t1291155347erdc6asw3pey1gt/1al5l1291155349.ps (opens in new window)

Click here to open pdf file.

http://www.freestatistics.org/blog/date/2010/Nov/30/t1291155347erdc6asw3pey1gt/2al5l1291155349.png (opens in new window)

http://www.freestatistics.org/blog/date/2010/Nov/30/t1291155347erdc6asw3pey1gt/2al5l1291155349.ps (opens in new window)

Click here to open pdf file.

http://www.freestatistics.org/blog/date/2010/Nov/30/t1291155347erdc6asw3pey1gt/33cm51291155349.png (opens in new window)

http://www.freestatistics.org/blog/date/2010/Nov/30/t1291155347erdc6asw3pey1gt/33cm51291155349.ps (opens in new window)

Click here to open pdf file.

http://www.freestatistics.org/blog/date/2010/Nov/30/t1291155347erdc6asw3pey1gt/43cm51291155349.png (opens in new window)

http://www.freestatistics.org/blog/date/2010/Nov/30/t1291155347erdc6asw3pey1gt/43cm51291155349.ps (opens in new window)

Click here to open pdf file.

http://www.freestatistics.org/blog/date/2010/Nov/30/t1291155347erdc6asw3pey1gt/53cm51291155349.png (opens in new window)

http://www.freestatistics.org/blog/date/2010/Nov/30/t1291155347erdc6asw3pey1gt/53cm51291155349.ps (opens in new window)

Click here to open pdf file.

http://www.freestatistics.org/blog/date/2010/Nov/30/t1291155347erdc6asw3pey1gt/63cm51291155349.png (opens in new window)

http://www.freestatistics.org/blog/date/2010/Nov/30/t1291155347erdc6asw3pey1gt/63cm51291155349.ps (opens in new window)

Click here to open pdf file.

http://www.freestatistics.org/blog/date/2010/Nov/30/t1291155347erdc6asw3pey1gt/7wllq1291155349.png (opens in new window)

http://www.freestatistics.org/blog/date/2010/Nov/30/t1291155347erdc6asw3pey1gt/7wllq1291155349.ps (opens in new window)

Click here to open pdf file.

Parameters (Session):

par1 = 50 ; par2 = 12 ;

Parameters (R input):

par1 = 50 ; par2 = 12 ;

R code (references can be found in the software module):

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)
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'))
bitmap(file='plot1.png')
plot(r$t[,1],type='p',ylab='simulated values',main='Simulation of Mean')
grid()
dev.off()
bitmap(file='plot2.png')
plot(r$t[,2],type='p',ylab='simulated values',main='Simulation of Median')
grid()
dev.off()
bitmap(file='plot3.png')
plot(r$t[,3],type='p',ylab='simulated values',main='Simulation of Midrange')
grid()
dev.off()
bitmap(file='plot4.png')
densityplot(~r$t[,1],col='black',main='Density Plot',xlab='mean')
dev.off()
bitmap(file='plot5.png')
densityplot(~r$t[,2],col='black',main='Density Plot',xlab='median')
dev.off()
bitmap(file='plot6.png')
densityplot(~r$t[,3],col='black',main='Density Plot',xlab='midrange')
dev.off()
z <- data.frame(cbind(r$t[,1],r$t[,2],r$t[,3]))
colnames(z) <- list('mean','median','midrange')
bitmap(file='plot7.png')
boxplot(z,notch=TRUE,ylab='simulated values',main='Bootstrap Simulation - Central Tendency')
grid()
dev.off()
load(file='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='mytable.tab')