Repository of Reproducible Computations

Free Statistics

of Irreproducible Research!

Author's title

Author

*The author of this computation has been verified*

R Software Module

rwasp_BoxplotAndScatterplotMatrix.wasp

Title produced by software

Boxplot and Trimmed Means

Date of computation

Mon, 01 Dec 2014 15:27:08 +0000

Cite this page as follows

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2014/Dec/01/t1417447652ktw3x4qck11z2qx.htm/, Retrieved Thu, 16 May 2024 16:36:55 +0000

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=261921, Retrieved Thu, 16 May 2024 16:36:55 +0000

QR Codes:

Paste this QR Code to cite your computation.

Original text written by user:

IsPrivate?

No (this computation is public)

User-defined keywords

Estimated Impact

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)

-       [Boxplot and Trimmed Means] [cognition in preg...] [2014-12-01 15:27:08] [59ea4bc8516da637c340908901d58948] [Current]

Feedback Forum

Post a new message

Dataseries X:

Download CSV

Histogram

Boxplots

4.8	5	4.4	4.6	4.8	4.6
5.4	5.4	4.2	5.6	4.2	4
6	5.6	5.2	4	4.6	4.8
6	5.6	5.2	5	4.8	4
4	5.4	4.4	4.4	4.6	5
6	6	5.5	6	6.2	5
4.2	4	4	4	4.4	4
4.2	4.2	4	4.4	4.2	4
4.2	4.6	4	4.6	4	4.4
4	4.2	3.8	4	4.2	4
4.6	4.2	4	4.8	4.4	4
5.666666667	4.666666667	5.25	5.333333333	5.2	5
3.8	3.6	3.6	4.4	3.8	3.8
4.8	5	4.4	5.2	5.2	4.6
4.2	4.6	4	5.2	4.6	4
4	4	4.4	4	4.2	4.6
5.4	5.8	4.333333333	5.2	5.4	4.666666667
4.6	5.8	4.6	4.8	4.6	4
4	4	4	4	4.4	4
4.4	5.4	4.4	4.4	5	4.4
5	4.8	5	5.6	4.8	4.8
5.4	5.8	5.2	6	6.2	5.8
5.2	5.4	5.2	4.8	5.2	4.6
4.4	4.8	4.6	4.2	4.6	4.75
4	4.6	4.6	4	4.6	5
3.4	4.6	2.333333333	2.8	4.2	2.75
4.4	4	4	4.8	4	4
4	4	3.8	4.2	4.6	3.8
4.4	4.8	4	4.8	4.6	4.4
4.4	5	4.6	5.2	5	4.8
4	4	4.75	4.6	4.6	4.25
4.8	5.2	4.4	4.4	5.6	4.4
4	4	4	4	4	4
4	4.4	4	3.6	3.8	4.5
5	5.6	5.2	4.6	5.8	5
4	4.6	3.8	4	4.2	3.5
3.4	2.2	2.5	2.2	3.4	2
4	4	3.2	4	3.8	4.4
5.4	4.6	4.2	5.6	4.8	4.4
4.2	4.2	4	3.6	4.2	4.25
4	4	3.25	4	4	3.25
5.6	5	4.2	5.4	5	4.2
4.8	5	4.75	5	5	4.75
4.4	4.8	4.4	4.4	4.6	4.4
4.6	5	4	4.6	5	4
4	3.8	4.2	3.8	3.8	4.2
5.2	4.6	4	4.25	5	5
4.8	4.2	4.2	4.6	4.8	4.6
4.6	5	4.4	5.2	5	4.4
4	4	4	4	4.2	4
4.2	5.2	4.6	4.2	4.8	4.4
4.6	4.8	4.4	4.2	4.6	4.5
4.4	4.2	4	4.6	5	5.2
5	5	4.8	5.6	5	4.4
5	4.2	2.25	5.4	4	3.25
4	4	4.8	4	4	4.6
3.8	3.2	3.75	3.6	3.2	3.75
5.2	6.2	5.6	4.4	5.2	5
4.2	5.6	4.2	4.6	4.6	3.8
4.4	4.6	4.2	4.4	3.6	4
4.2	5.2	5	4	4.8	5.4
6	6	5.8	5	4.8	5
4.8	4.8	5	4.8	4.8	5
4	4	4	4	4	4.2
4.4	4.4	4.4	4.6	5	4.4
4	4.4	4.4	4	4.4	4
4.8	4.8	4.4	5	5.2	4.6
4.8	4.6	4.4	4.2	4	4.2
5.6	6.4	5.8	4.8	5.6	5.2
6	6	6	6.8	7	6.4
4	4.2	4	4	4.2	4
4.2	4	4	4.2	4	4
5	5.2	5	5	5.2	5
4.4	4.4	4	5	4.2	4
4.2	4.6	4.2	4.6	4.6	4.2

