Skip to main content

A comparison of methods to calculate a total merit index using stochastic simulation

Abstract

Background

Modern dairy cattle breeding goals include several production and more and more functional traits. Estimated breeding values (EBV) that are combined in the total merit index usually come from single-trait models or from multivariate models for groups of traits. In most cases, a multivariate animal model based on phenotypic data for all traits is not feasible and approximate methods based on selection index theory are applied to derive the total merit index. Therefore, the objective of this study was to compare a full multitrait animal model with two approximate multitrait models and a selection index approach based on simulated data.

Methods

Three production and two functional traits were simulated to mimic the national Austrian Brown Swiss population. The reference method for derivation of the total merit index was a multitrait evaluation based on all phenotypic data. Two of the approximate methods were variations of an approximate multitrait model that used either yield deviations or de-regressed breeding values. The final method was an adaptation of the selection index method that is used in routine evaluations in Austria and Germany. Three scenarios with respect to residual covariances were set up: residual covariances were equal to zero, or half of or equal to the genetic covariances.

Results

Results of both approximate multitrait models were very close to those of the reference method, with rank correlations of 1. Both methods were nearly unbiased. Rank correlations for the selection index method showed good results when residual covariances were zero but correlations with the reference method decreased when residual covariances were large. Furthermore, EBV were biased when residual covariances were high.

Conclusions

We applied an approximate multitrait two-step procedure to yield deviations and de-regressed breeding values, which led to nearly unbiased results. De-regressed breeding values gave even slightly better results. Our results confirmed that ignoring residual covariances when a selection index approach is applied leads to remarkable bias. This could be relevant in terms of selection accuracy. Our findings suggest that the approximate multitrait approach applied to de-regressed breeding values can be used in routine genetic evaluation.

Background

In dairy cattle breeding programs, selection is focused on different production traits and more and more on functional traits. Usually, estimated breeding values (EBV) or phenotypes [1] are combined into a total merit index (TMI) [2]. Traits or groups of traits are usually evaluated separately based on different statistical models [3]. This is also the case in the joint Austrian-German genetic evaluation of one dairy and several dual-purpose cattle breeds. The TMI and several sub-indices are based on a selection index method [1] for all cattle breeds, which was proposed by [4]. Currently, the TMI of Fleckvieh (dual-purpose Simmental) and Brown Swiss populations consists of up to 30 production and functional traits. EBV for the TMI as well as for several sub-indices are estimated either using univariate (e.g. protein yield) or multivariate (e.g. calving ease and stillbirth) methods by applying animal or sire-maternal-grandsire models (the latter for functional longevity only). Some of these models include repeated measures, such as somatic cell count [2]. Subsequently, EBV for individual traits are combined to form TMI or other sub-indices by assuming that residual covariances between traits or groups of traits are equal to zero. Due to the large number of traits involved and to methodological constraints, the additive genetic (co)variance matrix to combine EBV into a TMI cannot be estimated multivariately for all traits. Additive genetic correlations between the traits in the TMI are mainly obtained from the literature [4]. Furthermore, we have observed that the TMI shows an upward trend in bias for animals with low reliabilities (r2) (unpublished results). We hypothesize that this is due to ignoring residual covariances when combining traits or trait groups into the TMI. Full multivariate estimation of all traits based on phenotypic data could represent the optimum methodology [5-7], but in routine genetic evaluations substantial restrictions (e.g. computer power and computational considerations) make this approach infeasible. Thus, an approximate multivariate model using a two-step procedure was proposed and validated using simulated data [3,8,9]. In the first step, genetic (co)variances were estimated based on yield deviations (YD). The associated weights for YD were calculated from univariate analyses and YD were adjusted for all environmental effects. In the second step, a multivariate animal model that included a random genetic effect for the animal and a fixed year effect was applied [3,9]. Due to improved connectedness of the data and simultaneous estimation of genetically correlated traits, accuracies of EBV can be increased compared to those from a univariate analysis. In particular, functional traits, which are characterized by low heritabilities (h2), benefit more from a multivariate animal model than moderately to highly heritable traits [6,10]. A set of characteristics using an approximate multivariate model based on YD exists and is well described [3,7,11,12]. However, YD cannot always be obtained in routine evaluation settings, e.g. for EBV from Interbull. Therefore, the objective of this study was to compare four methods to calculate TMI using: a full multitrait model; approximate multitrait models applied to YD and de-regressed estimated breeding values (drEBV); and the current selection index method which is used in the joint Austrian-German genetic evaluations. The consequences of ignoring residual covariances when a TMI is computed were investigated. Hence, three scenarios with respect to residual covariances were set up using data that were generated by stochastic simulation of data that represent a simplified breeding scheme of the Austrian Brown Swiss cattle population.

