Home » date » 2009 » Oct » 21 »

WS3P2

*The author of this computation has been verified*

R Software Module: rwasp_centraltendency.wasp (opens new window with default values)

Title produced by software: Central Tendency

Date of computation: Wed, 21 Oct 2009 09:47:30 -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/21/t1256140202upg6wb2gfkdcveb.htm/, Retrieved Wed, 21 Oct 2009 17:50:04 +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/21/t1256140202upg6wb2gfkdcveb.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 «

10689 10773 6718 6519 9354 9224 13355 12006 8888 5602 1913 4426 5848 3683 2094 2893 1957 4288 9278 8979 2578 4292 3510 5027 8047 6601 4550 7264 2227 2081 4189 6650 2579 2388 3056 1938 3700 4365 4161 4264 4606 6012 5586 7659 10776 5846 9714 10095 10433 10942 3646 7629 8813 8516 9939 10231 5604 10880 13297 7982 13076 14680 2716 7200 3236 4846 5788 8372

Output produced by software:

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

Central Tendency - Ungrouped Data
Measure	Value	S.E.	Value/S.E.
Arithmetic Mean	6618.73529411765	405.03506748194	16.3411413615755
Geometric Mean	5740.53269949587
Harmonic Mean	4884.26038411593
Quadratic Mean	7402.650544787
Winsorized Mean ( 1 / 22 )	6599.61764705882	399.616002045776	16.5148983355847
Winsorized Mean ( 2 / 22 )	6598.4705882353	399.089972533424	16.5337919826654
Winsorized Mean ( 3 / 22 )	6594.19117647059	395.736387667169	16.6630903348129
Winsorized Mean ( 4 / 22 )	6532.01470588235	381.147747117645	17.1377497447628
Winsorized Mean ( 5 / 22 )	6463.55882352941	363.848999050293	17.7643990787398
Winsorized Mean ( 6 / 22 )	6472.29411764706	360.411265827029	17.9580793702306
Winsorized Mean ( 7 / 22 )	6481.14705882353	355.209781877783	18.2459700984628
Winsorized Mean ( 8 / 22 )	6480.91176470588	355.12680823573	18.2495706164878
Winsorized Mean ( 9 / 22 )	6487.92647058824	350.18416534149	18.5271840154777
Winsorized Mean ( 10 / 22 )	6476.30882352941	339.393034412084	19.0820322366016
Winsorized Mean ( 11 / 22 )	6470	329.651146520237	19.6268087288536
Winsorized Mean ( 12 / 22 )	6477.76470588235	320.735763229788	20.1965775211710
Winsorized Mean ( 13 / 22 )	6500.32352941176	307.990651147963	21.1055871507246
Winsorized Mean ( 14 / 22 )	6482	296.329682867242	21.8742852126089
Winsorized Mean ( 15 / 22 )	6410.75	282.51425284735	22.6917754958859
Winsorized Mean ( 16 / 22 )	6396.86764705882	279.170043234744	22.9138756183804
Winsorized Mean ( 17 / 22 )	6498.61764705882	261.048225514976	24.8943184127716
Winsorized Mean ( 18 / 22 )	6441.17647058824	250.087773828715	25.7556631896759
Winsorized Mean ( 19 / 22 )	6436.70588235294	243.44719530172	26.4398440671108
Winsorized Mean ( 20 / 22 )	6421.70588235294	239.203840267629	26.8461654928622
Winsorized Mean ( 21 / 22 )	6331.22058823529	225.559438004022	28.0689677375522
Winsorized Mean ( 22 / 22 )	6308.25	215.682704771155	29.2478249783320
Trimmed Mean ( 1 / 22 )	6567.89393939394	392.270536092904	16.7432762215881
Trimmed Mean ( 2 / 22 )	6534.1875	383.35770683059	17.0446227728703
Trimmed Mean ( 3 / 22 )	6498.93548387097	372.860375736376	17.4299440401410
Trimmed Mean ( 4 / 22 )	6462.95	361.626266658394	17.8719042168061
Trimmed Mean ( 5 / 22 )	6442.70689655172	353.489859069515	18.2260020514047
Trimmed Mean ( 6 / 22 )	6437.64285714286	348.970767231568	18.4475132636855
Trimmed Mean ( 7 / 22 )	6430.37037037037	344.136893495767	18.6855012987716
Trimmed Mean ( 8 / 22 )	6420.88461538462	339.248405869878	18.9267937720167
Trimmed Mean ( 9 / 22 )	6410.68	332.897674344044	19.2572087282733
Trimmed Mean ( 10 / 22 )	6398.52083333333	325.857659921894	19.6359380806547
Trimmed Mean ( 11 / 22 )	6387.02173913043	319.251503068850	20.0062385853607
Trimmed Mean ( 12 / 22 )	6375.36363636364	312.782385041232	20.3827451329243
Trimmed Mean ( 13 / 22 )	6361.54761904762	306.151014868813	20.7791165473469
Trimmed Mean ( 14 / 22 )	6343.4	300.061657766430	21.1403217832574
Trimmed Mean ( 15 / 22 )	6325.68421052632	294.338325632873	21.4912013137437
Trimmed Mean ( 16 / 22 )	6314.97222222222	289.567648389504	21.8082795413935
Trimmed Mean ( 17 / 22 )	6304.73529411765	283.266370325914	22.2572672035290
Trimmed Mean ( 18 / 22 )	6280.5	278.411019114754	22.5583743774571
Trimmed Mean ( 19 / 22 )	6260.26666666667	273.806790396865	22.8638108557966
Trimmed Mean ( 20 / 22 )	6237.71428571429	268.051300355121	23.2705988646592
Trimmed Mean ( 21 / 22 )	6213.65384615385	259.796369621746	23.9174005980172
Trimmed Mean ( 22 / 22 )	6197.79166666667	251.342806633204	24.6587191003695
Median	5930
Midrange	8296.5
Midmean - Weighted Average at Xnp	6230.31428571429
Midmean - Weighted Average at X(n+1)p	6304.73529411765
Midmean - Empirical Distribution Function	6230.31428571429
Midmean - Empirical Distribution Function - Averaging	6304.73529411765
Midmean - Empirical Distribution Function - Interpolation	6304.73529411765
Midmean - Closest Observation	6230.31428571429
Midmean - True Basic - Statistics Graphics Toolkit	6304.73529411765
Midmean - MS Excel (old versions)	6314.97222222222
Number of observations	68

