Repository of Reproducible Computations

Free Statistics

of Irreproducible Research!

Author's title

Author

*Unverified author*

R Software Module

rwasp_BoxplotAndScatterplotMatrix.wasp

Title produced by software

Boxplot and Trimmed Means

Date of computation

Mon, 21 Oct 2013 07:28:33 -0400

Cite this page as follows

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2013/Oct/21/t138235493688lzk0dwiucjeg6.htm/, Retrieved Mon, 29 Apr 2024 02:34:13 +0000

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=217368, Retrieved Mon, 29 Apr 2024 02:34:13 +0000

QR Codes:

Paste this QR Code to cite your computation.

Original text written by user:

IsPrivate?

No (this computation is public)

User-defined keywords

Estimated Impact

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 data] [2013-10-21 11:28:33] [d41d8cd98f00b204e9800998ecf8427e] [Current]

Feedback Forum

Post a new message

Dataseries X:

Download CSV

Histogram

Boxplots

88	95	95		
98	92	87	
89	95	91	
81	97	84	
88	92	92	
95	104	84	
79	79	84	
91	78	69	
80	79	84	
91	86	84	
85	92	85	
94	78	89	
86	86	83	
82	71	86	
107	90	91	
97	93	90	
78	80	77	
91	90	88	
74	77	88	
95	86	88	
97	89	102
112	88	103	
80	86	84	
98	97	87	
102	81	88	
100	88	93	
97	82	101	
92	82	88
114	100	113	
88	91	95	
96	111	99	
91	102	85	
114	105	105	
91	92	105	
114	100	96	
89	82	88	
103	91	93	
83	71	93	
74	71	69	
88	89	83	
88	82	83	
98	105	93	
74	71	69	
86	91	86	
93	105	90	
85	82	87	
91	86	92	
75	72	75	
97	98	89	
88	80	88	
79	84	87	
106	109	113	
102	80	102	
98	112	101	
98	86	87	
114	112	105	
95	90	101	
98	96	79	
82	84	86	
97	94	95	
94	112	104	
84	81	79	
88	105	86	
83	77	81	
98	105	92	
106	112	95	
94	91	97	
91	92	87	
89	72	94	
102	100	105	
103	103	94	
96	89	88	
89	86	105	
112	112	113	
74	71	69	
74	73	69	
102	99	100	
106	85	113	
91	80	88	
94	88	104	
98	95	94	
85	86	91	
108	92	103	
90	96	95	
74	71	69	
81	88	86	
104	98	108	
90	108	83	
95	103	102	
85	81	94	
114	96	113	
114	112	113	
93	95	93	
104	92	105	
92	103	88	
100	95	90	
92	94	92	
77	79	81	
86	91	86	
87	87	85	
91	82	86	
102	97	97	
100	92	103	
113	112	102

Summary of computational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	0 seconds
R Server	'Sir Ronald Aylmer Fisher' @ fisher.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.
R Engine error message	Error in list(88, 95, 95, 98, 92, 87, 89, 95, 91, 81, 97, 84, 88, 92, : argument 313 is empty Calls: array Execution halted

\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 & 0 seconds \tabularnewline
R Server & 'Sir Ronald Aylmer Fisher' @ fisher.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
R Engine error message & Error in list(88, 95, 95, 98, 92, 87, 89, 95, 91, 81, 97, 84, 88, 92,  : 
  argument 313 is empty
Calls: array
Execution halted
 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=217368&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]0 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Sir Ronald Aylmer Fisher' @ fisher.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]
[ROW][C]R Engine error message[/C][C]Error in list(88, 95, 95, 98, 92, 87, 89, 95, 91, 81, 97, 84, 88, 92,  : 
  argument 313 is empty
Calls: array
Execution halted
[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=217368&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=217368&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 Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	0 seconds
R Server	'Sir Ronald Aylmer Fisher' @ fisher.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.
R Engine error message	Error in list(88, 95, 95, 98, 92, 87, 89, 95, 91, 81, 97, 84, 88, 92, : argument 313 is empty Calls: array Execution halted

Parameters (Session):

par1 = 2 ; par2 = TRUE ; par3 = 10 ;

Parameters (R input):

par1 = 2 ; par2 = TRUE ; par3 = 10 ;

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

Free Statistics

Description of Statistical Computation

Tree of Dependent Computations

Dataset

Tables (Output of Computation)

Figures (Output of Computation)

Input Parameters & R Code