Methods

Data simulation

A population structure that roughly reflected the Austrian breeding program of Brown Swiss cattle was simulated with the stochastic simulation program ADAM [13]. The simulated population comprised approximately 51 300 cows in 1710 herds. Five normally distributed traits were chosen to represent dairy, beef and fitness traits. Additional requirements were a wide range of heritabilities and genetic correlations between traits and economic importance. Fat yield (FY), protein yield (PY), somatic cell count (SCC) and non-return rate of cows (NRR) were measured on female animals. Net daily gain (NDG) was observed on approximately 60% of all male animals. No repeated records were assumed. The assumed heritabilities and genetic correlations for the five traits are in Table 1, which were obtained from estimates for the national Brown Swiss population [14]. Around 25% young bulls and 75% proven bulls were used for matings in the selection scheme. Breeding values and phenotypes for the five traits were simulated for animals from a base population. Animals were selected on a TMI based on multivariately EBV over 30 years. Based on values used for routine genetic evaluation, relative economic weights per additive genetic standard deviation were 5.4, 53.6, 4.3, 19.7 and 17% for FY, PY, NDG, SCC and NRR, respectively [14]. Three scenarios with respect to residual covariances between traits were simulated. In scenario 0, residual correlations between traits were assumed to be zero. For scenarios 1 and 2, residual correlations corresponded to 50% or 100% of the genetic correlations, respectively. These scenarios were chosen to specifically test the currently used selection index method for possible bias due to ignoring correlations between residual effects, as proposed by [4]. For each scenario, ten replicates were simulated. On average, breeding values and reliabilities of EBV were estimated for 755 567 animals for each scenario.

Table 1 Genetic parameters used for simulation (heritabilities on the diagonal, genetic correlations above the diagonal)

Methods to calculate the TMI

The first method (MULTI) was a full multivariate animal model based on phenotypic data using the true genetic and phenotypic parameters, assuming the following linear model:

$$ \mathbf{y} = \mathbf{X}\mathbf{b} + \mathbf{Z}\mathbf{a} + \mathbf{e}, $$
(1)

where y is the vector representing simulated phenotypic observations of FY, PY, NDG, SCC and NRR, respectively; X and Z represent the incidence matrices for fixed and random effects, respectively, b is the vector of the fixed herd-year effects, a is the vector of the random animal additive genetic effects and e is the vector of the random residual effects. Vector a was assumed to have a multivariate normal distribution, with MVN(0, G = G 0 ⊗ A), where G 0 is a 5 × 5 additive genetic variance-covariance matrix,⊗ is the Kronecker product of matrices, and A represents the numerator relationship matrix. Residuals e were assumed to be MVN(0, R = R 0 ⊗ I), where R 0 is the 5 × 5 residual variance-covariance matrix and I represents the identity matrix. Subsequently, the TMI was calculated as:

$$ {\mathrm{TMI}}_{\mathrm{MULTI}}={\widehat{\mathrm{a}}}_{\mathrm{FY}}{\upomega}_{\mathrm{FY}}+{\widehat{\mathrm{a}}}_{\mathrm{PY}}{\upomega}_{\mathrm{PY}}+{\widehat{\mathrm{a}}}_{\mathrm{NDG}}{\upomega}_{\mathrm{NDG}}+{\widehat{\mathrm{a}}}_{\mathrm{SCC}}{\upomega}_{\mathrm{SCC}}+{\widehat{\mathrm{a}}}_{\mathrm{NRR}}{\upomega}_{\mathrm{NRR}}, $$
(2)

where â represents the EBV for the traits and ω denotes the relative economic values. As the full multitrait estimation, MULTI, represents the optimum method, this was considered to be the reference method.

The second (YD) and third (DRP) methods were based on the approximate multitrait two-step procedure proposed by [3]. These methods use either YD or drEBV. For both methods, univariate genetic evaluations were used to calculate YD and drEBV for each trait. For the YD method, phenotypic observations were corrected for the fixed herd-year effect using the following model:

