Comparison of homogeneity and heterogeneity of residual variance using random regression test-day models for first lactation Japanese Holstein cows

Using a large‐scale data set that included first lactation test day records from 1975 to 2000 for Japanese Holsteins, genetic parameters for milk yield were estimated by using random regression (RR) test‐day models (TDM) with heterogeneous and homogeneous residual variances. It is necessary for the RR‐TDM to include a function that explains the shape of the lactation curve. The RR‐TDM with the LW curve, which combined Wilmink's curve and a Legendre polynomial, was used for fitting the model for milk yield. In recent years, increases in residual variance have been noted for Japanese dairy cattle. Thus, three kinds of heterogeneous residual variance over the calving year were considered: H1, H2 and HG. Linear and quadratic exponential functions for the calving year were used in H1 and H2, respectively. Residual variance of HG was divided into five groups according to calving year. Homogeneous residual variance was HO. All heterogeneous residual variances increased with calving year in an almost linear fashion. Residual variance increased over the study period. However, there is no need to consider heterogeneous residual variances in genetic evaluations, because the heterogeneity of residual variance over the years did not affect the ranking of top sires and cows.