library(drfit)
data(XY)
rXY <- drfit(XY)
rXY