maturity | R Documentation |
Determines whether an organism is sexually mature, based on qualitative or morphometric observational data.
is.mature(x, ...)
## S3 method for class 'scsbio'
is.mature(x, probability = FALSE, ...)
is.pubescent(x, ...)
## S3 method for class 'scsbio'
is.pubescent(x, ...)
is.adolescent(x, ...)
## S3 method for class 'scsbio'
is.adolescent(x, ...)
is.primiparous(x, ...)
## S3 method for class 'scsbio'
is.primiparous(x, ...)
is.multiparous(x, ...)
## S3 method for class 'scsbio'
is.multiparous(x, ...)
is.senile(x, ...)
## S3 method for class 'scsbio'
is.senile(x, ...)
maturity(x, ...)
## S3 method for class 'scsbio'
maturity(x, ...)
## S3 method for class 'numeric'
maturity(x, year, species, ...)
x |
Object or numeric maturity code. |
... |
Further arguments passed onto internal functions. |
probability |
Logical value specifying whether values are to be returned as probabilities of being mature. For individuals with missing morphometric data or observational sexual characters, a GAM binary regression is performed on the observations by year and the predicted probabilties are returned. |
year |
Survey year. |
species |
Species name or code. |
Logical vector, probability (when probability = TRUE
), or character string.
is.mature
Generic is.mature
function.
is.mature.scsbio
Snow crab biological data (scsbio
) maturity method.
is.primiparous
Generic is.primiparous
function.
is.primiparous.scsbio
Returns whether a female is primiparous, i.e. a first-time spawner.
is.multiparous
Generic is.multiparous
function.
is.multiparous.scsbio
Returns whether a female is multiparous, i.e. has previously spawned.
is.senile
Generic is.senile
function.
is.senile.scsbio
Returns whether an animal is senile, i.e. has decreased spawning capability because of old age.
maturity
Generic weight
method.
maturity.default
Return sexual maturity status.
maturity.numeric
Convert numeric codes to descriptive maturity character strings.
x <- read.scsbio(2010)
index <- is.mature(x)
plot(x$carapace.width[index], x$chela.height[index], xlim = c(0, 140), ylim = c(0, 40),
xlab = "Carapace width (mm)", ylab = "Chela height (mm)", type = "n")
points(x$carapace.width[index], x$chela.height[index], pch = 21, bg = "grey", cex = 0.4)
points(x$carapace.width[!index], x$chela.height[!index], pch = 21, bg = "red", cex = 0.4)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.