J. Dairy Sci. 2008. 91:4006-4012. doi:10.3168/jds.2007-0816
© 2008 American Dairy Science Association ®
Genetic Relationship Between Culling, Milk Production, Fertility, and Health Traits in Norwegian Red Cows
M. Holtsmark*,1,
B. Heringstad*,
,
P. Madsen
and
J. Ødegård*,
* Department of Animal and Aquacultural Sciences, Norwegian University of Life Sciences, PO Box 5003, N-1432 Ås, Norway
Geno Breeding and AI Association, PO Box 5003, N-1432 Ås, Norway
Department of Breeding and Genetics, Faculty of Agricultural Sciences, University of Århus, Research Centre Foulum, DK-8830 Tjele, Denmark
AKVAFORSK, Institute of Aquaculture Research, PO Box 5010, N-1432 Ås, Norway
1 Corresponding author: marte.holtsmark{at}umb.no
 |
ABSTRACT
|
|---|
First-lactation records on 836,452 daughters of 3,064 Norwegian Red sires were used to examine associations between culling in first lactation and 305-d protein yield, susceptibility to clinical mastitis, lactation mean somatic cell score (SCS), nonreturn rate within 56 d in heifers and primiparous cows, and interval from calving to first insemination. A Bayesian multivariate threshold-linear model was used for analysis. Posterior mean of heritability of liability to culling of primiparous cows was 0.04. The posterior means of the genetic correlations between culling and the other traits were –0.41 to 305-d protein yield, 0.20 to lactation mean SCS, 0.36 to clinical mastitis, 0.15 to interval from calving to first insemination, –0.11 to 56-d nonreturn as heifer, and –0.04 to 56-d nonreturn as primiparous cow. As much as 66% of the genetic variation in culling was explained by genetic variation in protein yield, clinical mastitis, interval of calving to first insemination, and 56-d nonreturn in heifers, whereas contribution from the SCS and 56-d nonreturn as primiparous cow was negligible, after taking the other traits into account. This implies that for breeds selected for a broad breeding goal, including functional traits such as health and fertility, most of the genetic variation in culling will probably be covered by other traits in the breeding goal. However, in populations where data on health and fertility is scarce or not available at all, selection against early culling may be useful in indirect selection for improved health and fertility. Regression of average sire posterior mean on birth-year of the sire indicate a genetic change equivalent to an annual decrease of the probability of culling in first-lactation Norwegian Red cattle by 0.2 percentage units. This genetic improvement is most likely a result of simultaneous selection for improved milk yield, health, and fertility over the last decades.
Key Words: genetic correlation culling fertility health
 |
INTRODUCTION
|
|---|
Globally, there is an increasing interest in including longevity in dairy cattle breeding programs. Different measures of longevity have been included into the routine breeding value evaluations in many countries [e.g., Australia, Canada, Denmark, France, and United States (Interbull, 2007)]. Longevity is an important economic factor in dairy farms. The length of the life of a cow is to a large extent decided by the individual farmer through voluntary culling of low-producing animals. However, diseases, injuries, and poor fertility may also force the farmer to cull a cow. This involuntary culling will often be related to decreased animal welfare and increased costs. Decreasing involuntary culling, especially in the first lactation(s), allows the farmer to increase the voluntary culling, such that profit is maximized (van Arendonk, 1985). Hence, the aim should be to reduce involuntary culling, rather than to increase herd life, per se. One approach is to examine the functional herd life of the cows, in which the herd life is adjusted for milk production, which is considered to be the main cause of voluntary culling (Short and Lawlor, 1992; Jairath et al., 1998; Chirinos et al., 2007). Alternatively, one may select directly to improve the performance for the underlying traits causing involuntary culling.
Studies have shown positive genetic correlation between risk of culling and clinical mastitis (CM; Heringstad et al., 2003a), SCC, interval between calving and first insemination, and number of inseminations (Roxström and Strandberg, 2002). Haile-Mariam et al. (2003, 2004) found negative genetic correlations between survival in early lactations and calving interval (also shown by Olori et al., 2002), SCC, and interval from calving to first service and positive genetic correlations to nonreturn after first insemination and 6-wk pregnancy rate. Although knowledge of the genetic parameters of early culling and its genetic relationship with important production and functional traits is crucial for successful implementation of survival traits in the breeding objective, as well as for prediction of expected genetic improvement, few studies have been done to quantify these relationships. The aims of this study were to estimate genetic parameters for culling in first lactation, defined as a binary trait, and its genetic correlations to 305-d protein yield, clinical mastitis, SCS, 56-d nonreturn rate in heifers and primiparous cows, and interval between calving and first insemination. Other objectives were to examine to what extent genetic variation in risk of early culling can be explained by other traits and how early culling has changed genetically over time.
 |
