Assessment of heterogeneity of residual variances using changepoint techniques

• Main Authors: Toro Miguel A, Carabaño Maria J, Rekaya Romdhane
• Format: Article
• Language: deu
• Published: BMC 2000-07-01
• Series: Genetics Selection Evolution
• Subjects: changepoint; heterogeneity; residual variance
• Online Access: http://www.gsejournal.org/content/32/4/383

<p>Abstract</p> <p>Several studies using test-day models show clear heterogeneity of residual variance along lactation. A changepoint technique to account for this heterogeneity is proposed. The data set included 100 744 test-day records of 10 869 Holstein-Friesian cows from northern Spain. A three-stage hierarchical model using the Wood lactation function was employed. Two unknown changepoints at times <it>T</it>₁and <it>T</it>₂, (0 <<it>T</it>₁<<it>T</it>₂<<it>t</it>_max), with continuity of residual variance at these points, were assumed. Also, a nonlinear relationship between residual variance and the number of days of milking <it>t </it>was postulated. The residual variance at a time <it>t</it>(<inline-formula><graphic file="1297-9686-32-4-383-i1.gif"/></inline-formula>) in the lactation phase <it>i </it>was modeled as: <inline-formula><graphic file="1297-9686-32-4-383-i2.gif"/></inline-formula> for (<it>i </it>= 1, 2, 3), where <it>λ</it>_{<it>ι </it>}is a phase-specific parameter. A Bayesian analysis using Gibbs sampling and the Metropolis-Hastings algorithm for marginalization was implemented. After a burn-in of 20 000 iterations, 40 000 samples were drawn to estimate posterior features. The posterior modes of <it>T</it>₁, <it>T</it>₂, <it>λ</it>₁, <it>λ</it>₂, <it>λ</it>₃, <inline-formula><graphic file="1297-9686-32-4-383-i3.gif"/></inline-formula>, <inline-formula><graphic file="1297-9686-32-4-383-i4.gif"/></inline-formula>, <inline-formula><graphic file="1297-9686-32-4-383-i5.gif"/></inline-formula> were 53.2 and 248.2 days; 0.575, -0.406, 0.797 and 0.702, 34.63 and 0.0455 kg², respectively. The residual variance predicted using these point estimates were 2.64, 6.88, 3.59 and 4.35 kg²at days of milking 10, 53, 248 and 305, respectively. This technique requires less restrictive assumptions and the model has fewer parameters than other methods proposed to account for the heterogeneity of residual variance during lactation.</p>