I\'m trying to fit a Richards function to this dataset:
structure(list(volume = c(50, 110, 183, 264, 357, 457, 562, 664, 761, 851,