Free Statistics

of Irreproducible Research!

Author's title

Author*Unverified author*
R Software Modulerwasp_bootstrapplot1.wasp
Title produced by softwareBootstrap Plot - Central Tendency
Date of computationMon, 02 May 2011 19:09:52 +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/2011/May/02/t1304363234ur6kqwiggfpv75z.htm/, Retrieved Sat, 11 May 2024 05:50:19 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=120823, Retrieved Sat, 11 May 2024 05:50:19 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywordsKDGP2W22
Estimated Impact22
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Bootstrap Plot - Central Tendency] [bootstrap plot ve...] [2011-05-02 19:09:52] [cd3c8dd726b701571412eb3280c696e7] [Current]
Feedback Forum

Post a new message
Dataseries X:
814
1150
1225
1691
1759
1754
2100
2062
2012
1897
1964
2186
966
1549
1538
1612
2078
2137
2907
2249
1883
1739
1828
1868
1138
1430
1809
1763
2200
2067
2503
2141
2103
1972
2181
2344
970
1199
1718
1683
2025
2051
2439
2353
2230
1852
2147
2286
1007
1665
1642
1518
1831
2207
2822
2393
2306
1785
2047
2171
1212
1335
2011
1860
1954
2152
2835
2224
2182
1992
2389
2724
891
1247
2017
2257
2255
2255
3057
3330
1896
2096
2374
2535
1041
1728
2201
2455
2204
2660
3670
2665
2639
2226
2586
2684
1185
1749
2459
2618
2585
3310
3923




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time4 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 & 4 seconds \tabularnewline
R Server & 'RServer@AstonUniversity' @ vre.aston.ac.uk \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=120823&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]4 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=120823&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=120823&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 time4 seconds
R Server'RServer@AstonUniversity' @ vre.aston.ac.uk







Estimation Results of Bootstrap
statisticQ1EstimateQ3S.D.IQR
mean2018.274271844662044.990291262142078.614077669950.73870885109460.3398058252428
median20472067209951.116689235193552
midrange2280.52368.52407103.670115052999126.5

\begin{tabular}{lllllllll}
\hline
Estimation Results of Bootstrap \tabularnewline
statistic & Q1 & Estimate & Q3 & S.D. & IQR \tabularnewline
mean & 2018.27427184466 & 2044.99029126214 & 2078.6140776699 & 50.738708851094 & 60.3398058252428 \tabularnewline
median & 2047 & 2067 & 2099 & 51.1166892351935 & 52 \tabularnewline
midrange & 2280.5 & 2368.5 & 2407 & 103.670115052999 & 126.5 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=120823&T=1

[TABLE]
[ROW][C]Estimation Results of 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]2018.27427184466[/C][C]2044.99029126214[/C][C]2078.6140776699[/C][C]50.738708851094[/C][C]60.3398058252428[/C][/ROW]
[ROW][C]median[/C][C]2047[/C][C]2067[/C][C]2099[/C][C]51.1166892351935[/C][C]52[/C][/ROW]
[ROW][C]midrange[/C][C]2280.5[/C][C]2368.5[/C][C]2407[/C][C]103.670115052999[/C][C]126.5[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=120823&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=120823&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 Bootstrap
statisticQ1EstimateQ3S.D.IQR
mean2018.274271844662044.990291262142078.614077669950.73870885109460.3398058252428
median20472067209951.116689235193552
midrange2280.52368.52407103.670115052999126.5



Parameters (Session):
par1 = 50 ;
Parameters (R input):
par1 = 50 ;
R code (references can be found in the software module):
par1 <- as.numeric(par1)
if (par1 < 10) par1 = 10
if (par1 > 5000) par1 = 5000
library(lattice)
library(boot)
boot.stat <- function(s,i)
{
s.mean <- mean(s[i])
s.median <- median(s[i])
s.midrange <- (max(s[i]) + min(s[i])) / 2
c(s.mean, s.median, s.midrange)
}
(r <- boot(x,boot.stat, R=par1, stype='i'))
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 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')