Repository of Reproducible Computations

Free Statistics

of Irreproducible Research!

Author's title

Author

*The author of this computation has been verified*

R Software Module

Ian.Hollidayrwasp_CARE Data Boxplot.wasp

Title produced by software

Boxplot and Trimmed Means

Date of computation

Tue, 27 Oct 2009 02:06:59 +0100

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/2009/Oct/27/t12566056724c2kggjo50yp3wa.htm/, Retrieved Tue, 07 May 2024 23:54:06 +0000

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=50825, Retrieved Tue, 07 May 2024 23:54:06 +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

134

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

-     [Boxplot and Trimmed Means] [Care Age 7 Data] [2009-10-18 00:49:58] [98fd0e87c3eb04e0cc2efde01dbafab6]
-   PD    [Boxplot and Trimmed Means] [comparison of var...] [2009-10-27 01:06:59] [04ab4c38fa9b627dffda5731b967a02a] [Current]

Feedback Forum

Post a new message

Dataseries X:

Download CSV

Histogram

Boxplots

113	111	107	117	114	119
113	111	107	117	114	119
113	111	107	117	114	119
113	111	107	117	114	119
113	111	107	114	113	114
113	111	107	113	113	113
108	110	107	113	113	113
105	110	107	113	111	110
105	107	105	112	110	110
105	107	104	112	110	109
105	107	104	112	110	108
105	106	102	110	110	108
105	106	101	107	109	108
104	105	100	106	109	108
104	105	100	105	108	108
103	105	100	105	106	108
103	105	100	104	106	106
103	105	99	104	105	105
102	104	99	104	105	104
102	103	98	102	105	104
102	102	98	102	105	103
102	101	98	102	105	103
101	100	98	102	105	103
101	100	97	102	104	103
101	98	97	102	104	102
100	98	97	99	104	102
99	98	96	99	104	102
97	97	96	99	102	102
97	97	96	99	102	102
96	96	95	98	102	102
95	96	95	98	102	102
95	96	94	98	102	101
95	96	94	98	102	101
95	96	94	98	102	101
95	95	93	98	102	101
94	95	93	98	102	100
94	94	93	98	100	100
94	94	93	98	100	100
94	94	93	98	100	100
93	94	93	98	99	100
93	94	92	98	98	99
93	94	92	97	98	99
93	92	92	97	97	99
93	92	92	97	97	99
92	91	92	96	97	98
92	91	92	96	97	98
92	91	92	95	97	98
92	90	91	95	97	98
91	90	91	95	97	98
91	90	91	94	97	98
91	90	90	94	97	98
90	90	90	94	97	98
90	90	90	92	96	97
90	90	90	92	95	97
89	89	90	92	95	97
89	89	89	92	95	97
88	89	89	92	95	97
88	89	89	92	95	97
88	89	89	92	95	97
88	89	89	91	95	96
88	89	89	91	95	96
88	88	88	91	95	96
88	88	88	91	95	96
88	88	88	91	95	96
88	87	87	90	95	95
88	87	86	90	93	95
87	87	86	90	93	95
87	87	86	90	93	95
87	86	86	90	93	94
87	86	85	90	93	94
87	85	84	90	93	94
87	84	84	90	92	94
86	84	84	90	92	94
86	84	84	89	92	93
86	84	84	89	92	93
86	84	83	89	92	93
86	84	83	89	92	93
86	84	83	89	92	92
86	84	83	89	91	92
85	83	83	89	91	92
85	83	83	89	91	92
85	83	82	88	91	92
84	82	82	88	91	92
84	81	82	88	91	91
84	81	82	88	91	91
84	80	82	87	91	91
84	80	82	87	91	91
84	79	81	87	90	91
83	78	81	87	90	91
83	77	80	87	90	90
83	75	79	87	90	90
83	73	78	86	90	90
81	73	77	86	89	89
81	73	77	86	89	89
79	73	77	86	89	89
79	72	75	86	88	89
77	72	74	86	88	89
75	69	73	86	88	89
69	67	73	86	88	89
69	67	73	86	88	89
69	67	73	86	88	88
69	67	73	86	87	88
69	67	73	86	87	88
69	67	73	86	87	88
			85	87	88
			85	86	88
			84	85	88
			84	85	88
			84	85	88
			84	84	87
			84	82	87
			84	82	87
			84	82	86
			84	82	86
			83	82	86
			83	82	86
			82	82	85
			82	82	84
			82	82	84
			82	82	84
			82	81	84
			82	81	83
			81	81	83
			81	80	83
			81	80	83
			80	79	83
			80	79	83
			80	79	83
			80	78	83
			80	78	82
			80	78	82
			80	78	82
			80	77	82
			80	77	82
			78	75	81
			78	75	81
			78	75	81
			77	75	80
			77	75	80
			77	75	79
			77	75	79
			75	74	78
			75	74	77
			75	73	77
			75	72	77
			74	72	75
			73	72	73
			73	72	73
			71	72	72
			71	72	72
			70	72	70
			70	72	70
			68	72	69
			68	72	69
			63	71	69
			62	71	68
			59	70	67
			58	68	64
			58	62	63
			58	62	63
			58	62	63
			58	62	63

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	'RServer@AstonUniversity' @ vre.aston.ac.uk
R Framework error message	Warning: there are blank lines in the 'Data X' field. Please, use NA for missing data - blank lines are simply deleted and are NOT treated as missing values.

