Open Access
Translator Disclaimer
1 March 2002 Age determination in leverets of European hare Lepus europaeus based on body measurements
Yves Bray, Stephane Champely, Dominique Soyez
Author Affiliations +

The objective of our work was to provide allometric relationships between different body measurements and age of young leverets of European hare Lepus europaeus (less than two months old) for use in live-trapping studies. Five morphometric variables were measured to predict the age of leverets: body mass, skull length, skull width, hind foot length and ear length. The hares came from two rearing centres near Paris (France) and were handled in 1995 and 1996. Measurement repeatability of each variable was compared in a pilot study (N = 51 leverets) and the effect of rearing centre on the growth curve was tested for the first 30 days after birth. We used a complementary sample (N = 168 leverets) to establish a model of age determination using I-splines regression. The skull length was the best candidate variable based on measurement repeatability, comparability between the two rearing centres and goodness of fit. I-spline modelling has minimal underlying statistical assumptions. In addition, it provides easy-to-use and powerful plots of the relation between age and measurements, including confidence limits estimated by a bootstrap procedure.

Understanding the population ecology of an animal species most often requires analyses of age-dependent processes (Steams 1992). As a special case, age-structured patterns in demographic parameters are particularly relevant to the analysis of population dynamics (Emlen & Pikitch 1989), and may interfere with population management, e.g. when considering game species and fragmented animal populations (Lebreton & Nichols 1998). Recent data about pre-weaning mortality and movements in the snowshoe hare Lepus americanus (O'Donoghue & Bergman 1992, O'Donoghue 1994) best illustrate the importance of analysing agerelated patterns in demographic parameters to identify key stages in the life cycles of animals. For European hares L. europaeus, the population growth rate seems to be highly dependent on breeding success (Marboutin & Péroux 1995), but the relative influence of changes in fecundity, juvenile survival and dispersal has never been documented. As part of an analysis of dispersal patterns in European hares from central France (Bray 1998), we investigated early movements and survival of wild new-born leverets using radio-tagged animals. It was therefore of great importance that leverets be aged accurately when first captured.

As in other mammal species (Morris 1972), the methods available for age determination in hares largely depend on whether animals are dead or alive. Some criteria have already been proposed for age determination of dead hares, and the most reliable criteria is eye-lens weight (Walhovd 1966, Pepin 1973, Broekhuizen & Maaskamp 1979, Suchentrunk, Willing & Hartl 1991). As far as live animals are concerned, it is possible to roughly distinguish between young (s7 months of age) and adults based on the disappearance of the epiphysial cartilage on the distal part of the ulna (Stroh 1931). Of other methods based on the ossification process (Lloyd 1968, Frylestam & von Schantz 1977, Pascal & Kovacs 1983) only Stroh's can be used under field conditions without laboratory equipment. Determining the age of very young specimens, such as leverets, has often been based on growth curves obtained under various circumstances. Sometimes the age has been estimated from repeated weighing of hand-reared hares (Hediger 1948, Pepin 1973, Biadi & Nouhaut 1975). Other times estimations have been based on recapture of young leverets found in the field and maintained in close captivity (Pilarska 1969, Broekhuizen 1971) or on leverets bom and maintained in large enclosures where natural conditions were imitated (Pielowski 1971). The resulting growth curves were highly variable (Fig. 1), suggesting that different local conditions and genetic backgrounds may induce changes in growth rates and, as a by-product, uncertainty in the prediction of the age. Last but not least, these curves only deal with body mass, and they do not provide any estimate of confidence limits. Other body measurements have already been used (e.g. length of hind foot, ear or skull), but the people who used these measurements were interested either in describing the species (Flux 1967), or in distinguishing between young and adults (Broekhuizen 1979).

Therefore, we investigated growth curves of five body measurements in hand-reared leverets, because recapture and repeated handling of wild animals would have been more invasive and would likely have induced more stress-related bias in growth rates. Furthermore, sample sizes could be increased when using hand- reared hares. The actual birth dates of the hares used were known, and the stability of results obtained for different rearing centres and observers would help in assessing the external validity of the method.

Figure 1.