Charts produced by software:

http://www.freestatistics.org/blog/date/2009/Oct/21/t1256140202upg6wb2gfkdcveb/1ny1t1256140048.png (open in new window)

http://www.freestatistics.org/blog/date/2009/Oct/21/t1256140202upg6wb2gfkdcveb/1ny1t1256140048.ps (open in new window)

http://www.freestatistics.org/blog/date/2009/Oct/21/t1256140202upg6wb2gfkdcveb/2ysl61256140048.png (open in new window)

http://www.freestatistics.org/blog/date/2009/Oct/21/t1256140202upg6wb2gfkdcveb/2ysl61256140048.ps (open in new window)

Parameters (Session):

Parameters (R input):

R code (references can be found in the software module):

geomean <- function(x) {

return(exp(mean(log(x))))

}

harmean <- function(x) {

return(1/mean(1/x))

}

quamean <- function(x) {

return(sqrt(mean(x*x)))

}

winmean <- function(x) {

x <-sort(x[!is.na(x)])

n<-length(x)

denom <- 3

nodenom <- n/denom

if (nodenom>40) denom <- n/40

sqrtn = sqrt(n)

roundnodenom = floor(nodenom)

win <- array(NA,dim=c(roundnodenom,2))

for (j in 1:roundnodenom) {

win[j,1] <- (j*x[j+1]+sum(x[(j+1):(n-j)])+j*x[n-j])/n

win[j,2] <- sd(c(rep(x[j+1],j),x[(j+1):(n-j)],rep(x[n-j],j)))/sqrtn

}

return(win)

}

trimean <- function(x) {

x <-sort(x[!is.na(x)])

n<-length(x)

denom <- 3

nodenom <- n/denom

if (nodenom>40) denom <- n/40

sqrtn = sqrt(n)

roundnodenom = floor(nodenom)

tri <- array(NA,dim=c(roundnodenom,2))

for (j in 1:roundnodenom) {

tri[j,1] <- mean(x,trim=j/n)

tri[j,2] <- sd(x[(j+1):(n-j)]) / sqrt(n-j*2)

}

return(tri)

}

midrange <- function(x) {

return((max(x)+min(x))/2)

}

q1 <- function(data,n,p,i,f) {

np <- n*p;

i <<- floor(np)

f <<- np - i

qvalue <- (1-f)*data[i] + f*data[i+1]

}

