Sexual dimorphism measures

Last updated

Although the subject of sexual dimorphism is not in itself controversial, the measures by which it is assessed differ widely. Most of the measures are used on the assumption that a random variable is considered so that probability distributions should be taken into account. In this review, a series of sexual dimorphism measures are discussed concerning both their definition and the probability law on which they are based. Most of them are sample functions, or statistics, which account for only partial characteristics, for example the mean or expected value, of the distribution involved. Further, the most widely used measure fails to incorporate an inferential support.

Contents

Introduction

It is widely known that sexual dimorphism is an important component of the morphological variation in biological populations (see, e.g., Klein and Cruz-Uribe, 1984; [1] Oxnard, 1987; [2] Kelley, 1993 [3] ). In higher Primates, sexual dimorphism is also related to some aspects of the social organization and behavior (Alexander et al., 1979; [4] Clutton-Brock, 1985 [5] ). Thus, it has been observed that the most dimorphic species tend to polygyny and a social organization based on male dominance, whereas in the less dimorphic species, monogamy and family groups are more common. Fleagle et al. (1980) [6] and Kay (1982), [7] on the other hand, have suggested that the behavior of extinct species can be inferred on the basis of sexual dimorphism and, e.g. Plavcan and van Schaick (1992) [8] think that sex differences in size among primate species reflect processes of an ecological and social nature. Some references on sexual dimorphism regarding human populations can be seen in Lovejoy (1981), [9] Borgognini Tarli and Repetto (1986) [10] and Kappelman (1996). [11]

These biological facts do not appear to be controversial. However, they are based on a series of different sexual dimorphism measures, or indices. Sexual dimorphism, in most works, is measured on the assumption that a random variable is being taken into account. This means that there is a law which accounts for the behavior of the whole set of values that compose the domain of the random variable, a law which is called distribution function. Because both studies of sexual dimorphism aim at establishing differences, in some random variable, between sexes and the behavior of the random variable is accounted for by its distribution function, it follows that a sexual dimorphism study should be equivalent to a study whose main purpose is to determine to what extent the two distribution functions - one per sex - overlap (see shaded area in Fig. 1, where two normal distributions are represented).

Fig. 1. Two normal distributions. Overlapnorm.jpg
Fig. 1. Two normal distributions.

Measures based on sample means

In Borgognini Tarli and Repetto (1986) an account of indices based on sample means can be seen. Perhaps, the most widely used is the quotient, [10]

where is the sample mean of one sex (e.g., male) and the corresponding mean of the other. Nonetheless, for instance,

have also been proposed.

Going over the works where these indices are used, the reader misses any reference to their parametric counterpart (see reference above). In other words, if we suppose that the quotient of two sample means is considered, no work can be found where, in order to make inferences, the way in which the quotient is used as a point estimate of

is discussed.

By assuming that differences between populations are the objective to analyze, when quotients of sample means are used it is important to point out that the only feature of these populations that seems to be interesting is the mean parameter. However, a population has also variance, as well as a shape which is defined by its distribution function (notice that, in general, this function depends on parameters such as means or variances).

Measures based on something more than sample means

Marini et al. (1999) [12] have illustrated that it is a good idea to consider something other than sample means when sexual dimorphism is analyzed. Possibly, the main reason is that the intrasexual variability influences both the manifestation of dimorphism and its interpretation.

Normal populations

Sample functions

It is likely that, within this type of indices, the one used the most is the well-known statistic with Student's t distribution see, for instance, Green, 1989. [13] Marini et al. (1999) [12] have observed that variability among females seems to be lower than among males, so that it appears advisable to use the form of the Student's t statistic with degrees of freedom given by the Welch-Satterthwaite approximation,

where are sample variances and sample sizes, respectively.

It is important to point out the following:

  • when this statistic is taken into account in sexual dimorphism studies, two normal populations are involved. From these populations two random samples are extracted, each one corresponding to a sex, and such samples are independent.
  • when inferences are considered, what we are testing by using this statistic is that the difference between two mathematical expectations is a hypothesized value, say

However, in sexual dimorphism analyses, it does not appear reasonably (see Ipiña and Durand, 2000 [14] ) to assume that two independent random samples have been selected. Rather on the contrary, when we sample we select some random observations - making up one sample - that sometimes correspond to one sex and sometimes to the other.

Taking parameters into account

