FreeStatistics.org

Free Statistics

of Irreproducible Research!

Description of Statistical Computation

Author's title

Author*Unverified author*
R Software Modulerwasp_BoxplotAndScatterplotMatrix.wasp
Title produced by softwareBoxplot and Trimmed Means
Date of computationSun, 20 Oct 2013 12:13:47 -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/2013/Oct/20/t1382285756fadm3rek23e23n8.htm/, Retrieved Sat, 27 Apr 2024 22:40:07 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=217095, Retrieved Sat, 27 Apr 2024 22:40:07 +0000
QR Codes:

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

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 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]
- R PD            [Boxplot and Trimmed Means] [boxplot computation] [2013-10-20 16:13:47] [d41d8cd98f00b204e9800998ecf8427e] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
3	6	10	13	25
5	9	29	23	35
10	11	32	33	45
15	15	41	53	60
56	30	44	64	63
59	38	5	67	63
62	38	69	68	66
63	38	70	69	67
68	45	70	69	67
68	45	71	69	69
72	72	72	70	70
72	72	72	70	70
72	72	73	72	71
72	72	74	72	71
73	72	74	73	73
73	72	75	73	73
76	72	77	75	74
76	72	77	77	75
76	73	77	77	75
76	74	78	77	75
76	74	78	78	77
77	75	78	79	78
77	75	78	79	78
77	75	78	80	79
77	75	78	80	79
78	75	78	81	79
79	75	79	81	79
79	75	79	81	79
79	77	79	82	80
79	77	79	82	81
80	78	80	82	81
80	78	81	82	81
80	78	81	82	82
80	78	81	83	82
80	79	82	83	82
81	79	82	83	82
81	79	82	83	82
81	80	82	83	82
81	80	83	83	82
81	81	83	83	82
82	81	83	83	82
82	81	83	84	83
82	81	84	84	83
82	81	84	84	83
82	81	84	84	84
82	81	84	85	84
83	81	84	86	84
83	81	84	86	84
84	81	84	86	85
84	81	85	86	85
84	81	86	87	85
85	81	86	87	85
85	81	86	87	85
85	81	86	88	85
85	85	86	88	85
85	85	86	88	85
85	86	87	88	85
85	87	87	88	85
86	87	87	88	85
86	87	87	88	86
86	87	87	88	86
86	88	87	88	86
86	88	87	89	86
86	88	88	89	87
86	88	88	89	87
86	88	88	89	87
86	88	88	89	87
86	89	88	89	87
86	89	89	89	87
86	89	89	89	87
86	90	89	90	88
87	90	90	90	88
88	90	90	90	88
88	90	90	91	88
88	91	90	91	88
88	91	90	91	89
88	91	90	91	89
88	91	91	91	89
88	91	91	91	89
88	91	92	92	90
88	91	92	92	90
89	91	92	92	90
89	91	92	92	91
89	91	93	92	91
89	92	94	92	92
89	92	95	93	92
89	92	95	93	92
89	92	95	93	92
89	92	95	93	92
90	92	95	94	92
90	92	95	94	92
90	93	95	94	92
90	93	96	94	93
90	93	96	94	93
90	93	96	95	93
90	93	96	95	93
90	93	96	95	93
90	95	97	95	94
91	95	97	96	94
91	95	97	96	94
91	95	98	96	95
91	95	98	96	95
91	95	98	96	95
92	95	98	97	95
92	95	99	97	95
92	95	99	97	95
92	95	99	97	95
92	95	99	97	95
92	95	100	97	97
92	96	100	97	97
93	97	100	98	97
93	97	101	98	97
94	97	102	98	97
95	97	102	98	97
95	97	102	98	97
95	97	102	98	97
96	107	103	98	97
96	108	103	98	98
110	109	103	99	99
140	124	104	99	99
150	150	104	99	99

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'Sir Maurice George Kendall' @ kendall.wessa.net
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 & 'Sir Maurice George Kendall' @ kendall.wessa.net \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=217095&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]
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=217095&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=217095&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'Sir Maurice George Kendall' @ kendall.wessa.net
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
WJ7ARD6880868995
WJ7AMA7278879297
WJ7AWL69818795102
WJ7AKN6882.5889398
WJ7ASK6782869297

\begin{tabular}{lllllllll}
\hline
Boxplot statistics \tabularnewline
Variable & lower whisker & lower hinge & median & upper hinge & upper whisker \tabularnewline
WJ7ARD & 68 & 80 & 86 & 89 & 95 \tabularnewline
WJ7AMA & 72 & 78 & 87 & 92 & 97 \tabularnewline
WJ7AWL & 69 & 81 & 87 & 95 & 102 \tabularnewline
WJ7AKN & 68 & 82.5 & 88 & 93 & 98 \tabularnewline
WJ7ASK & 67 & 82 & 86 & 92 & 97 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=217095&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]WJ7ARD[/C][C]68[/C][C]80[/C][C]86[/C][C]89[/C][C]95[/C][/ROW]
[ROW][C]WJ7AMA[/C][C]72[/C][C]78[/C][C]87[/C][C]92[/C][C]97[/C][/ROW]
[ROW][C]WJ7AWL[/C][C]69[/C][C]81[/C][C]87[/C][C]95[/C][C]102[/C][/ROW]
[ROW][C]WJ7AKN[/C][C]68[/C][C]82.5[/C][C]88[/C][C]93[/C][C]98[/C][/ROW]
[ROW][C]WJ7ASK[/C][C]67[/C][C]82[/C][C]86[/C][C]92[/C][C]97[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=217095&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=217095&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
WJ7ARD6880868995
WJ7AMA7278879297
WJ7AWL69818795102
WJ7AKN6882.5889398
WJ7ASK6782869297






Boxplot Means
Variabletrimmed meanunbiased SD
WJ7ARD84.496.75
WJ7AMA85.0911.38
WJ7AWL87.248.35
WJ7AKN87.367.77
WJ7ASK85.797.62

\begin{tabular}{lllllllll}
\hline
Boxplot Means \tabularnewline
Variable & trimmed mean & unbiased SD \tabularnewline
WJ7ARD & 84.49 & 6.75 \tabularnewline
WJ7AMA & 85.09 & 11.38 \tabularnewline
WJ7AWL & 87.24 & 8.35 \tabularnewline
WJ7AKN & 87.36 & 7.77 \tabularnewline
WJ7ASK & 85.79 & 7.62 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=217095&T=2

[TABLE]
[ROW][C]Boxplot Means[/C][/ROW]
[ROW][C]Variable[/C][C]trimmed mean[/C][C]unbiased SD[/C][/ROW]
[ROW][C]WJ7ARD[/C][C]84.49[/C][C]6.75[/C][/ROW]
[ROW][C]WJ7AMA[/C][C]85.09[/C][C]11.38[/C][/ROW]
[ROW][C]WJ7AWL[/C][C]87.24[/C][C]8.35[/C][/ROW]
[ROW][C]WJ7AKN[/C][C]87.36[/C][C]7.77[/C][/ROW]
[ROW][C]WJ7ASK[/C][C]85.79[/C][C]7.62[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=217095&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=217095&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
WJ7ARD84.496.75
WJ7AMA85.0911.38
WJ7AWL87.248.35
WJ7AKN87.367.77
WJ7ASK85.797.62


Figures (Output of Computation)

Input Parameters & R Code

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