R version 2.8.0 (2008-10-20) Copyright (C) 2008 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(-183923544514060.0000 + ,-177072609114252.0000 + ,-228635109114254.0000 + ,-237447609114242.0000 + ,-127760109114251.0000 + ,-193010109114250.0000 + ,-220635109114248.0000 + ,-164510109114256.0000 + ,-268322609114246.0000 + ,-333697609114249.0000 + ,-342601091142503.0000 + ,-154885109114249.0000 + ,-977452805255622.0000 + ,101105654874272.0000 + ,254315487427317.0000 + ,-432693451257275.0000 + ,-163581845125727.0000 + ,-162831845125727.0000 + ,465431548742728.0000 + ,266681548742733.0000 + ,-107144345125727.0000 + ,424806548742729.0000 + ,76918154874273.0000 + ,196293154874273.0000 + ,20143298346296.0000 + ,122839188627955.0000 + ,-0.278581137204542 + ,429089188627948.0000 + ,875964188627952.0000 + ,843464188627952.0000 + ,57721418862795.0000 + ,173846418862796.0000 + ,-185966081137205.0000 + ,476589188627952.0000 + ,890964188627953.0000 + ,-685285811372049.0000 + ,272611247451482.0000 + ,146462182851318.0000 + ,162899682851318.0000 + ,100871828513171.0000 + ,279774682851318.0000 + ,212524682851317.0000 + ,248899682851317.0000 + ,-419753171486821.0000 + ,-578781714868267.0000 + ,528371828513174.0000 + ,274274682851318.0000 + ,414649682851317.0000 + ,310789511440004.0000 + ,36264044683984.0000 + ,260779468398400.0000 + ,403265446839839.0000 + ,32795294683984.0000 + ,193702946839840.0000 + ,31707794683984.0000 + ,20220294683984.0000 + ,32139044683984.0000 + ,178015446839840.0000 + ,164529468398398.0000 + ,-681720531601603.0000 + ,-157032224571473.0000 + ,-761812891716376.0000 + ,-817437891716377.0000 + ,-134556289171638.0000 + ,771312108283621.0000 + ,199881210828362.0000 + ,105256210828362.0000 + ,198381210828362.0000 + ,262568710828362.0000 + ,196193710828362.0000 + ,116312108283621.0000 + ,-145993789171638.0000 + ,-166853960582951.0000 + ,-202003025183115.0000 + ,434344748168846.0000 + ,-113378025183116.0000 + ,-113690525183116.0000 + ,-155940525183116.0000 + ,-210565525183116.0000 + ,-124440525183115.0000 + ,-642530251831158.0000 + ,-298628025183116.0000 + ,-154190525183116.0000 + ,231844748168843.0000 + ,-249675696594429.0000 + ,118175238805407.0000 + ,-180387261194593.0000 + ,-791997611945937.0000 + ,-815122611945933.0000 + ,-246762261194593.0000 + ,-105387261194593.0000 + ,-319262261194593.0000 + ,-720747611945935.0000 + ,-904497611945934.0000 + ,-800122611945933.0000 + ,119362738805407.0000 + ,-534974326059064.0000 + ,-114646497206071.0000 + ,-155208997206071.0000 + ,-500214972060714.0000 + ,-196333997206071.0000 + ,-145839972060711.0000 + ,-822089972060712.0000 + ,179160027939294.0000 + ,-162896497206071.0000 + ,-132271497206071.0000 + ,-168339972060710.0000 + ,815410027939288.0000 + ,275680831382616.0000 + ,-324682332175485.0000 + ,179692667824515.0000 + ,271567667824508.0000 + ,-123155733217549.0000 + ,108594266782451.0000 + ,679692667824511.0000 + ,340942667824517.0000 + ,-137182332175489.0000 + ,-113093233217549.0000 + ,543442667824512.0000 + ,149719266782451.0000 + ,153859095371138.0000 + ,-282899692290262.0000 + ,238147530770974.0000 + ,503350307709731.0000 + ,802253077097358.0000 + ,-612274692290265.0000 + ,-140852469229027.0000 + ,-287274692290261.0000 + ,946003077097336.0000 + ,-121914969229026.0000 + ,415225307709736.0000 + ,115897530770973.0000 + ,270373593596605.0000 + ,-91111705240504.0000 + ,332579475949605.0000 + ,-294867052405046.0000 + ,-737992052405041.0000 + ,509507947594957.0000 + ,-866742052405044.0000 + ,-954920524050376.0000 + ,-663617052405043.0000 + ,732632947594957.0000 + ,-216299205240504.0000 + ,-128924205240504.0000 + ,-142784376651817.0000 + ,270665587480184.0000 + ,605040587480183.0000 + ,356915587480177.0000 + ,163790587480182.0000 + ,-64870941251982.0000 + ,115504058748018.0000 + ,-303709412519815.0000 + ,87816558748018.0000 + ,205441558748018.0000 + ,-641209412519819.0000 + ,-322745941251982.0000 + ,-139606112663295.0000 + ,352448227365406.0000 + ,-431767726345942.0000 + ,178698227365400.0000 + ,255732273654045.0000 + ,129307322736540.0000 + ,-163176772634598.0000 + ,164807322736541.0000 + ,219948227365402.0000 + ,138619822736540.0000 + ,870573227365404.0000 + ,514323227365403.0000 + ,-804278486747727.0000 + ,-105191879672279.0000 + ,524562032772056.0000 + ,684331203277199.0000 + ,-0.879379672279542 + ,-105129379672280.0000 + ,-827543796722797.0000 + ,-132629379672279.0000 + ,102558120327720.0000 + ,231831203277203.0000 + ,-180379379672280.0000 + ,-26700437967228.0000 + ,301354489164073.0000 + ,23986384316243.0000 + ,904238843162428.0000 + ,316113843162422.0000 + ,812988843162427.0000 + ,250488843162426.0000 + ,-135761156837575.0000 + ,33548884316243.0000 + ,140736384316242.0000 + ,132361384316243.0000 + ,947988843162427.0000 + ,417388431624253.0000) > 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] 9.553976e+12 > sqrtn <- sqrt(length(x)) > (armse <- sd(x) / sqrtn) [1] 2.963775e+13 > (armose <- arm / armse) [1] 0.3223583 > (geo <- geomean(x)) [1] NaN Warning message: In log(x) : NaNs produced > (har <- harmean(x)) [1] -40.61959 > (qua <- quamean(x)) [1] 4.097133e+14 > (win <- winmean(x)) [,1] [,2] [1,] 9.660989e+12 2.961580e+13 [2,] 9.751184e+12 2.945775e+13 [3,] 1.013370e+13 2.933085e+13 [4,] 1.063783e+13 2.915663e+13 [5,] 1.063947e+13 2.911355e+13 [6,] 9.937685e+12 2.896246e+13 [7,] 8.999289e+12 2.879778e+13 [8,] 9.350245e+12 2.871563e+13 [9,] 9.041812e+12 2.861337e+13 [10,] 7.853480e+12 2.832006e+13 [11,] 7.366820e+12 2.775874e+13 [12,] 5.836763e+12 2.713977e+13 [13,] 6.690288e+12 2.693233e+13 [14,] 3.832664e+12 2.588164e+13 [15,] -7.011365e+11 2.527851e+13 [16,] -4.484165e+11 2.492651e+13 [17,] 1.081318e+12 2.463125e+13 [18,] 2.452590e+11 2.450818e+13 [19,] 2.632079e+12 2.406802e+13 [20,] 5.479510e+12 2.353701e+13 [21,] 7.343945e+12 2.260330e+13 [22,] 1.004831e+13 2.197205e+13 [23,] 1.441288e+13 2.061945e+13 [24,] 1.387165e+13 2.053651e+13 [25,] 1.487841e+13 2.030097e+13 [26,] 2.432156e+13 1.908039e+13 [27,] 2.526942e+13 1.892236e+13 [28,] 2.650020e+13 1.878365e+13 [29,] 2.507318e+13 1.857031e+13 [30,] 1.837534e+13 1.776869e+13 [31,] 2.016519e+13 1.744981e+13 [32,] 1.909449e+13 1.717800e+13 [33,] 1.830348e+13 1.697531e+13 [34,] 1.673218e+13 1.656650e+13 [35,] 1.655912e+13 1.639855e+13 [36,] 1.752326e+13 1.597444e+13 [37,] 1.695806e+13 1.526062e+13 [38,] 1.672444e+13 1.513515e+13 [39,] 1.833086e+13 1.493850e+13 [40,] 1.982024e+13 1.474412e+13 [41,] 2.130573e+13 1.457450e+13 [42,] 2.205688e+13 1.447465e+13 [43,] 2.327559e+13 1.435866e+13 [44,] 2.439174e+13 1.411648e+13 [45,] 2.475647e+13 1.391902e+13 [46,] 2.512414e+13 1.381467e+13 [47,] 2.561295e+13 1.356292e+13 [48,] 2.576939e+13 1.348985e+13 [49,] 2.569529e+13 1.333017e+13 [50,] 2.528349e+13 1.329351e+13 [51,] 2.330581e+13 1.297281e+13 [52,] 2.396390e+13 1.263541e+13 [53,] 2.437036e+13 1.260245e+13 [54,] 2.168748e+13 1.226458e+13 [55,] 1.982686e+13 1.206451e+13 [56,] 1.787909e+13 1.188335e+13 [57,] 1.631157e+13 1.174093e+13 [58,] 1.587798e+13 1.170233e+13 [59,] 1.701851e+13 1.150815e+13 [60,] 1.732859e+13 1.147863e+13 [61,] 1.676966e+13 1.139582e+13 [62,] 1.235010e+13 1.102380e+13 [63,] 1.240324e+13 1.099182e+13 [64,] 1.498155e+13 1.076420e+13 > (tri <- trimean(x)) [,1] [,2] [1,] 9.809618e+12 2.907547e+13 [2,] 9.961409e+12 2.849754e+13 [3,] 1.006991e+13 2.796735e+13 [4,] 1.004773e+13 2.744884e+13 [5,] 9.892094e+12 2.694581e+13 [6,] 9.732655e+12 2.641727e+13 [7,] 9.695795e+12 2.588258e+13 [8,] 9.804342e+12 2.533965e+13 [9,] 9.866976e+12 2.476974e+13 [10,] 9.969321e+12 2.417187e+13 [11,] 1.020829e+13 2.357090e+13 [12,] 1.050350e+13 2.300191e+13 [13,] 1.095331e+13 2.247043e+13 [14,] 1.133722e+13 2.192237e+13 [15,] 1.197253e+13 2.145797e+13 [16,] 1.298642e+13 2.102473e+13 [17,] 1.400679e+13 2.059532e+13 [18,] 1.494257e+13 2.016275e+13 [19,] 1.596057e+13 1.970625e+13 [20,] 1.684667e+13 1.925728e+13 [21,] 1.757417e+13 1.882458e+13 [22,] 1.820615e+13 1.844774e+13 [23,] 1.869379e+13 1.809958e+13 [24,] 1.894196e+13 1.784671e+13 [25,] 1.922761e+13 1.757909e+13 [26,] 1.946620e+13 1.730875e+13 [27,] 1.920638e+13 1.711906e+13 [28,] 1.888936e+13 1.692487e+13 [29,] 1.849989e+13 1.672381e+13 [30,] 1.817020e+13 1.652155e+13 [31,] 1.816010e+13 1.636539e+13 [32,] 1.806308e+13 1.621873e+13 [33,] 1.801396e+13 1.607828e+13 [34,] 1.800038e+13 1.593881e+13 [35,] 1.805908e+13 1.581597e+13 [36,] 1.812765e+13 1.569217e+13 [37,] 1.815496e+13 1.558669e+13 [38,] 1.820851e+13 1.552171e+13 [39,] 1.827428e+13 1.545616e+13 [40,] 1.827180e+13 1.539508e+13 [41,] 1.820423e+13 1.533837e+13 [42,] 1.806975e+13 1.528430e+13 [43,] 1.789779e+13 1.522771e+13 [44,] 1.766690e+13 1.516930e+13 [45,] 1.737921e+13 1.511826e+13 [46,] 1.706445e+13 1.507169e+13 [47,] 1.672118e+13 1.502261e+13 [48,] 1.634281e+13 1.498201e+13 [49,] 1.594167e+13 1.493648e+13 [50,] 1.552626e+13 1.489233e+13 [51,] 1.510995e+13 1.483953e+13 [52,] 1.475932e+13 1.480104e+13 [53,] 1.436414e+13 1.477853e+13 [54,] 1.393260e+13 1.474780e+13 [55,] 1.359635e+13 1.473493e+13 [56,] 1.332447e+13 1.472925e+13 [57,] 1.312427e+13 1.472958e+13 [58,] 1.298300e+13 1.473271e+13 [59,] 1.285350e+13 1.472893e+13 [60,] 1.266525e+13 1.473158e+13 [61,] 1.245207e+13 1.472515e+13 [62,] 1.225222e+13 1.471357e+13 [63,] 1.224762e+13 1.472635e+13 [64,] 1.224021e+13 1.472964e+13 > (midr <- midrange(x)) [1] -1.473198e+13 > midm <- array(NA,dim=8) > for (j in 1:8) midm[j] <- midmean(x,j) > midm [1] 1.425715e+13 1.634281e+13 1.425715e+13 1.634281e+13 1.634281e+13 [6] 1.425715e+13 1.634281e+13 1.672118e+13 > postscript(file="/var/www/html/rcomp/tmp/13hkh1227966225.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/297b01227966225.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/34miz1227966225.tab") > > system("convert tmp/13hkh1227966225.ps tmp/13hkh1227966225.png") > system("convert tmp/297b01227966225.ps tmp/297b01227966225.png") > > > proc.time() user system elapsed 1.071 0.350 1.172