Genetic evaluation using random regression models with different covariance functions for test‐day milk yield in an admixture population of Thailand goats

The objectives of this study were to compare covariance functions (CF) and estimate the heritability of milk yield from test‐day records among exotic (Saanen, Anglo‐Nubian, Toggenburg and Alpine) and crossbred goats (Thai native and exotic breed), using a random regression model. A total of 1472 records of test‐day milk yield were used, collected from 112 does between 2003 and 2006. CF of the study were Wilmink function, second‐ and third‐order Legendre polynomials, and linear splines 4 knots located at 5, 25, 90 and 155 days in milk (SP25–90) and 5, 35, 95 and 155 of days in milk (SP35–95). Variance components were estimated by restricted maximum likelihood method (REML). Goodness of fit, Akaike information criterion (AIC), percentage of squared bias (PSB), mean square error (MSE), and empirical correlation (RHO) between the observed and predicted values were used to compare models. The results showed that CF had an impact on (co)variance estimation in random regression models (RRM). The RRM with splines 4 knots located at 5, 25, 90 and 155 of days in milk had the lowest AIC, PSB and MSE, and the highest RHO. The heritability estimated throughout lactation obtained with this model ranged from 0.13 to 0.23.