[R] conditioning variable in panel.xyplot?

Felix Andrews felix at nfrac.org
Tue Mar 9 01:33:25 CET 2010


Alternatively

library(latticeExtra)

xyplot(y ~ x | site, d) +
  xyplot(y ~ x | site, q, col = "red")



(which is a shortcut for:)

xyplot(y ~ x | site, d) +
  as.layer(xyplot(y ~ x | site, q, col = "red"))


On 9 March 2010 11:17, Seth W Bigelow <sbigelow at fs.fed.us> wrote:
> Ah, wonderful, thank you for the code Deepayan. To recap for posterity: I
> have two datafiles, d and q: each has x-y coordinates that are conditioned
> by site (The actual data, for me, is maps of parent trees and their
> seedlings). I wanted to superimpose the xy plots of d and q, by site,
> without going to the trouble of merging the d & q datasets into a single
> dataset. The solution is to use the "which.packet" statement is
>
>
> d <- data.frame(site  = c(rep("A",12), rep("B",12)),
> x=rnorm(24),y=rnorm(24))                # Create the main xy dataset
> q <- data.frame(site  = c(rep("A",7), rep("B",7)),
> x=rnorm(14),y=rnorm(14))                # Create the alternate xy dataset
>
>
> q.split <- split(q, q$site)                     # Split up the alternate
> dataset by site
>
> mypanel <- function(..., alt.data) {
>    with(alt.data[[ which.packet()[1] ]],                       #
> which.packet passes index of the relevant data subset...
>         panel.xyplot(x = x, y = y, col="red"))                 # ... to
> panel.xyplot()
>    panel.xyplot(...)
> }
>
> xyplot(y ~ x | site, d, alt.data = q.split,             # After providing
> the alternative dataset and the panel...
>       panel = mypanel)                         # ...everything prints out
> properly, like magic!
>
>
>
> Dr. Seth  W. Bigelow
> Biologist, USDA-FS Pacific Southwest Research Station
>        [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at 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.
>



-- 
Felix Andrews / 安福立
Postdoctoral Fellow
Integrated Catchment Assessment and Management (iCAM) Centre
Fenner School of Environment and Society [Bldg 48a]
The Australian National University
Canberra ACT 0200 Australia
M: +61 410 400 963
T: + 61 2 6125 4670
E: felix.andrews at anu.edu.au
CRICOS Provider No. 00120C
-- 
http://www.neurofractal.org/felix/



More information about the R-help mailing list