In these models, lactation number and DIM were fitted as fixed effects, and cow number was fitted as a random effect. The coefficient of determination values were estimated from pseudo coefficient of determination values using the square of the correlation between fitted values and observed values.