[Rd] Re: [R] How to vectorize

[Gabor Grothendieck]

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 at gmail.com> wrote:
> On 6/6/05, PANTERA Laurent <laurent.pantera at 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))
>