dat <- ratkowsky.onions
# Model inverse yield as a quadratic. Could be better...
library(lattice)
dat <- transform(dat, iyield = 1/yield)
xyplot(iyield ~ density, data=dat, group=loc, auto.key=TRUE)
m1 <- lm(iyield ~ I(density^2)*loc, dat)
dat$pred <- predict(m1)
require(latticeExtra)
xyplot(iyield ~ density, data=dat, group=loc, auto.key=TRUE) +
xyplot(pred ~ density, data=dat, group=loc, type='l')
# To do: semiparametric / gam modelRun the code above in your browser using DataLab