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]]