$$ \mathbf{y}* = \mathbf{y}\ \hbox{--}\ \mathbf{X}\mathbf{b}, $$
(3)

where y* is the vector of yield deviations, y is the vector of phenotypic observations of the traits FY, PY, NDG, SCC and NRR, and b is the vector of herd-year effects. X is the incidence matrix of proper order. After correcting for fixed effects, all five traits were analysed together in the following multivariate animal model:

$$ \mathbf{y}* = \mathbf{X}\mathbf{b} + \mathbf{Z}\mathbf{a} + \mathbf{e}, $$
(4)

where y* indicates yield deviations of animals of a trait; X and Z represent the incidence matrices for fixed and random effects, respectively, b is the vector of year of birth fixed effects, a is the vector of random animal additive genetic effects and e denotes the vector of random residual effects. The year of birth effect was included in order to account for over- or underestimation of the genetic trend [9]. The use of YD requires weights to consider the different amount of information for each animal. Reliabilities were calculated with the program package ApaX [15], using the approximate Interbull method of [15]. Based on these reliabilities, effective own performances (EOP) were calculated and used as weighting factors for yield deviations in the multivariate estimation of breeding values for method YD. The following formula was applied to calculate EOP for trait i:

$$ {\mathrm{EOP}}_{\mathrm{i},\mathrm{j}}=\frac{\upalpha_{\mathrm{i}}}{1\hbox{-} {\mathrm{r}}_{\mathrm{i},\mathrm{j}}^2}-{\upalpha}_{\mathrm{i},} $$
(5)

where αi is the ratio of residual and additive genetic variances of trait i; and r2 represents the reliability of the own performance of animal j for trait i [16]. The assumption and size of the (co)variance matrices G 0 and R 0 for estimating breeding values and the equation for the TMI were based on the same formulas as in method MULTI. For method DRP, EBV were de-regressed using a univariate de-regression based on the approach of [17] and [18], which is implemented in the program package MiX99 [19]. The de-regression procedure uses the EBV and their respective effective daughter contributions as weights, with the general mean as the only fixed effect. Based on model Equation (4), EBV were estimated in a 5-trait animal model, with drEBV as the response variables and an overall mean, a year effect, additive genetic effects of animals for each trait and residual effects of drEBV as explanatory variables. As described for method YD, EOP were used as weighting factors to estimate breeding values.

The fourth method SI is the approach that is currently used in routine genetic evaluations in Austria and Germany to calculate the total merit index based on selection index theory. The EBV of the five traits were estimated univariately using the same model Equation (1) but assuming that G = A \( {\upsigma}_{\mathrm{a}}^2 \) and R = I \( {\upsigma}_{\mathrm{e}}^2 \), where \( {\upsigma}_{\mathrm{a}}^2 \) is the additive genetic variance and \( {\upsigma}_{\mathrm{e}}^2 \) is the residual variance. In order to obtain the TMI with method SI, covariances between EBV (σij) were calculated as:

$$ {\upsigma}_{\mathrm{i}\mathrm{j}} = {\mathrm{r}}_{\mathrm{gij}}{{\mathrm{r}}^2}_{\mathrm{i}}{{\mathrm{r}}^2}_{\mathrm{j}}{\upsigma}_{\mathrm{ai}}{\upsigma}_{\mathrm{aj}}, $$
(6)

where rgij is the genetic correlation between traits i and j; r2 i,j are the reliabilities of EBV of traits i and j, and σai,aj are the additive genetic standard deviations of traits i and j. In this equation, residual correlations are neglected, i.e. assumed to be zero.

YD, drEBV and EBV were computed using the program package MiX99 [19]. Reliabilities of EBV were calculated with the program package ApaX [15]. Reliabilities of the TMI of methods MULTI, YD and DRP were calculated using the approach of [20]. Reliabilities of the TMI of method SI were calculated using the formula described by [4]. For all methods, genetic parameters were not re-estimated, but the true (simulated) parameters were used. All EBV were standardized to an additive genetic standard deviation of 12 and 100 for the mean of cohorts from birth year 18 to 22. Spearman rank correlations between TMI of different methods were calculated using SAS 9.2. [21]. To estimate breeding values and reliabilities by calculating genetic trends and for all further analyses, the last 20 years of the simulated period were used.

Results

Correlations

