Home » date » 2007 » Oct » 28 »

Hypothesis Testing with EDA (part 1) Q4

*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: Sun, 28 Oct 2007 08:25:15 -0700

Cite this page as follows:

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL http://www.freestatistics.org/blog/date/2007/Oct/28/9iiy06ztuf18qbx1193584708.htm/, Retrieved Sun, 26 May 2013 00:27:19 +0000

Original text written by user:

IsPrivate?

No (this computation is public)

User-defined keywords:

Hypothesis Testing with EDA (part 1) Q4

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 compuational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	3 seconds
R Server	'Gwilym Jenkins' @ 72.249.127.135

Estimation Results of Blocked Bootstrap
statistic	Q1	Estimate	Q3	S.D.	IQR
mean	85.65	86.8934426229508	87.9040983606557	1.64820091260918	2.25409836065573
median	86.4	87.3	88	1.87490628757129	1.59999999999999
midrange	87.85	88.1	88.85	1.05908182338157	1

Charts produced by software:

http://www.freestatistics.org/blog/date/2007/Oct/28/9iiy06ztuf18qbx1193584708/1pmgv1193585111.png (opens in new window)

http://www.freestatistics.org/blog/date/2007/Oct/28/9iiy06ztuf18qbx1193584708/1pmgv1193585111.ps (opens in new window)

Click here to open pdf file.

http://www.freestatistics.org/blog/date/2007/Oct/28/9iiy06ztuf18qbx1193584708/2n1vu1193585111.png (opens in new window)

http://www.freestatistics.org/blog/date/2007/Oct/28/9iiy06ztuf18qbx1193584708/2n1vu1193585111.ps (opens in new window)

Click here to open pdf file.

http://www.freestatistics.org/blog/date/2007/Oct/28/9iiy06ztuf18qbx1193584708/35t8i1193585111.png (opens in new window)

http://www.freestatistics.org/blog/date/2007/Oct/28/9iiy06ztuf18qbx1193584708/35t8i1193585111.ps (opens in new window)

Click here to open pdf file.

http://www.freestatistics.org/blog/date/2007/Oct/28/9iiy06ztuf18qbx1193584708/46cwi1193585111.png (opens in new window)

http://www.freestatistics.org/blog/date/2007/Oct/28/9iiy06ztuf18qbx1193584708/46cwi1193585111.ps (opens in new window)

Click here to open pdf file.

http://www.freestatistics.org/blog/date/2007/Oct/28/9iiy06ztuf18qbx1193584708/5gol81193585111.png (opens in new window)

http://www.freestatistics.org/blog/date/2007/Oct/28/9iiy06ztuf18qbx1193584708/5gol81193585111.ps (opens in new window)

Click here to open pdf file.

http://www.freestatistics.org/blog/date/2007/Oct/28/9iiy06ztuf18qbx1193584708/6hdsj1193585111.png (opens in new window)

http://www.freestatistics.org/blog/date/2007/Oct/28/9iiy06ztuf18qbx1193584708/6hdsj1193585111.ps (opens in new window)

Click here to open pdf file.

http://www.freestatistics.org/blog/date/2007/Oct/28/9iiy06ztuf18qbx1193584708/7psi01193585111.png (opens in new window)

http://www.freestatistics.org/blog/date/2007/Oct/28/9iiy06ztuf18qbx1193584708/7psi01193585111.ps (opens in new window)

Click here to open pdf file.

Parameters (Session):

par1 = 61 ;

Parameters (R input):

par1 = 61 ;

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')