Body weight growth curves of European hare leverets during their first 100 days of life. Data were obtained from the following sources: Pepin 1979 (on wild hares in France), Hedinger 1948 (on hand-reared hares in Germany), Pepin 1973 and our data (on hand-reared hares in France), Pielowski 1971 (on five leverets bom and kept in semicaptivity in Poland), Pilarska 1969 in Poland and Broekhuizen 1971 in the Netherlands on leverets bom in the field and kept in close captivity.


Material and methods

We used leverets up to two months of age from two rearing centres; measurements were done in a first step by two observers to assess repeatability where some error measurement of variables could occur. Autocorrelation in data was limited by subsampling when repeated measurements of the same individuals occurred. Changes in body measurements as a function of age were modelled by means of semi-parametric non-linear functions (I-splines; Ramsay 1988). These functions assume limited underlying statistical assumptions, and provide easy-to-use and powerful plots of the relation between age and measurements, including confidence limits estimated by a bootstrap procedure.

Sample collection

Data on age and body measurements were collected during 1995–1996, in two rearing centres located close to Paris (France). For one rearing centre, data were only available for the first month of age. As measuring accuracy may partly depend on an observer-effect, a pilot study of measurement repeatability was conducted, based on the comparison of paired data from two observers. A total of 51 leverets were used, from which both observers recorded independently the same body measurements.

Figure 2.

Skull length (SL; condylobasal length) corresponded to the distance between the back of the skull and the tip of nose (A); maximum skull width (SW: larger value of bizigomatic bredth) (B); hind foot length (FL; including nail) was measured as the distance between the foremost nail point and the end of the heel bone (C); ear length (EL; notch to tip) was measured excluding hairs (D).


Then, a complementary sample was obtained, to achieve a total sample size of N1 = 219. Measurements (N2 = 522) were performed from birth to the age of 65 days, every third day during the first month, and every week during the second month.

Variables and measurements

We considered five variables that could be related to body growth and age: body mass (BM), skull length (SL) and width (SW), hind foot length (FL; including nail), and ear length (EL; notch-to-tip excluding hairs). Mass was measured to the nearest gram using an electronic weighing machine, and further transformed using cube root. Skull length and width were measured to the nearest 0.1 millimetre using a calliper rule. Hind foot and ear lengths were measured to the nearest millimetre using a steel scale. Details on measurements of SL, SW, FL, and EL are presented in Figure 2.

Statistical analyses

We performed systematic subsampling within data from every leveret, to limit time-related autocorrelation in measurements. Serial autocorrelation was no longer significant when consecutive data were spaced using a 6-day period. For example, when computed from long series, the first-order coefficient of serial autocorrelation (time-step = three days) ranged within 0.70–0.75 in SL (upper 95% confidence limits: 0.57), whereas the second-order coefficient of correlation (time-step = six days) ranged within 0.42–0.48 (upper 95% confidence limits ≥0.8). As a result, the statistical analysis was based on N3 = 388 measurements in any of the five variables analysed.

Figure 3.

Bootstrap procedure for I-spline modelling of (X,Y). I-spline fitting as Equation 3 in A), histogram of adjusted errors as Equation 5 in B), 1,000 bootstrap I-spline models as Equation 6 in C), and 95% prediction intervals for skull length in D). The original data were superposed on the plot.


Covariance analysis of growth curves

In order to compare the growth curves and to test the equality of the growth curves from the two rearing centres, we used an extension of the classical covariance analysis suitable for growth curves.

Let Xij be the age of leveret i from centre j (j = 1,2) and yj(xij) be the response variable. In order to gain some flexibility in the modelling, we used a semi-parametric approach based on B-splines (de Boor 1978). This can be done using two nested linear models of the form:

where B() and Bj() were (cubic) B-splines with v degrees of freedom, and ϵij was the conventional error term assumed independently normal with zero mean and constant variance.

The model in Equation 1 consists of one growth curve common to the two centres, whereas the model in Equation 2 represents two different growth curves.

To compare the two models, the corresponding test of the centre effect is an F-test. An example of such computations using the S-Plus® software was described in Chambers & Hastie (1991: 280).

Modelling with I-splines

Let (xi,yi) be the data at hand where x is the age of the leveret and y a growth variable. A monotone spline regression (we used the notations from Ramsay 1988) was used to model the growth curves:

where ϵi was the error term and Ij() a set of monotonous transformation functions called I-splines. They corresponded here to a piecewise polynomial of degree k defined on an interval [L,U] with m interior knots positioned at appropriate x-quantiles.