Table 2 shows the Spearman rank correlations between TMI obtained using the reference method MULTI and from methods YD, DRP and SI categorized by TMI reliability for all animals for scenarios 0 and 2. Correlations for YD and especially DRP were almost 1. Correlations between the reference method and method SI were slightly lower, particularly for scenario 2, in which the residual (and thus the phenotypic) correlation was equal to the genetic correlation. Results for scenario 1 are not shown since values were always between those of scenarios 0 and 2. The simulated selection program resulted in a strong genetic trend, therefore results were also analyzed within year groups. Across all animals and in the last 20 years of the simulation (Table 3), correlations between TMI from the different methods were high and ranged from 0.983 to 1 for all scenarios. As expected, correlations decreased slightly when split into year groups. Correlations for method SI were more strongly affected, especially when non-zero residual covariances were simulated. In addition, the population was subdivided into the following four groups: bulls with progeny information (BP), bulls without progeny information (BNP), females with progeny information (FP) and females without progeny information (FNP). Correlations between TMI from the reference method MULTI and the three other methods for BP, BNP, FP and FNP across all 20 years and scenarios ranged from 0.984 to 1 (results not shown). For bulls and dams with progeny, a tendency for higher correlations was observed. Method SI showed the lowest correlations and method DRP showed the highest correlations with method MULTI. With regard to scenarios, rank correlations were highest for scenario 0 followed by scenarios 1 and 2. When grouped by reliability, rank correlations of TMI from MULTI with those from YD and DRP were always close to 1. However, correlations between MULTI and SI grouped by reliability were lower, especially when non-zero residual correlations were simulated (scenarios 1 and 2, results not shown). For BNP bulls, correlations even declined to 0.907.

Table 2 Rank correlations of TMI obtained using alternate methods with TMI obtained using the multivariate method within TMI reliability groups
Table 3 Rank correlations with multivariate TMI (MULTI) within year groups for different TMI methods for scenarios 0 and 2

Bias

Bias was defined as the difference between the estimated TMI of each method and the MULTI TMI (e.g. TMIYD-TMIMULTI). This was done for all animals and scenarios. Since the reference method MULTI essentially had no bias, the presented bias is equivalent to the bias from true TMI. Results for scenarios 0 and 2 for all animals grouped by reliability are in Table 4. For all scenarios, methods YD and DRP showed almost no bias. However, method SI showed a relevant bias in both scenarios. The TMI of animals with a reliability lower than 50% were overestimated, whereas that of animals with a reliability greater than 50% were underestimated. When residual correlations were non-zero (scenario 2), this bias was more pronounced. Method SI led to a markedly overestimated genetic trend, which was expressed as a downwards bias during the first years and an upwards bias during the last years (Table 5). This was especially evident for scenario 2 and even more for the top 10% animals based on TMI (Figure 1). When non-zero residual covariances were simulated (scenario 2), the bias was similar to that for scenario 0 but slightly more pronounced. Figure 2 shows the bias for all bulls with progeny performance (BP) and for bulls without progeny performance (BNP) for method SI for scenario 0 (currently used method). A difference of about half a genetic standard deviation was observed between the average TMI of BP and BNP at the end of the simulated period. Selecting the best 10% bulls each year from groups BP and BNP resulted in an even greater difference (Figure 3) for method SI. Differences between groups BP and BNP were even higher when the simulated residual covariances were greater (results are not shown). The difference between FP and FNP was close to 0 for method SI for all scenarios.

Table 4 Bias of approximate TMI methods relative to TMI from multivariate analysis within TMI reliability groups for scenarios 0 and 2
Table 5 Bias with different TMI methods from multivariate TMI within year groups for scenarios 0 and 2
Figure 1
figure 1

Time trend of bias (TMI-TMI MULTI ) with different TMI methods for the top 10% animals per year in scenario 0.

Figure 2
figure 2

Time trend of bias (TMI-TMI MULTI ) for bulls with progeny (BP) and bulls without progeny (BNP) (method SI, scenario 0).

Figure 3
figure 3

Time trend of bias (TMI-TMI MULTI ) for the top 10% bulls with progeny (BP) and bulls without progeny (BNP) by year (method SI, scenario 0).

Reliabilities and standard deviations

