• QR Code (size=2)
• QR Code (size=3)
• QR Code (size=4)
• QR Code (size=5)
• QR Code (size=6)

-0.0447135254057231
-0.0681917505477611
0.197918269607264
0.364977510288303
151.337.996.933.282
-0.359461693682944
0.457275231630557
-0.576200433829858
-0.35806144578699
126.000.911.969.899
-134.266.139.996.902
-0.35328957148349
0.161395641520833
-0.857823002973426
-0.555144912300313
-0.727027913044969
-0.369042261654297
-0.0578971497923536
-0.769036818583236
-0.853703554579139
0.700854566671402
-0.676058766811168
0.868578036660101
-0.108845889757913
-157.887.204.236.239
-112.268.950.514.506
0.396923851905675
0.0148495697982022
0.140259353790972
0.36986509106495
-0.279776101924479
0.304546129479045
0.892466302712004
0.0893800831529647
0.135269394216119
-120.736.758.947.528
-0.230372958786352
-0.0877971490505129
-0.333232967805547
0.601966275119466
0.489790331921813
-0.300400372795234
0.101351037939037
0.596346950586006
-0.47683886806596
-0.418379377480445
-0.117113003681212
-0.14718075434479
0.532293994236984
-0.976528360213378
0.331179682756783
0.869619341824002
-0.452578488121401
-0.421389494638729
-0.0683600581132898
0.326102644105208
0.848949464326408
0.455827855734715
0.823473300306631
0.933736135720925
123.343.911.513.023
0.408339888885549
0.287976014213179
-0.211203400023215
-0.267468886934195
-111.645.429.924.937
-0.0317247398908312
0.547972212918
0.120802505224662
-0.868321571498795
-0.311155680966178
-0.380689102119736
-0.305618643653553
-0.409940938228127
-0.00967793776262078
0.495426344923938
-0.704347377621235
-0.063469879208485
-0.51389881333356
0.58951146076434
-0.349423192929551
0.641305316981306
0.0343658031275834
-0.052752170661352
-0.880138087903201
-0.110617811300259
0.725017567873533
-0.505210289918464
101.506.155.749.305
0.432165488939686
-0.605966485103475
-0.999657590804615
-0.273898494723388
0.276038635783788
100.336.898.137.961
0.0129978236300762
-0.560227111581836
-0.549704221644895
-0.264710628773928
0.279628858544481
0.678638849882516
0.661069134033666
-0.704450399639958
-0.201260460274066
0.350125398121403
0.512074040089319
0.855206906548568
0.194069541456819
0.723929039883118
118.131.067.430.535
0.142881527617021
0.017925799932583
-0.497956664933812
-0.299719686679989
0.135788266531892
-0.170331262015063
-103.309.848.477.483
-0.173819270986012
-0.963024089117913
0.698379739279324
-0.367089110910225
-0.263848476436583
-0.418819019840626
-112.067.684.130.117
0.050058520577361
0.608903578527243
0.134690293435698
0.330041347787194
0.127838692569412
0.581319132516513
-0.0387079435373926
-0.870097194938493
-0.462113796849748
-0.775456514875374
1.395.869.467.623
-0.373271731449615
-0.268080176510708
108.656.232.048.462
-0.325879955921928
0.306618547884575
-0.552774192453413
0.928735518612731
0.0930554068876903
0.793725485236845
-0.0655196055780801
0.319389737259029
-0.466224915112025
-0.262872803122903
0.0101377300968685
0.166709542388814
-0.279488697025563
-0.416437753013679
-0.221155653982297
-0.180257931962942
-0.698512272287848
-0.0309661506328094
-0.218027303452013
-0.373584442457897
0.0529514273142018
0.16272689702866
0.827781138510815
-0.724325872700113
-0.470004756981956
104.570.491.743.386
-0.192863384714999
-0.57118024241685
0.531291140789378
-0.305718227359383
0.338050048474844
0.476242122330853
-0.813898476355091
-0.0347526256916884
0.3072752726873
0.298265453953271
-0.356205547514185
-0.356604247457695
0.1619143745199
0.157622107243565
0.301928819166098
-0.561981939800843
-0.054344435964539
-0.242625453728938
-0.0291188126122691
0.228470615616782
-0.510675391591998
112.387.561.916.988
-112.598.799.910.622
0.290323760058784
0.190186117613264
0.0853119358914465
-0.70803898396401
0.0821216895898958
-0.307060856617796
0.525053351110503
-0.601608778527328
0.53818181553833
-0.306512070731858
0.654542044064766
-0.537374859392012
-0.186945154315006
-0.0415286501896012
-0.0882630415428376
-0.256492309570936
-0.413059227096436
-0.270151685278362
-0.435615700218065
0.546649355950973
0.323884533274175
0.661932000134738
0.402774353219003
-0.428179650529287
-0.185363732952725
-0.146275810216839
0.177901302311314
-0.370503387950714
0.205238677066994
-0.0894107718995766
-0.0207802658241447
-0.00240155924277742
0.027237259338887
-0.304012117517099
150.787.083.654.866
0.259049823633338
-0.546357645795879
0.581469268528633
0.330213798462246
-0.983366580775191
-0.736262986208729
-0.338744406461511
0.76053719760488
-0.261709569948618
-0.476048053943579
-0.110367750063661
169.080.994.427.217
0.149871963217241
-0.894443959356019
0.0854707859389636
-0.117916015475136
-0.215892867114885
-0.394092028282386
0.0924263195771785
0.058591669911839
0.253245871086277
0.370568028643805
-0.386523945873864
0.44711088562219
0.623154288482276
-0.184203126128357
0.705478214339145
-0.317030724633186
-0.927475247660286
-0.0393758868189332
100.274.756.958.968
0.812619294578401
0.298293415694001
0.152551035780635
-0.150639291567251
0.110741928054607
0.552786542794576
0.043495317643748
0.363209762554804
-0.022797008470584
0.504138678148394
0.138779289305952
-0.0853035308417428
-0.496484911315823
-0.0852251911408867
-0.247830972313322
-0.120849214187901
-0.451770225462982
0.612150462832856
0.350625301752204
-0.359290488889893
-0.484436529031112
0.339430852270173
-0.0417944819768093
-0.0100068969567184
-0.403912575239055
0.151831207953865
-0.229701003816159
-0.217361643531165
-0.332188937105677
0.264417994027107
0.176869376138016
-0.176677800263361
-0.135670367154218
-0.827944200296059
-0.0722734399204414
-0.117607845832689
0.314635053724339
-0.154386181086592
0.163002765901339
-0.244739215669379
-0.204481562449311
0.0360259309576101
0.00417168114425691
0.266111866388402
-0.651003114649863
0.589994534676515
0.312322955620105
0.552669928748546
-0.0962287489408798
-0.467083613141848
-0.198196785314342
0.395575343631614
0.175735234004858
0.315653622056923
-0.161149078505793
0.882040919470625
0.0637848764428182
0.858633851994241
0.822761597216032
17.726.688.669.086
-0.566981098222408
0.153261675817818
-0.279232930441738
0.049572050471916
-116.973.330.350.582
-0.0821094023296669
0.0681021355632607
-0.201544972117859
0.0764370043960761
-0.52636014292932
-0.0804139922908793
-0.39060373403853
-0.366201530340365
-0.236325627413792
-0.0197748907207889
-0.400593230734952
0.273305313443463
0.571670833577912
0.355107022563698
-0.598053247795535
0.0680124314430453
0.0825157705042455
-0.236688909377822
-0.720993035857933
0.445855527021645
-0.278184399175499
-0.813695706192672
0.0382223634028053
0.287981462997034
-0.397401294088841
0.453441800201533
-0.336993631715903
0.0350016901299891
-0.12907406621208
-0.92948724921503
0.112538313206602
-0.266214839380531
0.318919353798501
-0.235823743806231
0.312056037491925
-0.607348899940247
0.821439552875221
-0.330769416907408
-0.0368470896547259
-0.22341151284797
-0.0162393530813855
0.494328753928939

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('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('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('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('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('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('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('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('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('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('midmean.htm', 'Midmean', 'click to view the definition of the Midmean')
mylabel <- paste(mymid,hyperlink('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('midmean.htm', 'Midmean', 'click to view the definition of the Midmean')
mylabel <- paste(mymid,hyperlink('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('midmean.htm', 'Midmean', 'click to view the definition of the Midmean')
mylabel <- paste(mymid,hyperlink('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('midmean.htm', 'Midmean', 'click to view the definition of the Midmean')
mylabel <- paste(mymid,hyperlink('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('midmean.htm', 'Midmean', 'click to view the definition of the Midmean')
mylabel <- paste(mymid,hyperlink('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('midmean.htm', 'Midmean', 'click to view the definition of the Midmean')
mylabel <- paste(mymid,hyperlink('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('midmean.htm', 'Midmean', 'click to view the definition of the Midmean')
mylabel <- paste(mymid,hyperlink('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('midmean.htm', 'Midmean', 'click to view the definition of the Midmean')
mylabel <- paste(mymid,hyperlink('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')