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 computationTue, 07 Dec 2010 20:34:46 +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/Dec/07/t1291753993j50i126celtgd8r.htm/, Retrieved Fri, 03 May 2024 23:04:09 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=106719, Retrieved Fri, 03 May 2024 23:04:09 +0000
QR Codes:

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

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] [notched boxplots ...] [2010-12-07 20:34:46] [dc77c696707133dea0955379c56a2acd] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
57619
57077
56161
55614
57029
57667
57275
57514
58519
58759
58170
58277
61050
63191
64035
63358
62149
61789
60973
61229
60566
61922
62945
64810
66641
67675
66585
66573
66193
65545
64835
67184
69421
72359
74730
75067
74927
74845
75335
78121
79944
82779
81484
79939
80611
79088
79992
77800
76277
75425
74524
72844
70768
70172
67726
68032
67096
68261
67363
67322
57465
57031
54046
47435
43623
48866
53075
53778
51237
50393
49947
50660
50307
52235
51783
51602
47498
45981
43588
42791
41325
38917
36971
33628
31445
28478
19247
12209
10467
10127
10344
10117
8142
6404
4519
3482
2373
1772
1178
746
503
293
161
91
55
20
14
8
2
0
1

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'Gwilym Jenkins' @ 72.249.127.135

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=106719&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'Gwilym Jenkins' @ 72.249.127.135






Estimation Results of Blocked Bootstrap
statisticQ1EstimateQ3S.D.IQR
mean42477.342342342348356.405405405454510.76576576588158.1803186346412033.4234234234
median5178357465612299609.300946607149446

\begin{tabular}{lllllllll}
\hline
Estimation Results of Blocked Bootstrap \tabularnewline
statistic & Q1 & Estimate & Q3 & S.D. & IQR \tabularnewline
mean & 42477.3423423423 & 48356.4054054054 & 54510.7657657658 & 8158.18031863464 & 12033.4234234234 \tabularnewline
median & 51783 & 57465 & 61229 & 9609.30094660714 & 9446 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=106719&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]42477.3423423423[/C][C]48356.4054054054[/C][C]54510.7657657658[/C][C]8158.18031863464[/C][C]12033.4234234234[/C][/ROW]
[ROW][C]median[/C][C]51783[/C][C]57465[/C][C]61229[/C][C]9609.30094660714[/C][C]9446[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=106719&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=106719&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
mean42477.342342342348356.405405405454510.76576576588158.1803186346412033.4234234234
median5178357465612299609.300946607149446






95% Confidence Intervals
MeanMedian
Lower Bound47772.138333433856797.5480195514
Upper Bound49473.762567467158132.4519804486

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=106719&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 Bound47772.138333433856797.5480195514
Upper Bound49473.762567467158132.4519804486


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