Home » date » 2009 » Oct » 26 »

Week Three, exercise 4

*The author of this computation has been verified*

R Software Module: Ian.Holliday/rwasp_CARE Data Boxplot.wasp (opens new window with default values)

Title produced by software: Boxplot and Trimmed Means

Date of computation: Mon, 26 Oct 2009 09:13:50 -0600

Cite this page as follows:

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL http://www.freestatistics.org/blog/date/2009/Oct/26/t1256570136bjimfh5w7kxj00t.htm/, Retrieved Wed, 22 May 2013 22:46:58 +0000

Original text written by user:

IsPrivate?

No (this computation is public)

User-defined keywords:

System-generated keywords (parent):

t1255827101xzhmyn1xdafjo0g (pk = 47220)

Estimated Impact

Dataseries X:

» Textfile « » CSV « » Correlation Matrix « » Notched Boxplots «

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

Output produced by software:

Enter (or paste) a matrix (table) containing all data (time) series. Every column represents a different variable and must be delimited by a space or Tab. Every row represents a period in time (or category) and must be delimited by hard returns. The easiest way to enter data is to copy and paste a block of spreadsheet cells. Please, do not use commas or spaces to seperate groups of digits!

Summary of computational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	1 seconds
R Server	'Gwilym Jenkins' @ 72.249.127.135
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
Variable	lower whisker	lower hinge	median	upper hinge	upper whisker
WJ10AFS	74	86	92	99	114
WJ10AVA	71	82	91	97.5	112
WJ10ARD	69	86	90	99.5	113

Boxplot Means
Variable	trimmed mean	unbiased SD
WJ10AFS	93.09	10.45
WJ10AVA	90.72	11.14
WJ10ARD	91.53	10.4

Charts produced by software:

http://www.freestatistics.org/blog/date/2009/Oct/26/t1256570136bjimfh5w7kxj00t/1xyj21256570028.png (opens in new window)

http://www.freestatistics.org/blog/date/2009/Oct/26/t1256570136bjimfh5w7kxj00t/1xyj21256570028.ps (opens in new window)

Click here to open pdf file.

http://www.freestatistics.org/blog/date/2009/Oct/26/t1256570136bjimfh5w7kxj00t/2g8gx1256570028.png (opens in new window)

http://www.freestatistics.org/blog/date/2009/Oct/26/t1256570136bjimfh5w7kxj00t/2g8gx1256570028.ps (opens in new window)

Click here to open pdf file.

Parameters (Session):

par1 = 3 ; par2 = FALSE ; par3 = 0 ;

Parameters (R input):

par1 = 3 ; par2 = FALSE ; par3 = 0 ;

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('http://www.xycoon.com/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('http://www.xycoon.com/lower_whisker.htm','lower whisker','definition of lower whisker'),1,TRUE)
a<-table.element(a,hyperlink('http://www.xycoon.com/lower_hinge.htm','lower hinge','definition of lower hinge'),1,TRUE)
a<-table.element(a,hyperlink('http://www.xycoon.com/central_tendency.htm','median','definitions about measures of central tendency'),1,TRUE)
a<-table.element(a,hyperlink('http://www.xycoon.com/upper_hinge.htm','upper hinge','definition of upper hinge'),1,TRUE)
a<-table.element(a,hyperlink('http://www.xycoon.com/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('http://www.xycoon.com/trimmed_mean.htm','trimmed mean','definition of trimmed mean'),1,TRUE)
a<-table.element(a,hyperlink('http://www.xycoon.com/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(z[i], trim=par3/100), digits=2))
a<-table.element(a,round(sd(z[i]), digits=2))
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable2.tab')