Chakraborty and Majumder (1982) [15] have proposed an index of sexual dimorphism that is the overlapping area - to be precise, its complement - of two normal density functions (see Fig. 1). Therefore, it is a function of four parameters (expected values and variances, respectively), and takes the shape of the two normals into account. Inman and Bradley (1989) [16] have discussed this overlapping area as a measure to assess the distance between two normal densities.

Regarding inferences, Chakraborty and Majumder proposed a sample function constructed by considering the Laplace-DeMoivre's theorem (an application to binomial laws of the central limit theorem). According to these authors, the variance of such a statistic is,

where is the statistic, and (male, female) stand for the estimate of the probability of observing the measurement of an individual of the sex in some interval of the real line, and the sample size of the i sex, respectively. Notice that this implies that two independent random variables with binomial distributions have to be regarded. One of such variables is number of individuals of the f sex in a sample of size composed of individuals of the f sex, which seems nonsensical.

Mixture models

Authors such as Josephson et al. (1996) [17] believe that the two sexes to be analyzed form a single population with a probabilistic behavior denominated a mixture of two normal populations. Thus, if is a random variable which is normally distributed among the females of a population and likewise this variable is normally distributed among the males of the population, then,

is the density of the mixture with two normal components, where are the normal densities and the mixing proportions of both sexes, respectively. See an example in Fig. 2 where the thicker curve represents the mixture whereas the thinner curves are the functions.

Fig. 2. A mixture of two normal components. Overlapmix.jpg
Fig. 2. A mixture of two normal components.

It is from a population modelled like this that a random sample with individuals of both sexes can be selected. Note that on this sample tests which are based on the normal assumption cannot be applied since, in a mixture of two normal components, is not a normal density.

Josephson et al. limited themselves to considering two normal mixtures with the same component variances and mixing proportions. [17] As a consequence, their proposal to measure sexual dimorphism is the difference between the mean parameters of the two normals involved. In estimating these central parameters, the procedure used by Josephson et al. is the one of Pearson's moments. [17] Nowadays, the EM expectation maximization algorithm (see McLachlan and Basford, 1988 [18] ) and the MCMC Markov chain Monte Carlo Bayesian procedure (see Gilks et al., 1996 [19] ) are the two competitors for estimating mixture parameters.

Possibly the main difference between considering two independent normal populations and a mixture model of two normal components is in the mixing proportions, which is the same as saying that in the two independent normal population model the interaction between sexes is ignored. This, in turn implies that probabilistic properties change (see Ipiña and Durand, 2000 [14] ).

The MI measure

Ipiña and Durand (2000, [14] 2004 [20] ) have proposed a measure of sexual dimorphism called . This proposal computes the overlapping area between the and functions, which represent the contribution of each sex to the two normal components mixture (see shaded area in Fig. 2). Thus, can be written,

being the real line.

The smaller the overlapping area the greater the gap between the two functions and , in which case the sexual dimorphism is greater. Obviously, this index is a function of the five parameters that characterize a mixture of two normal components . Its range is in the interval , and the interested reader can see, in the work of the authors who proposed the index, the way in which an interval estimate is constructed.

Measures based on non-parametric methods

Marini et al. (1999) [12] have suggested the Kolmogorov-Smirnov distance as a measure of sexual dimorphism. The authors use the following form of the statistic,

with being sample cumulative distributions corresponding to two independent random samples.

Such a distance has the advantage of being applicable whatever the form of the random variable distributions concerned, yet they should be continuous. The use of this distance assumes that two populations are involved. Further, the Kolmogorov-Smirnov distance is a sample function whose aim is to test that the two samples under analysis have been selected from a single distribution. If one accepts the null hypothesis, then there is not sexual dimorphism; otherwise, there is.

See also

Related Research Articles

<span class="mw-page-title-main">Normal distribution</span> Probability distribution

In statistics, a normal distribution or Gaussian distribution is a type of continuous probability distribution for a real-valued random variable. The general form of its probability density function is

The weighted arithmetic mean is similar to an ordinary arithmetic mean, except that instead of each of the data points contributing equally to the final average, some data points contribute more than others. The notion of weighted mean plays a role in descriptive statistics and also occurs in a more general form in several other areas of mathematics.

<span class="mw-page-title-main">Multivariate normal distribution</span> Generalization of the one-dimensional normal distribution to higher dimensions

