Hi all, I suppose this is a very simple question, but as I've lost already a bit of time with it, without being able to get what I wanted, I'm addressing the question to the group in the hope someone can help me.
I pretend to plot the richness of herbaceous species (RichHN) as a function of time since remnant isolation (Isol) conditioned to the area of the remnant (fArea - this is a factor with 3 classes). I also want to differentiate the sites plotted using a grouping binary variable (Urb). So far, I've written this: extinction<-xyplot(RichHN~Isol|fArea, data = Data, type = "p", groups=Urb, aspect="xy", layout=c(3,1), between=list(x=0.5, y=0.5), cex=1, pch = c(21, 1), auto.key=list(space="bottom"), main="Influence of time since remnant isolation and remnant area on species diversity", xlab="Time since remnant isolation (rank)", ylab="Richness of native herbaceous") print(extinction) What do I need to add, in order I can have: - the graph dots labeled with the variable SiteName - an overall trend line for each panel - the dots with Urb=1 filled in black; and the dots with Urb=0 filled in white I've been back and forwards, trying to use panel function, but I didn't manage to do it right... Thank you for your help in advance. Regards, C. ------------------------------------ Cristina Estima Ramalho PhD Candidate MSc in GIS Ecosystem Restoration Lab School of Plant Biology (M090) The University of Western Australia 35 Stirling Highway, Crawley WA 6009 Perth, Australia Phone: +61 (08) 6488 4655; Fax: +61 (08) 6488 7461 www.plants.uwa.edu.au/research/ecosystem_restoration --------------------------------------------------------------------------------- [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.