The parameters αj were estimated in a least square sense with respect to positive constraints on the parameters for j > 0. These positivity constraints were sufficient to ensure monotonicity.

The resulting fit was summarised by the classical coefficient of determination (percentage of variance explained by the model). We developed S-Plus® routines based on the gradient projection algorithm to fit this model (see the Appendix in Ramsay 1988).

Calibration by bootstrapping residuals

First, our aim was to obtain prediction intervals for the response variable y computed on a grid of x values. Secondly, the resulting plot was inverted to provide calibration intervals for x.

The I-spline model (Equation 3) was expressed as

The errors 8 were assumed to be a random sample from an unknown distribution having null expectation and constant variance.

In order to define a prediction interval at x = xnew, we used a bootstrap approach. As a matter of fact, a direct mathematical solution seemed difficult for two reasons. First, the residuals were clearly non-normal (Fig. 3B) and second, the model was fitted by a complex optimisation routine.

To calculate a bootstrap data set y*(xnew) for y(xnew) and thus a prediction interval by using appropriate quantiles, one needed to sample independently bootstrap data sets μ*(xnew) for μ(xnew) and ϵ* for ϵ, and then to compute:

For bootstrapping ϵ* in Equation 4, let ri be the residuals of the I-spline fitting (Fig. 3A):
The 'adjusted errors'
were closer to errors (Efron & Tibshirani 1993: 122). The bootstrap dataset ϵ* was obtained by sampling with replacement of the adjusted errors ei (see Fig. 3B).

For bootstrapping μ *(xnew) in Equation 4, we used the approach described as 'bootstrapping residuals' in Efron & Tibshirani (1993: 111).

First, bootstrap responses yi* were generated using

Second, the I-spline model was fitted to the new data (xi, yi*) giving bootstrap parameters αj* (see Fig. 3C). Third, one could predict


Collinearity between body measurements

The five predictors were strongly correlated (all crosscorrelation coefficients within the range of 0.83–0.99), and this collinearity prevented consideration of the potential of multivariate calibration (Brown 1982). Therefore, the problem was to select the best univariate predictor from the set of body measurements, based on measurement repeatability, comparability between rearing centres and quality of the modelling.

Reliability of body measurements

Measurement repeatability

In the pilot study which was done on 51 leverets, a paired t-test was performed on the data measured by each of the two observers. Skull length (SL) and skull width (SW) were the most promising variables (Table 1). Despite the fact that the two observers were experienced and had discussed the protocol for a long time, they failed to reproduce comparable measurements for hind foot (FL) and ear lengths (EL).

Rearing centre variability

The two rearing centres were compared using hares less than 30 days old. The differences in the growth curves obtained from the two rearing centres depended on each body measurement (Fig. 4). The covariance analysis (using cubic B-splines with df = 4) statistically evidenced this centre effect (Table 2). Nevertheless, the most comparable variable was skull length.

Table 1.

Repeatability of the body measurements obtained by the two observers (N = 51 leverets) assessed using paired t-test.


Figure 4.

Relationships between age of European hare leverets and the five body measurements (BM: body mass (A); SL: skull length (B); SW: maximum skull width (C); FL: hind foot length (D); EL: ear length; E)) and a cube root transformation of the body mass (F). Data were indicated by o and +, and B-spline fitting by continuous and dashed lines for the first and second rearing centre, respectively.


Growth curves modelling

The I-spline model was fitted to the five measurements using polynomial with k = 3 degrees and m = 2 interior knots. Every variable was well fitted by the versatile I-spline model, but the coefficient of determination was largest using skull length (R2SL = 0.93; R23√bm = 0-93; R2fl = 0.91; R2el = 0.88; R2SW = 0.79). Moreover, its residuals were more homoscedastic.

Considering the three criteria (measurement repeatability, comparability between the two rearing centres and goodness of fit by the I-spline model), the skull length was the variable of choice. So, the calibration intervals were defined only for this measurement.

Age calibration based on skull length

To compute the prediction intervals for skull length, we used a grid of 65 age values regularly spaced from day 1 to day 65. Each 95% prediction interval was defined using 1,000 bootstrap predicted values (Equation 4).