MATERIALS AND METHODS
|
|---|
Data
Health, fertility, production, and culling data were extracted from the Norwegian Dairy Herd Recording System. The data set included 836,452 daughters of 3,064 Norwegian Red sires. Only first-crop daughters were included (difference between birth year of sire and birth year of daughter <6 yr). First calving had to be between 1980 and 2004 and age at first calving between 20 and 40 mo. The data set was further restricted to cows from herd-5-yr periods with at least 5 primiparous cows.
Nonreturn rate within 56 d in heifers (NR0) and primiparous cows (NR1) was scored as 0 = return or 1 = nonreturn based on whether or not the cow had a second insemination from 6 to 56 d after the first insemination. A new insemination 0 to 5 d after an initial one was considered to be a double insemination. Clinical mastitis in the interval from 30 d before to 300 d after first calving was scored as 1 = treated at least once and 0 = no treatment. Interval from calving to first insemination (CFI) had to be between 20 and 210 d. Lactation mean SCS (LSCS) and 305-d protein yield (PY305) records were calculated for cows having at least 2 test-day records. For cows with incomplete lactations, PY305 was predicted using the same procedure as used for genetic evaluation of Norwegian Red (E. Sehested, Geno Breeding and AI Association, Norway; personal communication). Culling from the herd in first lactation (CULL) was scored as 0 for cows with a recorded second calving and for cows culled later than 365 d after first calving, and 1 otherwise. Data were extracted more than 1 yr after the last calving, such that all cows had the opportunity to survive up to 365 d after their first calving (98% of the cows in the data set had a culling record). Using these definitions, a total of 28.5% of the cows were culled during first lactation, 19% of the cows were treated for clinical mastitis, the overall nonreturn rate was 75% for heifers and 67% for primiparous cows, and 8% of the heifers and 10% of the cows had double inseminations. The number of records for each trait and pair-wise trait combinations are given in Table 1
. The sire pedigree file had information on sires and maternal grand sires of the 3,064 sires with daughters in the data set and included a total of 3,756 males.
View this table:
[in this window]
[in a new window]
|
Table 1. Number of cows with records of the different trait combinations: 305-d protein yield (PY305), interval from calving to first insemination (CFI), lactation mean SCS (LSCS), culling in first lactation (CULL), clinical mastitis (CM), and 56-d nonreturn for heifers (NR0) and primiparous cows (NR1)
|
|
Model
Binary traits were analyzed using threshold models (Falconer and Mackay, 1996) postulating an underlying continuous variable, the liability (
), such that the observed binary response (y) takes value 0 if
is smaller than a threshold value (T) and 1 otherwise:

