| 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.matureGeneric is.mature function.
is.mature.scsbioSnow crab biological data (scsbio) maturity method.
is.primiparousGeneric is.primiparous function.
is.primiparous.scsbioReturns whether a female is primiparous, i.e. a first-time spawner.
is.multiparousGeneric is.multiparous function.
is.multiparous.scsbioReturns whether a female is multiparous, i.e. has previously spawned.
is.senileGeneric is.senile function.
is.senile.scsbioReturns whether an animal is senile, i.e. has decreased spawning capability because of old age.
maturityGeneric weight method.
maturity.defaultReturn sexual maturity status.
maturity.numericConvert 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.