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

Tue, 28 Oct 2014 11:07:58 +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/Oct/28/t1414494497e51dmtgdxn91k42.htm/, Retrieved Sat, 11 May 2024 21:13:31 +0000

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=247888, Retrieved Sat, 11 May 2024 21:13:31 +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] [Care Age 10 Data] [2009-10-26 09:01:50] [98fd0e87c3eb04e0cc2efde01dbafab6]
-   PD  [Boxplot and Trimmed Means] [Care Age 7 Data] [2009-10-26 18:36:29] [98fd0e87c3eb04e0cc2efde01dbafab6]
-   P     [CARE Data - Boxplots and Scatterplot Matrix] [CARE Data] [2010-10-19 14:16:27] [3fdd735c61ad38cbc9b3393dc997cdb7]
- RM        [CARE Data - Boxplots and Scatterplot Matrix] [CARE data - works...] [2011-10-17 10:23:12] [98fd0e87c3eb04e0cc2efde01dbafab6]
- RMP         [Boxplot and Trimmed Means] [CARE Study Age 7 ] [2013-10-17 12:59:45] [34296d8f7657c52ed60d5bff9133afec]
- RM D            [Boxplot and Trimmed Means] [year 10] [2014-10-28 11:07:58] [bb19d6e3e90acb3dc02677e1de219dc1] [Current]
-  M                [Boxplot and Trimmed Means] [year 10] [2014-10-31 13:06:56] [893db485fc2886352518738b0ddc43e2] 

Feedback Forum

Post a new message

Dataseries X:

Download CSV

Histogram

Boxplots

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

Summary of computational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	1 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 & 1 seconds \tabularnewline
R Server & 'Sir Ronald Aylmer Fisher' @ fisher.wessa.net \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=247888&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]1 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=247888&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=247888&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	1 seconds
R Server	'Sir Ronald Aylmer Fisher' @ fisher.wessa.net

Boxplot statistics
Variable	lower whisker	lower hinge	median	upper hinge	upper whisker
WJ7ARD	73	85	91	98	113
WJ7AMA	72	82	91	97	112
WJ7AWL	69	86	92	102	114
WJ7AKN	69	84.5	91	98	110
WJ7ASK	73	84	91	97	114

\begin{tabular}{lllllllll}
\hline
Boxplot statistics \tabularnewline
Variable & lower whisker & lower hinge & median & upper hinge & upper whisker \tabularnewline
WJ7ARD & 73 & 85 & 91 & 98 & 113 \tabularnewline
WJ7AMA & 72 & 82 & 91 & 97 & 112 \tabularnewline
WJ7AWL & 69 & 86 & 92 & 102 & 114 \tabularnewline
WJ7AKN & 69 & 84.5 & 91 & 98 & 110 \tabularnewline
WJ7ASK & 73 & 84 & 91 & 97 & 114 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=247888&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]WJ7ARD[/C][C]73[/C][C]85[/C][C]91[/C][C]98[/C][C]113[/C][/ROW]
[ROW][C]WJ7AMA[/C][C]72[/C][C]82[/C][C]91[/C][C]97[/C][C]112[/C][/ROW]
[ROW][C]WJ7AWL[/C][C]69[/C][C]86[/C][C]92[/C][C]102[/C][C]114[/C][/ROW]
[ROW][C]WJ7AKN[/C][C]69[/C][C]84.5[/C][C]91[/C][C]98[/C][C]110[/C][/ROW]
[ROW][C]WJ7ASK[/C][C]73[/C][C]84[/C][C]91[/C][C]97[/C][C]114[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=247888&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=247888&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
WJ7ARD	73	85	91	98	113
WJ7AMA	72	82	91	97	112
WJ7AWL	69	86	92	102	114
WJ7AKN	69	84.5	91	98	110
WJ7ASK	73	84	91	97	114

Boxplot Notches
Variable	lower bound	median	upper bound
WJ7ARD	89.62	91	92.38
WJ7AMA	89.41	91	92.59
WJ7AWL	90.31	92	93.69
WJ7AKN	89.57	91	92.43
WJ7ASK	89.62	91	92.38

\begin{tabular}{lllllllll}
\hline
Boxplot Notches \tabularnewline
Variable & lower bound & median & upper bound \tabularnewline
WJ7ARD & 89.62 & 91 & 92.38 \tabularnewline
WJ7AMA & 89.41 & 91 & 92.59 \tabularnewline
WJ7AWL & 90.31 & 92 & 93.69 \tabularnewline
WJ7AKN & 89.57 & 91 & 92.43 \tabularnewline
WJ7ASK & 89.62 & 91 & 92.38 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=247888&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]WJ7ARD[/C][C]89.62[/C][C]91[/C][C]92.38[/C][/ROW]
[ROW][C]WJ7AMA[/C][C]89.41[/C][C]91[/C][C]92.59[/C][/ROW]
[ROW][C]WJ7AWL[/C][C]90.31[/C][C]92[/C][C]93.69[/C][/ROW]
[ROW][C]WJ7AKN[/C][C]89.57[/C][C]91[/C][C]92.43[/C][/ROW]
[ROW][C]WJ7ASK[/C][C]89.62[/C][C]91[/C][C]92.38[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=247888&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=247888&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
WJ7ARD	89.62	91	92.38
WJ7AMA	89.41	91	92.59
WJ7AWL	90.31	92	93.69
WJ7AKN	89.57	91	92.43
WJ7ASK	89.62	91	92.38

Boxplot Means
Variable	trimmed mean	unbiased SD
WJ7ARD	90.93	9.24
WJ7AMA	90.86	9.98
WJ7AWL	92.96	10.49
WJ7AKN	91.85	9.45
WJ7ASK	90.96	9.67

\begin{tabular}{lllllllll}
\hline
Boxplot Means \tabularnewline
Variable & trimmed mean & unbiased SD \tabularnewline
WJ7ARD & 90.93 & 9.24 \tabularnewline
WJ7AMA & 90.86 & 9.98 \tabularnewline
WJ7AWL & 92.96 & 10.49 \tabularnewline
WJ7AKN & 91.85 & 9.45 \tabularnewline
WJ7ASK & 90.96 & 9.67 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=247888&T=3

[TABLE]
[ROW][C]Boxplot Means[/C][/ROW]
[ROW][C]Variable[/C][C]trimmed mean[/C][C]unbiased SD[/C][/ROW]
[ROW][C]WJ7ARD[/C][C]90.93[/C][C]9.24[/C][/ROW]
[ROW][C]WJ7AMA[/C][C]90.86[/C][C]9.98[/C][/ROW]
[ROW][C]WJ7AWL[/C][C]92.96[/C][C]10.49[/C][/ROW]
[ROW][C]WJ7AKN[/C][C]91.85[/C][C]9.45[/C][/ROW]
[ROW][C]WJ7ASK[/C][C]90.96[/C][C]9.67[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=247888&T=3

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=247888&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
WJ7ARD	90.93	9.24
WJ7AMA	90.86	9.98
WJ7AWL	92.96	10.49
WJ7AKN	91.85	9.45
WJ7ASK	90.96	9.67

Figure 1

PNG link

Postscript link

PDF link

Figure 2

PNG link

Postscript link

PDF link

Parameters (Session):

par1 = 2 ; par2 = TRUE ; par3 = 5 ;

Parameters (R input):

par1 = 2 ; par2 = TRUE ; par3 = 5 ;

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(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