The mean reliabilities (r2) and standard deviations (SD) for the TMI of methods MULTI, YD and DRP were almost equal (numerically ±0.1 to ±0.3). Thus, only the results of the reference MULTI and SI methods are presented in Table 6. Estimated reliabilities were slightly higher for method SI than for the other methods for scenario 0, when grouped by years. For scenario 2, for which non-zero residual covariances were simulated, reliabilities for method SI were substantially higher than for MULTI for each year group. Since the estimated reliabilities depend on the method used, realised (true) reliabilities obtained from the squared correlation of estimated true breeding values for TMI are also presented. Realised reliabilities for SI were on average 4% lower than for MULTI.

Table 6 Reliabilities and standard deviations (SD) of TMI from methods MULTI and SI for scenarios 0 and 2 grouped by year and loss in selection response (SR) with SI compared to MULTI

To quantify the presumed inflation of TMI from method SI, the expected SD of TMI was calculated as, \( \sqrt{{\mathrm{r}}^2\upsigma {\scriptscriptstyle \frac{2}{\mathrm{T}}}} \) using realised reliabilities and the true genetic standard deviation\( {\upsigma}_{\mathrm{T}}^2 \)within year (Table 6). The variance for TMI for method MULTI was very close to the expected values. For method SI, the observed SD of the TMI was notably higher than the expected SD for scenario 0 (from 0.4 to 1) and especially for scenario 2 (from 1.6 to 2.1). These findings were similar for all scenarios and groups (results not shown).

Selection response

To demonstrate the practical implications of the lower realised reliabilities and the inflated SD of EBV with method SI, the expected loss in selection response using SI compared to MULTI was analysed. If we consider the breeder’s equation, we can assume that selection intensity, genetic standard deviation and generation interval do not depend on the TMI method. Thus, differences in selection response between the analysed methods depend only on the reliability of the TMI. This led to an expected loss in response of 3 to 5% with method SI compared to MULTI (Table 6). Since the realised reliabilities for methods YD and DRP were almost identical to those for MULTI, no loss in response is expected for these methods.

Discussion

In this study, four methods were used to combine several traits into a total merit index and were compared with different assumptions regarding residual covariances. Both approximate multitrait two-step procedures that used either YD or drEBV led to results that were comparable to those of the full multivariate animal model (MULTI). Our findings agree with those of [8,9] who compared a full multitrait model, an approximate two-step procedure applied to YD and a combination of single trait models. The approximate two-step procedure was not as efficient as the full multitrait model but superior to the single-trait approach in terms of genetic response. Results of this study substantiate some drawbacks of the method that is currently used to calculate a TMI in routine genetic evaluations in Germany and Austria. Spearman rank correlations of TMI between the reference MULTI and the YD and DRP methods were close to 1 for all scenarios for any categories of reliability or year. Correlations for method SI were lower than for methods YD and DRP even when zero residual covariances were simulated. Rank correlations between MULTI and SI decreased with increasing residual covariances. Applying selection index theory (e.g. method SI) to calculate a TMI is valid when traits are not or only slightly correlated [3,12]. However, this is not the case for the routinely calculated TMI in Austria and Germany due to the wide range of dairy, beef and functional traits which are correlated to a certain degree (e.g. genetic correlations between fat kg and fertility, dressing percentage or milkability are -0.20, -0.15 and 0.25, respectively [14]). In particular, some of the lowly heritable functional traits are correlated with production or conformation traits to a considerable degree [12,22]. Furthermore, it is well known that ignoring residual covariances when animals are recorded for different traits in the same environment is not valid [3]. This was confirmed by our results with method SI. As residual covariances increased, correlations of TMI with the reference method decreased and deviated more from the full multivariate model. In addition, a downwards bias was observed for animals with high reliabilities, e.g. bulls with progeny, and an upwards bias for animals with low reliabilities, e.g. young animals without progeny. Hence, the bias for the top 10% animals can be relevant in terms of selection decisions across birth cohorts. The average difference between TMI from the reference MULTI and SI methods can be up to half a genetic SD, which leads to substantial re-ranking of bulls. This can cause selection bias, in particular when early selection decisions are made on young bulls. It should be noted that, in this simulation study, young bulls without progeny cannot be compared with young bulls with genomic EBV from routine evaluations. In the simulation, young bulls had own data only for NDG, while for all other traits a pedigree index was used. This means that the problem of bias may be less severe in routine evaluations.

