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 07:37:27 -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/t1382441928is7zb68f25wx412.htm/, Retrieved Mon, 29 Apr 2024 04:56:34 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=218229, Retrieved Mon, 29 Apr 2024 04:56:34 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact64
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  D            [Boxplot and Trimmed Means] [predicting outcom...] [2013-10-22 11:37:27] [d41d8cd98f00b204e9800998ecf8427e] [Current]
Feedback Forum

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





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
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 & 3 seconds \tabularnewline
R Server & 'George Udny Yule' @ yule.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=218229&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]
[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=218229&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=218229&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
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
WJ7ASK67849097114
WJ10AFS788894100114
WJ10AVA7385.591.598.5112
WJ10ARD7986.591.5101113
WJ10AMA72849099111
WJ10AKN 75849197107

\begin{tabular}{lllllllll}
\hline
Boxplot statistics \tabularnewline
Variable & lower whisker & lower hinge & median & upper hinge & upper whisker \tabularnewline
WJ7ASK & 67 & 84 & 90 & 97 & 114 \tabularnewline
WJ10AFS & 78 & 88 & 94 & 100 & 114 \tabularnewline
WJ10AVA & 73 & 85.5 & 91.5 & 98.5 & 112 \tabularnewline
WJ10ARD & 79 & 86.5 & 91.5 & 101 & 113 \tabularnewline
WJ10AMA & 72 & 84 & 90 & 99 & 111 \tabularnewline
WJ10AKN
 & 75 & 84 & 91 & 97 & 107 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=218229&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]WJ7ASK[/C][C]67[/C][C]84[/C][C]90[/C][C]97[/C][C]114[/C][/ROW]
[ROW][C]WJ10AFS[/C][C]78[/C][C]88[/C][C]94[/C][C]100[/C][C]114[/C][/ROW]
[ROW][C]WJ10AVA[/C][C]73[/C][C]85.5[/C][C]91.5[/C][C]98.5[/C][C]112[/C][/ROW]
[ROW][C]WJ10ARD[/C][C]79[/C][C]86.5[/C][C]91.5[/C][C]101[/C][C]113[/C][/ROW]
[ROW][C]WJ10AMA[/C][C]72[/C][C]84[/C][C]90[/C][C]99[/C][C]111[/C][/ROW]
[ROW][C]WJ10AKN
[/C][C]75[/C][C]84[/C][C]91[/C][C]97[/C][C]107[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=218229&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=218229&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
WJ7ASK67849097114
WJ10AFS788894100114
WJ10AVA7385.591.598.5112
WJ10ARD7986.591.5101113
WJ10AMA72849099111
WJ10AKN 75849197107







Boxplot Notches
Variablelower boundmedianupper bound
WJ7ASK88.299091.71
WJ10AFS92.069495.94
WJ10AVA89.491.593.6
WJ10ARD89.1691.593.84
WJ10AMA87.589092.42
WJ10AKN 88.899193.11

\begin{tabular}{lllllllll}
\hline
Boxplot Notches \tabularnewline
Variable & lower bound & median & upper bound \tabularnewline
WJ7ASK & 88.29 & 90 & 91.71 \tabularnewline
WJ10AFS & 92.06 & 94 & 95.94 \tabularnewline
WJ10AVA & 89.4 & 91.5 & 93.6 \tabularnewline
WJ10ARD & 89.16 & 91.5 & 93.84 \tabularnewline
WJ10AMA & 87.58 & 90 & 92.42 \tabularnewline
WJ10AKN
 & 88.89 & 91 & 93.11 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=218229&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]WJ7ASK[/C][C]88.29[/C][C]90[/C][C]91.71[/C][/ROW]
[ROW][C]WJ10AFS[/C][C]92.06[/C][C]94[/C][C]95.94[/C][/ROW]
[ROW][C]WJ10AVA[/C][C]89.4[/C][C]91.5[/C][C]93.6[/C][/ROW]
[ROW][C]WJ10ARD[/C][C]89.16[/C][C]91.5[/C][C]93.84[/C][/ROW]
[ROW][C]WJ10AMA[/C][C]87.58[/C][C]90[/C][C]92.42[/C][/ROW]
[ROW][C]WJ10AKN
[/C][C]88.89[/C][C]91[/C][C]93.11[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=218229&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=218229&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
WJ7ASK88.299091.71
WJ10AFS92.069495.94
WJ10AVA89.491.593.6
WJ10ARD89.1691.593.84
WJ10AMA87.589092.42
WJ10AKN 88.899193.11







Boxplot Means
Variabletrimmed meanunbiased SD
WJ7ASK90.9610.96
WJ10AFS94.479.32
WJ10AVA92.189.99
WJ10ARD92.978.79
WJ10AMA92.2510.28
WJ10AKN 91.18.16

\begin{tabular}{lllllllll}
\hline
Boxplot Means \tabularnewline
Variable & trimmed mean & unbiased SD \tabularnewline
WJ7ASK & 90.96 & 10.96 \tabularnewline
WJ10AFS & 94.47 & 9.32 \tabularnewline
WJ10AVA & 92.18 & 9.99 \tabularnewline
WJ10ARD & 92.97 & 8.79 \tabularnewline
WJ10AMA & 92.25 & 10.28 \tabularnewline
WJ10AKN
 & 91.1 & 8.16 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=218229&T=3

[TABLE]
[ROW][C]Boxplot Means[/C][/ROW]
[ROW][C]Variable[/C][C]trimmed mean[/C][C]unbiased SD[/C][/ROW]
[ROW][C]WJ7ASK[/C][C]90.96[/C][C]10.96[/C][/ROW]
[ROW][C]WJ10AFS[/C][C]94.47[/C][C]9.32[/C][/ROW]
[ROW][C]WJ10AVA[/C][C]92.18[/C][C]9.99[/C][/ROW]
[ROW][C]WJ10ARD[/C][C]92.97[/C][C]8.79[/C][/ROW]
[ROW][C]WJ10AMA[/C][C]92.25[/C][C]10.28[/C][/ROW]
[ROW][C]WJ10AKN
[/C][C]91.1[/C][C]8.16[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=218229&T=3

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=218229&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
WJ7ASK90.9610.96
WJ10AFS94.479.32
WJ10AVA92.189.99
WJ10ARD92.978.79
WJ10AMA92.2510.28
WJ10AKN 91.18.16



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