- Original article
- Open Access
Mortality in over 350,000 Insured Swedish Dogs from 1995–2000: II. Breed-Specific Age and Survival Patterns and Relative Risk for Causes of Death
Acta Veterinaria Scandinavica volume 46, Article number: 121 (2005)
This study continues analysis from a companion paper on over 350,000 insured Swedish dogs up to 10 years of age contributing to more than one million dog-years at risk during 1995–2000. The age patterns for total and diagnostic mortality and for general causes of death (trauma, tumour, locomotor, heart and neurological) are presented for numerous breeds. Survival estimates at five, eight and 10 years of age are calculated. Survival to 10 years of age was 75% or more in Labrador and golden retrievers, miniature and toy poodles and miniature dachshunds and lowest in Irish wolfhounds (91% dead by 10 years). Multivariable analysis was used to estimate the relative risk for general and more specific causes of death between breeds accounting for gender and age effects, including two-way interactions. Older females had tumour as a designated cause of death more often than males in most breeds, but not in the Bernese mountain dog. Information presented in this and the companion paper inform our understanding of the population level burden of disease, and support decision-making at the population and individual level about health promotion efforts and treatment and prognosis of disease events.
Mortalitet hos hund baserad på data från över 350 000 svenska hundar försäkrade under 1995 till 2000: II. Rasspecifika ålders och överlevnadsmönster och relativ risk för olika dödsorsaker.
Detta är del 2 i en studie där data på över 350 000 försäkrade svenska hundar under 10 års ålder analyserats. Resultaten grundas på mer än en miljon hundår av risk under 1995 till 2000. Åldersmönstret för den totala och diagnostiska mortaliteten, och också för generella dödsorsaker (trauma, tumörer, problem från rörelseapparaten, hjärta och neurologiska problem) presenteras för ett flertal hundraser. Beräknad överlevnad vid 5, 8 och 10 års ålder presenteras. Överlevnaden vid 10 års ålder var 75% eller högre hos labrador och golden retriever, dvärg- och mellanpudel samt dvärgtaxar och lägst hos irländska varghundar (91% var döda vid 10 års ålder). Multivariabel analys användes för att beräkna relativ risk för generella och mer specifika dödsorsaker för olika raser, med hänsyn taget till könsoch ålderseffekter vilket också inkluderar två-vägs interaktioner. I de flesta hundraser hade äldre tikar tumör som dödsorsak oftare än äldre hanhundar, med undantag för berner sennen hundar. Information som presenteras i denna del (2) och del 1 ger insikt om dödsorsakspanoramat hos hund. Detta kan ge stöd för beslut på både populations och individuell nivå om vilka åtgärder somkan vidtas för att förbättra hundars hälsa, såväl behandlingsmässigt som prognostiskt.
The burden of disease within the purebred dog population has become the focus of discussions on welfare considerations and with regard to the ethics of breeding practices . There is a risk that breeding to eliminate specific genetic defects, without clear understanding of the inter-relationships between risk factors and disease expression could do more harm than good . Improving the overall health and welfare of the dog population requires that we quantify and monitor disease and death across breeds and over time.
Identifying those breeds with significantly increased risk of certain health problems is important for dog breeders, researchers and the veterinary community. In the accompanying paper, we presented mortality rates for various breeds of insured dogs under 10 years of age . Rates were presented using the exact time at risk for over 350,000 dogs contributing to over a million dog-years at risk (DYAR) during six years for total mortality (all deaths), diagnostic mortality (those deaths for which an insurance claim was processed), for deaths within various diagnostic categories and for some specific diagnoses. Breed rates for total mortality ranged from 230 to 1,574 and diagnostic mortality from 168 to 1,319 deaths per 10,000 DYAR.
In addition to mortality rates (MRs), proportional mortality estimates were presented. At the population level it was informative to see that, for example, German shepherd dogs accounted for 7% of the insured population and almost 12% of deaths. Bernese mountain dogs were over-represented in mortality at twice their occurrence in the population. Golden retrievers, on the other hand, accounted for less than 4% of the deaths while being more than 6% of the population.
Owner-pet mismatch has been suggested as an important cause of pet relinquishment. Specific information about expected longevity  and disease risks should be part of an informed choice for prospective dog owners, as individuals may have differing preferences and concerns. Decisions regarding highly invasive and costly veterinary care in older dogs should include consideration of both the likely duration and quality of life to be expected. Although many anecdotal reports on average age at death are available, and there are some better-quantified statistics, good estimates on rates of disease and expected longevity are not available for most breeds.
In addition to breed differences it is well recognized that the occurrence of many causes of death vary by age and other host factors. Differences between males and females for the rate of death due to certain conditions have been shown (e.g. ). In order to examine the combined effects of breed, gender and age, multivariable modelling techniques are useful.
This study is a further analysis of the database presented in the companion paper. The objectives of this paper are to present the age pattern for various causes of death for numerous breeds of Swedish dogs insured between 1995 and 2000, to calculate survival estimates at five, eight and 10 years of age and to present the relative risk for general and more specific causes of death between breeds accounting for gender and age effects.
Materials and methods
Details on the database and analysis, including calculation of rates and proportional mortalities have been presented in the companion paper . Briefly, dogs covered by life insurance from 1995 to 2000 in a Swedish animal insurance company (Agria Insurance, PO 70306, SE-107 23 Stockholm, Sweden) were included in this study. Most dogs were enrolled as puppies and life insurance was, because of the terms of the insurance policy, terminated by 10 years of age. Deaths may be claimed, in which case there is a recorded diagnosis for the cause of death or euthanasia, usually made by a veterinarian, but which might be 'dead-no diagnosis'. These deaths are included in 'diagnostic' MRs. Deaths are often recorded, even when not claimed and there is no diagnosis. Deaths with or without claims are included in 'total' MR calculations. Variables downloaded from the insurance company database were breed, date of birth, date of death, gender (male/female, not neuter status), diagnostic code for death, date when the dog entered or left the insurance program and information on the type of insurance for which the dog was enrolled.
A hierarchical diagnostic registry is used by veterinarians and at the insurance company to assign causes which are classified within organ system and process . For this study, the registered causes of death were partitioned into six diagnostic categories – tumours, trauma, locomotor disorders, heart, neurological and other problems. Tumours were those diagnoses listed under process neoplastic. The diagnostic category trauma included all diagnoses that were listed as traumatic processes; locomotor disorders were all those in systems skeletal, muscles and joints, except for the traumatic or neoplastic processes. Similarly, the diagnostic category heart consisted of all heart diagnoses, and neurological disorders included all diagnoses said to emanate from the nervous system, except for those said to be neoplastic or traumatic in origin. Cases with the diagnosis epileptiforme seizures were included under the diagnostic category neurological (from system unspecified). The specific diagnosis cruciate rupture was included under the diagnostic category locomotor and gastric dilatation/volvulus was included under the diagnostic category other (both from the traumatic process). All diagnoses that did not belong in any of the preceding categories were said to be in the diagnostic category other.
Breed-specific survival, for total and diagnostic mortality, to ages five, eight and 10 years has been estimated using Cox regression, using the baseline survival statement. There were no independent variables; each dog was entered the 1st of January 1995 or the day of enrolment. If not a case, nor censored during the year, it was censored the 31st of December 2000. For example, a dog that died at five years of age was coded as "failed" at five years of age. In addition, the survival to eight years for those alive at five, and to 10 years for those alive at eight was calculated.
The age- and breed-specific and age-, breed-and diagnostic category-specific MRs (hazards) were constructed for common and high-risk breeds (for breed groupings see table 1) using the SMOOTH macro , which creates age-specific hazards from the survival function computed by the SAS (SAS Institute Inc., Cary, NC, 27513, USA) procedure PHREG. The macro provides a smoothed estimate of the hazard curve using a kernel smoothing method. This involves arbitrarily setting the WIDTH parameter, which influences the degree of smoothing, to achieve a reasonable curve; in this case one-fifth of the range of event times was chosen. The mortality rates were plotted against age using different scales to adjust for marked differences across breeds. Poisson regression was used regressing total deaths, all deaths with a diagnosis and death due to specific diagnostic categories on age, gender and breed. In these analyses 20 dummy variables were created for the 20 breeds (common and high-risk) using 'other' breeds as the initial baseline. Age was entered as a continuous variable, after verifying that age could be considered adequately linearly related to the outcome. Two-way interactions were tested, one at a time, upon the unreduced main effect models. Interactions that were significant at P ≤ 0.05 were kept for further modelling. Models were reduced based on the type 3 criterion. A P-value of 0.001 was considered significant. For the gender-breed and age-breed interactions only those terms that were statistically significant were kept in the models. Model fit was inspected by plotting standardised deviance residuals against covariates and predicted values. Over-dispersion was not evaluated as each dog belonged to multiple records in the dataset. The statistical software program SAS (SAS Institute Inc., Cary, NC, 27513, USA) was used to analyse the data. The PHREG procedure was used for Cox regression and the GENMOD procedure for Poisson regression.
For the presentation of the Poisson models' breed-specific mortality rate ratios (MRR), estimates and SE's have been reported for all the included effects. Single and multiple effect MRR's were constructed based on the formula (for a single effect):
which shows the rate ratio for a one-unit change in xi and where E is the size of the exposure. The 99.9% confidence intervals (99.9% CI) were constructed taking βi ± 3.291 × SE, and then exponentiating these limits. Given the size of the data base, 99.9% limits were chosen to yield conservative results. The point estimate for a given age, for example a 9 < 10 year old female of a specific breed is: exp(estimate for breed + estimate for gender + estimate for age × 9 + the gender-age interaction for age × 9 + the breed-age interaction for age × 9 + the gender-breed interaction). This gives the MRR with respect to breed, age and gender simultaneously and expresses the risk relative to the baseline.
Figures 1 and 2 present the MRs from Cox regression (hazards) for the total and diagnostic mortality, as well as diagnostic category-specific mortality in six selected common breeds of dogs and six selected high-risk breeds respectively. Note that the scales on the y-axis vary from 7,000 in the Irish wolfhound and St Bernard to 800 in the mongrel and golden retriever. It is evident that in different breeds, different cause-specific mortality curves depart more clearly from the cluster of causes with lower MRs.
The probability of death (1-survival) to five, eight and 10 years of age for total and diagnostic mortality is presented for all breeds and groups, and for all dogs by gender in Table 1. For example, 91% of wolfhounds were dead by 10 years of age and given that an Irish wolfhound reached the age of eight years, the probability that it would be dead by 10 years of age was 76%.
Table 2 shows MRRs, with 99.9% CIs, the estimates and standard errors from the Poisson regressions for total and diagnostic mortality and tables 3 and 4 show the same for the models by diagnostic category. Model validation plots of standardised deviance residuals against predicted values demonstrated in general a satisfactory fit. However, some deaths with low predicted values had large residuals and did not fit the models. There was no pattern in the plots of deviance residuals against different covariates. In all models all main effects were statistically significant, although the insignificant terms within main effects were successively removed. For the trauma model younger animals were at increased risk; for all other models mortality rates increased with age. In all models the main effect for gender was negative indicating that, overall, females had lower mortality compared to males.
However, positive age-gender interactions were found for total mortality, diagnostic mortality, and tumours, indicating that for these conditions older females had a significantly increased risk compared to males. The age-gender interaction term for locomotor was negative; i.e. older males had higher relative mortality.
There were several significant age-breed interaction terms in several models, namely for total and diagnostic mortality, for locomotor, heart and neurological. In all but the neurological model these interactions were both positive and negative. For example in the heart model, the age-breed interactions for Cavalier King Charles spaniels and dachshunds were positive, while the one for German shepherds was negative (further illustrated below).
Statistically significant terms for certain breed-gender interactions were found in the diagnostic mortality and the tumour model (see below). For trauma, only main effect variables were significant, and, therefore, the interpretation of the breed-specific MRR is the effect of breed after adjusting for age and gender differences.
To be able to interpret the multivariable results from the Poisson regression more completely with respect to all of the variables, including the interactions, figures 3, 4, 5 have been constructed. These exemplify how the MRRs vary across breeds, genders and ages for the models with respect to diagnostic mortality, tumours and heart.
In figure 3, lines are shown for three breeds and the baseline for diagnostic mortality. One of the lines, for the German shepherd, is shown because there was interaction between breed and gender. By comparing the St Bernards and the baseline breeds to the breed with the breed-gender interaction, German shepherds, it is evident that the MR in females is consistently lower in male German shepherds, while in the other breeds there is no difference between genders in the older age categories.
In figure 4, comparing the Bernese mountain dogs with the Irish wolfhounds and the other breeds (boxer and baseline), it can be seen that the male Bernese mountain dogs have a higher MRR for tumours at all ages, whereas for the other breeds the tumour rate is actually higher in the females at higher ages.
Figure 5 shows results for the heart category. The positive age-breed interaction for the Cavalier King Charles spaniel clearly makes the curves bend upwards, while the negative breed interaction in the German shepherd visually has a smaller effect. There are large MRRs by breed and gender for a few breeds (e. g. Cavalier King Charles spaniels and Irish wolfhounds) compared to baseline, especially for males at older ages.
Age patterns of death are informative both within and across breeds. This paper presents comprehensive information on mortality by including individual breed patterns of death for overall mortality and by category-specific causes. Survival estimates are clinically useful information within a breed or for contrasting across breeds. Finally, multivariable models quantify the relative rate between breeds adjusted for effects of age and gender for different causes of death. Most of the age effects are, presumably, largely biological, however the model also adjusts for potential differences in age distributions between breeds that might be due to differences in enrolment in insurance or other factors. Although this study presents appropriate quantitative measures, the output can be somewhat overwhelming. Exhaustive presentation of relative rates of death for every breed, both genders and at all ages would be beyond the scope of a research publication. Therefore, estimates and standard errors from the modelling procedures have been presented to allow readers to calculate specific MRRs for any gender-age-breed combination in which they are interested. Graphic figures have been included to provide a more obvious display of some relationships across breeds.
In addition to differences between breeds, it was shown in the companion paper that for some conditions there was a statistically significant effect of gender . For example, males tended to have higher MRs for trauma and females had higher rates of cancer, except in the Bernese mountain dog, where males had twice the rate of tumours. Such complex relationships can be examined using statistical modelling. As the goal of this study was to produce clinically relevant statistics on death in dogs, a practical and parsimonious approach to modelling was taken.
Cox regression was used for the survival analysis as an appropriate technique to model time to death. Multivariable Cox regression models were initially attempted but proved unmanageable in terms of computer and time constraints, therefore Poisson regression was used for the calculation of MRRs, allowing inclusion of a primary baseline category of all other breeds.
Exploring the fit of the models several large residuals were found, mainly for low predicted values, however this was not related to low age categories. Also, the model fit was worst in the total and diagnostic mortality models (data not shown), likely partly because these contained larger numerators, and accordingly by chance some more deaths will be ill-fitting (and more deaths will be well-fitting). Further, all types of deaths were included in the outcomes and, as only a few variables were available to model, it is likely that important predictors of deaths were not included, thus limiting the fit of the model.
Unfortunately over-dispersion could not be evaluated because multiple records belonged to each dog, as the dogs moved through the age categories. However, in previous work over-dispersion has never been a problem working with the Agria insurance database (, in press, unpublished observations).
Age patterns and survival
Different breeds age at different speed, however in many breeds the mortality starts to increase after about six to seven years of age, similar to findings of increased morbidity and mortality seen in previous work by these authors . Superficially, this may be seen as supporting anecdotal suggestions that 'geriatric profiling' be started on all breeds at seven years of age. However, this increase must be considered together with the magnitude of mortality at specific ages in individual breeds of dogs. There is no accepted mortality level at which we designate dogs to be geriatric, or even senior. However, considering that for some breeds over 50% of dogs are dead by eight years of age (e.g. Irish wolfhounds, St. Bernard, great Dane) whereas in other breeds less than 25% have died (Labrador and golden retrievers, poodle, miniature dachshund), it is inappropriate to consider them as having equivalent biological age or similar risk of age-related diseases that might be detected by screening programs.
For tumours, the mortality rate starts to increase after six or seven years in Labrador retrievers and mongrels, although it is still below 200 deaths per 10,000 DYAR. In Bernese mountain dogs tumour deaths are common starting by four years of age and account for over 2,000 deaths per 10,000 DYAR by eight years of age. For boxers there is a sharp increase to very high levels after the age of six years. These findings highlight the probability that some previously published breed risk estimates are misleading. Proportional estimates of death from data sources with no population at risk (e.g. postmortem registries) have likely made erroneous conclusions about cancer risks across breeds .
The age patterns are reflected in the survival statistics and highlight differences across breeds. For Cavalier King Charles spaniels there is a very low early mortality with a sharp rise after five years of age. Only 7% of Cavalier King Charles spaniels are dead by five years of age, but 48% are dead by 10 years. This is in contrast to a breed like the German shepherd, where the early mortality is high (20% mortality by five years) but by 10 years of age the survival (49%) is quite similar to the Cavalier King Charles spaniel. Boxers have relatively low early mortality – one in nine will die before five years of age. One in five will die between five and eight and of those alive at eight, one-third will die by 10.
Information on the average survival pattern for different breeds is useful clinical information. Prospective owners should understand that in getting, for example, either a Bernese mountain dogs or a golden retriever, the likelihood of the dog living past 10 years of age is very different. Most owners, once their dog becomes seriously ill or dies, are keen to know if the condition is rare or common, in general or in their breed. Many owners are especially concerned if they feel their dog has died prematurely, that is at a younger than average age. Veterinarians and their clients can also use survival statistics as part of informed decision-making regarding expensive veterinary care and for conditions in which quality of life is an issue. For example, an Irish wolfhound that has survived to eight years has a 76% chance of being dead before 10, whereas an eight-year old golden retriever has a 90% chance of surviving to 10.
Prospective owners may wish to consider breed-specific patterns of death or disease before choosing a breed. Owner preferences and expectations about, for example, longevity, may influence their satisfaction with a chosen pet . For health management at the individual dog level, it is also useful for veterinarians to have quantified estimates of the increased risk certain breeds have for certain conditions. These findings inform diagnostic decision-making and facilitate communication with owners.
The MRRs highlight both the relative magnitude and breed differences for total and diagnostic mortality. Differences between the total and diagnostic mortality rates within breeds were discussed in the companion paper . Elective euthanasia for behaviour problems or dog-owner mismatch does not get reimbursed. Total mortality statistics probably best reflect overall mortality across breeds. However, for examining disease conditions as causes of death, the diagnostic mortality statistics are more informative, even though these are only available for the subset of dogs for which claims have been submitted.
Risk of tumour, in general was higher in males (MRR < 1), however the interaction term (age-gender, > 1) tells us that older females were at higher risk of tumours. This is likely due to the high rate of mammary tumours in females ( in press). Most interesting in the tumour model was the interaction term for Bernese mountain dogs. Overall this breed is at a much greater risk of dying from tumour than the baseline breeds and the interaction term tells us that the effect at older ages is the opposite in the Bernese mountain dog compared to other breeds, i.e. at older ages males have a higher risk of death from tumours than females. Note that confidence intervals are not included in figures 3, 4, 5.
Dachshunds (both breed groups) and poodles were at decreased risk of death due to tumours, as were Cavalier King Charles spaniels. One explanation might be a lower biological age, for example, in poodles. As only one quarter of poodles die before 10, poodles in this study population must be considered to have not even reached 'middle age'.
Traumatic deaths were more common in males and at younger ages (although there was no age-gender interaction). As was discussed in the companion paper, drevers had a high rate of death due to traumatic causes, and that was reflected in the multivariable models, where they were at 4.5 times increased risk than the baseline group. The risk for dachshunds and greyhounds were also high, although from different causes, i.e. mainly car accidents for the former and a great deal of fractures in the latter (data not shown). Interestingly German shepherds, Cavalier King Charles spaniels, Labradors, springer spaniels and golden retrievers had a significantly reduced risk of death due to trauma compared to baseline breeds. Reasons for this are unknown, but could relate to training and use of these animals and behavioural characteristics of the dogs or their owners.
Golden retrievers were at low risk for mortality in this study – only 22% died before 10 years. Golden retrievers were significantly less likely to die of trauma and heart disease and were in the baseline (average) risk group for neurological and tumour causes of death. They were at increased risk in the first age category for locomotor problems, but this effect waned with age as demonstrated by a negative age-breed interaction.
Overall, males were at increased risk and risk increased with age for locomotor causes of death. The interaction term indicates that the difference between males and females was more pronounced at higher ages. Many breeds were at increased risk of death from locomotor problems, but springer spaniels, mongrels and poodles had a significantly decreased risk. In five breeds (table 3) there were negative age-breed interactions. In four of these the risk in these breeds were pronounced at early age, but the risk waned with age. However, in the miniature dachshund, the risk was similar to baseline at low ages and less than baseline in older ages. Females were at reduced risk of dying from heart diseases compared to males. Risk increased quite sharply with increasing age. Several breeds had markedly increased risk of death due to heart causes, for example, Irish wolfhounds and great Danes were at 29 and 21 times increased risk, respectively, compared to baseline breeds. The Cavalier King Charles spaniel had a low MRR of 3, but also a positive age-breed interaction. This made the MRR in the first age category low, 1.7 and 3.0 for females and males respectively, but high, 446 and 816 in the last age category (fig 5). The miniature breeds (poodle and dachshund), the two retriever breeds and mongrels were at decreased risk.
The risk of death due to neurologic causes increased with age and males were 1.4 times more likely to die from neurologic causes compared to females. There were significant age-breed interactions in two breeds, boxers and German shepherds. Both were positive. For example, in boxers the risk in the lowest age category was not significantly different from baseline, but with age the MRR increased to 13.5 and 19 in females and males respectively (using the formulas from the part of materials and methods). In baseline breeds these MRR were 1.6 and 2.3 and females and males respectively. Further, St. Bernards and great Danes were 3.5 to five times more likely to die from neurological causes, and both regular and miniature dachshunds were over three times less likely to die of neurologic causes than other breeds (baseline). Labradors were also significantly less likely to die of neurologic causes than baseline breeds. As was seen and discussed in the companion paper, a large proportion of neurologic deaths were due to epileptic disease.
Information presented in this and the companion paper inform our understanding of the population level burden of disease and support decision-making at the population and individual level about health promotion efforts and treatment and prognosis of disease events. Findings such as these provide some new insights and, even where they confirm clinical impressions, offer the advantage of both relative and absolute quantification of the direction and magnitude of breed differences and the effects of age and gender. Population statistics on health and disease of dog populations will be a crucial part of breeding strategies and monitoring as we enter into a new era with knowledge of the canine genome. Increased disease in purebred populations has been identified as a welfare issue and as stewards of animal health the veterinary profession should work with other stakeholders in animal health to reduce the occurrence of preventable disease.
Allison P: Survival analysis using the SAS® system: a practical guide. 1995, Cary
Bonnett BN, Egenvall A, Olson P, Hedhammar Å: Mortality in over 350,000 insured Swedish dogs from 1995–2000: I. Breed-, gender-, age- and cause-specific rates. Acta vet scand. 2005, 46: 105-120. 10.1186/1751-0147-46-105.
Craig LE: Cause of death in dogs according to breed: a necropsy survey of five breeds. J Am Anim Hosp Assoc. 2001, 37: 438-443.
Egenvall A, Bonnett BN, Olson P, Hedhammar Å: Gender, age, breed and geographic pattern of morbidity and mortality in insured dogs during 1995 and 1996. Vet Rec. 2000, 146: 519-525.
Egenvall A, Bonnett BN, Öhagen P, Olson P, Hed-hammar Å, von Euler H: Incidence of and survival after mammary tumours in a population of over 80,000 insured female dogs in Sweden from 1995–2002. Prev Vet Med.
Meyers-Wallen VN: Ethics and genetic selection in purebred dogs. Reprod Domest Anim. 2003, 38: 73-76. 10.1046/j.1439-0531.2003.00384.x.
Michell A: Longevity of British breeds of dogs and relationships with sex, size, cardiovascular variables and disease. Vet Rec. 1999, 145: 625-629.
Ott RS: Animal selection and breeding techniques that create diseased populations and compromise welfare. J Amer Vet Med Ass. 1996, 12: 1969-1974.
Svenska Djursjukhusföreningen (Swedish animal hospital association): Diagnosregister för häst, hund och katt. Diagnostic registry for the horse, the dog and the cat. Taberg. 1993, In Swedish
This work has been supported by grants from the Foundation for Research, Agria Insurance and the Swedish Kennel Club.
Reprints may be obtained from: A. Egenvall, Department of Clinical Sciences, Faculty of Veterinary Medicine and Animal Science, Swedish University of Agricultural Sciences, SE-750 07 Uppsala, Sweden.
About this article
Cite this article
Egenvall, A., Bonnett, B., Hedhammar, Å. et al. Mortality in over 350,000 Insured Swedish Dogs from 1995–2000: II. Breed-Specific Age and Survival Patterns and Relative Risk for Causes of Death. Acta Vet Scand 46, 121 (2005) doi:10.1186/1751-0147-46-121
- competing risk analysis
- survival analysis.