Free Statistics

of Irreproducible Research!

Author's title

Author*The author of this computation has been verified*
R Software Modulerwasp_CARE Data Boxplot.wasp
Title produced by softwareCARE Data - Boxplots and Scatterplot Matrix
Date of computationSun, 14 Oct 2012 08:24:02 -0400
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2012/Oct/14/t1350217470e3e3zmhezje5xcp.htm/, Retrieved Sun, 28 Apr 2024 06:36:29 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=176110, Retrieved Sun, 28 Apr 2024 06:36:29 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact56
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]
- R  D          [CARE Data - Boxplots and Scatterplot Matrix] [year 10- 3 variables] [2012-10-14 12:24:02] [4ab5704ae4f56c211390cb70e8a77bb1] [Current]
Feedback Forum

Post a new message
Dataseries X:
71	69	67
71	69	67
71	69	67
71	69	67
71	69	67
71	69	67
72	75	69
72	77	72
73	79	72
77	79	73
77	81	73
78	81	73
78	83	73
79	83	75
79	83	77
79	83	78
80	84	79
80	84	80
80	84	80
80	84	81
81	84	81
81	84	82
81	85	83
82	85	83
82	85	83
82	86	84
82	86	84
82	86	84
82	86	84
84	86	84
84	86	84
85	86	84
86	87	84
86	87	85
86	87	86
86	87	86
86	87	87
86	87	87
86	88	87
86	88	87
87	88	88
88	88	88
88	88	88
88	88	89
88	88	89
89	88	89
89	88	89
89	88	89
90	89	89
90	89	89
90	90	90
91	90	90
91	90	90
91	91	90
91	91	90
91	91	90
92	92	90
92	92	91
92	92	91
92	92	91
92	93	92
92	93	92
92	93	94
92	93	94
93	93	94
94	94	94
94	94	94
95	94	94
95	94	95
95	95	95
95	95	96
95	95	96
96	95	96
96	95	96
96	96	96
97	97	97
97	97	97
97	99	98
98	100	98
98	101	98
99	101	100
100	101	100
100	102	101
100	102	102
102	102	103
103	102	104
103	103	105
103	103	105
104	103	105
105	104	105
105	104	105
105	105	106
105	105	106
105	105	107
108	105	107
109	105	107
111	105	110
112	108	110
112	113	111
112	113	111
112	113	111
112	113	111
112	113	111
112	113	111




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
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 & 2 seconds \tabularnewline
R Server & 'Sir Maurice George Kendall' @ kendall.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=176110&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]
[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=176110&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=176110&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
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
Variablelower whiskerlower hingemedianupper hingeupper whisker
WJ10AVA71829197.5112
WJ10ARD69869099.5113
WJ10AMA 67849098111

\begin{tabular}{lllllllll}
\hline
Boxplot statistics \tabularnewline
Variable & lower whisker & lower hinge & median & upper hinge & upper whisker \tabularnewline
WJ10AVA & 71 & 82 & 91 & 97.5 & 112 \tabularnewline
WJ10ARD & 69 & 86 & 90 & 99.5 & 113 \tabularnewline
WJ10AMA
 & 67 & 84 & 90 & 98 & 111 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=176110&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]WJ10AVA[/C][C]71[/C][C]82[/C][C]91[/C][C]97.5[/C][C]112[/C][/ROW]
[ROW][C]WJ10ARD[/C][C]69[/C][C]86[/C][C]90[/C][C]99.5[/C][C]113[/C][/ROW]
[ROW][C]WJ10AMA
[/C][C]67[/C][C]84[/C][C]90[/C][C]98[/C][C]111[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=176110&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=176110&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
WJ10AVA71829197.5112
WJ10ARD69869099.5113
WJ10AMA 67849098111







Boxplot Means
Variabletrimmed meanunbiased SD
WJ10AVA90.7211.14
WJ10ARD91.5310.4
WJ10AMA 90.3511.84

\begin{tabular}{lllllllll}
\hline
Boxplot Means \tabularnewline
Variable & trimmed mean & unbiased SD \tabularnewline
WJ10AVA & 90.72 & 11.14 \tabularnewline
WJ10ARD & 91.53 & 10.4 \tabularnewline
WJ10AMA
 & 90.35 & 11.84 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=176110&T=2

[TABLE]
[ROW][C]Boxplot Means[/C][/ROW]
[ROW][C]Variable[/C][C]trimmed mean[/C][C]unbiased SD[/C][/ROW]
[ROW][C]WJ10AVA[/C][C]90.72[/C][C]11.14[/C][/ROW]
[ROW][C]WJ10ARD[/C][C]91.53[/C][C]10.4[/C][/ROW]
[ROW][C]WJ10AMA
[/C][C]90.35[/C][C]11.84[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=176110&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=176110&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
Variabletrimmed meanunbiased SD
WJ10AVA90.7211.14
WJ10ARD91.5310.4
WJ10AMA 90.3511.84



Parameters (Session):
par1 = 3 ; par2 = FALSE ; par3 = 0 ;
Parameters (R input):
par1 = 3 ; par2 = FALSE ; par3 = 0 ;
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')