Repository of Reproducible Computations

Free Statistics

of Irreproducible Research!

Author's title

Author

*Unverified author*

R Software Module

rwasp_BoxplotAndScatterplotMatrix.wasp

Title produced by software

Boxplot and Trimmed Means

Date of computation

Mon, 21 Oct 2013 17:43:29 -0400

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/2013/Oct/21/t1382391833tiy7yg8o0rabl09.htm/, Retrieved Sun, 28 Apr 2024 20:40:16 +0000

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=217843, Retrieved Sun, 28 Apr 2024 20:40:16 +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]
- R PD            [Boxplot and Trimmed Means] [Year 10- Trim 20%] [2013-10-21 21:43:29] [d41d8cd98f00b204e9800998ecf8427e] [Current]

Feedback Forum

Post a new message

Dataseries X:

Download CSV

Histogram

Boxplots

95	81	84
87	90	97
91	91	88
84	87	87
92	89	92
84	85	91
84	81	80
69	80	73
84	75	73
84	84	77
85	91	86
89	87	82
83	90	82
86	72	73
91	96	99
90	88	90
77	73	78
88	86	100
88	67	73
88	89	95
102	94	93
103	88	91
84	77	81
87	89	98
88	104	88
93	89	90
101	107	85
88	96	79
113	111	107
95	84	82
99	83	98
85	82	84
105	105	99
105	100	89
96	111	105
88	90	94
93	87	92
93	90	81
69	67	NA
83	84	86
83	84	90
93	96	101
69	67	73
86	90	88
90	83	90
87	78	84
92	100	89
75	67	73
89	94	98
88	89	82
87	89	77
113	111	107
102	88	84
101	94	93
87	111	96
105	111	107
101	105	95
79	97	89
86	83	82
95	105	92
104	69	100
79	80	75
86	87	94
81	67	83
92	92	92
95	102	102
97	90	98
87	98	94
94	94	83
105	101	107
94	96	107
88	94	91
105	105	83
113	105	100
69	84	77
69	72	74
100	96	92
113	94	100
88	107	92
104	103	93
94	95	89
91	73	86
103	110	97
95	84	82
69	67	73
86	73	93
108	84	97
83	79	93
102	91	89
94	92	86
113	107	89
113	111	107
93	89	90
105	110	104
88	98	96
90	95	93
92	84	104
81	73	83
86	86	92
85	98	83
86	90	83
97	97	96
103	106	84
102	106	107

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	'George Udny Yule' @ yule.wessa.net
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 & 'George Udny Yule' @ yule.wessa.net \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=217843&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]'George Udny Yule' @ yule.wessa.net[/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=217843&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=217843&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	'George Udny Yule' @ yule.wessa.net
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	84	87	90	94	102
WJ10AMA	81	86.5	90	94	101
WJ10AKN	82	86	90	93	98

\begin{tabular}{lllllllll}
\hline
Boxplot statistics \tabularnewline
Variable & lower whisker & lower hinge & median & upper hinge & upper whisker \tabularnewline
WJ10ARD & 84 & 87 & 90 & 94 & 102 \tabularnewline
WJ10AMA & 81 & 86.5 & 90 & 94 & 101 \tabularnewline
WJ10AKN & 82 & 86 & 90 & 93 & 98 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=217843&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]84[/C][C]87[/C][C]90[/C][C]94[/C][C]102[/C][/ROW]
[ROW][C]WJ10AMA[/C][C]81[/C][C]86.5[/C][C]90[/C][C]94[/C][C]101[/C][/ROW]
[ROW][C]WJ10AKN[/C][C]82[/C][C]86[/C][C]90[/C][C]93[/C][C]98[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=217843&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=217843&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	84	87	90	94	102
WJ10AMA	81	86.5	90	94	101
WJ10AKN	82	86	90	93	98

Boxplot Means
Variable	trimmed mean	unbiased SD
WJ10ARD	90.31	4.72
WJ10AMA	90.08	5.14
WJ10AKN	89.77	4.65

\begin{tabular}{lllllllll}
\hline
Boxplot Means \tabularnewline
Variable & trimmed mean & unbiased SD \tabularnewline
WJ10ARD & 90.31 & 4.72 \tabularnewline
WJ10AMA & 90.08 & 5.14 \tabularnewline
WJ10AKN & 89.77 & 4.65 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=217843&T=2

[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]90.31[/C][C]4.72[/C][/ROW]
[ROW][C]WJ10AMA[/C][C]90.08[/C][C]5.14[/C][/ROW]
[ROW][C]WJ10AKN[/C][C]89.77[/C][C]4.65[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=217843&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=217843&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
WJ10ARD	90.31	4.72
WJ10AMA	90.08	5.14
WJ10AKN	89.77	4.65

Figure 1

PNG link

Postscript link

PDF link

Figure 2

PNG link

Postscript link

PDF link

Parameters (Session):

par1 = 2 ; par2 = FALSE ; par3 = 20 ;

Parameters (R input):

par1 = 2 ; par2 = FALSE ; par3 = 20 ;

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