Fitting milk production curves through nonlinear mixed models