FreeStatistics.org

Free Statistics

of Irreproducible Research!

Description of Statistical Computation

Author's title

Author*Unverified author*
R Software Modulerwasp_bootstrapplot1.wasp
Title produced by softwareBootstrap Plot - Central Tendency
Date of computationMon, 08 Sep 2008 06:54:32 -0600
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2008/Sep/08/t12208785411s9pghpez1ehc7d.htm/, Retrieved Sat, 18 May 2024 03:26:04 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=14335, Retrieved Sat, 18 May 2024 03:26:04 +0000
QR Codes:

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

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Bootstrap Plot - Central Tendency] [] [2008-09-04 13:43:55] [74be16979710d4c4e7c6647856088456]
-         [Bootstrap Plot - Central Tendency] [parent child test...] [2008-09-08 12:54:32] [d41d8cd98f00b204e9800998ecf8427e] [Current]
-           [Bootstrap Plot - Central Tendency] [hc7d?] [2008-09-09 08:30:45] [b98453cac15ba1066b407e146608df68]
-             [Bootstrap Plot - Central Tendency] [8jy5] [2008-09-09 08:32:39] [b98453cac15ba1066b407e146608df68]
- R PD          [Bootstrap Plot - Central Tendency] [parent child test...] [2008-09-22 11:23:17] [74be16979710d4c4e7c6647856088456]
-             [Bootstrap Plot - Central Tendency] [hc7d again] [2008-09-09 08:37:05] [b98453cac15ba1066b407e146608df68]
-           [Bootstrap Plot - Central Tendency] [hc7d] [2008-09-09 08:40:26] [b98453cac15ba1066b407e146608df68]
- R PD      [Bootstrap Plot - Central Tendency] [parent child test...] [2008-09-22 11:00:21] [74be16979710d4c4e7c6647856088456]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
        67
        35
        33
        37
        33
        42
        37
        65
        31
        44
        37
        58
        64
        38
        38
        36
        31
        68
        41
        62
        42
        63
        69
        32
        32
        59
        65
        36
        56
        64
        64
        61
        34
        40
        57
        47
        67
        38
        70
        31
        62
        34
        62
        69
        52
        40
        61
        34
        32
        70
        47
        70
        70
        63
        38
        64
        42
        66
        70
        31
        60
        45
        48
        33
        40
        35
        62
        32
        41
        69
        67
        37
        33
        61
        34
        50
        34
        66
        32
        68
        46
        46
        41
        31
        36
        61
        33
        39
        32
        67
        67
        66
        70
        37
        69
        66
        70
        42
        36
        45
        59
        37
        70
        70
        69
        53
        62
        68
        35
        61
        31
        60
        65
        33
        66
        66
        33
        35
        68
        60
        70
        45
        54
        34
        41
        39
        70
        41
        68
        66
        39
        65
        63
        35
        65
        53
        32
        42
        67
        37
        62
        70
        61
        52
        40
        66
        69
        42
        70
        57

Tables (Output of Computation)





Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time3 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 & 3 seconds \tabularnewline
R Server & 'Gwilym Jenkins' @ 72.249.127.135 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=14335&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]3 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=14335&T=0

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






Estimation Results of Bootstrap
statisticQ1EstimateQ3S.D.IQR
mean50.34551.093333333333351.84666666666671.117541604808051.50166666666667
median4752554.756680356722758
midrange50.550.550.50.02236067977499790

\begin{tabular}{lllllllll}
\hline
Estimation Results of Bootstrap \tabularnewline
statistic & Q1 & Estimate & Q3 & S.D. & IQR \tabularnewline
mean & 50.345 & 51.0933333333333 & 51.8466666666667 & 1.11754160480805 & 1.50166666666667 \tabularnewline
median & 47 & 52 & 55 & 4.75668035672275 & 8 \tabularnewline
midrange & 50.5 & 50.5 & 50.5 & 0.0223606797749979 & 0 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=14335&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]50.345[/C][C]51.0933333333333[/C][C]51.8466666666667[/C][C]1.11754160480805[/C][C]1.50166666666667[/C][/ROW]
[ROW][C]median[/C][C]47[/C][C]52[/C][C]55[/C][C]4.75668035672275[/C][C]8[/C][/ROW]
[ROW][C]midrange[/C][C]50.5[/C][C]50.5[/C][C]50.5[/C][C]0.0223606797749979[/C][C]0[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=14335&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=14335&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 Bootstrap
statisticQ1EstimateQ3S.D.IQR
mean50.34551.093333333333351.84666666666671.117541604808051.50166666666667
median4752554.756680356722758
midrange50.550.550.50.02236067977499790


Figures (Output of Computation)

Input Parameters & R Code

Parameters (Session):
par1 = 500 ;
Parameters (R input):
par1 = 500 ; par2 = ; par3 = ; par4 = ; par5 = ; par6 = ; par7 = ; par8 = ; par9 = ; par10 = ; par11 = ; par12 = ; par13 = ; par14 = ; par15 = ; par16 = ; par17 = ; par18 = ; par19 = ; par20 = ;
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')