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 softwareBoxplot and Trimmed Means
Date of computationMon, 26 Oct 2009 06:30:56 -0600
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2009/Oct/26/t1256560303vopzlz8bhei7z7h.htm/, Retrieved Thu, 02 May 2024 14:11:38 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=50506, Retrieved Thu, 02 May 2024 14:11:38 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact138
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Boxplot and Trimmed Means] [Care Age 7 Data] [2009-10-18 00:49:58] [98fd0e87c3eb04e0cc2efde01dbafab6]
- R PD  [Boxplot and Trimmed Means] [Scatter/Boxplot W...] [2009-10-22 16:56:04] [cb1b5098fbe5e6779d164349741d91ea]
-    D      [Boxplot and Trimmed Means] [Scatter/Boxplot W...] [2009-10-26 12:30:56] [fd4f510ca22152ae67cc1e51e700ce02] [Current]
Feedback Forum

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




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time2 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 & 2 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=50506&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]'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=50506&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=50506&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'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
WJ10ARD67848996113
WJ10AMA678389.596113
WJ10AKN588188.597117
WJ7ARD58809097117
WJ7AMA63829198119
WJ7AKN67848897113

\begin{tabular}{lllllllll}
\hline
Boxplot statistics \tabularnewline
Variable & lower whisker & lower hinge & median & upper hinge & upper whisker \tabularnewline
WJ10ARD & 67 & 84 & 89 & 96 & 113 \tabularnewline
WJ10AMA & 67 & 83 & 89.5 & 96 & 113 \tabularnewline
WJ10AKN & 58 & 81 & 88.5 & 97 & 117 \tabularnewline
WJ7ARD & 58 & 80 & 90 & 97 & 117 \tabularnewline
WJ7AMA & 63 & 82 & 91 & 98 & 119 \tabularnewline
WJ7AKN & 67 & 84 & 88 & 97 & 113 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=50506&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]67[/C][C]84[/C][C]89[/C][C]96[/C][C]113[/C][/ROW]
[ROW][C]WJ10AMA[/C][C]67[/C][C]83[/C][C]89.5[/C][C]96[/C][C]113[/C][/ROW]
[ROW][C]WJ10AKN[/C][C]58[/C][C]81[/C][C]88.5[/C][C]97[/C][C]117[/C][/ROW]
[ROW][C]WJ7ARD[/C][C]58[/C][C]80[/C][C]90[/C][C]97[/C][C]117[/C][/ROW]
[ROW][C]WJ7AMA[/C][C]63[/C][C]82[/C][C]91[/C][C]98[/C][C]119[/C][/ROW]
[ROW][C]WJ7AKN[/C][C]67[/C][C]84[/C][C]88[/C][C]97[/C][C]113[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=50506&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=50506&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
WJ10ARD67848996113
WJ10AMA678389.596113
WJ10AKN588188.597117
WJ7ARD58809097117
WJ7AMA63829198119
WJ7AKN67848897113







Boxplot Means
Variabletrimmed meanunbiased SD
WJ10ARD90.0410.99
WJ10AMA89.8610.01
WJ10AKN88.5512.28
WJ7ARD89.2212.37
WJ7AMA90.312.55
WJ7AKN89.7110.87

\begin{tabular}{lllllllll}
\hline
Boxplot Means \tabularnewline
Variable & trimmed mean & unbiased SD \tabularnewline
WJ10ARD & 90.04 & 10.99 \tabularnewline
WJ10AMA & 89.86 & 10.01 \tabularnewline
WJ10AKN & 88.55 & 12.28 \tabularnewline
WJ7ARD & 89.22 & 12.37 \tabularnewline
WJ7AMA & 90.3 & 12.55 \tabularnewline
WJ7AKN & 89.71 & 10.87 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=50506&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.04[/C][C]10.99[/C][/ROW]
[ROW][C]WJ10AMA[/C][C]89.86[/C][C]10.01[/C][/ROW]
[ROW][C]WJ10AKN[/C][C]88.55[/C][C]12.28[/C][/ROW]
[ROW][C]WJ7ARD[/C][C]89.22[/C][C]12.37[/C][/ROW]
[ROW][C]WJ7AMA[/C][C]90.3[/C][C]12.55[/C][/ROW]
[ROW][C]WJ7AKN[/C][C]89.71[/C][C]10.87[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=50506&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=50506&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
WJ10ARD90.0410.99
WJ10AMA89.8610.01
WJ10AKN88.5512.28
WJ7ARD89.2212.37
WJ7AMA90.312.55
WJ7AKN89.7110.87



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), digits=2))
a<-table.element(a,round(sd(z[i]), digits=2))
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable2.tab')