Free Statistics

of Irreproducible Research!

Author's title

Author*The author of this computation has been verified*
R Software Modulerwasp_BoxplotAndScatterplotMatrix.wasp
Title produced by softwareBoxplot and Trimmed Means
Date of computationMon, 27 Oct 2014 23:27:22 +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/2014/Oct/27/t1414452478wgbwm7hx48ajo8e.htm/, Retrieved Fri, 10 May 2024 16:30:19 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=247673, Retrieved Fri, 10 May 2024 16:30:19 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact59
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's stats] [2014-10-27 23:27:22] [a053a3e71b53d87369bc27bb638db423] [Current]
-   P               [Boxplot and Trimmed Means] [30% trimming on t...] [2014-10-27 23:38:25] [0af727c7d6fddd5894090204e2de4756]
Feedback Forum

Post a new message
Dataseries X:
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




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time2 seconds
R Server'Sir Maurice George Kendall' @ kendall.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 & 2 seconds \tabularnewline
R Server & 'Sir Maurice George Kendall' @ kendall.wessa.net \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=247673&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]'Sir Maurice George Kendall' @ kendall.wessa.net[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=247673&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=247673&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 time2 seconds
R Server'Sir Maurice George Kendall' @ kendall.wessa.net







Boxplot statistics
Variablelower whiskerlower hingemedianupper hingeupper whisker
WJ10AFS78859197107
WJ10AVA7781879297
WJ10ARD77848891101
WJ10AMA7281889096
WJ10AKN738086.59199

\begin{tabular}{lllllllll}
\hline
Boxplot statistics \tabularnewline
Variable & lower whisker & lower hinge & median & upper hinge & upper whisker \tabularnewline
WJ10AFS & 78 & 85 & 91 & 97 & 107 \tabularnewline
WJ10AVA & 77 & 81 & 87 & 92 & 97 \tabularnewline
WJ10ARD & 77 & 84 & 88 & 91 & 101 \tabularnewline
WJ10AMA & 72 & 81 & 88 & 90 & 96 \tabularnewline
WJ10AKN & 73 & 80 & 86.5 & 91 & 99 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=247673&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]WJ10AFS[/C][C]78[/C][C]85[/C][C]91[/C][C]97[/C][C]107[/C][/ROW]
[ROW][C]WJ10AVA[/C][C]77[/C][C]81[/C][C]87[/C][C]92[/C][C]97[/C][/ROW]
[ROW][C]WJ10ARD[/C][C]77[/C][C]84[/C][C]88[/C][C]91[/C][C]101[/C][/ROW]
[ROW][C]WJ10AMA[/C][C]72[/C][C]81[/C][C]88[/C][C]90[/C][C]96[/C][/ROW]
[ROW][C]WJ10AKN[/C][C]73[/C][C]80[/C][C]86.5[/C][C]91[/C][C]99[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=247673&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=247673&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
WJ10AFS78859197107
WJ10AVA7781879297
WJ10ARD77848891101
WJ10AMA7281889096
WJ10AKN738086.59199







Boxplot Notches
Variablelower boundmedianupper bound
WJ10AFS87.289194.72
WJ10AVA83.598790.41
WJ10ARD85.838890.17
WJ10AMA85.218890.79
WJ10AKN83.0986.589.91

\begin{tabular}{lllllllll}
\hline
Boxplot Notches \tabularnewline
Variable & lower bound & median & upper bound \tabularnewline
WJ10AFS & 87.28 & 91 & 94.72 \tabularnewline
WJ10AVA & 83.59 & 87 & 90.41 \tabularnewline
WJ10ARD & 85.83 & 88 & 90.17 \tabularnewline
WJ10AMA & 85.21 & 88 & 90.79 \tabularnewline
WJ10AKN & 83.09 & 86.5 & 89.91 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=247673&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]WJ10AFS[/C][C]87.28[/C][C]91[/C][C]94.72[/C][/ROW]
[ROW][C]WJ10AVA[/C][C]83.59[/C][C]87[/C][C]90.41[/C][/ROW]
[ROW][C]WJ10ARD[/C][C]85.83[/C][C]88[/C][C]90.17[/C][/ROW]
[ROW][C]WJ10AMA[/C][C]85.21[/C][C]88[/C][C]90.79[/C][/ROW]
[ROW][C]WJ10AKN[/C][C]83.09[/C][C]86.5[/C][C]89.91[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=247673&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=247673&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
WJ10AFS87.289194.72
WJ10AVA83.598790.41
WJ10ARD85.838890.17
WJ10AMA85.218890.79
WJ10AKN83.0986.589.91







Boxplot Means
Variabletrimmed meanunbiased SD
WJ10AFS90.677.81
WJ10AVA86.836.3
WJ10ARD88.085.42
WJ10AMA86.57.39
WJ10AKN85.837.72

\begin{tabular}{lllllllll}
\hline
Boxplot Means \tabularnewline
Variable & trimmed mean & unbiased SD \tabularnewline
WJ10AFS & 90.67 & 7.81 \tabularnewline
WJ10AVA & 86.83 & 6.3 \tabularnewline
WJ10ARD & 88.08 & 5.42 \tabularnewline
WJ10AMA & 86.5 & 7.39 \tabularnewline
WJ10AKN & 85.83 & 7.72 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=247673&T=3

[TABLE]
[ROW][C]Boxplot Means[/C][/ROW]
[ROW][C]Variable[/C][C]trimmed mean[/C][C]unbiased SD[/C][/ROW]
[ROW][C]WJ10AFS[/C][C]90.67[/C][C]7.81[/C][/ROW]
[ROW][C]WJ10AVA[/C][C]86.83[/C][C]6.3[/C][/ROW]
[ROW][C]WJ10ARD[/C][C]88.08[/C][C]5.42[/C][/ROW]
[ROW][C]WJ10AMA[/C][C]86.5[/C][C]7.39[/C][/ROW]
[ROW][C]WJ10AKN[/C][C]85.83[/C][C]7.72[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=247673&T=3

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=247673&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
WJ10AFS90.677.81
WJ10AVA86.836.3
WJ10ARD88.085.42
WJ10AMA86.57.39
WJ10AKN85.837.72



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')