The 4 binary traits (NR0, NR1, CM, and CULL) and the 3 continuous traits (PY305, CFI, and LSCS), which were assumed to follow a Gaussian distribution, were analyzed using a multivariate linear-threshold sire model. In matrix notation, the model can be written as:
where Y is a vector of liabilities for the binary traits and phenotypic observations for the continuous traits;
is a vector of trait-specific fixed effects, where βNR0 includes effects of double insemination (2 classes), age at first insemination (21 classes), and month-year of first insemination (316 classes) for heifers; βNR1 contains effects of double insemination, age at first calving (21 classes), and month-year of first insemination after first calving (297 classes) for cows; whereas βCFI, βCM, βLSCS, βPY305, and βCULL include effects of age at first calving (21 classes) and month-year of first calving (288 classes); h is a vector of random herd-5-yr effects
where hNR0 contains the random effects of the herd-5-yr period of first insemination as heifer, where hNR1, hCFI, hCM, hLSCS, hPY305, and hCULL contain the random herd-5-yr effects of first calving;
is a vector of sire effects; e is a vector of random residuals; and X, Z1, and Z2 are the corresponding incidence matrices.
Residuals were assumed independent between cows and correlated across traits within cow, following a multivariate normal distribution, e ~N(0, RI), where R is the residual (co)variance matrix for the 7 traits. For all binary traits, residual variances (in the liability scale) were set equal to 1 (Sorensen and Gianola, 2002). The residual (co)variance matrix was sampled using the approach of Korsgaard et al. (2005), combining the findings of Korsgaard et al. (2003) and González (2004), imposing restrictions on the residual variances of the liabilities for the binary traits.
Prior Distributions.
Independent bounded uniform priors were assumed for each of the fixed effects. Multivariate normal prior distributions were assigned for sire and herd-5-yr effects:
respectively, where G0 and H0 are the 7 x 7 (co)variance matrices for sire and herd-5-yr effects, A is the additive genetic relationship matrix of the 3,756 sires, and I is an identity matrix with rank equal to number of herd-5-yr classes (69,798). The matrices G0, H0, and R were required to be positive definite.
Sampling and Convergence Diagnostics.
The data were analyzed using Gibbs sampling in the RJMC routine of the DMU software (Madsen and Jensen, 2006). One chain with a total of 40,000 rounds, with samples stored at every 10th round, was evaluated with the Raftery and Lewis convergence diagnostics (Raftery and Lewis, 1992) using the BOA software package (Smith, 2005). To estimate the 2.5% quantile of each parameter to within ± 0.005 with 95% probability, a chain of 290,000 iterations (saving every 10th round) was required. To achieve this, the original chain was prolonged to 85,000 samples (after discarding the first 20,000 as burn-in). In addition, 3 independent chains, each consisting of 76,000 samples, were branched out from the original chain. Independence of the chains was assured by changing random seed numbers for the random number generators and by excluding the first 1,000 rounds of each chain as an additional burn-in. Posterior means were calculated across all chains. Within each round of sampling, heritabilities were calculated as h2 = 4
2S / (
2S +
2E), where
2S and
2E are the estimated sire and the residual variances, respectively.
Sire Evaluation and Genetic Change
The posterior means of sire-transmitting abilities for CULL were transformed from the liability scale to the probability scale using:

where pi is the expected probability of culling in first lactation for daughters of sire i, µ is the overall mean liability of CULL, and
i is the posterior mean of sire i. Genetic change for CULL was evaluated by plotting average sire posterior means against birth-year of the sire. The average annual change was estimated as the slope of the corresponding linear regression.
Amount of Genetic Variance for Culling Explained by Other Traits
Health, fertility, and production traits influence culling decisions. The relative amount of genetic variation in CULL that can be explained by genetic variation in the other traits included in our study was examined through Cholesky decomposition of the posterior mean of the genetic correlation matrix (with CULL as the last trait in the matrix). The proportion of free genetic variation in liability of CULL is then the square of the last diagonal element of the resulting matrix L (L772) (i.e., the proportion of genetic variance not explained by other traits).
 |