In probability theory and statistics, the multivariate normal distribution, multivariate Gaussian distribution, or joint normal distribution is a generalization of the one-dimensional (univariate) normal distribution to higher dimensions. One definition is that a random vector is said to be k-variate normally distributed if every linear combination of its k components has a univariate normal distribution. Its importance derives mainly from the multivariate central limit theorem. The multivariate normal distribution is often used to describe, at least approximately, any set of (possibly) correlated real-valued random variables each of which clusters around a mean value.

<span class="mw-page-title-main">Log-normal distribution</span> Probability distribution

In probability theory, a log-normal (or lognormal) distribution is a continuous probability distribution of a random variable whose logarithm is normally distributed. Thus, if the random variable X is log-normally distributed, then Y = ln(X) has a normal distribution. Equivalently, if Y has a normal distribution, then the exponential function of Y, X = exp(Y), has a log-normal distribution. A random variable which is log-normally distributed takes only positive real values. It is a convenient and useful model for measurements in exact and engineering sciences, as well as medicine, economics and other topics (e.g., energies, concentrations, lengths, prices of financial instruments, and other metrics).

Students <i>t</i>-distribution Probability distribution

In probability and statistics, Student's t-distribution is a continuous probability distribution that generalizes the standard normal distribution. Like the latter, it is symmetric around zero and bell-shaped.

In statistics, maximum likelihood estimation (MLE) is a method of estimating the parameters of an assumed probability distribution, given some observed data. This is achieved by maximizing a likelihood function so that, under the assumed statistical model, the observed data is most probable. The point in the parameter space that maximizes the likelihood function is called the maximum likelihood estimate. The logic of maximum likelihood is both intuitive and flexible, and as such the method has become a dominant means of statistical inference.

<span class="mw-page-title-main">Rayleigh distribution</span> Probability distribution

In probability theory and statistics, the Rayleigh distribution is a continuous probability distribution for nonnegative-valued random variables. Up to rescaling, it coincides with the chi distribution with two degrees of freedom. The distribution is named after Lord Rayleigh.

<span class="mw-page-title-main">Logistic distribution</span> Continuous probability distribution

In probability theory and statistics, the logistic distribution is a continuous probability distribution. Its cumulative distribution function is the logistic function, which appears in logistic regression and feedforward neural networks. It resembles the normal distribution in shape but has heavier tails. The logistic distribution is a special case of the Tukey lambda distribution.

In probability and statistics, a mixture distribution is the probability distribution of a random variable that is derived from a collection of other random variables as follows: first, a random variable is selected by chance from the collection according to given probabilities of selection, and then the value of the selected random variable is realized. The underlying random variables may be random real numbers, or they may be random vectors, in which case the mixture distribution is a multivariate distribution.

<span class="mw-page-title-main">Lévy distribution</span>

In probability theory and statistics, the Lévy distribution, named after Paul Lévy, is a continuous probability distribution for a non-negative random variable. In spectroscopy, this distribution, with frequency as the dependent variable, is known as a van der Waals profile. It is a special case of the inverse-gamma distribution. It is a stable distribution.

von Mises distribution Probability distribution on the circle

In probability theory and directional statistics, the von Mises distribution is a continuous probability distribution on the circle. It is a close approximation to the wrapped normal distribution, which is the circular analogue of the normal distribution. A freely diffusing angle on a circle is a wrapped normally distributed random variable with an unwrapped variance that grows linearly in time. On the other hand, the von Mises distribution is the stationary distribution of a drift and diffusion process on the circle in a harmonic potential, i.e. with a preferred orientation. The von Mises distribution is the maximum entropy distribution for circular data when the real and imaginary parts of the first circular moment are specified. The von Mises distribution is a special case of the von Mises–Fisher distribution on the N-dimensional sphere.

<span class="mw-page-title-main">Chi distribution</span> Probability distribution

In probability theory and statistics, the chi distribution is a continuous probability distribution over the non-negative real line. It is the distribution of the positive square root of a sum of squared independent Gaussian random variables. Equivalently, it is the distribution of the Euclidean distance between a multivariate Gaussian random variable and the origin. It is thus related to the chi-squared distribution by describing the distribution of the positive square roots of a variable obeying a chi-squared distribution.

In statistics, the method of moments is a method of estimation of population parameters. The same principle is used to derive higher moments like skewness and kurtosis.

