library(drfit) data(antifoul) rantifoul <- drfit(antifoul) rantifoul