R version 2.10.0 (2009-10-26) 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.0192508842078761 + ,-0.0243384070148590 + ,0.0158952358252125 + ,-0.0260722519603377 + ,0.00884040098416066 + ,-0.0212744425298722 + ,0.0374645687523337 + ,0.000287491611075654 + ,-0.00339379103432857 + ,0.0171818334669544 + ,-0.033426965424278 + ,0.00217934347065340 + ,-0.0274115854113506 + ,-0.0569170628639589 + ,-0.000209124152662078 + ,-0.00400646079770227 + ,-0.00748529583959452 + ,-0.00585685997080524 + ,0.0752980051262616 + ,0.00262906294324769 + ,-0.00639835830921532 + ,0.0257097451976159 + ,0.0327139661508631 + ,0.0363171189336884 + ,-0.00129610688362125 + ,-0.0219713395125436 + ,-0.0176410998712386 + ,0.0411130767250018 + ,0.0124816435524092 + ,0.020793906880278 + ,-0.0484360434735618 + ,0.0263730615102121 + ,-0.0200050532987077 + ,-0.0134631900868672 + ,-0.013113310509618 + ,0.00219519903716346 + ,0.00349361484962575 + ,0.0314382751789365 + ,0.0140064877187480 + ,-0.0167139079425391 + ,0.0122645988561186 + ,0.00114204009422627 + ,-0.0409497518101819 + ,-0.00376200146549827 + ,0.000868011842158073 + ,-0.0322586305360591 + ,0.0682477804173895 + ,-0.0118187390548836 + ,-0.0359119197592944 + ,-0.00236630329867407 + ,0.0668542802302875 + ,0.000748549879677433 + ,-0.0097707458515968 + ,-0.0168718298288839 + ,-0.00438484896246825 + ,0.00552586381473103 + ,0.0219829902429555 + ,-0.00573500441714712 + ,0.0238920416800163 + ,-0.0221259220432080 + ,-0.0310115591186089 + ,0.00225903855869600 + ,-0.00289136404543944 + ,0.00486404360736588 + ,-0.0379857416078679 + ,-0.0192308004977295 + ,-0.0604836676986188 + ,-0.0523951106616251 + ,-0.0264459462048099 + ,0.0190059857944224 + ,0.051579422018824 + ,0.0105776720459108 + ,-0.00355602540245227 + ,-0.0163280881061046 + ,-0.0107869871318581 + ,-0.00460914546299842 + ,-0.0694566072240325 + ,0.0275743671281808 + ,0.103706651222640 + ,-0.0682976491656814 + ,0.00956417905378418 + ,0.00777995668698742 + ,-0.000471725526274496 + ,-0.0170756133168365 + ,-0.0219398274893341 + ,0.0266415497232234 + ,-0.00718744004864717 + ,0.0255235333060244 + ,-0.0689794396176606 + ,-0.00269863966055686 + ,-0.00134991577133056 + ,-0.05643065977133 + ,-0.0254379905104139 + ,0.165898547230958 + ,0.0537206816137991 + ,0.0642377469045096 + ,-0.00592693533019978 + ,0.0194027109676582 + ,-0.0100689167955129 + ,0.0582650246376592 + ,0.0165272380190866 + ,0.0724558629310153 + ,0.0306429617300824 + ,-0.0387863822572402 + ,0.0310497861315557 + ,-0.0857565364396057 + ,0.00579522547072407 + ,-0.0384897916145452 + ,0.0360585980557015 + ,-0.00320594887898824 + ,0.0137113324902744 + ,-0.0240598741976394 + ,0.0184468748372574 + ,0.00325157937184180 + ,0.0254779586910709 + ,-0.0283919147746814 + ,-0.0227883523855816 + ,-0.0430964564458684 + ,-0.0216370051000690 + ,-0.0140707505733256 + ,-0.0154464375494437 + ,-0.0406067597697643 + ,-0.0137719609823467 + ,-0.0143160597857540 + ,-0.0155991144608775 + ,-0.0225547231065048 + ,0.0228793700809542 + ,-0.0265463902544472 + ,0.0287954700277272 + ,-0.0300888958159422 + ,0.00711730170513459 + ,-0.0388479619658106 + ,-0.0210517741393705 + ,0.0209016660704892 + ,-0.00736398444160426 + ,-0.00695634652202435 + ,-0.0215118032211888 + ,0.0130815818198251 + ,-0.0294445694816839 + ,-0.0117484040652358 + ,0.0237400494796378 + ,-0.000850108661828754 + ,-0.0633816370082633 + ,0.0161387564701593 + ,0.0039191450328204 + ,-0.0300291383913910 + ,-0.00743127215616451 + ,-0.0087841632926546 + ,0.0442649141203364 + ,-0.0331697190770205 + ,-0.0340606317626467 + ,-0.0469800100804107 + ,0.00632880283139771 + ,-0.0186025106792085 + ,0.0183753561468092 + ,0.0199570007781874 + ,-0.0224484565035172 + ,0.0120439983018481 + ,0.0276616630681987 + ,0.0595111245787349 + ,-0.0756624215189736 + ,0.0273649927990675 + ,0.0215180661735496 + ,0.0274133738668157 + ,-0.00463098036463181 + ,0.0393431417763206 + ,-0.00890522610843237 + ,0.0664244228906462 + ,-0.0242474922714065 + ,0.0613863253342411 + ,0.0261613141046088 + ,-0.0436635790987991 + ,0.0076741977290358 + ,0.0201218132900302) > 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.0006930345 > sqrtn <- sqrt(length(x)) > (armse <- sd(x) / sqrtn) [1] 0.002578591 > (armose <- arm / armse) [1] -0.2687648 > (geo <- geomean(x)) [1] NaN Warning message: In log(x) : NaNs produced > (har <- harmean(x)) [1] -0.0410296 > (qua <- quamean(x)) [1] 0.03392316 > (win <- winmean(x)) [,1] [,2] [1,] -0.000992447 0.0024566386 [2,] -0.001247652 0.0023736183 [3,] -0.001288427 0.0023632670 [4,] -0.001369492 0.0023436447 [5,] -0.001268270 0.0023141778 [6,] -0.001183163 0.0022964445 [7,] -0.001127648 0.0022606607 [8,] -0.001236385 0.0022359004 [9,] -0.001124643 0.0021912756 [10,] -0.000968725 0.0021499316 [11,] -0.001163963 0.0020935774 [12,] -0.001082917 0.0020431340 [13,] -0.001587032 0.0019602567 [14,] -0.001667904 0.0019057859 [15,] -0.001790917 0.0018827271 [16,] -0.001801931 0.0018421350 [17,] -0.001908021 0.0018277512 [18,] -0.001904083 0.0018209524 [19,] -0.002214262 0.0017717133 [20,] -0.002122522 0.0017277363 [21,] -0.001945978 0.0016976215 [22,] -0.001917296 0.0016831405 [23,] -0.002127501 0.0016525900 [24,] -0.002158221 0.0016222607 [25,] -0.001991587 0.0016019341 [26,] -0.001877774 0.0015850550 [27,] -0.001876009 0.0015832998 [28,] -0.001898356 0.0015613038 [29,] -0.001767662 0.0015388545 [30,] -0.001635148 0.0015182467 [31,] -0.001561456 0.0014947360 [32,] -0.001577229 0.0014893452 [33,] -0.001514999 0.0014816201 [34,] -0.001700955 0.0014374993 [35,] -0.001510348 0.0014134658 [36,] -0.001669609 0.0013934451 [37,] -0.001820323 0.0013704305 [38,] -0.001644170 0.0013344956 [39,] -0.001729963 0.0013159042 [40,] -0.001730306 0.0013112100 [41,] -0.001812673 0.0012886226 [42,] -0.001815143 0.0012813154 [43,] -0.001944335 0.0012672097 [44,] -0.001968081 0.0012504886 [45,] -0.002080298 0.0012335785 [46,] -0.002036455 0.0012260568 [47,] -0.002298697 0.0011898181 [48,] -0.002190525 0.0011466305 [49,] -0.002087544 0.0011173447 [50,] -0.002151750 0.0011102452 [51,] -0.002521195 0.0010424526 [52,] -0.002322084 0.0010089381 [53,] -0.002341658 0.0009762967 [54,] -0.002464603 0.0009538436 [55,] -0.002483291 0.0009432808 [56,] -0.002430117 0.0009260652 [57,] -0.002671664 0.0008623176 [58,] -0.002958602 0.0008281231 > (tri <- trimean(x)) [,1] [,2] [1,] -0.001167035 0.0023694316 [2,] -0.001345731 0.0022740500 [3,] -0.001396522 0.0022182860 [4,] -0.001434290 0.0021620706 [5,] -0.001451477 0.0021073036 [6,] -0.001490832 0.0020555908 [7,] -0.001546597 0.0020034870 [8,] -0.001612508 0.0019540908 [9,] -0.001664948 0.0019049658 [10,] -0.001732779 0.0018590279 [11,] -0.001820243 0.0018154652 [12,] -0.001889451 0.0017763409 [13,] -0.001968469 0.0017406538 [14,] -0.002003438 0.0017125610 [15,] -0.002032398 0.0016885623 [16,] -0.002052124 0.0016651099 [17,] -0.002071559 0.0016440316 [18,] -0.002083688 0.0016225882 [19,] -0.002096454 0.0015998930 [20,] -0.002088403 0.0015800827 [21,] -0.002086154 0.0015626967 [22,] -0.002095089 0.0015464785 [23,] -0.002106074 0.0015299981 [24,] -0.002104788 0.0015146632 [25,] -0.002101664 0.0015004682 [26,] -0.002107943 0.0014865602 [27,] -0.002120780 0.0014726144 [28,] -0.002134148 0.0014572771 [29,] -0.002146779 0.0014421707 [30,] -0.002166733 0.0014273008 [31,] -0.002194261 0.0014124829 [32,] -0.002226551 0.0013979435 [33,] -0.002259243 0.0013821116 [34,] -0.002296264 0.0013650130 [35,] -0.002325558 0.0013497930 [36,] -0.002365290 0.0013347095 [37,] -0.002398915 0.0013194832 [38,] -0.002426680 0.0013043368 [39,] -0.002464003 0.0012903328 [40,] -0.002498843 0.0012760539 [41,] -0.002535182 0.0012601265 [42,] -0.002569251 0.0012440528 [43,] -0.002604753 0.0012262813 [44,] -0.002635827 0.0012072838 [45,] -0.002667264 0.0011871056 [46,] -0.002694942 0.0011656075 [47,] -0.002726077 0.0011414833 [48,] -0.002746361 0.0011176744 [49,] -0.002772873 0.0010949072 [50,] -0.002805760 0.0010716665 [51,] -0.002805760 0.0010451887 [52,] -0.002852781 0.0010227411 [53,] -0.002878895 0.0010003364 [54,] -0.002905619 0.0009779307 [55,] -0.002927823 0.0009541681 [56,] -0.002950506 0.0009270001 [57,] -0.002977455 0.0008966581 [58,] -0.002993549 0.0008701339 > (midr <- midrange(x)) [1] 0.04007101 > midm <- array(NA,dim=8) > for (j in 1:8) midm[j] <- midmean(x,j) > midm [1] -0.002857713 -0.002604753 -0.002604753 -0.002604753 -0.002635827 [6] -0.002604753 -0.002604753 -0.002604753 > postscript(file="/var/www/rcomp/tmp/1xh3q1261081372.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/rcomp/tmp/2460w1261081372.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/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/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/rcomp/tmp/39o9v1261081372.tab") > > try(system("convert tmp/1xh3q1261081372.ps tmp/1xh3q1261081372.png",intern=TRUE)) character(0) > try(system("convert tmp/2460w1261081372.ps tmp/2460w1261081372.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.260 0.700 1.809