<span class="mw-page-title-main">Inverse Gaussian distribution</span> Family of continuous probability distributions

In probability theory, the inverse Gaussian distribution is a two-parameter family of continuous probability distributions with support on (0,∞).

In estimation theory and decision theory, a Bayes estimator or a Bayes action is an estimator or decision rule that minimizes the posterior expected value of a loss function. Equivalently, it maximizes the posterior expectation of a utility function. An alternative way of formulating an estimator within Bayesian statistics is maximum a posteriori estimation.

<span class="mw-page-title-main">Dirichlet process</span> Family of stochastic processes

In probability theory, Dirichlet processes are a family of stochastic processes whose realizations are probability distributions. In other words, a Dirichlet process is a probability distribution whose range is itself a set of probability distributions. It is often used in Bayesian inference to describe the prior knowledge about the distribution of random variables—how likely it is that the random variables are distributed according to one or another particular distribution.

In probability theory and statistics, the normal-gamma distribution is a bivariate four-parameter family of continuous probability distributions. It is the conjugate prior of a normal distribution with unknown mean and precision.

The auxiliary particle filter is a particle filtering algorithm introduced by Pitt and Shephard in 1999 to improve some deficiencies of the sequential importance resampling (SIR) algorithm when dealing with tailed observation densities.

A product distribution is a probability distribution constructed as the distribution of the product of random variables having two other known distributions. Given two statistically independent random variables X and Y, the distribution of the random variable Z that is formed as the product is a product distribution.

In machine learning, the kernel embedding of distributions comprises a class of nonparametric methods in which a probability distribution is represented as an element of a reproducing kernel Hilbert space (RKHS). A generalization of the individual data-point feature mapping done in classical kernel methods, the embedding of distributions into infinite-dimensional feature spaces can preserve all of the statistical features of arbitrary distributions, while allowing one to compare and manipulate distributions using Hilbert space operations such as inner products, distances, projections, linear transformations, and spectral analysis. This learning framework is very general and can be applied to distributions over any space on which a sensible kernel function may be defined. For example, various kernels have been proposed for learning from data which are: vectors in , discrete classes/categories, strings, graphs/networks, images, time series, manifolds, dynamical systems, and other structured objects. The theory behind kernel embeddings of distributions has been primarily developed by Alex Smola, Le Song , Arthur Gretton, and Bernhard Schölkopf. A review of recent works on kernel embedding of distributions can be found in.

