R version 2.9.0 (2009-04-17) Copyright (C) 2009 The R Foundation for Statistical Computing ISBN 3-900051-07-0 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > x <- c(-0.0042361726044575 + ,-4.3964831775156E-14 + ,0.0042361726044557 + ,-0.0042361726044575 + ,0.0042361726044557 + ,4.2188474935756E-14 + ,4.2188474935756E-14 + ,4.2188474935756E-14 + ,4.2188474935756E-14 + ,-0.0042361726044575 + ,0.0042361726044557 + ,4.2188474935756E-14 + ,4.2188474935756E-14 + ,-0.0042361726044575 + ,-4.3964831775156E-14 + ,-4.3964831775156E-14 + ,-4.3964831775156E-14 + ,-4.3964831775156E-14 + ,-4.3964831775156E-14 + ,0.0042361726044557 + ,4.2188474935756E-14 + ,4.2188474935756E-14 + ,-0.0042361726044575 + ,-4.3964831775156E-14 + ,0.0042361726044557 + ,-0.0042361726044575 + ,-4.3964831775156E-14 + ,-4.3964831775156E-14 + ,-4.3964831775156E-14 + ,0.0042361726044557 + ,-0.0042361726044575 + ,-4.3964831775156E-14 + ,-4.3964831775156E-14 + ,-4.3964831775156E-14 + ,-4.3964831775156E-14 + ,0.0042361726044557 + ,4.2188474935756E-14 + ,-0.0042361726044575 + ,-4.3964831775156E-14 + ,-4.3964831775156E-14 + ,0.0052924177514559 + ,-0.0052924177515092 + ,0.0052924177514559 + ,-9.3258734068513E-15 + ,-9.3258734068513E-15 + ,-0.0052924177515092 + ,-4.3964831775156E-14 + ,-4.3964831775156E-14 + ,-4.3964831775156E-14 + ,-4.3964831775156E-14 + ,0.0052924177514559 + ,-9.3258734068513E-15 + ,-9.3258734068513E-15 + ,-9.3258734068513E-15 + ,-9.3258734068513E-15 + ,-0.0052924177515092 + ,0.0052924177514559 + ,-9.3258734068513E-15 + ,-9.3258734068513E-15 + ,-9.3258734068513E-15 + ,-0.0052924177515092 + ,0.0052924177514559 + ,-0.0052924177515092 + ,-4.3964831775156E-14 + ,-4.3964831775156E-14 + ,0.0052924177514559 + ,-9.3258734068513E-15 + ,-0.0052924177515092 + ,-4.3964831775156E-14 + ,0.0052924177514559 + ,-9.3258734068513E-15 + ,-9.3258734068513E-15 + ,-0.0052924177515092 + ,-4.3964831775156E-14 + ,0.0015906684159561 + ,-3.3306690738755E-14 + ,-0.0053120974849334 + ,0.0053120974848566 + ,-3.3306690738755E-14 + ,-3.3306690738755E-14 + ,-3.3306690738755E-14 + ,-3.3306690738755E-14 + ,-3.3306690738755E-14 + ,-3.3306690738755E-14 + ,-0.0053120974849334 + ,-4.3520742565306E-14 + ,0.0053120974848566 + ,-3.3306690738755E-14 + ,-3.3306690738755E-14 + ,-3.3306690738755E-14 + ,-3.3306690738755E-14 + ,-3.3306690738755E-14 + ,-3.3306690738755E-14 + ,-3.3306690738755E-14 + ,-3.3306690738755E-14 + ,-0.0053120974849334 + ,-4.3520742565306E-14 + ,-4.3520742565306E-14 + ,-4.3520742565306E-14 + ,-0.0093646169310437 + ,0.014676714415899 + ,-0.0053120974849334 + ,0.0053120974848566 + ,-3.3306690738755E-14 + ,-3.3306690738755E-14 + ,-3.3306690738755E-14 + ,-0.0053120974849334 + ,0.0013306721856563 + ,-0.91895384929365 + ,0.91895384929359 + ,-4.8849813083507E-14 + ,-4.8849813083507E-14 + ,-4.8849813083507E-14 + ,-0.0047987295499485 + ,-4.3520742565306E-14 + ,-4.3520742565306E-14 + ,-4.3520742565306E-14 + ,-4.3520742565306E-14 + ,-4.3520742565306E-14 + ,0.0047987295498562 + ,-4.8849813083507E-14 + ,-0.0026631174194489 + ,-0.013423020332165 + ,-0.0029774010356243 + ,-0.0048913140998841 + ,-4.3076653355456E-14 + ,-4.3076653355456E-14 + ,0.004891314099857 + ,1.5987211554602E-14 + ,1.5987211554602E-14 + ,-0.0048913140998841 + ,-4.3076653355456E-14 + ,-4.3076653355456E-14 + ,-4.3076653355456E-14 + ,-4.3076653355456E-14 + ,-4.3076653355456E-14 + ,0.004891314099857 + ,1.5987211554602E-14 + ,-0.0048913140998841 + ,0.004891314099857 + ,1.5987211554602E-14 + ,1.5987211554602E-14 + ,1.5987211554602E-14 + ,1.5987211554602E-14 + ,-0.0048913140998841 + ,0.004891314099857 + ,1.5987211554602E-14 + ,1.5987211554602E-14 + ,-0.0048913140998841 + ,0.004891314099857 + ,-0.0048913140998841 + ,-4.3076653355456E-14 + ,-4.3076653355456E-14 + ,-4.3076653355456E-14 + ,-4.3076653355456E-14 + ,-4.3076653355456E-14 + ,-4.3076653355456E-14 + ,-4.3076653355456E-14 + ,-4.3076653355456E-14 + ,-4.3076653355456E-14 + ,-4.3076653355456E-14 + ,-4.3076653355456E-14 + ,-4.3076653355456E-14 + ,-4.3076653355456E-14 + ,-4.3076653355456E-14 + ,-4.3076653355456E-14 + ,-4.3076653355456E-14 + ,-4.3076653355456E-14 + ,0.004891314099857 + ,-0.0048913140998841 + ,-4.3076653355456E-14 + ,0.004891314099857 + ,1.5987211554602E-14 + ,-0.0048913140998841 + ,-4.3076653355456E-14 + ,0.004891314099857 + ,-0.0048913140998841 + ,-4.3076653355456E-14 + ,-4.3076653355456E-14 + ,-4.3076653355456E-14 + ,-4.3076653355456E-14 + ,0.004891314099857 + ,1.5987211554602E-14 + ,-0.0048913140998841 + ,-4.3076653355456E-14 + ,0.004891314099857 + ,1.5987211554602E-14 + ,1.5987211554602E-14 + ,1.5987211554602E-14 + ,1.5987211554602E-14 + ,-0.0048913140998841 + ,-0.018974857760943 + ,-1.0214051826551E-14 + ,-1.0214051826551E-14 + ,-1.0214051826551E-14 + ,0.02684357289639 + ,-2.4424906541753E-14 + ,-0.026843572896424 + ,-1.0214051826551E-14 + ,0.02684357289639) > ylimmax = '' > ylimmin = '' > main = 'Robustness of Central Tendency' > #'GNU S' R Code compiled by R2WASP v. 1.0.44 () > #Author: Prof. Dr. P. Wessa > #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/ > #Source of accompanying publication: Office for Research, Development, and Education > #Technical description: Write here your technical program description (don't use hard returns!) > 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)) [1] -0.0001135653 > sqrtn <- sqrt(length(x)) > (armse <- sd(x) / sqrtn) [1] 0.006522883 > (armose <- arm / armse) [1] -0.01741030 > (geo <- geomean(x)) [1] NaN Warning message: In log(x) : NaNs produced > (har <- harmean(x)) [1] -7.512089e-14 > (qua <- quamean(x)) [1] 0.09201666 > (win <- winmean(x)) [,1] [,2] [1,] -1.135653e-04 3.850024e-04 [2,] -3.487818e-05 3.607886e-04 [3,] -1.341035e-04 2.800857e-04 [4,] -2.402278e-04 2.239217e-04 [5,] -1.389148e-04 2.070312e-04 [6,] -1.389148e-04 2.070312e-04 [7,] -1.396036e-04 2.069402e-04 [8,] -1.396036e-04 2.069402e-04 [9,] -1.396036e-04 2.069402e-04 [10,] -1.386196e-04 2.068168e-04 [11,] -1.386196e-04 2.068168e-04 [12,] -1.386196e-04 2.068168e-04 [13,] -1.386196e-04 2.068168e-04 [14,] -1.666968e-04 2.032074e-04 [15,] -1.666968e-04 2.032074e-04 [16,] -1.666968e-04 2.032074e-04 [17,] -1.326030e-04 1.989970e-04 [18,] -1.326030e-04 1.989970e-04 [19,] -1.326030e-04 1.989970e-04 [20,] -1.326030e-04 1.989970e-04 [21,] -1.326030e-04 1.989970e-04 [22,] -1.326030e-04 1.989970e-04 [23,] -1.326030e-04 1.989970e-04 [24,] -1.437132e-04 1.975940e-04 [25,] -2.140328e-04 1.890088e-04 [26,] -2.140328e-04 1.890088e-04 [27,] -2.140328e-04 1.890088e-04 [28,] -2.010710e-04 1.874039e-04 [29,] -1.195002e-04 1.776181e-04 [30,] -1.195002e-04 1.776181e-04 [31,] -1.195002e-04 1.776181e-04 [32,] -5.427809e-04 1.332125e-04 [33,] -5.856803e-04 1.298943e-04 [34,] -8.118945e-04 1.173507e-04 [35,] -8.118945e-04 1.173507e-04 [36,] -8.118945e-04 1.173507e-04 [37,] -5.790218e-04 8.341122e-05 [38,] -5.193079e-04 7.479611e-05 [39,] -1.802780e-14 2.560643e-15 [40,] -1.802780e-14 2.560643e-15 [41,] -1.802780e-14 2.560643e-15 [42,] -1.802780e-14 2.560643e-15 [43,] -2.261080e-14 1.866627e-15 [44,] -2.261080e-14 1.866627e-15 [45,] -2.261080e-14 1.866627e-15 [46,] -2.261080e-14 1.866627e-15 [47,] -2.261080e-14 1.866627e-15 [48,] -2.261080e-14 1.866627e-15 [49,] -2.261080e-14 1.866627e-15 [50,] -2.261080e-14 1.866627e-15 [51,] -2.261080e-14 1.866627e-15 [52,] -2.261080e-14 1.866627e-15 [53,] -2.261080e-14 1.866627e-15 [54,] -2.261080e-14 1.866627e-15 [55,] -2.261080e-14 1.866627e-15 [56,] -2.261080e-14 1.866627e-15 [57,] -2.261080e-14 1.866627e-15 [58,] -2.995160e-14 1.140011e-15 [59,] -2.995160e-14 1.140011e-15 [60,] -2.995160e-14 1.140011e-15 [61,] -2.995160e-14 1.140011e-15 [62,] -2.995160e-14 1.140011e-15 [63,] -2.995160e-14 1.140011e-15 [64,] -2.995160e-14 1.140011e-15 [65,] -2.995160e-14 1.140011e-15 [66,] -2.995160e-14 1.140011e-15 > (tri <- trimean(x)) [,1] [,2] [1,] -1.147125e-04 3.380776e-04 [2,] -1.158830e-04 2.809617e-04 [3,] -1.576381e-04 2.272885e-04 [4,] -1.658098e-04 2.048053e-04 [5,] -1.462261e-04 1.991011e-04 [6,] -1.477817e-04 1.971925e-04 [7,] -1.493707e-04 1.951576e-04 [8,] -1.508874e-04 1.930031e-04 [9,] -1.524373e-04 1.907028e-04 [10,] -1.540218e-04 1.882448e-04 [11,] -1.557523e-04 1.856330e-04 [12,] -1.575223e-04 1.828362e-04 [13,] -1.593329e-04 1.798376e-04 [14,] -1.611856e-04 1.766180e-04 [15,] -1.607224e-04 1.735783e-04 [16,] -1.602483e-04 1.703082e-04 [17,] -1.597627e-04 1.667838e-04 [18,] -1.617110e-04 1.634308e-04 [19,] -1.637074e-04 1.598081e-04 [20,] -1.657538e-04 1.558849e-04 [21,] -1.678519e-04 1.516244e-04 [22,] -1.700039e-04 1.469826e-04 [23,] -1.722117e-04 1.419057e-04 [24,] -1.744777e-04 1.363269e-04 [25,] -1.761868e-04 1.303209e-04 [26,] -1.741411e-04 1.246211e-04 [27,] -1.720393e-04 1.182810e-04 [28,] -1.698792e-04 1.111657e-04 [29,] -1.698792e-04 1.032951e-04 [30,] -1.707146e-04 9.555495e-05 [31,] -1.731887e-04 8.658419e-05 [32,] -1.757356e-04 7.590547e-05 [33,] -1.586159e-04 6.952115e-05 [34,] -1.390079e-04 6.235876e-05 [35,] -1.085605e-04 5.477516e-05 [36,] -7.716165e-05 4.515133e-05 [37,] -4.476602e-05 3.157717e-05 [38,] -2.147675e-05 2.147675e-05 [39,] -2.742433e-14 2.204313e-15 [40,] -2.782589e-14 2.156421e-15 [41,] -2.824106e-14 2.102677e-15 [42,] -2.867055e-14 2.042160e-15 [43,] -2.911511e-14 1.973717e-15 [44,] -2.938522e-14 1.963300e-15 [45,] -2.966516e-14 1.951024e-15 [46,] -2.995546e-14 1.936651e-15 [47,] -3.025672e-14 1.919904e-15 [48,] -3.056956e-14 1.900457e-15 [49,] -3.089468e-14 1.877928e-15 [50,] -3.123279e-14 1.851863e-15 [51,] -3.158471e-14 1.821715e-15 [52,] -3.195129e-14 1.786824e-15 [53,] -3.233347e-14 1.746377e-15 [54,] -3.273227e-14 1.699359e-15 [55,] -3.314879e-14 1.644476e-15 [56,] -3.358425e-14 1.580041e-15 [57,] -3.358425e-14 1.503779e-15 [58,] -3.403995e-14 1.412495e-15 [59,] -3.470936e-14 1.408840e-15 [60,] -3.491096e-14 1.403544e-15 [61,] -3.512290e-14 1.396318e-15 [62,] -3.534600e-14 1.386810e-15 [63,] -3.558115e-14 1.374589e-15 [64,] -3.582936e-14 1.359125e-15 [65,] -3.609176e-14 1.339753e-15 [66,] -3.636960e-14 1.315630e-15 > (midr <- midrange(x)) [1] -2.997602e-14 > midm <- array(NA,dim=8) > for (j in 1:8) midm[j] <- midmean(x,j) > midm [1] -2.911511e-14 -2.911511e-14 -2.911511e-14 -2.911511e-14 -2.911511e-14 [6] -2.911511e-14 -2.911511e-14 -2.911511e-14 > postscript(file="/var/www/html/rcomp/tmp/1ibw01290523387.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > 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() null device 1 > postscript(file="/var/www/html/rcomp/tmp/2s2v31290523387.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > 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() null device 1 > > #Note: the /var/www/html/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/html/rcomp/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="/var/www/html/rcomp/tmp/3vnm21290523387.tab") > > try(system("convert tmp/1ibw01290523387.ps tmp/1ibw01290523387.png",intern=TRUE)) character(0) > try(system("convert tmp/2s2v31290523387.ps tmp/2s2v31290523387.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.057 0.356 2.161