RESULTS AND DISCUSSION
|
|---|
Heritability
The posterior distributions of the heritabilities were all sharp and symmetric (results not shown). The posterior mean (SD) of heritability of liability to CULL was 0.04 (0.002; Table 2
). This point estimate was somewhat greater than the threshold model estimate of Heringstad et al. (2003a) of 0.02. Applying a linear sire model, Nielsen et al. (1999) found a heritability of 0.04 to 0.07 for survival from first calving to the end of the second lactation (binary trait). Also applying linear models, Jairath et al. (1998) estimated a heritability of 0.03 for functional survival in the first lactation, whereas Olori et al. (2002) estimated the heritability for survival in first lactation, preadjusted for milk yield, to be 0.01. Because heritability estimates from linear models are frequency dependent and may vary with the trait definition, the latter estimates are not directly comparable to the estimates obtained in our study. The heritabilities for the other traits (Table 2
) showed good agreement with estimates from several previous studies on Norwegian Red (with partly overlapping data sets, e.g., Heringstad et al. (2003b); Ødegård et al. (2004); Andersen-Ranberg et al. (2005a, b); Andersen-Ranberg and Heringstad (2006).
View this table:
[in this window]
[in a new window]
|
Table 2. Posterior mean (standard deviation) of within-herd heritability (h2), sire variance (on the diagonal), sire covariance (above the diagonal), and genetic correlations (below the diagonal) between 305-d protein yield (PY305), interval from calving to first insemination (CFI), lactation mean SCS (LSCS), culling in first lactation (CULL), clinical mastitis (CM), and 56-d nonreturn for heifers (NR0) and primiparous cows (NR1)
|
|
Genetic Correlations
As shown in Table 2
, the posterior means of the genetic correlations between CULL and the other traits were greatest for PY305 and CM, –0.41 and 0.36, respectively, followed by LSCS, CFI, NR0, and NR1. The genetic correlation between NR1 and CULL was not significantly different from 0.0 (23.7% of the samples were equal or greater than 0.0). Table 2
shows a weak negative correlation between CULL and NR0 and a weak positive genetic correlation between CULL and CFI. The latter is in agreement with Haile-Mariam et al. (2004), who found a genetic correlation of –0.10 between survival in first lactation (the opposite of culling) and interval from calving to first service. The weak correlation between CULL and CFI may be explained by farmers deciding to cull cows that come very late into heat after first calving, rather than inseminating them. This would result in missing data for CFI and NR1 and cause selection in data for interval from calving to first heat, which is not accounted for as long as CFI is missing. Hence, the cows with the longest interval from calving to first heat (not recorded) are likely to be overrepresented among the cows with missing observations and may potentially cause some bias in estimates of some genetic correlations.
Table 2
shows a moderate, negative genetic correlation between CULL and PY305 of –0.41 and moderate positive genetic correlations of CULL with CM of 0.36 and LSCS of 0.20. Heringstad et al. (2003a) also found a moderate genetic correlation of 0.53 between CULL and CM when analyzing a smaller data set from the same breed. Nielsen et al. (1999) reported genetic correlations ranging from –0.37 to –0.75 between survival (the opposite of culling) and udder disease in first lactation for Danish Red, Danish Holstein, and Danish Jersey. Short and Lawlor (1992) found a genetic correlation between 305-d milk yield in first lactation (a trait strongly correlated with PY305) and occurrence of a second calving (the opposite of CULL) of 0.33, which is similar to our estimated correlation between CULL and PY305.
The moderate to small, unfavorable genetic correlations between PY305 and the health and fertility traits (Table 2
) were in agreement with previous studies (Carlén et al., 2004; Heringstad et al., 2005; Andersen-Ranberg and Heringstad, 2006). The estimated genetic correlations within and between the fertility and health traits (Table 2
) were all similar to estimates from previous studies in Norwegian Red (Ødegård et al., 2004; Andersen-Ranberg et al., 2005a,b; Andersen-Ranberg and Heringstad, 2006).
Herd 5-Yr and Residual Correlations
The herd-5-yr and residual correlations are given in Tables 3
and 4
, respectively. Most herd-5-yr correlations were small, except for NR0 and NR1, CFI and NR1, and CFI and NR0, which were moderate. These correlations show that, as expected, herd-management effect fertility. Residual correlations between traits were generally close to 0, except the correlations between PY305 and CULL and CM and CULL, which were small.
View this table:
[in this window]
[in a new window]
|
Table 3. Posterior mean (standard deviation) of herd-5-yr variance (on the diagonal), herd-5-yr covariance (above the diagonal), and herd-5-yr correlations (below the diagonal) between 305-d protein yield (PY305), interval from calving to first insemination (CFI), lactation mean SCS (LSCS), culling in the first lactation (CULL), clinical mastitis (CM), and 56-d nonreturn for heifers (NR0) and primiparous cows (NR1)
|
|
View this table:
[in this window]
[in a new window]
|
Table 4. Posterior mean (standard deviation) of residual variance (on the diagonal), residual covariance (above the diagonal), and residual correlations (below the diagonal) between 305-d protein yield (PY305), interval from calving to first insemination (CFI), lactation mean SCS (LSCS), culling in first lactation (CULL), clinical mastitis (CM), and 56-d nonreturn for heifers (NR0) and primiparous cows (NR1)
|
|
Genetic Change in Culling of Primiparous Cows
Average sire posterior mean for probability of culling by birth-year of sire is given in Figure 1
. The figure shows a steady decrease of the sire posterior means for probability of culling of primiparous cows by 0.2% units per year, which is likely to be a result of genetic improvement of milk yield, health, and fertility in Norwegian Red as a result of selection for a broad breeding goal (Geno, 2007). However, as observed by van Arendonk (1985), the reduction of involuntary culling may not increase the productive herd life of the cows, because the farmer may increase the rate of voluntary culling. This may especially occur when the carcass value approaches the cost of raising a replacement heifer (Allaire and Cunningham, 1980; van Arendonk, 1985). Norwegian Red is selected for both milk and beef production, and because most Norwegian dairy farmers use the Norwegian Red both for milk and beef production, this affects the rate of voluntary culling.
Genetic Variation in Culling Can Be Explained by Other Traits
A Cholesky decomposition of the genetic correlation matrix (Table 2
) revealed that 66% of the genetic variation for CULL was explained by genetic variation in PY305, CM, CFI, and NR0. The traits NR1 and LSCS had no significant effect on genetic variation in liability of CULL that was not already explained by the other traits. Examining first-lactation health and survival data in Danish dairy breeds, Nielsen et al. (1999) found that survival had moderate to high genetic correlations with udder diseases, ranging from –0.37 to –0.75, and feet and leg diseases, ranging from –0.42 to –0.43. It was concluded that although the frequency of feet and leg diseases in the material was low, these diseases tended to cripple the cows and therefore often resulted in culling. Although the major part of the genetic variation in liability to culling in Norwegian Red is explained by traits that are included in the current study, there are also several other traits that may be of some importance for genetic variation in liability to early culling (e.g., other diseases, conformation traits, and temperament) that are currently included in the breeding program for Norwegian Red (Geno, 2007).
Generally, including information on a correlated trait will usually improve accuracy of selection. However, this is not the case if the new trait is a function of the formerly included traits (e.g., protein percentage in addition to milk and protein yield). For early culling, phenotypes of the underlying traits are often direct causes of both voluntary and involuntary culling (e.g., poor milk yield, health, fertility, or a combination of those). Thus, for breeding programs in which the most important of these factors are directly recorded and selected for, there is most likely little to gain by including data on early culling in genetic evaluation. However, all costs associated with each underlying trait, including the risk of involuntary culling, should be taken into account in calculation of their economical values. For dairy cattle populations in which registration of health or fertility traits, or both, are less comprehensive, information on early culling, which may be rather easy to obtain, could be useful in indirect selection for improvement of the underlying traits.
 |
CONCLUSIONS
|
|---|
As much as 66% of the genetic variation in culling in primiparous cows can be explained by genetic variation in protein yield, clinical mastitis, and fertility (CFI, and nonreturn in heifers). Further research is needed to examine whether there will be any additional value of including early culling in selection programs that already have direct selection on health, fertility, and production traits. However, for populations in which direct information on health and fertility traits is scarce, early culling may be a useful tool for indirect selection for improved health and fertility.
 |
ACKNOWLEDGEMENTS
|
|---|
Access to the data was granted by the Norwegian Dairy Herd Recording System in agreement no. 005/2006. This study is a part of the project "Selection for reduced risk of early culling in dairy cattle: Comparison of alternative approaches to improve animal welfare" (project no. 173276/I10), financed by the Research Council of Norway.
Received for publication October 31, 2007.
Accepted for publication April 28, 2008.
 |
REFERENCES
|
|---|
Allaire, F. R., and E. P. Cunningham. 1980. Culling on low milk yield and its economic consequences for the dairy herd. Livest. Prod. Sci. 7:349–359.[CrossRef]
Andersen-Ranberg, I. M., and B. Heringstad. 2006. Genetic associations between female fertility, mastitis and protein yield in Norwegian Red. Proc. 8th World Congr. Genet. Appl. Livest. Prod., Belo Horizonte, Brazil. Commun. No. 1–20.
Andersen-Ranberg, I. M., B. Heringstad, D. Gianola, Y. M. Chang, and G. Klemetsdal. 2005a. Comparison between bivariate models for 56-day nonreturn and interval from calving to first lactation in Norwegian Red. J. Dairy Sci. 88:2190–2198.[Abstract/Free Full Text]
Andersen-Ranberg, I. M., G. Klemetsdal, B. Heringstad, and T. Steine. 2005b. Heritabilities, genetic correlations, and genetic change for female fertility and protein yield in Norwegian dairy cattle. J. Dairy Sci. 88:348–355.[Abstract/Free Full Text]
Carlén, E., E. Strandberg, and A. Roth. 2004. Genetic parameters for clinical mastitis, SCS, and production in the first three lactations of Swedish Holstein cows. J. Dairy Sci. 87:3062–3070.[Abstract/Free Full Text]
Chirinos, Z., M. J. Carabaño, and D. Hernández. 2007. Genetic evaluation of length of productive life in the Spanish Holstein Friesian population. Model validation and genetic parameters estimation. Livest. Prod. Sci. 106:120–131.[CrossRef]
Falconer, D. S., and T. F. C. Mackay. 1996. Introduction to Quantitative Genetics. 4th ed. Addison Wesley Longman Ltd., Essex, UK.
Geno. 2007. Net merit index. http://www.geno.no/ Accessed July 2007.
González, L. G. 2004. Data augmentation in the Bayesian multivariate probit model. Discussion paper SERP2004001 (PDF), University of Sheffield. http://www.shef.ac.uk/content/1/c6/06/31/58/SERP2004001.pdf Accessed July 2007.
Haile-Mariam, M., P. J. Bowman, and M. E. Goddard. 2003. Genetic and environmental relationship among calving interval, survival, persistency of milk yield and SCC in dairy cattle. Livest. Prod. Sci. 80:189–200.[CrossRef]
Haile-Mariam, M., P. J. Bowman, and M. E. Goddard. 2004. Genetic parameters of fertility traits and their correlations with production, type, workability, liveweight, survival index and cell count. Aust. J. Agric. Res. 55:77–87.[CrossRef]
Heringstad, B., Y. M. Chang, D. Gianola, and G. Klemetsdal. 2005. Genetic associations between susceptibility to clinical mastitis and protein yield in Norwegian dairy cattle. J. Dairy Sci. 88:1509–1514.[Abstract/Free Full Text]
Heringstad, B., R. Rekaya, D. Gianola, G. Klemetsdal, and K. A. Weigel. 2003a. Bivariate analysis of liability to clinical mastitis and to culling in first lactation cows. J. Dairy Sci. 86:653–660.[Abstract/Free Full Text]
Heringstad, B., R. Rekaya, D. Gianola, G. Klemetsdal, and K. A. Weigel. 2003b. Genetic change for clinical mastitis in Norwegian Cattle: A threshold model analysis. J. Dairy Sci. 86:369–375.[Abstract/Free Full Text]
Interbull. 2007. Description of National Genetic Evaluation Systems for dairy cattle traits as applied in different Interbull member countries. http://www-interbull.slu.se/national_ges_info2/framesida-ges.htm Accessed July 2007.
Jairath, L., J. C. M. Dekkers, L. R. Schaeffer, Z. Liu, E. B. Burnside, and B. Kolstad. 1998. Genetic evaluation for herd life in Canada. J. Dairy Sci. 81:550–562.[Abstract]
Korsgaard, I. R., A. H. Andersen, P. Madsen, and J. Ødegård. 2005. Another useful reparameterisation to obtain samples from conditional inverse Wishart distributions. I: Proc. EAAP, Book of abstracts 11. Wageningen Pers, Wageningen, the Netherlands.
Korsgaard, I. R., M. Sandø Lund, D. Sorensen, D. Gianola, P. Madsen, and J. Jensen. 2003. Multivariate Bayesian analysis of Gaussian, right censored Gaussian, ordered categorical and binary traits using Gibbs sampling. Genet. Sel. Evol. 35:159–183.[CrossRef][Medline]
Madsen, P., and J. Jensen. 2006. A Users Guide to DMU. A Package for Analysing Multivariate Mixed Models. Version 6, release 4.6. Danish Institute of Agricultural Sciences, Tjele, Denmark.
Nielsen, U. S., G. A. Pedersen, J. Pedersen, and J. Jensen. 1999. Genetic variation in disease traits and their relationship with survival in Danish dairy cattle. Interbull Bull. 21:170–178.
Ødegård, J., B. Heringstad, and G. Klemetsdal. 2004. Short communication: Bivariate genetic analysis of clinical mastitis and SCC in Norwegian dairy cattle. J. Dairy Sci. 87:3515–3517.[Abstract/Free Full Text]
Olori, V. E., T. H. E. Meuwissen, and R. F. Veerkamp. 2002. Calving interval and survival breeding values as measure of cow fertility in a pasture-based production system with seasonal calving. J. Dairy Sci. 85:689–696.[Abstract]
Raftery, A. L., and S. Lewis. 1992. How many iterations in the Gibbs sampler? Pages 763–774 in Bayesian Statistics 4. J. M. Bernardo, J. O. Berger, A. P. Dawid, and A. F. M. Smith, ed. Oxford University Press, Oxford, UK.
Roxström, A., and E. Strandberg. 2002. Genetic analysis of functional fertility-, mastitis-, and production-determined length of productive life in Swedish dairy cattle. Livest. Prod. Sci. 74:125–135.[CrossRef]
Short, T. H., and T. J. Lawlor. 1992. Genetic parameters of conformation traits, milk yield, and herd life in Holsteins. J. Dairy Sci. 75:1987–1998.[Abstract]
Smith, B. J. 2005. Bayesian output analysis program (BOA), version 1.1. Users manual. http://www.public-health.uiowa.edu/boa/BOA.pdf Accessed July 2007.
Sorensen, D., and D. Gianola. 2002. Likelihood, Bayesian, and MCMC methods in quantitative genetics. Springer-Verlag, New York, NY.
van Arendonk, J. A. M. 1985. Studies on the replacement policies in dairy cattle. II. Optimum policy and influence of changes in production and prices. Livest. Prod. Sci. 13:101–121.[CrossRef]
This article has been cited by other articles:

|
 |

|
 |
 
M. Holtsmark, B. Heringstad, and J. Odegard
Predictive abilities of different statistical models for analysis of survival data in dairy cattle
J Dairy Sci,
November 1, 2009;
92(11):
5730 - 5738.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
M. Holtsmark, B. Heringstad, P. Madsen, and J. Odegard
Short Communication: Use of Culling Information in Genetic Evaluation of Fertility and Mastitis in Norwegian Red Cows
J Dairy Sci,
November 1, 2008;
91(11):
4429 - 4432.
[Abstract]
[Full Text]
[PDF]
|
 |
|