(Co)variance components for milk, fat, and protein yields during first and second lactations were estimated from test-day data from 23,029 Holstein cows from 37 herds in Pennsylvania and Wisconsin using a multitrait test-day model. Canonical transformation was used with an expectation-maximization algorithm. To allow description of (co)variances within and across yield traits and parities, four...