上の図は例によって library(lattice)
によるもの．
欠測している「年・樹種」をぬくには，
xyplot() の
layout
と
skip
オプションを使う．
概要はこんなかんじで．
v.sp.year <- sort(apply(
expand.grid(levels(d$sp), 2001:2004),
1, function(x) paste(x[1], x[2])
))
print(xyplot(
flower ~ log(dbh) | sp.year,
data = d,
ylim = c(-0.1, 1.1),
skip = !(v.sp.year %in% levels(d$sp.year)),
layout = c(4, length(v.sp)),
panel = function(x, y, subscripts, ...) {
dsub <- d[subscripts,]
panel.points(
x,
y + (dsub$light.level - 2.5) * 0.05,
pch = 3,
col = col.ll[dsub$light.level]
)
}
))
drop.unused.levels
なるオプション
(help(xyplot))
を使ってもよさそうなんだけど，
これはいまいち使いかたがよくわからない．