These intervals were inverted to provide calibration intervals (Fig. 5). It was clear from the plot that a good estimation of the age may be given for values of skull length within the range of 40–75 mm. For skull lengths larger than 80 mm, the calibration intervals were not superiorly bounded and the median curve became almost vertical. So estimating the age of leverets older than 40 days could be somewhat hazardous.

Table 2.

F-test of rearing centre effect in covariance analysis of B- splines curves. This test compares the sum of squared residuals between the common growth curve model (Equation 1) and the two growth curves model (Equation 2) expressing the differences between the two rearing centres.


Figure 5.

Allometric relationship between skull length and age with 95% calibration intervals.



Until now, age determination in young European hare leverets in the field has been based on the relationship between body mass and age. However, growth curves described in the literature differ greatly (see Bray 1998 for a review). According to the chosen growth curve the weight of a weanling (aged 1 month) ranged (see Fig. 1) from 250 g (Broekhuizen 1971) to 1,000 g (Hediger 1948). Moreover, individual variations were not considered as the growth curve only provided a mean value. Considering this variability, it seems difficult to select the best-fitting curve. Indeed, body mass was influenced by the animal's condition, and may vary with geographical situation (Cabon-Raczynska & Raczynski 1974, Broekhuizen & Maaskamp 1979, Morris 1972, Pepin 1979), as shown for other species (see Myers 1958, Dudzinski & Mykytowycz 1960 for wild rabbits Oryctolagus cuniculus). The birth weight and growth rates during the first month of life both varied with litter size (Biadi & Nouhaud 1975), and individual variations among leverets of the same litter were significant (Pielowski 1971).

Our results provide different relationships between body measurement and age together with confidence limits. Skull length was unambiguously the best variable as it lowered the effect of both observer and centre and resulted in a better fit.

Morris (1972) suggested that variables using rigid linear dimensions vary less according to season or nutritional status than other volumetric variables such as body mass. In the European hare, the relationships between linear measurement and age have been poorly documented. Broekhuizen & Maaskamp (1979) have provided a relationship between hind foot length and age for leverets less than five months old. The growth of the hind foot decreased rapidly before the age of two months, so the authors argued that it was only possible to distinguish adults from leverets less than two months old. Previously, Cabon-Raczynska & Raczynski (1964) showed allometric growth of hind feet and ears and emphasised that the growth of these dimensions decreased rapidly before the age of six months. Keith (1968) provided more detailed information on hind-foot/age relationship in captive-reared snowshoe hares, and Flux (1967) outlined collinearity between body mass and these variables in young hares. These studies have shown a strong growth of hind feet and ears during the first month of life, but a greater variability of curves using soft organs (Flux 1970). Our results also showed a lower measurement repeatability of the ear length.

In most of the five variables studied, we found significant differences between the two rearing centres: growth curves of body mass, skull length and hind foot length had steeper slopes in the second centre, contrary to the growth curve of the ear. Skull width varied according to rearing centre more than other variables, and the individual variation for a given age was greater. Although the growth of skull length was less influenced by the rearing centre effect than other variables, some variability due to the animals' origin, health status and rearing conditions still remained, and should be documented more completely. Some complementary sets of body measurements from other rearing centres would help us to assess more accurately the actual level of residual variability and, as a result, to define meaningful and robust age classes to account for such residual variance. Further experiments would confirm the applicability of our findings to wild leverets.

Field data on showshoe hares (Keith 1968) showed that wild hares grew more slowly during their first three months of life than individuals in captivity. Conversely, data collected by Pepin (1974) in France, showed that after the age of 40 days wild European hares were heavier than captive hares. This result is in accordance with Petrusewicz's (1970) findings in Poland. Because growth rates in wild young hares seem to be higher, models established in rearing centres could likely overestimate the actual age when used in the wild. Nevertheless, the growth curves using handreared leverets (symbolised with fi01_31.gif, fi02_31.gif and fi03_31.gif in Figure 1) were higher and, despite the possible sampling variations, seemed more stable in particular before the age of 30 days than curves based on wild leverets kept in captivity (Pilarska 1969 (fi04_31.gif), Broekhuizen 1971(fi05_31.gif)) or living under field conditions in large enclosures in Poland (Pielowski 1971(fi06_31.gif)). In this way, the first curves (circles) could be more accurate than the last (squares and triangular symbols) in age determination of wild leverets. Considering the variability in growth rates, our non-parametric model provides at least a robust relative ranking of leverets (according to meaningful age classes). Relative age determination is, indeed, biologically meaningful for a great number of studies, e.g. changes in survival or spacing behaviour of growing juveniles.

