O objetivo deste trabalho foi estimar funções de covariância para a produção de leite no dia do controle (PLDC). Foram analisados 27 mil registros de PLDC de 3.362 primeiras lactações de vacas da raça Gir leiteira, paridas entre 1990 e 2007. As PLDC foram agrupadas em vinte classes quinzenais, analisadas por modelos de regressão aleatória, cujos efeitos aleatórios, genético-aditivo e de ambiente permanente foram modelados pelas funções de Wilmink (W) ou Ali & Schaeffer (AS). A modelagem da variância residual (VR) foi ajustada por meio de 1, 4, 6 ou 10 classes. As estimativas de herdabilidade para a PLDC variaram de 0,12 a 0,32, para a função AS, e de 0,09 a 0,33, para a função W, e foram maiores ao início da lactação. As correlações genéticas entre as PLDC decresceram de valores próximos à unidade, entre controles adjacentes, para valores negativos entre as PLDC da primeira e duas últimas quinzenas da lactação. O modelo que empregou a função AS com quatro classes de VR é uma opção parcimoniosa para o ajuste das PLDC de vacas Gir leiteira no Brasil.
bovinos de leite; função paramétrica; inferência bayesiana; parâmetros genéticos