q2 <- function(data,n,p,i,f) {

np <- (n+1)*p

i <<- floor(np)

f <<- np - i

qvalue <- (1-f)*data[i] + f*data[i+1]

}

q3 <- function(data,n,p,i,f) {

np <- n*p

i <<- floor(np)

f <<- np - i

if (f==0) {

qvalue <- data[i]

} else {

qvalue <- data[i+1]

}

}

q4 <- function(data,n,p,i,f) {

np <- n*p

i <<- floor(np)

f <<- np - i

if (f==0) {

qvalue <- (data[i]+data[i+1])/2

} else {

qvalue <- data[i+1]

}

}

q5 <- function(data,n,p,i,f) {

np <- (n-1)*p

i <<- floor(np)

f <<- np - i

if (f==0) {

qvalue <- data[i+1]

} else {

qvalue <- data[i+1] + f*(data[i+2]-data[i+1])

}

}

q6 <- function(data,n,p,i,f) {

np <- n*p+0.5

i <<- floor(np)

f <<- np - i

qvalue <- data[i]

}

q7 <- function(data,n,p,i,f) {

np <- (n+1)*p

i <<- floor(np)

f <<- np - i

if (f==0) {

qvalue <- data[i]

} else {

qvalue <- f*data[i] + (1-f)*data[i+1]

}

}

q8 <- function(data,n,p,i,f) {

np <- (n+1)*p

i <<- floor(np)

f <<- np - i

if (f==0) {

qvalue <- data[i]

} else {

if (f == 0.5) {

qvalue <- (data[i]+data[i+1])/2

} else {

if (f < 0.5) {

qvalue <- data[i]

} else {

qvalue <- data[i+1]

}

}

}

}

midmean <- function(x,def) {

x <-sort(x[!is.na(x)])

n<-length(x)

if (def==1) {

qvalue1 <- q1(x,n,0.25,i,f)

qvalue3 <- q1(x,n,0.75,i,f)

}

if (def==2) {

qvalue1 <- q2(x,n,0.25,i,f)

qvalue3 <- q2(x,n,0.75,i,f)

}

if (def==3) {

qvalue1 <- q3(x,n,0.25,i,f)

qvalue3 <- q3(x,n,0.75,i,f)

}

if (def==4) {

qvalue1 <- q4(x,n,0.25,i,f)

qvalue3 <- q4(x,n,0.75,i,f)

}

if (def==5) {

qvalue1 <- q5(x,n,0.25,i,f)

qvalue3 <- q5(x,n,0.75,i,f)

}

if (def==6) {

qvalue1 <- q6(x,n,0.25,i,f)

qvalue3 <- q6(x,n,0.75,i,f)

}

if (def==7) {

qvalue1 <- q7(x,n,0.25,i,f)

qvalue3 <- q7(x,n,0.75,i,f)

}

if (def==8) {

qvalue1 <- q8(x,n,0.25,i,f)

qvalue3 <- q8(x,n,0.75,i,f)

}

midm <- 0

myn <- 0

roundno4 <- round(n/4)

round3no4 <- round(3*n/4)

for (i in 1:n) {

if ((x[i]>=qvalue1) & (x[i]<=qvalue3)){

midm = midm + x[i]

myn = myn + 1

}

}

midm = midm / myn

return(midm)

}

(arm <- mean(x))

sqrtn <- sqrt(length(x))

(armse <- sd(x) / sqrtn)

(armose <- arm / armse)

(geo <- geomean(x))

(har <- harmean(x))

(qua <- quamean(x))

(win <- winmean(x))

(tri <- trimean(x))

(midr <- midrange(x))

midm <- array(NA,dim=8)

for (j in 1:8) midm[j] <- midmean(x,j)

midm

bitmap(file='test1.png')

lb <- win[,1] - 2*win[,2]

ub <- win[,1] + 2*win[,2]

if ((ylimmin == '') | (ylimmax == '')) plot(win[,1],type='b',main=main, xlab='j', pch=19, ylab='Winsorized Mean(j/n)', ylim=c(min(lb),max(ub))) else plot(win[,1],type='l',main=main, xlab='j', pch=19, ylab='Winsorized Mean(j/n)', ylim=c(ylimmin,ylimmax))

lines(ub,lty=3)

lines(lb,lty=3)

grid()

dev.off()