\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 & 'RServer@AstonUniversity' @ vre.aston.ac.uk \tabularnewline
R Framework error message & Warning: there are blank lines in the 'Data X' field.
Please, use NA for missing data - blank lines are simply
 deleted and are NOT treated as missing values. \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=50825&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]'RServer@AstonUniversity' @ vre.aston.ac.uk[/C][/ROW]
[ROW][C]R Framework error message[/C][C]Warning: there are blank lines in the 'Data X' field.
Please, use NA for missing data - blank lines are simply
 deleted and are NOT treated as missing values.[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=50825&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=50825&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	'RServer@AstonUniversity' @ vre.aston.ac.uk
R Framework error message	Warning: there are blank lines in the 'Data X' field. Please, use NA for missing data - blank lines are simply deleted and are NOT treated as missing values.

Boxplot statistics
Variable	lower whisker	lower hinge	median	upper hinge	upper whisker
WJ10ARD	58	84	90	102	113
WJ10AMA	62	81	90	102	111
WJ10AKN	63	83	90	98	107
WJ7ARD	62	86	93	102	117
WJ7AMA	68	89	96.5	105	114
WJ7AKN	68	89	97.5	103	119

\begin{tabular}{lllllllll}
\hline
Boxplot statistics \tabularnewline
Variable & lower whisker & lower hinge & median & upper hinge & upper whisker \tabularnewline
WJ10ARD & 58 & 84 & 90 & 102 & 113 \tabularnewline
WJ10AMA & 62 & 81 & 90 & 102 & 111 \tabularnewline
WJ10AKN & 63 & 83 & 90 & 98 & 107 \tabularnewline
WJ7ARD & 62 & 86 & 93 & 102 & 117 \tabularnewline
WJ7AMA & 68 & 89 & 96.5 & 105 & 114 \tabularnewline
WJ7AKN & 68 & 89 & 97.5 & 103 & 119 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=50825&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]WJ10ARD[/C][C]58[/C][C]84[/C][C]90[/C][C]102[/C][C]113[/C][/ROW]
[ROW][C]WJ10AMA[/C][C]62[/C][C]81[/C][C]90[/C][C]102[/C][C]111[/C][/ROW]
[ROW][C]WJ10AKN[/C][C]63[/C][C]83[/C][C]90[/C][C]98[/C][C]107[/C][/ROW]
[ROW][C]WJ7ARD[/C][C]62[/C][C]86[/C][C]93[/C][C]102[/C][C]117[/C][/ROW]
[ROW][C]WJ7AMA[/C][C]68[/C][C]89[/C][C]96.5[/C][C]105[/C][C]114[/C][/ROW]
[ROW][C]WJ7AKN[/C][C]68[/C][C]89[/C][C]97.5[/C][C]103[/C][C]119[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=50825&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=50825&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
WJ10ARD	58	84	90	102	113
WJ10AMA	62	81	90	102	111
WJ10AKN	63	83	90	98	107
WJ7ARD	62	86	93	102	117
WJ7AMA	68	89	96.5	105	114
WJ7AKN	68	89	97.5	103	119

Boxplot Notches
Variable	lower bound	median	upper bound
WJ10ARD	87.77	90	92.23
WJ10AMA	87.39	90	92.61
WJ10AKN	88.14	90	91.86
WJ7ARD	91.01	93	94.99
WJ7AMA	94.51	96.5	98.49
WJ7AKN	95.76	97.5	99.24

\begin{tabular}{lllllllll}
\hline
Boxplot Notches \tabularnewline
Variable & lower bound & median & upper bound \tabularnewline
WJ10ARD & 87.77 & 90 & 92.23 \tabularnewline
WJ10AMA & 87.39 & 90 & 92.61 \tabularnewline
WJ10AKN & 88.14 & 90 & 91.86 \tabularnewline
WJ7ARD & 91.01 & 93 & 94.99 \tabularnewline
WJ7AMA & 94.51 & 96.5 & 98.49 \tabularnewline
WJ7AKN & 95.76 & 97.5 & 99.24 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=50825&T=2

[TABLE]
[ROW][C]Boxplot Notches[/C][/ROW]
[ROW][C]Variable[/C][C]lower bound[/C][C]median[/C][C]upper bound[/C][/ROW]
[ROW][C]WJ10ARD[/C][C]87.77[/C][C]90[/C][C]92.23[/C][/ROW]
[ROW][C]WJ10AMA[/C][C]87.39[/C][C]90[/C][C]92.61[/C][/ROW]
[ROW][C]WJ10AKN[/C][C]88.14[/C][C]90[/C][C]91.86[/C][/ROW]
[ROW][C]WJ7ARD[/C][C]91.01[/C][C]93[/C][C]94.99[/C][/ROW]
[ROW][C]WJ7AMA[/C][C]94.51[/C][C]96.5[/C][C]98.49[/C][/ROW]
[ROW][C]WJ7AKN[/C][C]95.76[/C][C]97.5[/C][C]99.24[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=50825&T=2

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

As an alternative you can also use a QR Code:

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

Boxplot Notches
Variable	lower bound	median	upper bound
WJ10ARD	87.77	90	92.23
WJ10AMA	87.39	90	92.61
WJ10AKN	88.14	90	91.86
WJ7ARD	91.01	93	94.99
WJ7AMA	94.51	96.5	98.49
WJ7AKN	95.76	97.5	99.24

Boxplot Means
Variable	trimmed mean	unbiased SD
WJ10ARD	91.16	12.58
WJ10AMA	90.48	13.14
WJ10AKN	90.04	10.66
WJ7ARD	94.22	12.59
WJ7AMA	95.81	12.05
WJ7AKN	96.6	11.7

\begin{tabular}{lllllllll}
\hline
Boxplot Means \tabularnewline
Variable & trimmed mean & unbiased SD \tabularnewline
WJ10ARD & 91.16 & 12.58 \tabularnewline
WJ10AMA & 90.48 & 13.14 \tabularnewline
WJ10AKN & 90.04 & 10.66 \tabularnewline
WJ7ARD & 94.22 & 12.59 \tabularnewline
WJ7AMA & 95.81 & 12.05 \tabularnewline
WJ7AKN & 96.6 & 11.7 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=50825&T=3

[TABLE]
[ROW][C]Boxplot Means[/C][/ROW]
[ROW][C]Variable[/C][C]trimmed mean[/C][C]unbiased SD[/C][/ROW]
[ROW][C]WJ10ARD[/C][C]91.16[/C][C]12.58[/C][/ROW]
[ROW][C]WJ10AMA[/C][C]90.48[/C][C]13.14[/C][/ROW]
[ROW][C]WJ10AKN[/C][C]90.04[/C][C]10.66[/C][/ROW]
[ROW][C]WJ7ARD[/C][C]94.22[/C][C]12.59[/C][/ROW]
[ROW][C]WJ7AMA[/C][C]95.81[/C][C]12.05[/C][/ROW]
[ROW][C]WJ7AKN[/C][C]96.6[/C][C]11.7[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=50825&T=3

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

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
WJ10ARD	91.16	12.58
WJ10AMA	90.48	13.14
WJ10AKN	90.04	10.66
WJ7ARD	94.22	12.59
WJ7AMA	95.81	12.05
WJ7AKN	96.6	11.7

Figure 1

PNG link

Postscript link

PDF link

Figure 2

PNG link

Postscript link

PDF link

Parameters (Session):

par1 = 3 ; par2 = TRUE ; par3 = 0 ;

Parameters (R input):

par1 = 3 ; par2 = TRUE ; par3 = 0 ;

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

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(z[i], trim=par3/100), digits=2))
a<-table.element(a,round(sd(z[i]), 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