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 ?
Thanks
Laurent
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))
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)) >