bitmap(file='test2.png')

lb <- tri[,1] - 2*tri[,2]

ub <- tri[,1] + 2*tri[,2]

if ((ylimmin == '') | (ylimmax == '')) plot(tri[,1],type='b',main=main, xlab='j', pch=19, ylab='Trimmed Mean(j/n)', ylim=c(min(lb),max(ub))) else plot(tri[,1],type='l',main=main, xlab='j', pch=19, ylab='Trimmed Mean(j/n)', ylim=c(ylimmin,ylimmax))

lines(ub,lty=3)

lines(lb,lty=3)

grid()

dev.off()

load(file='createtable')

a<-table.start()

a<-table.row.start(a)

a<-table.element(a,'Central Tendency - Ungrouped Data',4,TRUE)

a<-table.row.end(a)

a<-table.row.start(a)

a<-table.element(a,'Measure',header=TRUE)

a<-table.element(a,'Value',header=TRUE)

a<-table.element(a,'S.E.',header=TRUE)

a<-table.element(a,'Value/S.E.',header=TRUE)

a<-table.row.end(a)

a<-table.row.start(a)

a<-table.element(a,hyperlink('http://www.xycoon.com/arithmetic_mean.htm', 'Arithmetic Mean', 'click to view the definition of the Arithmetic Mean'),header=TRUE)

a<-table.element(a,arm)

a<-table.element(a,hyperlink('http://www.xycoon.com/arithmetic_mean_standard_error.htm', armse, 'click to view the definition of the Standard Error of the Arithmetic Mean'))

a<-table.element(a,armose)

a<-table.row.end(a)

a<-table.row.start(a)

a<-table.element(a,hyperlink('http://www.xycoon.com/geometric_mean.htm', 'Geometric Mean', 'click to view the definition of the Geometric Mean'),header=TRUE)

a<-table.element(a,geo)

a<-table.element(a,'')

a<-table.element(a,'')

a<-table.row.end(a)

a<-table.row.start(a)

a<-table.element(a,hyperlink('http://www.xycoon.com/harmonic_mean.htm', 'Harmonic Mean', 'click to view the definition of the Harmonic Mean'),header=TRUE)

a<-table.element(a,har)

a<-table.element(a,'')

a<-table.element(a,'')

a<-table.row.end(a)

a<-table.row.start(a)

a<-table.element(a,hyperlink('http://www.xycoon.com/quadratic_mean.htm', 'Quadratic Mean', 'click to view the definition of the Quadratic Mean'),header=TRUE)

a<-table.element(a,qua)

a<-table.element(a,'')

a<-table.element(a,'')

a<-table.row.end(a)

for (j in 1:length(win[,1])) {

a<-table.row.start(a)

mylabel <- paste('Winsorized Mean (',j)

mylabel <- paste(mylabel,'/')

mylabel <- paste(mylabel,length(win[,1]))

mylabel <- paste(mylabel,')')

a<-table.element(a,hyperlink('http://www.xycoon.com/winsorized_mean.htm', mylabel, 'click to view the definition of the Winsorized Mean'),header=TRUE)

a<-table.element(a,win[j,1])

a<-table.element(a,win[j,2])

a<-table.element(a,win[j,1]/win[j,2])

a<-table.row.end(a)

}

for (j in 1:length(tri[,1])) {

a<-table.row.start(a)

mylabel <- paste('Trimmed Mean (',j)

mylabel <- paste(mylabel,'/')

mylabel <- paste(mylabel,length(tri[,1]))

mylabel <- paste(mylabel,')')

a<-table.element(a,hyperlink('http://www.xycoon.com/arithmetic_mean.htm', mylabel, 'click to view the definition of the Trimmed Mean'),header=TRUE)

a<-table.element(a,tri[j,1])

a<-table.element(a,tri[j,2])

a<-table.element(a,tri[j,1]/tri[j,2])

a<-table.row.end(a)

}

a<-table.row.start(a)

a<-table.element(a,hyperlink('http://www.xycoon.com/median_1.htm', 'Median', 'click to view the definition of the Median'),header=TRUE)

a<-table.element(a,median(x))

a<-table.element(a,'')

a<-table.element(a,'')

a<-table.row.end(a)

a<-table.row.start(a)

a<-table.element(a,hyperlink('http://www.xycoon.com/midrange.htm', 'Midrange', 'click to view the definition of the Midrange'),header=TRUE)

