Foram utilizados 10.238 registros semanais de produção de leite no dia do controle leiteiro provenientes de 388 primeiras lactações de cabras da raça Saanen visando comparar diferentes modelos de regressão aleatória (MRA). Primeiramente, foram comparados cinco modelos, cujos termos exponenciais da função de Wilmink assumiram os seguintes valores -0,0350; -0,0500; -0,0565; -0,0680 e -0,1000 (W0350, W0500, W0565, W0680 e W1000, respectivamente), considerando-se homogeneidade de variância residual ao longo da lactação. No modelo W0500, o valor -0,0500 foi mantido, enquanto nos modelos W0350, W0565, W0680 e W1000 foram usados os valores -0,0350; -0,0565; -0,0680; e -0,1000, respectivamente, em substituição ao valor -0,0500, proposto no modelo original utilizado para bovinos de leite. Escolhido o melhor modelo, segundo o ln L, foram avaliadas, pelos critérios AIC, BIC e ln L, a homogeneidade e heterogeneidade da variância residual: homogeneidade, duas, três, quatro, cinco e seis classes ao longo da lactação. De acordo com o critério usado, o modelo W0350 apresenta o melhor ajuste dentre os avaliados. Com relação à variância residual, a utilização de seis classes ao longo da lactação é indicada pelos critérios AIC, BIC, ln L e teste de razão de verossimilhança. As estimativas de herdabilidade ao longo da lactação, para o melhor modelo, variam de 0,07 (2ª semana de lactação) a 0,25 (20ª semana de lactação).