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

-0.0447135254057231
-0.0681917505478084
0.197918269607118
0.364977510288305
151.337.996.933.159
-0.359461693683337
0.457275231643101
-0.576200433829823
-0.358061445786517
126.000.911.969.509
-134.266.139.991.407
-0.353289571497946
0.161395641494175
-0.85782300297493
-0.555144912298359
-0.727027913044961
-0.369042261655499
-0.0578971497992717
-0.769036818576891
-0.853703554581848
0.700854566672244
-0.676058766832693
0.868578036694034
-0.10884588975301
-157.887.204.236.819
-112.268.950.514.398
0.39692385190415
0.0148495697973627
0.140259353790028
0.369865091064737
-0.279776101924426
0.304546129475365
0.892466302715145
0.0893800831442478
0.135269394231309
-120.736.758.947.938
-0.230372958786898
-0.0877971490509319
-0.33323296780618
0.60196627511901
0.48979033192109
-0.300400372794546
0.101351037939174
0.596346950582894
-0.476838868063972
-0.418379377485807
-0.117113003673828
-0.147180754345273
0.532293994236285
-0.976528360213732
0.331179682756605
0.869619341823574
-0.452578488121829
-0.42138949463861
-0.0683600581129025
0.326102644103438
0.848949464327804
0.455827855730794
0.823473300311693
0.933736135720844
123.343.911.512.994
0.408339888885158
0.287976014212968
-0.211203400023637
-0.267468886934765
-111.645.429.924.972
-0.0317247398904611
0.547972212916532
0.120802505225608
-0.868321571501203
-0.311155680962737
-0.380689102119501
-0.305618643654202
-0.409940938228271
-0.00967793776263202
0.495426344923855
-0.704347377621387
-0.0634698792084991
-0.513898813333212
0.589511460763565
-0.349423192928886
0.641305316979721
0.0343658031299827
-0.0527521706610295
-0.880138087903553
-0.110617811300362
0.725017567873349
-0.50521028991854
1.015.061.557.493
0.432165488939549
-0.605966485103154
-0.999657590805216
-0.273898494723043
0.276038635782484
100.336.898.138.124
0.0129978236302341
-0.560227111581804
-0.549704221644982
-0.264710628774072
0.279628858544352
0.678638849882419
0.661069134033749
-0.704450399639775
-0.201260460274452
0.350125398121578
0.512074040088484
0.855206906549666
0.194069541456855
0.723929039883086
118.131.067.430.518
0.142881527616885
0.0179257999324864
-0.497956664934098
-0.299719686680144
0.135788266531753
-0.170331262015486
-103.309.848.477.463
-0.173819270986569
-0.963024089117263
0.698379739279439
-0.367089110910285
-0.263848476436372
-0.418819019840622
-112.067.684.130.113
0.0500585205773671
0.608903578527248
0.134690293435901
0.330041347787177
0.127838692569612
0.581319132516197
-0.038707943536881
-0.870097194938441
-0.462113796849836
-0.775456514875537
139.586.946.762.297
-0.373271731449701
-0.268080176510487
108.656.232.048.461
-0.325879955921929
0.30661854788457
-0.552774192453423
0.92873551861249
0.0930554068879337
0.793725485236928
-0.0655196055781843
0.31938973725902
-0.466224915112164
-0.262872803122968
0.0101377300966978
0.166709542388761
-0.27948869702557
-0.416437753013764
-0.221155653982262
-0.180257931963139
-0.698512272287555
-0.030966150632744
-0.218027303452022
-0.373584442457798
0.052951427314271
0.16272689702871
0.827781138510899
-0.724325872700048
-0.470004756981797
104.570.491.743.371
-0.192863384714994
-0.571180242416833
0.531291140789517
-0.305718227359467
0.338050048474922
0.476242122330785
-0.81389847635506
-0.0347526256916431
0.307275272687164
0.298265453953272
-0.356205547514144
-0.356604247457693
0.16191437451988
0.157622107243432
0.301928819166272
-0.561981939800824
-0.0543444359644914
-0.242625453729012
-0.0291188126122587
0.228470615616776
-0.510675391591978
112.387.561.916.997
-112.598.799.910.629
0.290323760058927
0.190186117613125
0.0853119358914609
-0.708038983963899
0.0821216895899207
-0.307060856617881
0.525053351110559
-0.601608778527341
0.538181815538459
-0.306512070731953
0.654542044064886
-0.537374859392079
-0.186945154314925
-0.0415286501896635
-0.0882630415428836
-0.256492309570866
-0.41305922709641
-0.270151685278353
-0.435615700218057
0.546649355951026
0.323884533274196
0.661932000134874
0.402774353219053
-0.428179650529226
-0.185363732952745
-0.146275810216944
0.177901302311226
-0.370503387950712
0.205238677067021
-0.0894107718996266
-0.0207802658241128
-0.00240155924277281
0.0272372593388978
-0.30401211751709
150.787.083.654.869
0.259049823633271
-0.546357645795726
0.581469268528577
0.330213798462058
-0.983366580775167
-0.736262986208755
-0.338744406461673
0.760537197604814
-0.261709569948611
-0.476048053943411
-0.110367750063647
169.080.994.427.218
0.149871963217245
-0.894443959355815
0.0854707859389192
-0.11791601547532
-0.215892867114871
-0.394092028282406
0.0924263195771732
0.0585916699118239
0.253245871086333
0.370568028643833
-0.386523945873816
0.447110885622225
0.623154288482199
-0.184203126128336
0.705478214339162
-0.317030724633314
-0.927475247660233
-0.0393758868190344
100.274.756.958.954
0.812619294578388
0.298293415694118
0.152551035780658
-0.150639291567325
0.110741928054486
0.552786542794449
0.0434953176436756
0.36320976255479
-0.0227970084706981
0.504138678148357
0.138779289305853
-0.0853035308417627
-0.496484911315902
-0.0852251911409686
-0.247830972313436
-0.120849214187909
-0.451770225463026
0.612150462832891
0.350625301752197
-0.359290488889771
-0.48443652903107
0.339430852270146
-0.0417944819768526
-0.0100068969566531
-0.403912575239062
0.151831207953882
-0.229701003816198
-0.217361643531106
-0.332188937105688
0.264417994027136
0.176869376138034
-0.176677800263279
-0.13567036715417
-0.827944200296052
-0.0722734399204132
-0.117607845832746
0.314635053724412
-0.154386181086533
0.163002765901451
-0.244739215669361
-0.204481562449259
0.0360259309576058
0.00417168114426594
0.266111866388443
-0.651003114649842
0.589994534676577
0.312322955620027
0.55266992874864
-0.0962287489408765
-0.467083613141822
-0.198196785314411
0.395575343631519
0.175735234004809
0.315653622056974
-0.16114907850581
0.882040919470631
0.0637848764427299
0.858633851994301
0.822761597215925
177.266.886.690.859
-0.56698109822252
0.153261675817821
-0.27923293044209
0.0495720504717388
-116.973.330.350.604
-0.0821094023297278
0.0681021355630733
-0.201544972117822
0.0764370043961384
-0.526360142929291
-0.080413992290807
-0.390603734038536
-0.366201530340293
-0.236325627413761
-0.0197748907207431
-0.400593230734885
0.273305313443572
0.571670833577946
0.355107022563815
-0.598053247795437
0.0680124314430932
0.0825157705041395
-0.236688909377818
-0.720993035857936
0.445855527021637
-0.278184399175569
-0.813695706192547
0.0382223634028161
0.287981462997023
-0.397401294088736
0.453441800201652
-0.336993631715909
0.0350016901300973
-0.129074066212103
-0.92948724921499
0.112538313206619
-0.266214839380579
0.318919353798613
-0.235823743806211
0.312056037492047
-0.607348899940245
0.821439552875318
-0.330769416907474
-0.0368470896545833
-0.223411512848023
-0.0162393530813711
0.494328753928909

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