dat <- minnesota.barley.yield
datw <- minnesota.barley.weather
# Weather trends over time
require("latticeExtra")
useOuterStrips(xyplot(cdd~mo|year*site, datw, groups=year,
subset=(mo > 3 & mo < 10), scales=list(alternating=FALSE),
type='l', auto.key=list(columns=5)))
# Total cooling/heating/precip in Apr-Aug for each site/yr
ww <- subset(datw, mo>=4 & mo<=8)
ww <- aggregate(cbind(cdd,hdd,precip)~site+year, data=ww, sum)
# Average yield per each site/env
yy <- aggregate(yield~site+year, dat, mean)
minn <- merge(ww, yy)
# Higher yields generally associated with cooler temps, more precip
splom(~minn[,c('cdd','hdd','precip','yield')]|site, data=minn,
groups=year, auto.key=list(columns=5), varname.cex=.8,
pscales=2, axis.text.cex=.5, axis.line.tck=0.2)Run the code above in your browser using DataLab