Free Statistics

of Irreproducible Research!

Author's title

Author*Unverified author*
R Software Modulerwasp_BoxplotAndScatterplotMatrix.wasp
Title produced by softwareBoxplot and Trimmed Means
Date of computationTue, 22 Oct 2013 08:03:32 -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/22/t1382443435ymnk59wkuw2tr4b.htm/, Retrieved Mon, 29 Apr 2024 07:53:26 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=218300, Retrieved Mon, 29 Apr 2024 07:53:26 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact78
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] [Compedium 3] [2013-10-22 12:03:32] [d41d8cd98f00b204e9800998ecf8427e] [Current]
Feedback Forum

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




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time3 seconds
R Server'George Udny Yule' @ yule.wessa.net

\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 & 3 seconds \tabularnewline
R Server & 'George Udny Yule' @ yule.wessa.net \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=218300&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]3 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'George Udny Yule' @ yule.wessa.net[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=218300&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=218300&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 time3 seconds
R Server'George Udny Yule' @ yule.wessa.net







Boxplot statistics
Variablelower whiskerlower hingemedianupper hingeupper whisker
WJ7ARD8085889298
WJ7AMA78869195102
WJ7AWL81879299108
WJ7AKN82889297101
WJ7ASK82859095103

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

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







Boxplot Notches
Variablelower boundmedianupper bound
WJ7ARD86.888889.12
WJ7AMA89.569192.44
WJ7AWL90.079293.93
WJ7AKN90.569293.44
WJ7ASK88.49091.6

\begin{tabular}{lllllllll}
\hline
Boxplot Notches \tabularnewline
Variable & lower bound & median & upper bound \tabularnewline
WJ7ARD & 86.88 & 88 & 89.12 \tabularnewline
WJ7AMA & 89.56 & 91 & 92.44 \tabularnewline
WJ7AWL & 90.07 & 92 & 93.93 \tabularnewline
WJ7AKN & 90.56 & 92 & 93.44 \tabularnewline
WJ7ASK & 88.4 & 90 & 91.6 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=218300&T=2

[TABLE]
[ROW][C]Boxplot Notches[/C][/ROW]
[ROW][C]Variable[/C][C]lower bound[/C][C]median[/C][C]upper bound[/C][/ROW]
[ROW][C]WJ7ARD[/C][C]86.88[/C][C]88[/C][C]89.12[/C][/ROW]
[ROW][C]WJ7AMA[/C][C]89.56[/C][C]91[/C][C]92.44[/C][/ROW]
[ROW][C]WJ7AWL[/C][C]90.07[/C][C]92[/C][C]93.93[/C][/ROW]
[ROW][C]WJ7AKN[/C][C]90.56[/C][C]92[/C][C]93.44[/C][/ROW]
[ROW][C]WJ7ASK[/C][C]88.4[/C][C]90[/C][C]91.6[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=218300&T=2

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

As an alternative you can also use a QR Code:  

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

Boxplot Notches
Variablelower boundmedianupper bound
WJ7ARD86.888889.12
WJ7AMA89.569192.44
WJ7AWL90.079293.93
WJ7AKN90.569293.44
WJ7ASK88.49091.6







Boxplot Means
Variabletrimmed meanunbiased SD
WJ7ARD88.195.14
WJ7AMA90.86.27
WJ7AWL92.327.47
WJ7AKN91.985.39
WJ7ASK90.145.83

\begin{tabular}{lllllllll}
\hline
Boxplot Means \tabularnewline
Variable & trimmed mean & unbiased SD \tabularnewline
WJ7ARD & 88.19 & 5.14 \tabularnewline
WJ7AMA & 90.8 & 6.27 \tabularnewline
WJ7AWL & 92.32 & 7.47 \tabularnewline
WJ7AKN & 91.98 & 5.39 \tabularnewline
WJ7ASK & 90.14 & 5.83 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=218300&T=3

[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]88.19[/C][C]5.14[/C][/ROW]
[ROW][C]WJ7AMA[/C][C]90.8[/C][C]6.27[/C][/ROW]
[ROW][C]WJ7AWL[/C][C]92.32[/C][C]7.47[/C][/ROW]
[ROW][C]WJ7AKN[/C][C]91.98[/C][C]5.39[/C][/ROW]
[ROW][C]WJ7ASK[/C][C]90.14[/C][C]5.83[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=218300&T=3

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

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
WJ7ARD88.195.14
WJ7AMA90.86.27
WJ7AWL92.327.47
WJ7AKN91.985.39
WJ7ASK90.145.83



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