O objetivo do presente estudo foi avaliar quatro modelos matemáticos quanto ao seu ajuste a curvas de lactação de vacas Holandesas pertencentes a rebanhos da região sudoeste do estado do Paraná, Brasil. Inicialmente 42.281 controles de produções leiteiras referentes aos anos de 2005 a 2011, foram da "Associação Paranaense de Criadores de Bovinos da Raça Holandesa (APCBRH)". Os dados que não possuíam a data de secagem e a produção total aos 305 dias foram excluídos, restando 15.142 controles referentes a 2.441 vacas Holandesas. Os dados foram classificados de acordo com a ordem de parição (variando de um a seis), e dentro de cada ordem de parição os animais foram separados em quartis (Q25%, Q50%,Q75% e Q100%) conforme a produção total aos 305 dias. Para cada ordem de parição, em cada quartil foram ajustados quatro modelos matemáticos, dois predominantemente empíricos (modelo de Brody e de Wood) e dois com características mais mecanicistas (modelo de Dijkstra e de Pollott). A qualidade de ajuste foi avaliada pelo critério de Akaike corrigido. O modelo de Wood apresentou o melhor ajuste em quase todas as situações avaliadas e, portanto, deve ser considerado como o melhor modelo para descrever, pelo menos empiricamente, as curvas de lactação de vacas Holandesas criadas no Sudoeste do Paraná.
gado leiteiro; produção de leite; seleção de modelos matemáticos; modelos não lineares