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

Tue, 22 Oct 2013 08:41:59 -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/22/t1382445760q2zh9p5wf7679re.htm/, Retrieved Sun, 28 Apr 2024 23:45:42 +0000

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=218430, Retrieved Sun, 28 Apr 2024 23:45:42 +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

124

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  D          [Boxplot and Trimmed Means] [Care study age 7 ...] [2013-10-22 00:40:18] [74be16979710d4c4e7c6647856088456]
-   PD            [Boxplot and Trimmed Means] [CARE Study age 10...] [2013-10-22 01:51:50] [74be16979710d4c4e7c6647856088456]
-   P                 [Boxplot and Trimmed Means] [Year 10 False] [2013-10-22 12:41:59] [d41d8cd98f00b204e9800998ecf8427e] [Current]

Feedback Forum

Post a new message

Dataseries X:

Download CSV

Histogram

Boxplots

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

Summary of computational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	3 seconds
R Server	'Sir Ronald Aylmer Fisher' @ fisher.wessa.net
R Framework error message	The field 'Names of X columns' contains a hard return which cannot be interpreted. Please, resubmit your request without hard returns in the 'Names of X columns'.

\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 & 'Sir Ronald Aylmer Fisher' @ fisher.wessa.net \tabularnewline
R Framework error message & The field 'Names of X columns' contains a hard return which cannot be interpreted.
Please, resubmit your request without hard returns in the 'Names of X columns'. \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=218430&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]'Sir Ronald Aylmer Fisher' @ fisher.wessa.net[/C][/ROW]
[ROW][C]R Framework error message[/C][C]The field 'Names of X columns' contains a hard return which cannot be interpreted.
Please, resubmit your request without hard returns in the 'Names of X columns'.[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=218430&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=218430&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	3 seconds
R Server	'Sir Ronald Aylmer Fisher' @ fisher.wessa.net
R Framework error message	The field 'Names of X columns' contains a hard return which cannot be interpreted. Please, resubmit your request without hard returns in the 'Names of X columns'.

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=218430&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=218430&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=218430&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=218430&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=218430&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=218430&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