Get growth curves giving weight as a function of age
getGrowthCurves(object, species = NULL, max_age = 20, percentage = FALSE)Arguments
- object
MizerSim or MizerParams object. If given a MizerSim object, uses the growth rates at the final time of a simulation to calculate the size at age. If given a MizerParams object, uses the initial growth rates instead.
- species
The species to be selected. Optional. By default all target species are selected. A vector of species names, or a numeric vector with the species indices, or a logical vector indicating for each species whether it is to be selected (TRUE) or not.
- max_age
The age up to which to run the growth curve. Default is 20.
- percentage
Boolean value. If TRUE, the size is given as a percentage of the maximal size.
Value
An array (species x age) containing the weight in grams.
See also
Other summary functions:
getBiomass(),
getDiet(),
getN(),
getSSB(),
getYieldGear(),
getYield()
Examples
growth_curves <- getGrowthCurves(NS_params, species = c("Cod", "Haddock"))
str(growth_curves)
#> num [1:2, 1:50] 0.001 0.001 0.87 2.546 4.499 ...
#> - attr(*, "dimnames")=List of 2
#> ..$ Species: chr [1:2] "Haddock" "Cod"
#> ..$ Age : chr [1:50] "0" "0.408163265306122" "0.816326530612245" "1.22448979591837" ...
library(ggplot2)
ggplot(melt(growth_curves)) +
geom_line(aes(Age, value)) +
facet_wrap(~ Species, scales = "free") +
ylab("Size[g]") + xlab("Age[years]")
