Skip to content

lattice to ggplot2 conversion help

2 messages · ashz, Ista Zahn

ashz
# 
Hi,

I am interested in ggplot2 and I found this lattice code very interesting
(http://addictedtor.free.fr/graphiques/graphcode.php?graph=48).

Code:
library(lattice) 
lattice.options(default.theme = canonical.theme(color = FALSE)) 


tmp <-
    expand.grid(geology = c("Sand","Clay","Silt","Rock"), 
                species = c("ArisDiff", "BracSera", "CynDact",
                            "ElioMuti", "EragCurS", "EragPseu"),
                dist = seq(1,9,1) ) 

tmp$height <- rnorm(216) 


sp <- list(superpose.symbol = list(pch = 1:6, cex = 1.2), 
           superpose.line = list(col = "grey", lty = 1)) 

# print is needed when you source() the file
print(xyplot(height ~ dist | geology, data = tmp, 
       groups = species,
       layout = c(2,2), 
       panel = function(x, y, type, ...) {
         panel.superpose(x, y, type="l", ...)
         lpoints(x, y, pch=16, col="white", cex=2) 
         panel.superpose(x, y, type="p",...)
       },
       par.settings = sp, 
       auto.key = list(columns = 2, lines = TRUE))) 



I will be very happy if someone can please explain me how to do it in
ggplot2 as it will be great help.

Cheers,
Ashz


--
View this message in context: http://r.789695.n4.nabble.com/lattice-to-ggplot2-conversion-help-tp3760001p3760001.html
Sent from the R help mailing list archive at Nabble.com.
Ista Zahn
# 
Hi Ashz,
On Mon, Aug 22, 2011 at 8:42 AM, ashz <ashz at walla.co.il> wrote:
The basic plot can be created with

ggplot(tmp, aes(x = dist, y = height)) +
  geom_point(aes(shape = species)) +
  geom_line(aes(group = species)) +
  facet_wrap(~geology)

and some of the formatting can be reproduced with

ggplot(tmp, aes(x = dist, y = height)) +
  geom_point(aes(shape = species), size = 4) +
  geom_line(aes(group = species), color = "gray60") +
  facet_wrap(~geology) +
  theme_bw() +
  opts(legend.position = "top", legend.direction = "horizontal")

I do have to say that I think this graph is  a mess though. Too many
jumbled points and lines to easily make sense of it. I would go with
small multiples all the way:

ggplot(tmp, aes(x = dist, y = height)) +
  geom_point() +
  geom_line() +
  facet_grid(species~geology) +
  theme_bw() +
  opts(legend.position = "top", legend.direction = "horizontal")

as this seems like a much clearer presentation of the data.

Best,
Ista