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 16:32:29 -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/t12565963986q7j1y36m5yuzco.htm/, Retrieved Thu, 02 May 2024 15:42:17 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=50785, Retrieved Thu, 02 May 2024 15:42:17 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact123
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]
-   PD    [Boxplot and Trimmed Means] [Compendium 3 year...] [2009-10-26 22:32:29] [1432877eae051137c16816da1447f292] [Current]
Feedback Forum

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




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

\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
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=50785&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]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=50785&T=0

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







Boxplot statistics
Variablelower whiskerlower hingemedianupper hingeupper whisker
WJ10ARD8891929598
WJ10AMA8588919295
WJ10AWL8688909395
WJ10AKN8488909296
WJ10ASK8488909294

\begin{tabular}{lllllllll}
\hline
Boxplot statistics \tabularnewline
Variable & lower whisker & lower hinge & median & upper hinge & upper whisker \tabularnewline
WJ10ARD & 88 & 91 & 92 & 95 & 98 \tabularnewline
WJ10AMA & 85 & 88 & 91 & 92 & 95 \tabularnewline
WJ10AWL & 86 & 88 & 90 & 93 & 95 \tabularnewline
WJ10AKN & 84 & 88 & 90 & 92 & 96 \tabularnewline
WJ10ASK & 84 & 88 & 90 & 92 & 94 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=50785&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]88[/C][C]91[/C][C]92[/C][C]95[/C][C]98[/C][/ROW]
[ROW][C]WJ10AMA[/C][C]85[/C][C]88[/C][C]91[/C][C]92[/C][C]95[/C][/ROW]
[ROW][C]WJ10AWL[/C][C]86[/C][C]88[/C][C]90[/C][C]93[/C][C]95[/C][/ROW]
[ROW][C]WJ10AKN[/C][C]84[/C][C]88[/C][C]90[/C][C]92[/C][C]96[/C][/ROW]
[ROW][C]WJ10ASK[/C][C]84[/C][C]88[/C][C]90[/C][C]92[/C][C]94[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=50785&T=1

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







Boxplot Means
Variabletrimmed meanunbiased SD
WJ10ARD92.353.1
WJ10AMA90.353.09
WJ10AWL89.942.79
WJ10AKN89.763.21
WJ10ASK90.182.91

\begin{tabular}{lllllllll}
\hline
Boxplot Means \tabularnewline
Variable & trimmed mean & unbiased SD \tabularnewline
WJ10ARD & 92.35 & 3.1 \tabularnewline
WJ10AMA & 90.35 & 3.09 \tabularnewline
WJ10AWL & 89.94 & 2.79 \tabularnewline
WJ10AKN & 89.76 & 3.21 \tabularnewline
WJ10ASK & 90.18 & 2.91 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=50785&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]92.35[/C][C]3.1[/C][/ROW]
[ROW][C]WJ10AMA[/C][C]90.35[/C][C]3.09[/C][/ROW]
[ROW][C]WJ10AWL[/C][C]89.94[/C][C]2.79[/C][/ROW]
[ROW][C]WJ10AKN[/C][C]89.76[/C][C]3.21[/C][/ROW]
[ROW][C]WJ10ASK[/C][C]90.18[/C][C]2.91[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=50785&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=50785&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
WJ10ARD92.353.1
WJ10AMA90.353.09
WJ10AWL89.942.79
WJ10AKN89.763.21
WJ10ASK90.182.91



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