Summary of computational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	2 seconds
R Server	'Sir Ronald Aylmer Fisher' @ fisher.wessa.net

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=261921&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 Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	2 seconds
R Server	'Sir Ronald Aylmer Fisher' @ fisher.wessa.net

Boxplot statistics
Variable	lower whisker	lower hinge	median	upper hinge	upper whisker
comp1	3.4	4	4.4	5	6
comp2	3.2	4.2	4.6	5.2	6.4
comp3	3.2	4	4.33	4.67	5.6
spec1	2.8	4	4.6	5	6
spec2	3.2	4.2	4.6	5	5.8
spec3	3.25	4	4.4	4.75	5.8

\begin{tabular}{lllllllll}
\hline
Boxplot statistics \tabularnewline
Variable & lower whisker & lower hinge & median & upper hinge & upper whisker \tabularnewline
comp1 & 3.4 & 4 & 4.4 & 5 & 6 \tabularnewline
comp2 & 3.2 & 4.2 & 4.6 & 5.2 & 6.4 \tabularnewline
comp3 & 3.2 & 4 & 4.33 & 4.67 & 5.6 \tabularnewline
spec1 & 2.8 & 4 & 4.6 & 5 & 6 \tabularnewline
spec2 & 3.2 & 4.2 & 4.6 & 5 & 5.8 \tabularnewline
spec3 & 3.25 & 4 & 4.4 & 4.75 & 5.8 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=261921&T=1

