[Rd] Re: [R] How to vectorize

From: Gabor Grothendieck <ggrothendieck_at_gmail.com>
Date: Wed 08 Jun 2005 - 13:20:34 GMT

To r-devel:

The discussion below on r-help brought out the following:

> plot(1:10)
> text(5,5,lab=expression(italic(22*"33")))

has the effect of italicizing 33 (which is a character string) but not 22 (which is not). Is this intended behavior? It seems strange.

I am using "R version 2.1.0, 2005-05-14" on Windows XP.

On 6/6/05, Gabor Grothendieck <ggrothendieck@gmail.com> wrote:
> On 6/6/05, PANTERA Laurent <laurent.pantera@irsn.fr> wrote:
> > Dear R-List,
> >
> > I would like to write nicely the names of some isotopes on a plot. The
> > code bellow works fine.
> >
> > plot(1:10,1:10)
> > text(c(2,4,8),c(2,4,8),labels=c(expression(italic(phantom(0)^{78}*Ge)),
> > expression(italic(phantom(0)^{137}*Cs)),
> > expression(italic(phantom(0)^{129*m}*Te))),
> > cex=3
> > )
> >
> > But, since I have a lot of isotopes to write on the plot, I would like
> > to construct automatically the labels. So I wrote the code below which
> > works fine.
> >
> > listenoms <- list(nom=c("Ge","Cs","Te"),num=c("78","137","129*m"))
> > n <- length(listenoms$nom)
> > resu <- "c("
> > for( i in 1:(n-1))
> > {
> > resu <- paste(resu,paste("expression(italic(phantom(0)^{",
> > listenoms$num[i],"}*",
> > listenoms$nom[i],")),",sep=""))
> > }
> > resu <- paste(resu,paste("expression(italic(phantom(0)^{",
> > listenoms$num[n],"}*",
> > listenoms$nom[n],")))",sep=""))
> > plot(1:10,1:10)
> > text(c(2,4,8),c(2,4,8),labels=eval(parse(text=resu)),cex=2)
> >
> > I assume there is a better way to do that using vectorization.
> > May you help me to find it ?
> >
>
> Assuming your setup:
>
> x <- c(2,4,8)
> nom <- c("Ge","Cs","Te")
> num <- c("78","137","129*m")
> plot(1:10)
>
> # the ith label can be generated via
> lab <- function(i)
> as.expression(bquote(italic(phantom(0)^{.(num[i])}*.(nom[i]))))
>
> # giving the following vectorized solution
> labs <- lapply(seq(x), lab)
> text(x,x,do.call(c,labs))
>
> # although a 'for' loop is arguably simpler
> plot(1:10)
> for(i in seq(x)) text(x[i], x[i], lab(i))
>



R-devel@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel Received on Wed Jun 08 23:22:47 2005

This archive was generated by hypermail 2.1.8 : Mon 20 Feb 2006 - 03:21:08 GMT