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 computationSun, 28 Nov 2010 13:09:51 +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/28/t1290949685404a3s8vzvj29z5.htm/, Retrieved Thu, 02 May 2024 19:15:39 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=102538, Retrieved Thu, 02 May 2024 19:15:39 +0000
QR Codes:

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

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]
- RM D    [Blocked Bootstrap Plot - Central Tendency] [Workshop 6 - Notc...] [2010-11-16 10:49:54] [8b017ffbf7b0eded54d8efebfb3e4cfa]
-    D      [Blocked Bootstrap Plot - Central Tendency] [Notched Boxplots ...] [2010-11-28 09:58:58] [8b017ffbf7b0eded54d8efebfb3e4cfa]
-    D          [Blocked Bootstrap Plot - Central Tendency] [Bootstrap Simulat...] [2010-11-28 13:09:51] [3de277db83c2673156e9464be2ef6f69] [Current]
-                 [Blocked Bootstrap Plot - Central Tendency] [Bootstrap Simulat...] [2010-11-28 13:10:40] [8b017ffbf7b0eded54d8efebfb3e4cfa]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
1579
2146
2462
3695
4831
5134
6250
5760
6249
2917
1741
2359
1511
2059
2635
2867
4403
5720
4502
5749
5627
2846
1762
2429
1169
2154
2249
2687
4359
5382
4459
6398
4596
3024
1887
2070
1351
2218
2461
3028
4784
4975
4607
6249
4809
3157
1910
2228
1594
2467
2222
3607
4685
4962
5770
5480
5000
3228
1993
2288
1580
2111
2192
3601
4665
4876
5813
5589
5331
3075
2002
2306
1507
1992
2487
3490
4647
5594
5611
5788
6204
3013
1931
2549
1504
2090
2702
2939
4500
6208
6415
5657
5964
3163
1997
2422

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'George Udny Yule' @ 72.249.76.132

\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 & 'George Udny Yule' @ 72.249.76.132 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=102538&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]'George Udny Yule' @ 72.249.76.132[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=102538&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=102538&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'George Udny Yule' @ 72.249.76.132






Estimation Results of Blocked Bootstrap
statisticQ1EstimateQ3S.D.IQR
mean3577.6406253606.822916666673643.5807291666748.251267419889365.9401041666665
median3018.53051.5311979.9462274881008100.5

\begin{tabular}{lllllllll}
\hline
Estimation Results of Blocked Bootstrap \tabularnewline
statistic & Q1 & Estimate & Q3 & S.D. & IQR \tabularnewline
mean & 3577.640625 & 3606.82291666667 & 3643.58072916667 & 48.2512674198893 & 65.9401041666665 \tabularnewline
median & 3018.5 & 3051.5 & 3119 & 79.9462274881008 & 100.5 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=102538&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]3577.640625[/C][C]3606.82291666667[/C][C]3643.58072916667[/C][C]48.2512674198893[/C][C]65.9401041666665[/C][/ROW]
[ROW][C]median[/C][C]3018.5[/C][C]3051.5[/C][C]3119[/C][C]79.9462274881008[/C][C]100.5[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=102538&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=102538&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
mean3577.6406253606.822916666673643.5807291666748.251267419889365.9401041666665
median3018.53051.5311979.9462274881008100.5






95% Confidence Intervals
MeanMedian
Lower Bound3610.46856111713044.39869531706
Upper Bound3619.802272216233058.60130468294

\begin{tabular}{lllllllll}
\hline
95% Confidence Intervals \tabularnewline
 & Mean & Median \tabularnewline
Lower Bound & 3610.4685611171 & 3044.39869531706 \tabularnewline
Upper Bound & 3619.80227221623 & 3058.60130468294 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=102538&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]3610.4685611171[/C][C]3044.39869531706[/C][/ROW]
[ROW][C]Upper Bound[/C][C]3619.80227221623[/C][C]3058.60130468294[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=102538&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=102538&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 Bound3610.46856111713044.39869531706
Upper Bound3619.802272216233058.60130468294


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