data(gomez.seedrate)
dat <- gomez.seedrate
require(lattice)
#xyplot(yield ~ rate, data=dat, group=rep)
bwplot(yield ~ rate, data=dat, horizontal=FALSE)
# Quadratic response
m1 <- lm(yield ~ poly(rate,2), dat)Run the code above in your browser using DataLab