Calculate an estimate of total sexual shape dimorphism using a linear model

Usage

SShDLM(Coords, Sex, zeroed = T, stat = "median", iter = 1000)

Arguments

Coords: Two-dimensional array of coordinates from geometric morphometric analysis (see two.d.array() in geomorph for correct formatting)
Sex: A character or factor vector recording sex for each individual as 'm' or 'f'
zeroed: a logical value stating whether you wish to correct the SShD for a natural zero (defaults to TRUE)
stat: A character value specifying the centering statistic to be used if zeroed is set to TRUE (currently supports the median (default) and minimum)
iter: a numeric value specifying the number of iterations the null SShD correction factor should be permuted for (default is 1000)

Examples

data(pupfish)
with(pupfish, SShDLM(coords, Sex, zeroed = F))
#> Warning: No data found for at least one sex. Sex variable must be formatted as 'm' and 'f' (see ConvertSexFormat)
#> [1] NA