FreeStatistics.org

Free Statistics

of Irreproducible Research!

Description of Statistical Computation

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 computationThu, 22 Oct 2009 13:43:48 -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/22/t12562406715awsupokemz3549.htm/, Retrieved Thu, 02 May 2024 16:47:27 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=49811, Retrieved Thu, 02 May 2024 16:47:27 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact112

Tree of Dependent Computations

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] [EXERCISE 4] [2009-10-22 19:43:48] [8b088281e1eac39fda503d484618b84c] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
69	67	73
69	67	73
69	67	73
69	67	73
69	67	73
69	67	73
75	69	73
77	72	74
79	72	75
79	73	77
81	73	77
81	73	77
83	73	78
83	75	79
83	77	80
83	78	81
84	79	81
84	80	82
84	80	82
84	81	82
84	81	82
84	82	82
85	83	82
85	83	83
85	83	83
86	84	83
86	84	83
86	84	83
86	84	83
86	84	84
86	84	84
86	84	84
87	84	84
87	85	84
87	86	85
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	89
88	89	89
88	89	89
88	89	89
88	89	89
89	89	89
89	89	90
90	90	90
90	90	90
90	90	90
91	90	90
91	90	91
91	90	91
92	90	91
92	91	92
92	91	92
92	91	92
93	92	92
93	92	92
93	94	92
93	94	92
93	94	93
94	94	93
94	94	93
94	94	93
94	95	93
95	95	93
95	96	94
95	96	94
95	96	94
95	96	95
96	96	95
97	97	96
97	97	96
99	98	96
100	98	97
101	98	97
101	100	97
101	100	98
102	101	98
102	102	98
102	103	98
102	104	99
103	105	99
103	105	100
103	105	100
104	105	100
104	105	100
105	106	101
105	106	102
105	107	104
105	107	104
105	107	105
105	110	107
108	110	107
113	111	107
113	111	107
113	111	107
113	111	107
113	111	107
113	111

Tables (Output of Computation)





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=49811&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=49811&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=49811&T=0

As an alternative you can also use a QR Code:  
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
WJ10ARD69869099.5113
WJ10AMA67849098111
WJ10AKN69839096107

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=49811&T=1

As an alternative you can also use a QR Code:  
QR Code


The GUIDs for individual cells are displayed in the table below:

Boxplot statistics
Variablelower whiskerlower hingemedianupper hingeupper whisker
WJ10ARD69869099.5113
WJ10AMA67849098111
WJ10AKN69839096107






Boxplot Means
Variabletrimmed meanunbiased SD
WJ10ARD91.5310.4
WJ10AMA90.3511.84
WJ10AKN89.569.4

\begin{tabular}{lllllllll}
\hline
Boxplot Means \tabularnewline
Variable & trimmed mean & unbiased SD \tabularnewline
WJ10ARD & 91.53 & 10.4 \tabularnewline
WJ10AMA & 90.35 & 11.84 \tabularnewline
WJ10AKN & 89.56 & 9.4 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=49811&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]91.53[/C][C]10.4[/C][/ROW]
[ROW][C]WJ10AMA[/C][C]90.35[/C][C]11.84[/C][/ROW]
[ROW][C]WJ10AKN[/C][C]89.56[/C][C]9.4[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=49811&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=49811&T=2

As an alternative you can also use a QR Code:  
QR Code


The GUIDs for individual cells are displayed in the table below:

Boxplot Means
Variabletrimmed meanunbiased SD
WJ10ARD91.5310.4
WJ10AMA90.3511.84
WJ10AKN89.569.4


Figures (Output of Computation)

Input Parameters & R Code

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