The ability to correctly sex marked birds is crucial to most behavioral or ecological studies (Greenwood, 1980; Andersson, 1994; Short and Balaban, 1994; Childress and Bennun, 2002; Barbraud et al., 2003) and for the management and conservation of species (Zavalaga and Paredes, 1997; Fernandez-Juricic, et al., 2009). However, sex determination in immature and adult individuals where the plumage is similar in both sexes based only on external morphological characters is often difficult (Cuthill et al., 1999).

To overcome this difficulty a range of techniques have been used such as laparoscopy (Petrides, 1950; Richter and Bourne, 1990; Richner, 1989), measuring the plasma testosterone levels during the breeding period (Czekala and Lasley, 1977; Bercovitz et al., 1978), vocalization analyses (Bourgeois et al., 2007), individual breeding or observation of territorial behavior (Castoro and Guhl, 1958; Flux and Innes, 2001; Fletcher and Hamer, 2003), and generalized molecular techniques (Griffiths et al., 1998; Bertault et al., 1999; Fridolfsson and Ellegren, 1999; Tomasulo et al., 2002; Dubiec and Zagalska-Neubauer, 2006; Balkız et al., 2007). However, despite the reliability and the large utilization of molecular techniques, these methods are time-consuming, intrusive and require extensive laboratory equipment, implying additional financial costs (Childress et al., 2005).

Sexing based on morphological measurements using multivariate statistical approaches (e.g., discriminant function analysis, DFA) is often a reasonable choice for quick and inexpensive but efficient sex identification in field studies on bird species presenting monomorphic plumage (Dechaume-Moncharmont et al., 2011). Indeed, the existence of a significant difference (even very small) between morphological measurements of males and females allows the discrimination between the sexes (Wilson, 1999; Donohue and Dufty, 2006; Moreno et al., 2007; Murphy, 2007; Ottvall and Gunnarsson, 2007; Cardoni et al., 2009). Diverse discriminant and logistic regression functions based on morphological measurements have been developed in many bird studies designed for widespread use by field researchers (Bosch, 1996; Balbontin et al., 2001; Bertellotti et al., 2002; Delvin et al., 2004; Figuerola et al., 2006, Hallgrimsson et al., 2011). However, different methods can be used to assess the reliability of these functions by examining the proportion of correctly classified males and females, such as, resubstitution (Zwarts et al., 1996; Copello et al., 2006), jackknife (O'Dwyer et al., 2006; Thorogood et al., 2009; Herring et al., 2008, 2010) or sample splitting methods (Setiawan et al., 2004; Meissner, 2005; Frey et al., 2008). Dechaume-Moncharmont et al. (2011) showed that the choice of the validation method may have a strong effect on the estimated discriminant rate, particularly for small sample sizes, and recommended the use of the jackknife method.

In this study, we used a large dataset (4013 birds measured across 11 years) of greater flamingo (Phoenicopterus roseus) chicks ringed and measured in the Camargue (southern France) and afterward sexed with molecular techniques (Bertault et al., 1999; Balkız et al., 2007). Greater flamingos lay a single egg and have a wide distribution, ranging from west Africa eastward throughout the Mediterranean to southwest and south Asia, and throughout sub-Saharan Africa (Johnson and Cézilly, 2007). Both sexes are similar in plumage but are sexually dimorphic in size when adults (Johnson and Cézilly, 2007). Although immature and adult male greater flamingos are significantly larger and heavier than females (Cramp and Simmons, 1977; Johnson et al., 1993), there remains a risk of error in sex determination based on size or weight due to the overlap of larger (and heavier) females with smaller (and lighter) males. This overlap is possibly wider in chicks, yet tarsus length of males appeared longer than that of females of exactly the same age in unfledged chicks from captivity (Studer-Thiersch, 1986).

The purposes of our study were to (1) examine whether DFA relying on morphological measurements of greater flamingo chicks predicts sex with a sufficient degree of accuracy, (2) test the applicability of one of the 11 single DFA formulas across cohorts and (3) evaluate through simulations the minimal sample size of chicks sexed by molecular methods to generate a DFA with a reliable discriminant rate.

As part of a long-term study on the reproductive biology of a colony of greater flamingos in Camargue, southern France (43∘25′ N, 4∘38′ E), a proportion of chicks were captured at the end of each breeding season (end of July–early August) and banded with a metal ring and a unique combination of plastic bands that allowed recognition of individuals. For each chick, only three external measurements were made (body weight, tarsus length and wing length; Childress et al., 2005) to avoid prolonged capture stress. Thus, between 1995 and 2008 (except for 2001, 2002 and 2007), a total of 4013 Flamingos chicks were captured. A total of 1828 males and 2187 females were genetically sexed through PCR (polymerase chain reaction) amplification of the CHD-Z and CHD-W genes, using DNA extracted from blood or feather samples (Bertault et al., 1999, Balkız et al., 2007).

In the greater flamingos, as in many other birds, adult males are larger than females. However, sexing flamingo chicks remains difficult because flamingos continue to grow after fledging, and captures generally include individuals of different ages. Sexual dimorphism occurs in flamingos chicks with females of 1.5 to 2.5 months being already smaller than males of the same age (Studer-Thiersch, 1986; but see Bertault et al., 2000). Our results demonstrate that the sex of greater flamingo chicks can be determined based on individual morphology using DFA and a subsample of individuals molecularly sexed. This method was also used for adult greater flamingos (Richter and Bourne, 1990) and other birds, where morphometric criteria may discriminate between males and females (Childress et al., 2005; Alarcos et al., 2007; Hurley et al., 2007; Ackerman et al., 2008; Herring et al., 2008; Lislevand et al., 2009; Herring et al., 2010).

However, we found that the reliability of one particular yearly discriminant function was reduced when attempting to apply it across cohorts. This result was consistent with the Evans et al. (1993) study on laughing gulls (Larus atricilla); they found a significant difference when DFA was applied at different years and localities. We also found that females were better discriminated than males in all years. The lower discriminant rate for males could result from an asymmetry in the distribution of the size of the males with an overrepresentation of smaller males. The discriminant rate of greater flamingo chicks varied from 67.1 % to 86.5 % and was lower than in other waterbird studies: Childress et al. (2005) on lesser flamingos Phœnicopterus minor (N=18 individuals, 94.0 %), Palomares et al. (1997) on black-headed gulls Larus ridibundus (N=143, 90.2 %) and Herring et al. (2008) on great egrets Ardea alba (N=76, 89.5 %). Furthermore, our simulation showed that for small samples, there is a wide variance in the proportion of misclassified birds. Indeed, a high or low discriminant rate can be obtained by chance (low or high misclassified); therefore the use of a formula constructed from a small dataset is problematic in the greater flamingos and should be avoided. Like previous studies, we recommend extreme caution when sexing birds based on DFA, particularly when discriminant equations are derived from small datasets (Brennan et al., 1991; Shealer and Cleary, 2007; Isaksson et al., 2008; Dechaume-Moncharmont et al., 2011). Our results show that sexing flamingo chicks based on DFA requires a large dataset and repeated sampling to include natural yearly variations on body condition during the breeding season. We therefore conclude that DFA is not an effective method to sex flamingo chicks.