In the method that is currently used to calculate a TMI, residual correlations are neglected for traits for which covariance is expected to occur, e.g. fat and protein yield, functional longevity and some fitness traits [22,23]. Results of this study based on simulated data imply that the variances of EBV with low reliability are in general inflated using method SI. This is confirmed by results obtained with real data on functional longevity and type traits in Austrian Fleckvieh cows [24]. Reliabilities and variances of the TMI were overestimated for functional longevity when combined with auxiliary type traits using method SI and ignoring residual covariances. Based on the results of this simulation, drEBV could be a good alternative to YD, since they can be easier to obtain in some cases and show equally good results. This could also help to include Interbull EBV in the national evaluation, since individual YD are not available at the international level. The current method for calculating TMI and other sub-indices could thus be replaced by an approximate two-step procedure using drEBV.

Two crucial points need to be clarified before the approximate multitrait approach can be implemented in routine genetic evaluations: (1) accurate genetic and residual (co)variance components must be estimated using a multivariate analysis of all traits in order to allow their inclusion in the two-step procedure; if this is not possible, it has been suggested by [8] to cluster traits that have genetic correlations above 0.10; and (2) furthermore, a genomic evaluation including international EBV (Interbull) must be implemented.

Conclusions

An approximate multitrait two-step procedure to compute TMI applied to drEBV led to nearly unbiased results. Fortunately, the outcomes for the multitrait method based on drEBV are equally as good as those on YD, which will facilitate its implementation, especially for specific traits for which it is difficult to obtain accurate YD. The advantages of these methods are greatest when residual covariances differ from zero. Although there are several crucial prerequisites before implementing an approximate multitrait two-step procedure in routine genetic evaluations, our results open up perspectives for the replacement of the current selection index method by this procedure based on drEBV.

