Free Statistics

of Irreproducible Research!

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 computationWed, 17 Nov 2010 00:21:07 +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/17/t1289953187vdmqzprb0czoz0l.htm/, Retrieved Thu, 28 Mar 2024 15:13:24 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=96498, Retrieved Thu, 28 Mar 2024 15:13:24 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact150
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] [Prijs van elektri...] [2010-11-17 00:21:07] [93ab421e12cd1017d2b38fdbcbdb62e0] [Current]
Feedback Forum

Post a new message
Dataseries X:
23,645
23,755
24,136
24,369
24,713
25,657
26,165
26,366
26,504
26,339
25,897
26,386
26,698
27,701
28,15
28,433
28,96
30,931
31,513
31,569
31,755
30,89
30,381
30,967
31,552
32,097
32,473
32,848
33,1
34,962
35,744
36,496
36,273
34,872
34,919
34,96
36,006
36,086
36,741
36,799
36,438
38,047
38,403
38,471
38,265
37,388
36,655
37,074
37,184
37,306
37,393
37,117
37,424
39,216
39,507
39,592
39,578
38,781
38,325
38,188
38,6
38,832
39,066
39,487
39,539
41,635
42,406
43,35
40,801
39,441
38,935
38,626
38,975
39,107
39,312
39,485
39,468
41,169
41,172
41,497
40,25
39,752
39,754




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time2 seconds
R Server'Sir Ronald Aylmer Fisher' @ 193.190.124.24

\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 & 'Sir Ronald Aylmer Fisher' @ 193.190.124.24 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=96498&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]'Sir Ronald Aylmer Fisher' @ 193.190.124.24[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=96498&T=0

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

As an alternative you can also use a 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'Sir Ronald Aylmer Fisher' @ 193.190.124.24







Estimation Results of Blocked Bootstrap
statisticQ1EstimateQ3S.D.IQR
mean33.801587349397634.925530120481935.98843674698801.693673968468482.18684939759036
median34.96236.79937.579752.383598489694972.61774999999999

\begin{tabular}{lllllllll}
\hline
Estimation Results of Blocked Bootstrap \tabularnewline
statistic & Q1 & Estimate & Q3 & S.D. & IQR \tabularnewline
mean & 33.8015873493976 & 34.9255301204819 & 35.9884367469880 & 1.69367396846848 & 2.18684939759036 \tabularnewline
median & 34.962 & 36.799 & 37.57975 & 2.38359848969497 & 2.61774999999999 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=96498&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]33.8015873493976[/C][C]34.9255301204819[/C][C]35.9884367469880[/C][C]1.69367396846848[/C][C]2.18684939759036[/C][/ROW]
[ROW][C]median[/C][C]34.962[/C][C]36.799[/C][C]37.57975[/C][C]2.38359848969497[/C][C]2.61774999999999[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=96498&T=1

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

As an alternative you can also use a QR Code:  

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

Estimation Results of Blocked Bootstrap
statisticQ1EstimateQ3S.D.IQR
mean33.801587349397634.925530120481935.98843674698801.693673968468482.18684939759036
median34.96236.79937.579752.383598489694972.61774999999999







95% Confidence Intervals
MeanMedian
Lower Bound34.809077760846636.6030251886752
Upper Bound35.118705371683536.9949748113248

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

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

As an alternative you can also use a QR Code:  

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

95% Confidence Intervals
MeanMedian
Lower Bound34.809077760846636.6030251886752
Upper Bound35.118705371683536.9949748113248



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