Home » date » 2009 » Oct » 24 »

workshop 3

*Unverified author*

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: Sat, 24 Oct 2009 10:53:15 -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/24/t12564032245odotx0600j76oe.htm/, Retrieved Sat, 24 Oct 2009 18:53:48 +0200

BibTeX entries for LaTeX users:

@Manual{KEY,
    author = {{YOUR NAME}},
    publisher = {Office for Research Development and Education},
    title = {Statistical Computations at FreeStatistics.org, URL http://www.freestatistics.org/blog/date/2009/Oct/24/t12564032245odotx0600j76oe.htm/},
    year = {2009},
}
@Manual{R,
    title = {R: A Language and Environment for Statistical Computing},
    author = {{R Development Core Team}},
    organization = {R Foundation for Statistical Computing},
    address = {Vienna, Austria},
    year = {2009},
    note = {{ISBN} 3-900051-07-0},
    url = {http://www.R-project.org},
}

Original text written by user:

IsPrivate?

No (this computation is public)

User-defined keywords:

Dataseries X:

» Textbox « » Textfile « » CSV «

95 81 84 89 81 83 87 90 97 90 91 102 91 91 88 82 96 75 84 87 87 82 85 92 92 89 92 71 74 94 84 85 91 80 78 88 84 81 80 62 71 77 69 80 73 77 72 63 84 75 73 80 87 73 84 84 77 87 72 85 85 91 86 96 90 80 89 87 82 80 87 77 83 90 82 85 89 67 86 72 73 89 110 98 91 96 99 89 95 89 90 88 90 73 72 63 77 73 78 88 70 97 88 86 100 86 79 77 88 67 73 83 92 96 88 89 95 91 93 84 102 94 93 98 97 94 103 88 91 80 75 70 84 77 81 84 88 100 87 89 98 90 105 88 88 104 88 92 82 97 93 89 90 86 89 78 101 107 85 117 114 108 88 96 79 88 81 81 113 111 107 86 73 96 95 84 82 78 82 84 99 83 98 104 110 104 85 82 84 89 93 87 105 105 99 99 114 98 105 100 89 84 91 89 96 111 105 89 97 93 88 90 94 95 97 89 93 87 92 58 68 63 93 90 81 77 72 84 69 67 NA 77 81 82 83 84 86 98 102 105 83 84 90 58 62 70 93 96 101 81 102 94 69 67 73 92 84 101 86 90 88 94 75 83 90 83 90 87 102 88 87 78 84 75 62 69 92 100 89 80 100 94 75 67 73 84 82 82 89 94 98 87 78 82 88 89 82 114 106 109 87 89 77 97 88 86 113 111 107 102 85 etc...

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	2 seconds
R Server	'Gwilym Jenkins' @ 72.249.127.135

Boxplot statistics
Variable	lower whisker	lower hinge	median	upper hinge	upper whisker
WJ10ARD	69	86	90	99.5	113
WJ10AMA	67	84	90	98	111
WJ10AKN	73	83	90	96	107
WJ7ARD	58	82	88.5	98	117
WJ7AMA	62	81	91	98	114
WJ7AKN	63	83	92	99	119

Boxplot Means
Variable	trimmed mean	unbiased SD
WJ10ARD	NA	NA
WJ10AMA	NA	NA
WJ10AKN	NA	NA
WJ7ARD	88.68	12.64
WJ7AMA	90.17	12.52
WJ7AKN	91.41	12.07

Charts produced by software:

http://www.freestatistics.org/blog/date/2009/Oct/24/t12564032245odotx0600j76oe/1vhta1256403193.png (open in new window)

http://www.freestatistics.org/blog/date/2009/Oct/24/t12564032245odotx0600j76oe/1vhta1256403193.ps (open in new window)

http://www.freestatistics.org/blog/date/2009/Oct/24/t12564032245odotx0600j76oe/2b2a51256403193.png (open in new window)

http://www.freestatistics.org/blog/date/2009/Oct/24/t12564032245odotx0600j76oe/2b2a51256403193.ps (open in new window)

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

This work is licensed under a Creative Commons Attribution-Noncommercial-Share Alike 3.0 License.

Software written by Ed van Stee & Patrick Wessa

Disclaimer

Information provided on this web site is provided "AS IS" without warranty of any kind, either express or implied, including, without limitation, warranties of merchantability, fitness for a particular purpose, and noninfringement. We use reasonable efforts to include accurate and timely information and periodically update the information, and software without notice. However, we make no warranties or representations as to the accuracy or completeness of such information (or software), and we assume no liability or responsibility for errors or omissions in the content of this web site, or any software bugs in online applications. Your use of this web site is AT YOUR OWN RISK. Under no circumstances and under no legal theory shall we be liable to you or any other person for any direct, indirect, special, incidental, exemplary, or consequential damages arising from your access to, or use of, this web site.

We may request personal information to be submitted to our servers in order to be able to:

personalize online software applications according to your needs
enforce strict security rules with respect to the data that you upload (e.g. statistical data)
manage user sessions of online applications
alert you about important changes or upgrades in resources or applications

We NEVER allow other companies to directly offer registered users information about their products and services. Banner references and hyperlinks of third parties NEVER contain any personal data of the visitor.

We do NOT sell, nor transmit by any means, personal information, nor statistical data series uploaded by you to third parties.

We carefully protect your data from loss, misuse, alteration, and destruction. However, at any time, and under any circumstance you are solely responsible for managing your passwords, and keeping them secret.

We store a unique ANONYMOUS USER ID in the form of a small 'Cookie' on your computer. This allows us to track your progress when using this website which is necessary to create state-dependent features. The cookie is used for NO OTHER PURPOSE. At any time you may opt to disallow cookies from this website - this will not affect other features of this website.

We examine cookies that are used by third-parties (banner and online ads) very closely: abuse from third-parties automatically results in termination of the advertising contract without refund. We have very good reason to believe that the cookies that are produced by third parties (banner ads) do NOT cause any privacy or security risk.

FreeStatistics.org is safe. There is no need to download any software to use the applications and services contained in this website. Hence, your system's security is not compromised by their use, and your personal data - other than data you submit in the account application form, and the user-agent information that is transmitted by your browser - is never transmitted to our servers.

As a general rule, we do not log on-line behavior of individuals (other than normal logging of webserver 'hits'). However, in cases of abuse, hacking, unauthorized access, Denial of Service attacks, illegal copying, hotlinking, non-compliance with international webstandards (such as robots.txt), or any other harmful behavior, our system engineers are empowered to log, track, identify, publish, and ban misbehaving individuals - even if this leads to ban entire blocks of IP addresses, or disclosing user's identity.

FreeStatistics.org is powered by