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 computationMon, 21 Oct 2013 18:18:49 -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/21/t1382393955nu48upqpr63sf5n.htm/, Retrieved Mon, 29 Apr 2024 04:38:06 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=217873, Retrieved Mon, 29 Apr 2024 04:38:06 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact84
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] [Year 10 and 7 data] [2013-10-21 22:18:49] [d41d8cd98f00b204e9800998ecf8427e] [Current]
Feedback Forum

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




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time5 seconds
R Server'Gertrude Mary Cox' @ cox.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 & 5 seconds \tabularnewline
R Server & 'Gertrude Mary Cox' @ cox.wessa.net \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=217873&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]5 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Gertrude Mary Cox' @ cox.wessa.net[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=217873&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=217873&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 time5 seconds
R Server'Gertrude Mary Cox' @ cox.wessa.net







Boxplot statistics
Variablelower whiskerlower hingemedianupper hingeupper whisker
WJ7ARD8085889298
WJ7AMA78869195102
WJ7AKN82889297101
WJ10ARD8789.594102.5113
WJ10AMA8489.594.5104.5111
WJ10AKN84909398.5107

\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
WJ7AKN & 82 & 88 & 92 & 97 & 101 \tabularnewline
WJ10ARD & 87 & 89.5 & 94 & 102.5 & 113 \tabularnewline
WJ10AMA & 84 & 89.5 & 94.5 & 104.5 & 111 \tabularnewline
WJ10AKN & 84 & 90 & 93 & 98.5 & 107 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=217873&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]WJ7AKN[/C][C]82[/C][C]88[/C][C]92[/C][C]97[/C][C]101[/C][/ROW]
[ROW][C]WJ10ARD[/C][C]87[/C][C]89.5[/C][C]94[/C][C]102.5[/C][C]113[/C][/ROW]
[ROW][C]WJ10AMA[/C][C]84[/C][C]89.5[/C][C]94.5[/C][C]104.5[/C][C]111[/C][/ROW]
[ROW][C]WJ10AKN[/C][C]84[/C][C]90[/C][C]93[/C][C]98.5[/C][C]107[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=217873&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=217873&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
WJ7AKN82889297101
WJ10ARD8789.594102.5113
WJ10AMA8489.594.5104.5111
WJ10AKN84909398.5107







Boxplot Notches
Variablelower boundmedianupper bound
WJ7ARD86.888889.12
WJ7AMA89.569192.44
WJ7AKN90.569293.44
WJ10ARD91.589496.42
WJ10AMA91.7194.597.29
WJ10AKN91.419394.59

\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
WJ7AKN & 90.56 & 92 & 93.44 \tabularnewline
WJ10ARD & 91.58 & 94 & 96.42 \tabularnewline
WJ10AMA & 91.71 & 94.5 & 97.29 \tabularnewline
WJ10AKN & 91.41 & 93 & 94.59 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=217873&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]WJ7AKN[/C][C]90.56[/C][C]92[/C][C]93.44[/C][/ROW]
[ROW][C]WJ10ARD[/C][C]91.58[/C][C]94[/C][C]96.42[/C][/ROW]
[ROW][C]WJ10AMA[/C][C]91.71[/C][C]94.5[/C][C]97.29[/C][/ROW]
[ROW][C]WJ10AKN[/C][C]91.41[/C][C]93[/C][C]94.59[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=217873&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=217873&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
WJ7AKN90.569293.44
WJ10ARD91.589496.42
WJ10AMA91.7194.597.29
WJ10AKN91.419394.59







Boxplot Means
Variabletrimmed meanunbiased SD
WJ7ARD88.195.14
WJ7AMA90.86.27
WJ7AKN91.985.39
WJ10ARD95.327.88
WJ10AMA95.348.07
WJ10AKN93.796.41

\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
WJ7AKN & 91.98 & 5.39 \tabularnewline
WJ10ARD & 95.32 & 7.88 \tabularnewline
WJ10AMA & 95.34 & 8.07 \tabularnewline
WJ10AKN & 93.79 & 6.41 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=217873&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]WJ7AKN[/C][C]91.98[/C][C]5.39[/C][/ROW]
[ROW][C]WJ10ARD[/C][C]95.32[/C][C]7.88[/C][/ROW]
[ROW][C]WJ10AMA[/C][C]95.34[/C][C]8.07[/C][/ROW]
[ROW][C]WJ10AKN[/C][C]93.79[/C][C]6.41[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=217873&T=3

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=217873&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
WJ7AKN91.985.39
WJ10ARD95.327.88
WJ10AMA95.348.07
WJ10AKN93.796.41



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