References

  1. Hazel LN, Lush JL. The efficiency of three methods of selection. J Hered. 1942;33:393–9.

    Google Scholar 

  2. Miesenberger J, Fuerst C. Experiences in selecting on a total merit index in the Austrian Fleckvieh breed. Biotech Anim Husbandry. 2006;22:17–27.

    Article  Google Scholar 

  3. Ducrocq V, Boichard D, Barbat A, Larroque H. Implementation of an approximate multitrait BLUP evaluation to combine production traits and functional traits into a total merit index. In Proceedings of the 52nd Annual Meeting of the European Association for Animal Production: 26-29 August 2001; Budapest; 2001.

  4. Miesenberger J. Zuchtzieldefinition und Indexselektion für die österreichische Rinderzucht, PhD thesis. Austria: University of Natural Resources and Life Sciences Vienna; 1997.

    Google Scholar 

  5. Besbes B, Ducrocq V, Protais M. An approximate total merit index combining linear traits, a survival trait and a categorical trait in laying hens. In Proceedings of the 7th World Congress on Genetics Applied to Livestock Production: 19-23 August 2002; Montpellier; CD-ROM communication n° 20-05.

  6. Mrode RA. Linear models for the prediction of animal breeding values. 3rd ed. Wallingford: CABI; 2014.

    Google Scholar 

  7. Tarrés J, Piedrafita J, Ducrocq V. Validation of an approximate approach to compute genetic correlations between longevity and linear traits. Genet Sel Evol. 2006;38:65–83.

    Article  PubMed Central  PubMed  Google Scholar 

  8. Lassen J, Sørensen MK, Madsen P, Ducrocq V. A stochastic simulation study on validation of an approximate multitrait model using preadjusted data for prediction of breeding values. J Dairy Sci. 2007;90:3002–11.

    Article  CAS  PubMed  Google Scholar 

  9. Lassen J, Sørensen MK, Madsen P, Ducrocq V. An approximate multitrait model for genetic evaluation in dairy cattle with a robust estimation of genetic trends. Genet Sel Evol. 2007;39:353–67.

    Article  PubMed Central  PubMed  Google Scholar 

  10. Thompson R, Meyer K. A review of theoretical aspects in the estimation of breeding values for multi-trait selection. Livest Prod Sci. 1986;15:299–313.

    Article  Google Scholar 

  11. Ducrocq V. Multiple trait prediction. Principles and problems. In Proceedings of the 5th World Congress on Genetics Applied to Livestock Production: 6-11 August 1994; Guelph. 1994;18:455-62.

  12. Lassen J. Inference in multivariate mixed models for prediction of breeding values, PhD thesis. Tjele, Denmark: The Royal Veterinary and Agricultural University, Frederiksberg, Denmark; Danish Institute of Agricultural Science; 2006.

    Google Scholar 

  13. Pedersen LD, Sørensen AC, Henryon M, Ansari-Mahyari S, Berg P. ADAM: A computer program to simulate selective breeding schemes for animals. Livest Sci. 2009;121:343–4.

    Article  Google Scholar 

  14. Fuerst C, Dodenhoff J, Egger-Danner C, Emmerling R, Hamann H, Krogmeier D, et al. Zuchtwertschätzung beim Rind - Grundlagen, Methoden und Interpretation. 2014. http://zar.at/Downloads/Zuchtwerte.html. Accessed 12 Jan 2015.

  15. Strandén I, Lidauer M, Mäntysaari EA, Pösö J. Calculation of Interbull weighting factors for the Finnish test day model. In Proceedings of the Interbull technical workshop;22-23 October 2000; Verden. Interbull Bulletin. 2000;26:78–9.

    Google Scholar 

  16. Edel C, Emmerling R, Götz KU. Optimized aggregation for MA-BLUP evaluation in German Fleckvieh. In Proceedings of the Interbull meeting, 21-24 August 2009; Barcelona. Interbull Bulletin. 2009;40:178–83.

    Google Scholar 

  17. Jairath J, Dekkers JCM, Schaeffer L, Liu Z, Burnside EB, Kolstad B. Genetic evaluation for herd life in Canada. J Dairy Sci. 1998;81:550–62.

    Article  CAS  PubMed  Google Scholar 

  18. Schaeffer LR. Multiple trait international bull comparison. Livest Prod Sci. 2001;69:145–53.

    Article  Google Scholar 

  19. Lidauer M, Matilainen K, Mäntysaari E, Pitkänen T, Taskinen M, Strandén I. MiX99. General program for solving large mixed model equations with preconditioned conjugate gradient method. Release VII/2013 beta; Jokioinen, Finland: MTT Agrifood Research Finland; 2013

  20. Tier B, Meyer K. Approximating prediction error covariance among additive genetic effects within animals in multiple-trait and random regression models. J Anim Breed Genet. 2004;121:77–89.

    Article  Google Scholar 

  21. SAS. User´s guide. Cary, NC, USA: SAS Institute Inc.; 2008

  22. Govignon-Gion A, Dassonneville R, Balloche G, Ducrocq V. Genetic evaluation of mastitis in dairy cattle in France. In Proceedings of the Interbull meeting, 28-31 May 2012; Cork. Interbull Bulletin. 2012;46:121–6.

    Google Scholar 

  23. Samoré AB, Schneider MDel P, Canavesi F, Bagnato A, Groen AF. Relationship between somatic cell count and functional longevity assessed using survival analysis in Italian Holstein-Friesian cow. Livest Prod Sci. 2003;80:211–20.

    Article  Google Scholar 

  24. Pfeiffer C, Fuerst-Waltl B, Ducrocq V, Fuerst C. Approximate multivariate genetic evaluation of functional longevity and type traits in Austrian Fleckvieh cattle. In Proceedings of the 10th World Congress on Genetics Applied to Livestock Production: 17-22 August 2014; Vancouver; 2014. https://asas.org/docs/default-source/wcgalp-posters/386_paper_8337_manuscript_96_0.pdf?sfvrsn=2.

Download references

Acknowledgements

The authors gratefully thank Anders Christian Sørensen and his colleagues from Aarhus University for granting access and support with the simulation program ADAM. We also gratefully acknowledge the funding of the project OptiGene (‘Optimization of long-term genetic progress of Austrian cattle breeds with emphasis on health and genomic selection’; project number 100808) by the Austrian Federal Ministry of Agriculture, Forestry, Environment and Water Management, the Federations of Austrian Fleckvieh, Brown-Swiss, Pinzgauer and Tyrolean Grey cattle and the Federation of Austrian Cattle Breeders (ZAR).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Christina Pfeiffer.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

CP, BFW, HS and CF designed the basic concept, including methods and scenarios, of the simulation study. CP, HS and FS worked on the simulation. CP and CF did the analyses. CP, BFW, HS and CF interpreted and discussed the results. CP drafted the manuscript. All authors read and approved the final version of the manuscript.

Rights and permissions

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Pfeiffer, C., Fuerst-Waltl, B., Schwarzenbacher, H. et al. A comparison of methods to calculate a total merit index using stochastic simulation. Genet Sel Evol 47, 36 (2015). https://doi.org/10.1186/s12711-015-0118-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12711-015-0118-4

Keywords