Free Statistics

of Irreproducible Research!

Author's title

Author*The author of this computation has been verified*
R Software ModuleIan.Hollidayrwasp_CARE Data Boxplot.wasp
Title produced by softwareCARE Data - Boxplots and Scatterplot Matrix
Date of computationSat, 23 Oct 2010 13:28:10 +0000
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2010/Oct/23/t1287840490ofzxetgig5z8ryl.htm/, Retrieved Sun, 05 May 2024 05:44:27 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=87685, Retrieved Sun, 05 May 2024 05:44:27 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact99
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]
-   PD        [CARE Data - Boxplots and Scatterplot Matrix] [Compendium 3] [2010-10-23 13:28:10] [2371c5d03e6d993afc5955d1ad229356] [Current]
Feedback Forum

Post a new message
Dataseries X:
82	81	83	84	82
82	82	83	84	83
82	82	84	84	83
82	82	84	84	83
82	82	84	85	84
83	82	84	86	84
83	82	84	86	84
84	82	84	86	84
84	82	84	86	85
84	82	85	87	85
84	82	86	87	85
84	84	86	87	85
84	85	86	88	85
84	85	86	88	85
84	85	86	88	85
85	86	86	88	85
85	87	87	88	85
86	87	87	88	85
86	87	87	88	85
86	87	87	88	86
86	88	87	88	86
86	88	87	89	86
86	88	87	89	86
86	88	88	89	87
86	88	88	89	87
86	88	88	89	87
86	89	88	89	87
86	89	88	89	87
86	89	89	89	87
86	90	89	90	87
87	90	89	90	88
87	90	90	90	88
87	90	90	91	88
87	90	90	91	88
87	91	90	91	88
87	91	90	91	89
88	91	90	91	89
88	91	91	91	89
88	91	91	92	89
88	91	92	92	90
89	91	92	92	90
89	91	92	92	90
89	91	92	92	91
89	92	93	92	91
89	92	94	93	92
89	92	95	93	92
89	92	95	93	92
89	92	95	93	92
90	92	95	94	92
90	92	95	94	92
90	93	95	94	92
90	93	95	94	92
90	93	96	94	93
90	93	96	95	93
90	93	96	95	93
90	93	96	95	93
90	95	96	95	93
91	95	97	96	94
91	95	97	96	94
91	95	97	96	94
91	95	98	96	95
91	95	98	96	95
92	95	98	97	95
92	95	98	97	95
92	95	99	97	95
92	95	99	97	95
92	95	99	97	95
92	95	99	97	97
92	96	100	97	97
94	97	100	98	97
94	97	100	98	97
94	97	101	98	97
95	97	102	98	97
95	97	102	98	97
95	97	102	98	97
96	97	102	98	97
96	97	103	98	98
97	97	103	99	99
97	97	103	99	99
97	98	104	99	99









Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time1 seconds
R Server'Gwilym Jenkins' @ 72.249.127.135
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 & 1 seconds \tabularnewline
R Server & 'Gwilym Jenkins' @ 72.249.127.135 \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=87685&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]'Gwilym Jenkins' @ 72.249.127.135[/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=87685&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=87685&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 Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time1 seconds
R Server'Gwilym Jenkins' @ 72.249.127.135
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
Variablelower whiskerlower hingemedianupper hingeupper whisker
WJ7ARD8686.588.59091
WJ7AMA88909192.595
WJ7AWL8789929597
WJ7AKN8890929496
WJ7ASK8687.5909294

\begin{tabular}{lllllllll}
\hline
Boxplot statistics \tabularnewline
Variable & lower whisker & lower hinge & median & upper hinge & upper whisker \tabularnewline
WJ7ARD & 86 & 86.5 & 88.5 & 90 & 91 \tabularnewline
WJ7AMA & 88 & 90 & 91 & 92.5 & 95 \tabularnewline
WJ7AWL & 87 & 89 & 92 & 95 & 97 \tabularnewline
WJ7AKN & 88 & 90 & 92 & 94 & 96 \tabularnewline
WJ7ASK & 86 & 87.5 & 90 & 92 & 94 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=87685&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]86[/C][C]86.5[/C][C]88.5[/C][C]90[/C][C]91[/C][/ROW]
[ROW][C]WJ7AMA[/C][C]88[/C][C]90[/C][C]91[/C][C]92.5[/C][C]95[/C][/ROW]
[ROW][C]WJ7AWL[/C][C]87[/C][C]89[/C][C]92[/C][C]95[/C][C]97[/C][/ROW]
[ROW][C]WJ7AKN[/C][C]88[/C][C]90[/C][C]92[/C][C]94[/C][C]96[/C][/ROW]
[ROW][C]WJ7ASK[/C][C]86[/C][C]87.5[/C][C]90[/C][C]92[/C][C]94[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=87685&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=87685&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
Variablelower whiskerlower hingemedianupper hingeupper whisker
WJ7ARD8686.588.59091
WJ7AMA88909192.595
WJ7AWL8789929597
WJ7AKN8890929496
WJ7ASK8687.5909294







Boxplot Notches
Variablelower boundmedianupper bound
WJ7ARD87.6388.589.37
WJ7AMA90.389191.62
WJ7AWL90.59293.5
WJ7AKN919293
WJ7ASK88.889091.12

\begin{tabular}{lllllllll}
\hline
Boxplot Notches \tabularnewline
Variable & lower bound & median & upper bound \tabularnewline
WJ7ARD & 87.63 & 88.5 & 89.37 \tabularnewline
WJ7AMA & 90.38 & 91 & 91.62 \tabularnewline
WJ7AWL & 90.5 & 92 & 93.5 \tabularnewline
WJ7AKN & 91 & 92 & 93 \tabularnewline
WJ7ASK & 88.88 & 90 & 91.12 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=87685&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]87.63[/C][C]88.5[/C][C]89.37[/C][/ROW]
[ROW][C]WJ7AMA[/C][C]90.38[/C][C]91[/C][C]91.62[/C][/ROW]
[ROW][C]WJ7AWL[/C][C]90.5[/C][C]92[/C][C]93.5[/C][/ROW]
[ROW][C]WJ7AKN[/C][C]91[/C][C]92[/C][C]93[/C][/ROW]
[ROW][C]WJ7ASK[/C][C]88.88[/C][C]90[/C][C]91.12[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=87685&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=87685&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
Variablelower boundmedianupper bound
WJ7ARD87.6388.589.37
WJ7AMA90.389191.62
WJ7AWL90.59293.5
WJ7AKN919293
WJ7ASK88.889091.12







Boxplot Means
Variabletrimmed meanunbiased SD
WJ7ARD88.31.72
WJ7AMA91.152.06
WJ7AWL92.053.32
WJ7AKN91.92.34
WJ7ASK902.59

\begin{tabular}{lllllllll}
\hline
Boxplot Means \tabularnewline
Variable & trimmed mean & unbiased SD \tabularnewline
WJ7ARD & 88.3 & 1.72 \tabularnewline
WJ7AMA & 91.15 & 2.06 \tabularnewline
WJ7AWL & 92.05 & 3.32 \tabularnewline
WJ7AKN & 91.9 & 2.34 \tabularnewline
WJ7ASK & 90 & 2.59 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=87685&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]88.3[/C][C]1.72[/C][/ROW]
[ROW][C]WJ7AMA[/C][C]91.15[/C][C]2.06[/C][/ROW]
[ROW][C]WJ7AWL[/C][C]92.05[/C][C]3.32[/C][/ROW]
[ROW][C]WJ7AKN[/C][C]91.9[/C][C]2.34[/C][/ROW]
[ROW][C]WJ7ASK[/C][C]90[/C][C]2.59[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=87685&T=3

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=87685&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
Variabletrimmed meanunbiased SD
WJ7ARD88.31.72
WJ7AMA91.152.06
WJ7AWL92.053.32
WJ7AKN91.92.34
WJ7ASK902.59



Parameters (Session):
par1 = 3 ; par2 = TRUE ; par3 = 25 ;
Parameters (R input):
par1 = 3 ; par2 = TRUE ; par3 = 25 ;
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, na.rm=TRUE), digits=2))
a<-table.element(a,round(sd(z[i], na.rm=TRUE), digits=2))
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable2.tab')