# NOT RUN {
data(holland.arthropods)
dat <- holland.arthropods
# use log count to make it possible to have same scale for insects
if(require(reshape2) & require(lattice)){
grays <- colorRampPalette(c("#d9d9d9","#252525"))
dat2 <- melt(dat, id.var=c('row','col'))
contourplot(log(value) ~ col*row|variable, dat2,
col.regions=grays(7), region=TRUE,
main="holland.arthropods - log counts in winter wheat")
}
# }
# NOT RUN {
# individual species
require(lattice)
grays <- colorRampPalette(c("#d9d9d9","#252525"))
contourplot(linyphiidae ~ col*row, dat, at=c(0,40,80,120,160,200), region=TRUE,
col.regions=grays(5),
main="holland.arthropods - linyphiidae counts in winter wheat")
contourplot(n.brevicollis ~ col*row, dat, region=TRUE)
contourplot(linyphiidae~ col*row, dat, region=TRUE)
contourplot(collembola ~ col*row, dat, region=TRUE)
contourplot(carabidae ~ col*row, dat, region=TRUE)
contourplot(lycosidae ~ col*row, dat, region=TRUE)
contourplot(weedcover ~ col*row, dat, region=TRUE)
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab