Estimated heterogeneity of phenotypic variance of test-day yield with a structural variance model

First-lactation test-day milk, fat, and protein yields from New York, Wisconsin, and California herds from 1990 through 2000 were adjusted additively for age and lactation stage. A random regression model with third-order Legendre polynomials for permanent environmental and genetic effects was used. The model included a random effect with the same polynomial regressions for 2 yr of calvings within herd (herd-time effect) to provide herd-specific lactation curves that can change every 2 yr. (Co)variance components were estimated using expectation-maximization REML simultaneously with phenotypic variances that were modeled using a structural variance model. Maximum heritability for test-day milk yield was estimated to be approximately 20% around 200 to 250 d in milk; heritabilities were slightly lower for test-day fat and protein yields. Herd-time effects explained 12 to 20% of phenotypic variance and had the greatest impact at start of lactation. Variances of test-day yields increased with time, subclass size, and milking frequency. Test month had limited influence on variance. Variance increased for cows in herds with low and high milk yields and for early and late lactation stages. Repeatabilities of variances observed for a given class of herd, test-day, and milking frequency were 14 to 17% across nested variance subclasses based on lactation stage.