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]]
______________________________________________
[email protected] 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.