Figure 6.

Relationship between body mass (in g) and skull length (in mm) with dashed lines indicating the 95% prediction intervals (computed using the linear and homoscedastic relation between the cube root of body mass and skull length, and then back-transformed to the original scale).


Considering linear variables, only cranial measurements did not differ significantly between observers in our study. Skull length was the best candidate variable when considering goodness-of-fit and external validity (centre effect and repeatability).

From a practical point of view, Figure 5 is very easy to use as it allows a quick estimation of age in leverets from skull measurement. For example, a skull length of 70 mm corresponds to an age of 23 days (95% confidence limits ranged within 17–32 days). Calibration limits widen with age. They become unbounded when the skull length is above 80 mm. Therefore, the method should be limited to leverets less than 40 days old; for older leverets no confidence limits can be computed.

Using only skull length to estimate the age of leverets may be risky as measurement errors always occur. In order to detect such errors, one could also use the strong relationship between skull length and body mass (r = 0.99) to check compatibility between these two measurements (Fig. 6).


we are especially grateful to Eric Marboutin for improving the first draft of this manuscript. Sincere thanks go to Regis Péroux and the editorial team for their helpful comments. Thanks also to Annie Valla for her manuscript preparation.



Biadi, F. & Nouhaud, B. 1975: Élevage du liévre: bilan de six années d'élevage au Centre expérimental de Saint-Benoist. - Bulletin Mensuel de l'Office National de la Chasse, Spécial scientifique et technique 4: 69–126. (In French). Google Scholar


Bray, Y. 1998: Vers une meilleure connaissance des flux démographiques chez le lièvre dEurope (Lepus europaeus). - Mémoire de diplôme supérieur d'étude et de recherche, Université de Bourgogne, Dijon, France, 85 pp. (In French). Google Scholar


Broekhuizen, S. 1971: Age determination and age composition of hare populations. - Proceedings of the Xth International Union of Game Biologists (IUGB) Congress, Paris, France, 3–7 May 1971, pp. 477–489. Google Scholar


Broekhuizen, S. 1979: Survival in adult European Hares. - Acta Theriologica 24(34): 465–473. Google Scholar


Broekhuizen, S. & Maaskamp, F. 1979: Age determination in the European hare (Lepus europaeus Pallas) in the Netherlands. - Zeitschrift fur Saugetierkunde 44: 162–175. Google Scholar


Brown, P.J. 1982: Multivariate calibration. - Journal of the Royal Statistical Society B 44: 287–321. Google Scholar


Caboń-Raczyńska, K. & Raczyński, J. 1964: Studies on the European hare -II- Variations in the weight and dimensions of the body weight of certain internal organs. - Acta Theriologica 9(16): 233–248. Google Scholar


Caboń-Raczyńska, K. & Raczynski, J. 1974: Studies on the European hare - XXX - Variability of the body weight of European hares. - Acta Theriologica 19(5): 169–180. Google Scholar


Chambers, J.M. & Hastie, T.J. 1991: Statistical models in S. - Chapman and Hall, London, 608 pp. Google Scholar


de Boor, C. 1978: A practical guide to splines. - Springer-Verlag, New York, 386 pp. Google Scholar


Dudzinski, M.L. & Mykytowycz, R. 1960: Analysis of weights and growth rates of an experimental colony of wild rabbits Oryctolagus cuniculus. - CSIRO Wildlife Research 5: 102–115. Google Scholar


Efron, B. & Tibshirani, R.J. 1993: An introduction to the bootstrap. - Chapman and Hall, New York, 436 pp. Google Scholar


Emlen, J.M. & Pikitch, E. 1989: Animal population dynamics: identification of critical components. - Ecology Modelling 44: 253–273. Google Scholar


Flux, J.E.C. 1967: Reproduction and body weights of the hare, Lepus europaeus, in New-Zealand. - New Zealand Journal of Science 10: 357–401. Google Scholar


Flux, J.E.C. 1970: Life history of the Mountain Hare (Lepustimidus scoticus) in north-east Scotland. - Journal of Zoology 161: 75–123. Google Scholar


