if (FALSE) {
library(agridat)
data(holland.arthropods)
dat <- holland.arthropods
# use log count to make it possible to have same scale for insects
libs(reshape2, 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")
if(0){
# individual species
libs(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)
}
}
Run the code above in your browser using DataLab