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]")