Frylestam, B. & von Schantz, T. 1977: Age determination of European hare based on periosteal growth lines. - Mammal Review 7: 151–154. Google Scholar


Hediger, H. 1948: Die zucht des Feldhasen (Lepus europaeus Pallas) in Gefangenschaft. - Physiologia Comparata Oecologia 1: 46–62. (In German). Google Scholar


Keith, L.B. 1968: Techniques for snowshoe hare population studies. - Journal of Wildlife Management 32(4): 801–834. Google Scholar


Lebreton, J.D. & Nichols, J.D. 1998: Ecology and management of subdivided animal populations: report on a workshop. - GibierFaune Sauvage Game Wildlife 15(1): 83–91. Google Scholar


Lloyd, H.G. 1968: Observation on breeding in the Brown hare (Lepus europaeus) during the first pregnancy of season. - Journal of Zoology 156: 521–528. Google Scholar


Marboutin, E. & Péroux, R. 1995: Survival pattern of European hare in a decreasing population. - Journal of Applied Ecology 32: 809–816. Google Scholar


Morris, P. 1972: A review of mammalian age determination methods. - Mammal Review 2: 69–104. Google Scholar


Myers, K. 1958: Further observations on the use of field enclosures for the study of the wild rabbits, Oryctolagus cuniculus (L.). - CSIRO Wildlife Research 3: 40–49. Google Scholar


O'Donoghue, M. 1994: Early survival of juvenile snowshoe hares. - Ecology 75(6): 1582–1592. Google Scholar


O'Donoghue, M. & Bergman, C.M. 1992: Early movements and dispersal of juvenile snowshoe hares. - Canadian Journal of Zoology 70: 1787–1791. Google Scholar


Pascal, M. & Kovacs, G. 1983: La détermination de 1′âge individuel chez le lièvre européen par la technique squelettochronologique. - Revue d'Ecologie (Terre Vie) 37:171–186. (In French). Google Scholar


Pépin, D. 1973: Recherche d'un critère d'âge chez le lièvre Lepus europaeus.- Annales de Zoologie - Écologie animale 5(2): 271–281. (In French). Google Scholar


Pépin, D. 1974: Mise au point de techniques pour l'étude de population de lièvre. - Bulletin Mensuel de l'Office National de la Chasse 2: 77–119. (In French). Google Scholar


Pépin, D. 1979: Body weight of hares in the Paris Basin (France). - In: Myers, K. & Mac Innes, C.D. (Eds.); Proceedings of a World Lagomorph Conference, Guelph, Ontario, Canada, pp. 229–238. Google Scholar


Petrusewicz, K. 1970: Dynamics and production of the hare population in Poland. - Acta Theriologica XV(26): 413–445. Google Scholar


Pielowski, Z. 1971: The individual growth curve of the hare. - Acta Theriologica XVI(5): 79–88. Google Scholar


Pilarska, J. 1969: Individual growth curve and food consumption by European hare in laboratory conditions. - Bulletin de l'académie polonaise des sciences Cl II, XVII(5): 299–305. Google Scholar


Ramsay, J.O. 1988: Monotone regression splines in action. - Statistical Sciences 3: 425–461. Google Scholar


Steams, S.C. 1992: The evolution of life histories. - Oxford University Press, Oxford, 249 pp. Google Scholar


Stroh, G. 1931: Zwei sichrer Altermerkmale beim Hasen. - Berliner Tierärztliche Wochenschrift 47(12): 180–181. (In German). Google Scholar


Suchentrunk, F., Willing, R. & Hartl, B.G. 1991: On eye lens weight and other age criteria of the Brown hare (Lepus europaeus Pallas, 1778). - Zeitschrift für Säugetierkunde 56: 365–374. Google Scholar


Walhovd, H. 1966: Reliability of age criteria for Danish hares (Lepus europaeus Pallas). - Danish Review of Game Biology 4: 106–128. Google Scholar
Yves Bray, Stephane Champely, and Dominique Soyez "Age determination in leverets of European hare Lepus europaeus based on body measurements," Wildlife Biology 8(1), 31-39, (1 March 2002).
Received: 16 November 2000; Accepted: 2 March 2001; Published: 1 March 2002

age criteria
European hare
Lepus europaeus
nonparametric modelling
Get copyright permission
Back to Top