if (FALSE) {
library(agridat)
data(ratkowsky.onions)
dat <- ratkowsky.onions
# Model inverse yield as a quadratic. Could be better...
libs(lattice)
dat <- transform(dat, iyield = 1/yield)
m1 <- lm(iyield ~ I(density^2)*loc, dat)
dat$pred <- predict(m1)
libs(latticeExtra)
foo <- xyplot(iyield ~ density, data=dat, group=loc, auto.key=TRUE,
main="ratkowski.onions",ylab="Inverse yield")
foo + xyplot(pred ~ density, data=dat, group=loc, type='l')
}
Run the code above in your browser using DataLab