FreeStatistics.org

Free Statistics

of Irreproducible Research!

Description of Statistical Computation

Author's title

Author*The author of this computation has been verified*
R Software Modulerwasp_bootstrapplot.wasp
Title produced by softwareBlocked Bootstrap Plot - Central Tendency
Date of computationFri, 26 Nov 2010 14:43:41 +0000
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2010/Nov/26/t129078252292ugjc0z5pezx7u.htm/, Retrieved Fri, 03 May 2024 21:47:28 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=101949, Retrieved Fri, 03 May 2024 21:47:28 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact95

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Bivariate Data Series] [Bivariate dataset] [2008-01-05 23:51:08] [74be16979710d4c4e7c6647856088456]
- RMPD  [Blocked Bootstrap Plot - Central Tendency] [Colombia Coffee] [2008-01-07 10:26:26] [74be16979710d4c4e7c6647856088456]
-  M D      [Blocked Bootstrap Plot - Central Tendency] [mini tutorial wor...] [2010-11-26 14:43:41] [0dfe009a651fec1e160584d659799586] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
4946
4490
4851
4591
4279
4191
4285
4516
4197
4404
4373
5307
5320
4356
4484
4210
4018
3912
3972
3886
3892
4242
4134
4743
5116
4823
5489
4263
4221
4076
3715
3715
3784
4141
3968
4767
5019
4343
4853
4154
4035
3996
4734
3778
3887
3953
3987
4436
4803
4672
4560
4289
3961
3943
3932
3816
3834
4130
4467
4447

Tables (Output of Computation)





Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time2 seconds
R Server'RServer@AstonUniversity' @ vre.aston.ac.uk

\begin{tabular}{lllllllll}
\hline
Summary of computational transaction \tabularnewline
Raw Input & view raw input (R code)  \tabularnewline
Raw Output & view raw output of R engine  \tabularnewline
Computing time & 2 seconds \tabularnewline
R Server & 'RServer@AstonUniversity' @ vre.aston.ac.uk \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=101949&T=0

[TABLE]
[ROW][C]Summary of computational transaction[/C][/ROW]
[ROW][C]Raw Input[/C][C]view raw input (R code) [/C][/ROW]
[ROW][C]Raw Output[/C][C]view raw output of R engine [/C][/ROW]
[ROW][C]Computing time[/C][C]2 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'RServer@AstonUniversity' @ vre.aston.ac.uk[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=101949&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=101949&T=0

As an alternative you can also use a QR Code:  
QR Code


The GUIDs for individual cells are displayed in the table below:

Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time2 seconds
R Server'RServer@AstonUniversity' @ vre.aston.ac.uk






Estimation Results of Blocked Bootstrap
statisticQ1EstimateQ3S.D.IQR
mean4296.754328.433333333334359.6583333333345.345275850222362.9083333333338
median42094252.5428566.167581463080176

\begin{tabular}{lllllllll}
\hline
Estimation Results of Blocked Bootstrap \tabularnewline
statistic & Q1 & Estimate & Q3 & S.D. & IQR \tabularnewline
mean & 4296.75 & 4328.43333333333 & 4359.65833333333 & 45.3452758502223 & 62.9083333333338 \tabularnewline
median & 4209 & 4252.5 & 4285 & 66.1675814630801 & 76 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=101949&T=1

[TABLE]
[ROW][C]Estimation Results of Blocked Bootstrap[/C][/ROW]
[ROW][C]statistic[/C][C]Q1[/C][C]Estimate[/C][C]Q3[/C][C]S.D.[/C][C]IQR[/C][/ROW]
[ROW][C]mean[/C][C]4296.75[/C][C]4328.43333333333[/C][C]4359.65833333333[/C][C]45.3452758502223[/C][C]62.9083333333338[/C][/ROW]
[ROW][C]median[/C][C]4209[/C][C]4252.5[/C][C]4285[/C][C]66.1675814630801[/C][C]76[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=101949&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=101949&T=1

As an alternative you can also use a QR Code:  
QR Code


The GUIDs for individual cells are displayed in the table below:

Estimation Results of Blocked Bootstrap
statisticQ1EstimateQ3S.D.IQR
mean4296.754328.433333333334359.6583333333345.345275850222362.9083333333338
median42094252.5428566.167581463080176






95% Confidence Intervals
MeanMedian
Lower Bound4324.441459591644239.05919939715
Upper Bound4333.375207075034249.94080060285

\begin{tabular}{lllllllll}
\hline
95% Confidence Intervals \tabularnewline
 & Mean & Median \tabularnewline
Lower Bound & 4324.44145959164 & 4239.05919939715 \tabularnewline
Upper Bound & 4333.37520707503 & 4249.94080060285 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=101949&T=2

[TABLE]
[ROW][C]95% Confidence Intervals[/C][/ROW]
[ROW][C][/C][C]Mean[/C][C]Median[/C][/ROW]
[ROW][C]Lower Bound[/C][C]4324.44145959164[/C][C]4239.05919939715[/C][/ROW]
[ROW][C]Upper Bound[/C][C]4333.37520707503[/C][C]4249.94080060285[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=101949&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=101949&T=2

As an alternative you can also use a QR Code:  
QR Code


The GUIDs for individual cells are displayed in the table below:

95% Confidence Intervals
MeanMedian
Lower Bound4324.441459591644239.05919939715
Upper Bound4333.375207075034249.94080060285


Figures (Output of Computation)

Input Parameters & R Code

Parameters (Session):
par1 = 500 ; par2 = 12 ;
Parameters (R input):
par1 = 500 ; 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)
c(s.mean, s.median)
}
(r <- tsboot(x, boot.stat, R=par1, l=12, sim='fixed'))
z <- data.frame(cbind(r$t[,1],r$t[,2]))
colnames(z) <- list('mean','median')
bitmap(file='plot7.png')
b <- boxplot(z,notch=TRUE,ylab='simulated values',main='Bootstrap Simulation - Central Tendency')
grid()
dev.off()
b
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.end(a)
table.save(a,file='mytable.tab')

a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'95% Confidence Intervals',3,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'',1,TRUE)
a<-table.element(a,'Mean',1,TRUE)
a<-table.element(a,'Median',1,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Lower Bound',1,TRUE)
a<-table.element(a,b$conf[1,1])
a<-table.element(a,b$conf[1,2])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Upper Bound',1,TRUE)
a<-table.element(a,b$conf[2,1])
a<-table.element(a,b$conf[2,2])
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable1.tab')