# NOT RUN {
data(ferritin)
temp <- ferritin[ferritin$period <4,]
plot(temp$old.lot, temp$new.lot, type='n', log='xy',
xlab="Old lot", ylab="New Lot")
text(temp$old.lot, temp$new.lot, temp$period,
col=temp$period)
# }
Run the code above in your browser using DataLab