Statistical analysis of somatic cell scores via mixed model methodology for longitudinal data