[TABLE]
[ROW][C]Boxplot statistics[/C][/ROW]
[ROW][C]Variable[/C][C]lower whisker[/C][C]lower hinge[/C][C]median[/C][C]upper hinge[/C][C]upper whisker[/C][/ROW]
[ROW][C]comp1[/C][C]3.4[/C][C]4[/C][C]4.4[/C][C]5[/C][C]6[/C][/ROW]
[ROW][C]comp2[/C][C]3.2[/C][C]4.2[/C][C]4.6[/C][C]5.2[/C][C]6.4[/C][/ROW]
[ROW][C]comp3[/C][C]3.2[/C][C]4[/C][C]4.33[/C][C]4.67[/C][C]5.6[/C][/ROW]
[ROW][C]spec1[/C][C]2.8[/C][C]4[/C][C]4.6[/C][C]5[/C][C]6[/C][/ROW]
[ROW][C]spec2[/C][C]3.2[/C][C]4.2[/C][C]4.6[/C][C]5[/C][C]5.8[/C][/ROW]
[ROW][C]spec3[/C][C]3.25[/C][C]4[/C][C]4.4[/C][C]4.75[/C][C]5.8[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=261921&T=1

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

As an alternative you can also use a QR Code:

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

Boxplot statistics
Variable	lower whisker	lower hinge	median	upper hinge	upper whisker
comp1	3.4	4	4.4	5	6
comp2	3.2	4.2	4.6	5.2	6.4
comp3	3.2	4	4.33	4.67	5.6
spec1	2.8	4	4.6	5	6
spec2	3.2	4.2	4.6	5	5.8
spec3	3.25	4	4.4	4.75	5.8

Boxplot Means
Variable	trimmed mean	unbiased SD
comp1	4.57	0.64
comp2	4.71	0.74
comp3	4.34	0.7
spec1	4.55	0.71
spec2	4.63	0.65
spec3	4.37	0.63

\begin{tabular}{lllllllll}
\hline
Boxplot Means \tabularnewline
Variable & trimmed mean & unbiased SD \tabularnewline
comp1 & 4.57 & 0.64 \tabularnewline
comp2 & 4.71 & 0.74 \tabularnewline
comp3 & 4.34 & 0.7 \tabularnewline
spec1 & 4.55 & 0.71 \tabularnewline
spec2 & 4.63 & 0.65 \tabularnewline
spec3 & 4.37 & 0.63 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=261921&T=2

[TABLE]
[ROW][C]Boxplot Means[/C][/ROW]
[ROW][C]Variable[/C][C]trimmed mean[/C][C]unbiased SD[/C][/ROW]
[ROW][C]comp1[/C][C]4.57[/C][C]0.64[/C][/ROW]
[ROW][C]comp2[/C][C]4.71[/C][C]0.74[/C][/ROW]
[ROW][C]comp3[/C][C]4.34[/C][C]0.7[/C][/ROW]
[ROW][C]spec1[/C][C]4.55[/C][C]0.71[/C][/ROW]
[ROW][C]spec2[/C][C]4.63[/C][C]0.65[/C][/ROW]
[ROW][C]spec3[/C][C]4.37[/C][C]0.63[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=261921&T=2

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

As an alternative you can also use a QR Code:

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

Boxplot Means
Variable	trimmed mean	unbiased SD
comp1	4.57	0.64
comp2	4.71	0.74
comp3	4.34	0.7
spec1	4.55	0.71
spec2	4.63	0.65
spec3	4.37	0.63

Figure 1

PNG link

Postscript link

PDF link

Figure 2

PNG link

Postscript link

PDF link

Parameters (Session):

par1 = 3 ; par2 = FALSE ; par3 = 0 ;

Parameters (R input):

par1 = 3 ; par2 = FALSE ; par3 = 0 ;

R code (references can be found in the software module):

par3 <- '0'
par2 <- 'TRUE'
par1 <- '3'
par1 <- as.numeric(par1) #colour
par2<- as.logical(par2) # Notches ?
par3<-as.numeric(par3) # % trim
if(par3>45){par3<-45;warning('trim limited to 45%')}
if(par3<0){par3<-0;warning('negative trim makes no sense. Trim is zero.')}
lotrm<-as.integer(length(y[1,])*par3/100)+1
hitrm<-as.integer(length(y[1,])*(100-par3)/100)
y1<-array(dim=c(dim(y)[1], hitrm-lotrm+1), dimnames=list(dimnames(y)[[1]], 1:(hitrm-lotrm+1) ))
for(i in 1:dim(y)[1]){
tmp<-order(y[i,])
y1[i,]<- y[i, tmp[lotrm:hitrm] ]
}
bitmap(file='test2.png')
pairs(t(y))
dev.off()
y<-y1
z <- as.data.frame(t(y))
bitmap(file='test1.png')
(r<-boxplot(z ,xlab=xlab,ylab=ylab,main=main,notch=par2,col=par1))
dev.off()
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,hyperlink('overview.htm','Boxplot statistics','Boxplot overview'),6,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Variable',1,TRUE)
a<-table.element(a,hyperlink('lower_whisker.htm','lower whisker','definition of lower whisker'),1,TRUE)
a<-table.element(a,hyperlink('lower_hinge.htm','lower hinge','definition of lower hinge'),1,TRUE)
a<-table.element(a,hyperlink('central_tendency.htm','median','definitions about measures of central tendency'),1,TRUE)
a<-table.element(a,hyperlink('upper_hinge.htm','upper hinge','definition of upper hinge'),1,TRUE)
a<-table.element(a,hyperlink('upper_whisker.htm','upper whisker','definition of upper whisker'),1,TRUE)
a<-table.row.end(a)
for (i in 1:length(y[,1]))
{
a<-table.row.start(a)
a<-table.element(a,dimnames(t(x))[[2]][i],1,TRUE)
for (j in 1:5)
{
a<-table.element(a,round(r$stats[j,i], digits=2))
}
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable.tab')
if (par2){
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Boxplot Notches',4,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Variable',1,TRUE)
a<-table.element(a,'lower bound',1,TRUE)
a<-table.element(a,'median',1,TRUE)
a<-table.element(a,'upper bound',1,TRUE)
a<-table.row.end(a)
for (i in 1:length(y[,1]))
{
a<-table.row.start(a)
a<-table.element(a,dimnames(t(x))[[2]][i],1,TRUE)
a<-table.element(a,round(r$conf[1,i], digits=2))
a<-table.element(a, round(r$stats[3,i], digits=2))
a<-table.element(a,round(r$conf[2,i], digits=2))
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable1.tab')
}
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Boxplot Means',4,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Variable',1,TRUE)
a<-table.element(a,hyperlink('trimmed_mean.htm','trimmed mean','definition of trimmed mean'),1,TRUE)
a<-table.element(a,hyperlink('unbiased1.htm','unbiased SD','definition of unbiased SD'),1,TRUE)
a<-table.row.end(a)
for (i in 1:length(y[,1]))
{
a<-table.row.start(a)
a<-table.element(a,dimnames(t(x))[[2]][i],1,TRUE)
a<-table.element(a,round(mean(unlist(z[i]), trim=par3/100, na.rm=TRUE), digits=2))
a<-table.element(a,round(sd(unlist(z[i]), na.rm=TRUE), digits=2))
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable2.tab')

Free Statistics

Description of Statistical Computation

Tree of Dependent Computations

Dataset

Tables (Output of Computation)

Figures (Output of Computation)

Input Parameters & R Code