Advances in Animal and Veterinary Sciences

Download PDF Download ePUB
 Research Article

Research Article

Advances in Animal and Veterinary Sciences 2 (7): 381 – 389

A Meta–analysis of the Impact of Parity on Dystocia and Stillbirth in Holstein Cattle

Firas Rashad Al–Samarai

    Department of Veterinary Public Health, University of Baghdad, Iraq

*Corresponding author:firas_rashad@yahoo.com

ARTICLE CITATION: Al–Samarai FR (2014). A meta–analysis of the impact of parity on dystocia and stillbirth in Holstein cattle. Adv. Anim. Vet. Sci. 2 (7): 381 – 389.
Received: 2014–04–06, Revised: 2014–08–05, Accepted: 2014–08–07
The electronic version of this article is the complete one and can be found online at ( http://dx.doi.org/10.14737/journal.aavs/2014/2.7.381.389 ) which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited

ABSTRACT

Dystocia and stillbirth are major factors reducing the productivity of dairy cattle. The objective of this study was to determine the effect of parity on the rates of dystocia and stillbirth. A meta–analysis was conducted to investigate the impact of first parity (primiparous) and later parities (multiparous) on dystocia and stillbirth in Holstein cattle. A total of 30 and 19 papers were analyzed for evaluation of two traits. Results revealed that primiparous cattle are more susceptible to dystocia [Odds Ratio (OR) = 2.68, 95% Confidence interval (CI) 2.51 to 2.85], stillbirth (OR = 2.18, 95% CI 1.84 to 2.58) as compared with multiparous. These results supported the opinion about the importance of considering primiparous and multiparous as different traits in genetic evaluation and shed light on the importance of improving genetics and environment of heifers to minimize the effect of dystocia and stillbirth in Holstein cattle.

INTRODUCTION

Meta–analyses can be defined as systematic reviews with pooled data (Ton et al., 2007). It was considered as a valuable method with unique properties: establishing whether scientific findings are consistent (Cook et al., 1998) and can be generalized across populations (Burrin and Britton, 1986), limit bias, improve reliability and accuracy of conclusions (Collett, 1994) and increase the power and precision of treatment effects (Bell and Bauman, 1997).

In recent years, breeders have shown increasing interest in selection of functional traits in dairy cattle (Mark, 2004), therefore they have focused to shift selection from traits that increasing–profit to reducing–costs traits (De Maturana, 2007). Health management has been emphasis in order to minimize losses due health disorders (Beaudea et al., 2000).

Dystocia and stillbirth are related terms; as dystocia associated with approximately 50% of calf mortality cases at birth (Mee, 2008). These two traits may result in direct losses due to calf and, dam mortality and premature culling, as well as indirect costs due to additional veterinary services, labor and treatment (Szucs et al., 2009).

Dystocia and stillbirth are generally scored on categorical or binary scales which make them sensitive to subjectivity (Dekkers, 1994).

Dystocia may be defined as delayed or difficult parturition. It’s an important problem in Holstein cattle since one birth of every 5 to first parity dams need assistance (Philipsson, 1996). Stillbirths are defined as a calf that dies just before, during, or within the first 24 to 48 h after birth with at least 260 days of gestation (Meyer et al., 2001; Chassagne et al., 1999).

Several studies revealed that primiparous and multiparous cows clearly differ in the rate of dystocia and stillbirths. Meyer et al., (2001) confirmed that statistical analysis of the two traits could be best when considering primiparous and multiparous cows as separated traits.

The aim of this study is to view an extract of estimations (Odds ratio) for the effect of primiparous and multiparous on dystocia and stillbirth in Holstein cattle.

MATERIALS AND METHODS

Extensive literature search of scientific electronic search engines (PubMed, Google Scholar, CAB, ISI Web of Knowledge, Science Direct, SciQuest, and Scirus) was conducted to identify primary studies carried out between 1980 and 2013. Following rigorous screening for appropriate subject matter, high quality of studies, and adequate statistical reporting, were extracted for meta–analysis. Several keyword combinations (dystocia, stillbirth, odds ratio, calving problems, Holstein cattle, meta–analysis) were used. Criteria examined included randomization of study, recording, statistical analysis. Analytic techniques described by Dohoo et al., (2003).

Articles were selected to meet the following criteria:

  1. Published in English

  2. Published as peer reviewed original articles

  3. Must had information about dystocia and stillbirth in first and later parity

  4. Articles of Holstein cows only were included in the analysis

  5. The non–peer–reviewed studies were assessed and included in the meta–analysis if they met the selection criteria

The articles selected were generally American or European as shown in Table 2 and 3.

Thirty papers were used to evaluate the impact of primiparous and multiparous cows on dystocia and nineteen papers for stillbirth. The scoring of dystocia was not constant in all papers, whereas stillbirth was recorded as dichotomous. The definition of dystocia was not standardized across studies (Table 1). Most studies classified dystocia within 5 categories including unassisted, easy, moderate, difficult, and very difficult. Some studies were recorded dystocia with four categories: = easy (non–assisted), 2 = moderate assistance (veterinarian called as precaution), 3 = difficult, 4 = very difficult with veterinary assistance. Some else recorded dystocia with three categories: no assist, easy and hard or two categories: unassisted and assisted.

Meta–analyses were conducted on dystocia and stillbirth using Comprehensive Meta–Analysis.V.2 software (2013), whereas forest plot was carried out by using MedCalc V.6 (2013). Guidelines for conducting appropriate meta–analysis were largely based on meta–

Data Analysis

Analysis of Potential Publication Bias
A funnel plot: two modes were available, one which plots a study’s effect size against its standard error and another which plots effect size against precision.

In the absence of bias the plot would be symmetric about the summary effect (Duval andTweedie2000):

Test the rank correlation (Kendall’s tau) between the standardized effect size and the variances (or standard errors) of these effects (Begg and Mazumdar, 1994).

Test the standardized effect (the regression of effect size divided by standard error on precision (inverse of standard error).

Analysis of Heterogeneity
Heterogeneity of the estimated OR was assessed using the Cochran’s Q statistic chi square test (Egger et al., 2001). If there was evidence of heterogeneity, then a random model (inverse variance) is preferred. The degree of heterogeneity was assessed by the I2 (I squared) statistic. This describes the percentage of the variability in effect estimates that is due to heterogeneity rather than sampling error (chance) (Higgins et al., 2003).

RESULTS AND DISCUSSION

Effect of Primiparous and Multiparous on Dystocia
Results shows that I2 = 99.72, and P = 0.000. An I2 value >50 may be considered indicative of substantial heterogeneity. In such case a random model is considered more suitable than fixed model (Rabiee et al., 2012). Estimated OR with a random model for dystocia is OR = 2.68, 95% CI 2.51 to 2.85 and the corresponding estimates in a fixed model is OR = 2.61, 95% CI 2.60 to 2.61 (Table 4).

It's obvious from Figure (1) the presence of bias as the OR values were distributed asymmetrically. Pooled OR was corrected according to fill and trim method "Durval and Tweedie".

Sex studies accounting for gender were needed in dystocia to be symmetrically distributed. The observed OR value of random effect was 2.68, 95% CI (2.51, 2.85), Q value 10615.92 whereas the adjusted OR was 2.42, 95% CI (2.26, 2.60), 13133.952.

Egger’s linear regression method, quantifies the bias captured by the funnel plot. Egger’s method uses the actual values of the effect sizes and their precision. The rank correlation test of Begg and Mazumdar (1994) showed that there was no significant correlation between effect and study size (P = 0.39). This was also confirmed by the regression test of Egger et al. (1997), which showed no significant association between study size and effect (intercept = 1.45, P = 0.36).

Heterogeneity in studies could be belonging to many reasons such as: studies conducted by different people, in different areas, with different definitions and at different times, which create a heterogeneous population of studies. Differences between studies in terms of the definition or measurement of outcomes, may lead to differences in observed effects (Lean et al., 2009). As I2 (I square) was significant (99.72), hence the sources of heterogeneity of response were investigated by meta–regression.

In our research, heterogeneity could be attributed to differences in definition of dystocia (categories). To investigate the validity of using this factor as predictor factor, data were analyzed using ANOVA. T–test was confirmed the significant (P < 0.05) differences between OR estimates. Hence, data were subjected to meta–regression. Two types of regression were used: fixed effect regression which shows that the slope is 0.13, 95%CI (0.12, 0.13), P = 0.000 with intercept 0.32, 95%CI (0.30, 0.34), P = 0.000 (Figure 2) and mixed effect regression which shows that the slope is 0.11, 95%CI (0.02, 0.20), P = 0.000 with intercept 0.52, 95%CI (0.12, 0.91), P = 0.000 (Figure 3). The significant effects of two regressions confirmed the effect of dystocia categories on the value of log OR. It was shown from the two Figures (2, 3) that the OR increased as category increasing.

Figure (4) illustrate each study represented by a circle proportional to its weight in the analysis. This view identifies which studies have the greatest impact on the slope of the regression line. Studies with five categories have more impact on slope.

Forest plots were used to provide illustration of the calculated OR per study as well as the overall pooled effect of all studies in the plot. The forest plot is a graphical presentation of the results that displays the point estimate and confidence interval of the effect observed in each study, along with the summary estimate and its confidence interval (Dohoo et al., 2003).

A forest plot of the studies of dystocia was shown in Figure (5).

Effect of Primiparous and Multiparous on Stillbirth
Results shows that I2 = 99.92, and P = 0.000. Estimated OR with random model for stillbirth is OR = 2.18, 95% CI 1.84 to 2.58 and the corresponding estimates in fixed model is OR = 2.06, 95% CI 2.05 to 2.07 (Table 5).

Duval and Tweedie (2000) reported that: when there was no missing study in the funnel plot, the observed and adjusted OR is identical.

Figure (6) shows that there was no bias as the studies were distributed symmetrically. The estimate of observed and adjusted OR are identical. (OR = 2.18, 95% CI 1.84 to 2.58).

Egger’s linear regression method was applied and results shows that the intercept (β0) is 1.76, 95% CI (21.70, 25.24), with t = 0.158, df = 17. The one–tailed p–value is 0.43, and the two–tailed p–value is 0.87. Begg and Masumdar rank correlation was used also and the Tau value was –0.29, P (one tailed) = 0.04 and P (two tailed) = 0.08.

It's obvious that the heterogeneity is detected in stillbirth but when there is heterogeneity that cannot readily be explained, one analytical approach is to incorporate it into a random effects model. In such case, we were unable to define the causing factors and then unable to apply meta–regression. A Forest plot for stillbirth was shown in Figure (7).

Results revealed that the test of heterogeneity confirmed the existance of a substantial heterogeneity in dystocia and stillbirth. So the estimation of OR by random model is more accurate in two mentioned traits. The OR of dystocia (2.68) is higher than stillbirth (2.18) which means that heifers is more likely to have dystocia as compared with stillbirth. Although the studies in dystocia were more as compared with stillbirth, results show that bias associated with estimation of OR was present in dystocia only. The current study confirmed that primiparous cows were most likely to have dystocia and stillbirth as compared with multiparous cows. Results also indicate that differentiation can be made among primiparous and multiparous cows in the risk of having dystocia and stillbirth. These differences among cows could be useful to aid the better management to minimize their harmful effects in the dairy herds; particularly both traits have a low heritability (Lin et al., 1989).

REFERENCES

Adamec V, Cassell BG, Smith EP, Pearson RE (2006). Effects of inbreeding in the dam on dystocia and stillbirths in US Holsteins. J. Dairy Sci. 89: 307 – 314.
http://dx.doi.org/10.3168/jds.S0022-0302(06)72095-1

Ansari–Lari M (2007).Study of perinatal mortality and dystocia in dairy cows in Fars province, south Iran. Inter. J. Dairy Sci. 2: 85 – 89.
http://dx.doi.org/10.3923/ijds.2007.85.89

Atashi H (2011). Factors affecting stillbirth and effects of stillbirth on subsequent lactation performance in a Holstein dairy herd in Isfahan. Iranian J. Vet. Res. Shiraz Uni. 12: 24 – 30.

Atashi H, Abdolmohammadi AR, Asaadi A, Akhlaghi A, Dadpasand M, Ahangari YJ (2012a). Using an incomplete gamma function to quantify the effect of dystocia on the lactation performance of Holstein dairy cows in Iran. J. Dairy Sci. 95:2718–2722. doi.org/10.3168/ jds. 2011–4954.

Atashi H, Zamiri MJ, Sayadnejad MB (2012b). The effect of maternal inbreeding on incidence of twinning, dystocia and stillbirth in Holstein cows of Iran. Iranian J. Vet. Res. Shiraz Uni. 13: 93–99.

Bar D, Ezra E (2005).Effects of common calving diseases on milk production in high yielding dairy cows. Israel J Vet Med. 60: 106 – 111.

Beaudea F, Seegers H, Ducrocq V, Fourichon C (2000). Effect of health disorders on culling in dairy cows: a review and critical discussion. Ann. Zootech. 49: 293 – 311.
http://dx.doi.org/10.1051/animres:2000102

Bell AW, Bauman DE (1997). Adaptations of glucose metabolism during pregnancy and lactation. J. Mamm. Biol. 2: 265 – 278.

Begg CB, Mazumdar M (1994). Operating characteristics of a rank correlation test for publication bias. Biometrics. 50: 1088 – 1101.
http://dx.doi.org/10.2307/2533446
PMid:7786990

Bergen W, Bates D (1984). Ionophores: Their effect on production efficiency and mode of action. J. Anim. Sci. 58: 1465 – 1483.
PMid:6378864

Berger PJ (1994). Genetic prediction for calving ease in the United States: data, models, and use by the dairy industry. J Dairy Sci. 77: 1146 – 1153.
http://dx.doi.org/10.3168/jds.S0022-0302(94)77051-X

Bicalho RC, Galvao KN, Cheong SH, Gilbert RO, Warnick LD, Guard CL (2007). Effect of stillbirth on dam's survival and reproduction performance in Holstein dairy cows. J. Dairy Sci. 90: 2797 – 2803. doi:10.3168/jds.2006–504
http://dx.doi.org/10.3168/jds.2006-504

Burrin D, Britton R (1986). Response to monensin in cattle during subacute acidosis. J. Anim. Sci. 63: 888 – 893.
PMid:3759716

Cady RA, Burnside EB (1982).Evaluation of dairy bulls in Ontario for calving ease of offspring. J. Dairy. Sci. 65: 2150 – 2156.
http://dx.doi.org/10.3168/jds.S0022-0302(82)82474-0

Chassagne M, Barnouin J, Chacornac JP (1999). Risk factors for stillbirth in Holstein heifers under field conditions in France: A prospective survey. Theriogenology. 51: 1477 – 1488
http://dx.doi.org/10.1016/S0093-691X(99)00091-6

Cleophas TJ, Zwinderman AH (2007). Meta–analysis. Circulation. 115: 2870 – 2875. http://circ.ahajournals.org/content/115/22/2870. doi: 10.1161/CIRCULATIONAHA.105.594960
http://dx.doi.org/10.1161/CIRCULATIONAHA.105.594960

Cole JB, Wiggans GR, Van Raden PM (2007). Genetic evaluation of stillbirth in United States Holsteins using a sire maternal grand sire threshold model. J. Dairy Sci. 90: 2480 – 2488. doi:10.3168/jds.2006–435
http://dx.doi.org/10.3168/jds.2006-435

Collett D (1994).Modelling Survival Data in Medical Research. Chapman and Hall, London, UK.
http://dx.doi.org/10.1007/978-1-4899-3115-3
PMCid:PMC1205911

Cook DJ, Mulrow CD, Haynes RB (1998). Systematic reviews: synthesis of the bestevidence for clinical decisions. Ann. Intern. Med. 317: 339 – 342.

Dakal K, Maltecca C, Cassady JP, Baloche G, Williams CM, and Washburn SP.(2013). Calf birth weight, gestation length, calving ease, and neonatal calf mortality in Holstein, Jersey, and crossbred cows in a pasture system. J. Dairy Sci. 96: 690 – 698. doi.org/10.3168/jds.2012-5817.
http://dx.doi.org/10.3168/jds.2012-5817
PMid:23084888

De La, Calle LML (2007). Integral study of calving ease in Spanish Holstein population. A thesis, Univ. of Zaragoza, Spain.

Dematawewa CMB, Berger PJ (1997). Effect of dystocia on yield, fertility, and cow losses and an economic evaluation of dystocia scores for Holsteins. J. Dairy Sci. 80: 754 – 761.
http://dx.doi.org/10.3168/jds.S0022-0302(97)75995-2

de Maturana, LE, Ugarte E, Gonza lez–Recio O (2007). Impact of calving ease on functional longevity and herd amortization costs in Basque Holsteins using survival analysis. J. Dairy Sci. 90: 4451 – 4457.
http://dx.doi.org/10.3168/jds.2006-734
PMid:17699066

Djemali M, Freeman AE, Berger PJ (1987). Reporting of dystocia scores and effects of dystocia on production, days open, and days dry from dairy herd improvement data. J. Dairy Sci.70: 2127 – 2131.
http://dx.doi.org/10.3168/jds.S0022-0302(87)80263-1

Dohoo I, Martin W, Stryhn H (2003). Meta–Analysis. Page 706 in Veterinary Epidemiologic Research. AVC Inc., Charlotte town, PEI, Canada.

Duval S, Tweedie R (2000). Trim and fill: A simple funnel–plot–based method of testing and adjusting for publication bias in meta–analysis.Biometrics,56:455–63.
http://dx.doi.org/10.1111/j.0006-341X.2000.00455.x
PMid:10877304

Eaglen SE, Coffey MP, Woolliams JA, Wall E (2012). Evaluating alternate models to estimate genetic parameters of calving traits in United Kingdom Holstein–Friesian dairy cattle. Genetics Selection Evolution. 44: 1 – 13. http://www.gsejournal.org/content/44/1/23.
http://dx.doi.org/10.1186/1297-9686-44-23 PMid:22839757 PMCid:PMC3468354

Egger M, Davey SG, Schneider M, Minder C (1997). Bias in meta–analysis detected by a simple graphical test. BMJ. 315: 629 – 634.
http://dx.doi.org/10.1136/bmj.315.7109.629
PMid:9310563 PMCid:PMC2127453

Egger M, Smith GD, Altman D (2001). Systematic reviews in health care: Meta–analysis in context. 2nd ed. BMJ Books, London, UK.
http://dx.doi.org/10.1002/9780470693926

Fiedlerova M, Rehak D, Vacek M, Volek J, Fiedler J, Simecek P, Masata O, Jilek F (2008).Analysis of non–genetic factors affecting calving difficulty in the Czech Holstein population. Czech J. Anim. Sci. 53: 284 – 291.

Gevrekci Y, Akbas Y, Kizilkaya K (2011). Comparison of different models in genetic analysis of dystocia. Kafkas Univ Vet Fak Derg.17: 387 – 392.

Gonzalez–Recio O, Lopez E, de Maturana, Gutierrez JP (2007). Inbreeding depression on female fertility and calving ease in Spanish dairy cattle. J. Dairy Sci. 90: 5744 – 5752. doi:10.3168/jds.2007–0203
http://dx.doi.org/10.3168/jds.2007-0203

Heins J, Hansen LB, Seykora AJ (2006).Calving difficulty and stillbirths of pure Holsteins versus crossbreds of Holstein with Normande, Montbeliarde, and Scandinavian Red. J. Dairy Sci. 89: 2805 – 2810.
http://dx.doi.org/10.3168/jds.S0022-0302(06)72357-8
http://dx.doi.org/10.3168/jds.S0022-0302(06)72356-6
http://dx.doi.org/10.3168/jds.S0022-0302(06)72545-0

Higgins J, Thompson S, Deeks J, Altman D (2003). Measuring inconsistency in meta–analyses. Br. Med. J. 327: 557 – 560.
http://dx.doi.org/10.1136/bmj.327.7414.557
PMid:12958120 PMCid:PMC192859

Johanson JM, Berger PJ (2003). Birth weight as a predictor of calving ease and perinatal mortality in Holstein cattle. J. Dairy Sci. 86: 3745 – 3755.
http://dx.doi.org/10.3168/jds.S0022-0302(03)73981-2

Lean IJ, Rabiee AR, Duffield TF, Dohoo IR (2009).Invited review: Use of meta–analysis in animal health and reproduction:Methods and applications. J. Dairy Sci. 92: 3545 – 3565.
http://dx.doi.org/10.3168/jds.2009-2140
PMid:19620636

Lin HK, Oitenacu PA, Van Vleck D, Ers HN, Smith RD (1989). Heritabilities of and genetic correlations among six health problems in Holstein cows. J. Dairy Sci. 72: 180 – 186.
http://dx.doi.org/10.3168/jds.S0022-0302(89)79095-0

Lombard JE, Garry FB, Tomlinson SM, Garber LP (2007). Impacts of dystocia on health and survival of dairy calves. J. Dairy Sci. 90: 1751 – 1760. doi:10.3168/jds.2006–295.
http://dx.doi.org/10.3168/jds.2006-295

Lopez ME, Legarra A, Varona L, Ugarte E (2007). Analysis of Fertility and Dystocia in Holsteins Using Recursive Models to Handle Censored and Categorical Data. J. Dairy Sci. 90: 2012 – 2024
http://dx.doi.org/10.3168/jds.2005-442
PMid:17369243

Mark T (2004). Applied genetic evaluations for production and functional traits in dairy cattle. J. Dairy Sci. 87: 2641 – 2652.
http://dx.doi.org/10.3168/jds.S0022-0302(04)73390-1

Martinez ML, Freeman AE, Berger PJ (1983). Genetic relationship between calf livability and calving difficulty of Holsteins. J. Dairy Sci. 66: 1494 – 1502.
http://dx.doi.org/10.3168/jds.S0022-0302(83)81964-X

Mee JF (2008). Newborn dairy calf management. Vet Clin North Am Food Anim Pract. 24(1): 1 – 17.
http://dx.doi.org/10.1016/j.cvfa.2007.10.002
PMid:18299029

Meyer CL, Berger PJ, Koehler KJ, Thompson JR, Sattler CG (2001). Phenotypic trends in incidence of stillbirth for Holsteins in the United States. J. Dairy Sci., 84, 515– 523.
http://dx.doi.org/10.3168/jds.S0022-0302(01)74502-X

Olson KM, Cassell BG, McAllister AJ, Washburn SP (2009). Dystocia, stillbirth, gestation length, and birth weight in Holstein, Jersey, and reciprocal crosses from a planned experiment J. Dairy Sci. 92: 6167 – 6175. doi: 10.3168/jds.2009–2260
http://dx.doi.org/10.3168/jds.2009-2260

Philipsson J (1996). Strategies to reduce problems in calving performance and stillbirths by selection and differential use of bulls. In: Proc. Int. Workshop Genet. Improvement of Functional Traits in Cattle. Gembloux, Belgium, Interbull Bulletin. p. 12: 65 – 71.

Rabiee AR, Breinhild K, Scott W, Golder HM, Block E, Lean IJ (2012). Effect of fat additions to diets of dairy cattle on milk production and components: A meta–analysis and meta–regression. J. Dairy Sci. 95: 3225 – 3247.
http://dx.doi.org/10.3168/jds.2011-4895
PMid:22612958

Steinbock L (2006). Comparative aspects on genetics of stillbirth and calving difficulty in Swedish dairy cattle breeds. Licentiate thesis, Swedish University of Agricultural Sciences, Uppsala. pub. epsilon. slu.se /1082/ 1/ SammanläggningLenaS.pdf‎

Steinbock L, Nasholm A, Berglund B, Johansson K, Philipsson J (2003). Genetic effects on stillbirth and calving difficulty in Swedish Holsteins at first and second calving. J. Dairy Sci. 86: 2228 – 2235.
http://dx.doi.org/10.3168/jds.S0022-0302(03)73813-2

Szucs E, Gulyas L, Cziszter LT, Demirkan I (2009).Stillbirth in dairy cattle: review. Zootehnieşi Biotehnologii. 42: 622 – 636.

Van Plet ML, de Jong GE, Roelfzma JE (2009). Analysis of calving traits with a multitrait animal model with a correlated direct and maternal effect and reciprocal crosses from a planned experiment. J. Dairy Sci. 92: 6167 – 6175.

Van Tassell CP, Wiggans GR, Misztal I (2003). Implementation of a sire–maternal grandsire model for evaluation of calving ease in the United States. J. Dairy Sci. 86: 3366 – 3373.
http://dx.doi.org/10.3168/jds.S0022-0302(03)73940-X

Wall E, Coffey M, Mrode R, Krzyzelewski T, Banos G (2008). Feasibility of producing calving ease evaluations for UK dairy cattle. Project funded by the Milk Development Council.
www.dairyco.org.uknon...download.aspx

Weller JI, Misztal I, Gianola D (1988).Genetic analysis of dystocia and calf mortality in Israeli–Holsteins by threshold and linear models. J. Dairy Sci. 71: 2491 – 2501.
http://dx.doi.org/10.3168/jds.S0022-0302(88)79836-7

Wiggans GR, Cole JB, Thornton LLM (2008). Multiparity evaluation of calving ease and stillbirth with separate genetic effects by parity. J. Dairy Sci. 91: 3173 – 3178. doi:10.3168/jds.2007–09.