library(faraway)
data(corrosion)
corrosion[c(1:3,12:13),]
library(lattice)
xyplot(loss~Fe, data=corrosion)