Fitting phenological curves with Generalized Linear Mixed Models (GLMMs)