a<-table.element(a,midr)

a<-table.element(a,'')

a<-table.element(a,'')

a<-table.row.end(a)

a<-table.row.start(a)

mymid <- hyperlink('http://www.xycoon.com/midmean.htm', 'Midmean', 'click to view the definition of the Midmean')

mylabel <- paste(mymid,hyperlink('http://www.xycoon.com/method_1.htm','Weighted Average at Xnp',''),sep=' - ')

a<-table.element(a,mylabel,header=TRUE)

a<-table.element(a,midm[1])

a<-table.element(a,'')

a<-table.element(a,'')

a<-table.row.end(a)

a<-table.row.start(a)

mymid <- hyperlink('http://www.xycoon.com/midmean.htm', 'Midmean', 'click to view the definition of the Midmean')

mylabel <- paste(mymid,hyperlink('http://www.xycoon.com/method_2.htm','Weighted Average at X(n+1)p',''),sep=' - ')

a<-table.element(a,mylabel,header=TRUE)

a<-table.element(a,midm[2])

a<-table.element(a,'')

a<-table.element(a,'')

a<-table.row.end(a)

a<-table.row.start(a)

mymid <- hyperlink('http://www.xycoon.com/midmean.htm', 'Midmean', 'click to view the definition of the Midmean')

mylabel <- paste(mymid,hyperlink('http://www.xycoon.com/method_3.htm','Empirical Distribution Function',''),sep=' - ')

a<-table.element(a,mylabel,header=TRUE)

a<-table.element(a,midm[3])

a<-table.element(a,'')

a<-table.element(a,'')

a<-table.row.end(a)

a<-table.row.start(a)

mymid <- hyperlink('http://www.xycoon.com/midmean.htm', 'Midmean', 'click to view the definition of the Midmean')

mylabel <- paste(mymid,hyperlink('http://www.xycoon.com/method_4.htm','Empirical Distribution Function - Averaging',''),sep=' - ')

a<-table.element(a,mylabel,header=TRUE)

a<-table.element(a,midm[4])

a<-table.element(a,'')

a<-table.element(a,'')

a<-table.row.end(a)

a<-table.row.start(a)

mymid <- hyperlink('http://www.xycoon.com/midmean.htm', 'Midmean', 'click to view the definition of the Midmean')

mylabel <- paste(mymid,hyperlink('http://www.xycoon.com/method_5.htm','Empirical Distribution Function - Interpolation',''),sep=' - ')

a<-table.element(a,mylabel,header=TRUE)

a<-table.element(a,midm[5])

a<-table.element(a,'')

a<-table.element(a,'')

a<-table.row.end(a)

a<-table.row.start(a)

mymid <- hyperlink('http://www.xycoon.com/midmean.htm', 'Midmean', 'click to view the definition of the Midmean')

mylabel <- paste(mymid,hyperlink('http://www.xycoon.com/method_6.htm','Closest Observation',''),sep=' - ')

a<-table.element(a,mylabel,header=TRUE)

a<-table.element(a,midm[6])

a<-table.element(a,'')

a<-table.element(a,'')

a<-table.row.end(a)

a<-table.row.start(a)

mymid <- hyperlink('http://www.xycoon.com/midmean.htm', 'Midmean', 'click to view the definition of the Midmean')

mylabel <- paste(mymid,hyperlink('http://www.xycoon.com/method_7.htm','True Basic - Statistics Graphics Toolkit',''),sep=' - ')

a<-table.element(a,mylabel,header=TRUE)

a<-table.element(a,midm[7])

a<-table.element(a,'')

a<-table.element(a,'')

a<-table.row.end(a)

a<-table.row.start(a)

mymid <- hyperlink('http://www.xycoon.com/midmean.htm', 'Midmean', 'click to view the definition of the Midmean')

mylabel <- paste(mymid,hyperlink('http://www.xycoon.com/method_8.htm','MS Excel (old versions)',''),sep=' - ')

a<-table.element(a,mylabel,header=TRUE)

a<-table.element(a,midm[8])

a<-table.element(a,'')

a<-table.element(a,'')

a<-table.row.end(a)

a<-table.row.start(a)

a<-table.element(a,'Number of observations',header=TRUE)

a<-table.element(a,length(x))

a<-table.element(a,'')

a<-table.element(a,'')

a<-table.row.end(a)

a<-table.end(a)

table.save(a,file='mytable.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