References

  1. Klei, Richard G.; Cruz-Uribe, Kathryn (1984). The Analysis of Animal Bones from Archaeological Sites. University of Chicago Press. ISBN   9780226439587.
  2. Oxnard, C.E. (1987). Fossils, Teeth and Sex: New Perspective in Human Evolution. University of Washington Press. ISBN   978-0295963891.
  3. Kelley, Jay (1993). "Taxonomic Implications of Sexual Dimorphism in Lufengpithecus". Species, Species Concepts and Primate Evolution. Boston, MA: Springer US. pp. 429–458. doi:10.1007/978-1-4899-3745-2_17. ISBN   978-1-4899-3747-6 . Retrieved 2023-07-27.
  4. Alexander, R. D.; Hoogland, J. L.; Howard, R. D.; Noonan, K. M.; Sherman, P. W. (1979). "Sexual dimorphism and breeding systems in pinnipeds, ungulates, primates and humans". In N. A. Chagnon; W. Irons; M. A. Scituate (eds.). Evolutionary Biology and Human Social Behavior: An Anthropological Perspective. Duxbury Press. pp. 402–435.
  5. Clutton-Brock, T. H. (1985). "Size, Sexual Dimorphism, and Polygyny in Primates". In Jungers, William L. (ed.). Size and Scaling in Primate Biology. Advances in Primatology. Boston, MA: Springer US. pp. 51–60. doi:10.1007/978-1-4899-3647-9_4. ISBN   978-1-4899-3647-9 . Retrieved 2023-07-27.
  6. Fleagle, John G.; Kay, Richard F.; Simons, Elwyn L. (Sep 1980). "Sexual dimorphism in early anthropoids". Nature . 287 (5780): 328–330. Bibcode:1980Natur.287..328F. doi:10.1038/287328a0. ISSN   1476-4687. PMID   6999362. S2CID   438852.
  7. Kay, Richard F. (1982-06-01). "Sivapithecus simonsi, a new species of miocene hominoid, with comments on the phylogenetic status of the ramapithecinae". International Journal of Primatology. 3 (2): 113–173. doi:10.1007/BF02693493. ISSN   1573-8604. S2CID   24825999.
  8. Plavcan, J. Michael; van Schaik, Carel P. (1992). "Intrasexual competition and canine dimorphism in anthropoid primates". American Journal of Physical Anthropology . 87 (4): 461–477. doi:10.1002/ajpa.1330870407. ISSN   0002-9483. PMID   1580353.
  9. Lovejoy, C. Owen (1981-01-23). "The Origin of Man". Science. 211 (4480): 341–350. Bibcode:1981Sci...211..341L. doi:10.1126/science.211.4480.341. ISSN   0036-8075. PMID   17748254.
  10. 1 2 Borgognini Tarli, S. M.; Repetto, E. (1986-02-01). "Methodological considerations on the study of sexual dimorphism in past human populations". Human Evolution. 1 (1): 51–66. doi:10.1007/BF02437285. ISSN   1824-310X. S2CID   85064651.
  11. Kappelman, John (1996). "The evolution of body mass and relative brain size in fossil hominids". Journal of Human Evolution. 30 (3): 243–276. doi:10.1006/jhev.1996.0021. ISSN   0047-2484.
  12. 1 2 3 Marini, Elisabetta; Racugno, Walter; Borgognini Tarli, Silvana M. (1999). <501::aid-ajpa6>3.0.co;2-7 "Univariate estimates of sexual dimorphism: The effects of intrasexual variability". American Journal of Physical Anthropology. 109 (4): 501–508. doi: 10.1002/(sici)1096-8644(199908)109:4<501::aid-ajpa6>3.0.co;2-7 . ISSN   0002-9483. PMID   10423265.
  13. Greene, David Lee (1989). "Comparison oft-tests for differences in sexual dimorphism between populations". American Journal of Physical Anthropology. 79 (1): 121–125. doi:10.1002/ajpa.1330790113. ISSN   0002-9483.
  14. 1 2 3 Ipiña, S (2000). "A Measure of Sexual Dimorphism in Populations which are Univariate Normal Mixtures". Bulletin of Mathematical Biology. 62 (5): 925–941. doi:10.1006/bulm.2000.0185. ISSN   0092-8240. PMID   11016090. S2CID   22840533.
  15. Chakraborty, Ranajit; Majumder, Partha P. (1982). "On Bennett's measure of sex dimorphism". American Journal of Physical Anthropology. 59 (3): 295–298. doi:10.1002/ajpa.1330590309. ISSN   0002-9483. PMID   7158663.
  16. Inman, Henry F.; Bradley, Edwin L. (1989). "The overlapping coefficient as a measure of agreement between probability distributions and point estimation of the overlap of two normal densities". Communications in Statistics - Theory and Methods. 18 (10): 3851–3874. doi:10.1080/03610928908830127. ISSN   0361-0926.
  17. 1 2 3 Josephson, Steven C.; Juell, Kenneth E.; Rogers, Alan R. (June 1996). <191::aid-ajpa3>3.0.co;2-0 "Estimating sexual dimorphism by method-of-moments". American Journal of Physical Anthropology. 100 (2): 191–206. doi:10.1002/(sici)1096-8644(199606)100:2<191::aid-ajpa3>3.0.co;2-0. ISSN   0002-9483. PMID   8771311.
  18. Lindsay, Bruce; McLachlan, G. L.; Basford, K. E.; Dekker, Marcel (Mar 1989). "Mixture Models: Inference and Applications to Clustering". Journal of the American Statistical Association. 84 (405): 337. doi:10.2307/2289892. ISSN   0162-1459. JSTOR   2289892. S2CID   119405289.
  19. "Full conditional distributions", Markov Chain Monte Carlo in Practice, Chapman and Hall/CRC, pp. 93–106, 1995-12-01, doi:10.1201/b14835-10, ISBN   978-0-429-17023-2 , retrieved 2023-07-27
  20. Ipiña, S (May 2004). "Inferential assessment of the MI index of sexual dimorphism: A comparative study with some other sexual dimorphism measures". Bulletin of Mathematical Biology. 66 (3): 505–522. doi:10.1016/j.bulm.2003.09.003. ISSN   0092-8240. PMID   15006446. S2CID   30298339.