R version 2.13.0 (2011-04-13) Copyright (C) 2011 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i486-pc-linux-gnu (32-bit) 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. > > RC.capture <- function (expression, collapse = NULL) { + resultConn <- textConnection('RC.resultText', open = 'w', local=TRUE) + sink(resultConn) + on.exit(function() { + sink() + close(resultConn) + }) + expression + on.exit(NULL) + sink() + close(resultConn) + return(paste(c(RC.resultText, ''), collapse = collapse, sep = '')) + } > RC.texteval <- function (sourceText, collapse = NULL, echo = TRUE) { + sourceConn <- textConnection(sourceText, open = 'r') + on.exit(close(sourceConn)) + result <- RC.capture(source(file = sourceConn, local = FALSE, echo = echo, print.eval = TRUE), collapse = collapse) + on.exit(NULL) + close(sourceConn) + res <- '' + for(i in 1:length(result)) { + if (result[i]!='') res <- paste(res,result[i],' + ',sep='') + } + return(res) + } > x <- c(38,32,35,33,37,29,31,36,35,38,31,34,35,38,37,33,32,38,38,32,33,31,38,39,32,32,35,37,33,33,28,32,31,37,30,33,31,33,31,33,32,33,32,33,28,35,39,34,38,32,38,30,33,38,32,32,34,34,36,34,28,34,35,35,31,37,35,27,40,37,36,38,39,41,27,30,37,31,31,27,36,38,37,33,34,31,39,34,32,33,36,32,41,28,30,36,35,31,34,36,36,35,37,28,39,32,35,39,35,42,34,33,41,33,34,32,40,40,35,36,37,27,39,38,31,33,32,39,36,33,33,32,37,30,38,29,22,35,35,34,35,34,34,35,23,31,27,36,31,32,39,37,38,39,34,31,32,37,36,32,35,36) > par3 = '1e-08' > par2 = '80' > par1 = '1' > par1 <- as.numeric(par1) > par2 <- as.numeric(par2) > par3 <- as.numeric(par3) > stem(x,par1,par2,par3) The decimal point is at the | 22 | 0 23 | 0 24 | 25 | 26 | 27 | 00000 28 | 00000 29 | 00 30 | 00000 31 | 000000000000000 32 | 00000000000000000000 33 | 000000000000000000 34 | 000000000000000 35 | 000000000000000000 36 | 0000000000000 37 | 0000000000000 38 | 00000000000000 39 | 0000000000 40 | 000 41 | 000 42 | 0 > > #Note: the /var/wessaorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/wessaorg/rcomp/createtable") > > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'Stem-and-leaf Plot',1,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,paste('
',RC.texteval('stem(x,par1,par2,par3)'),'
',sep='')) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/wessaorg/rcomp/tmp/1bmjl1322680281.tab") > > > > proc.time() user system elapsed 0.226 0.024 0.243