Data used to demonstrate size-structured modelling including bottom-up
linkages.
data(gulf_of_alaska)
loads a list that includes four components:
years
the set of years modeledtaxa
the set of taxa (including each stanza for age-structured populations)type
is a character vector indicating for each taxon (as vector names) whether it is a heterotroph, autotroph, or detritus poolbiomass_data
is a long-form data-frame with three columns, providing the Year, Mass (in relative units for copepods and euphausiids, and million metric tons for Pollock and sablefish), and Taxon for each year with available datacatch_data
is a long-form data-frame with three columns, providing the Year, Mass (in million metric tons), and Taxon for each year with available dataP_over_B
is a numeric vector with the unitless ratio of biomass production to population biomass for each taxon (as vector names)Q_over_B
is a numeric vector with the unitless ratio of biomass consumption to population biomass for each taxon (as vector names)EE
is a numeric vector with the ecotrophic efficiency for eachtaxa
, in this case provided only for the detritus pool and determining its turn-over rateDiet_proportions
is a numeric matrix where each column lists the proportion of biomass consumed for each predator (column) that is provided by each prey (row)agecomp_data
is a named list with age-composition samples for any age-structured population (named list element), where each list element is a matrix with named rows indicating years, and columns indicating integer agesstanza_groups
is a character vector with names matching taxa, and elements indicating multi-stanza groups (i.e., age-structured populations)Amax
is a numeric vector indicating the maximum age (in years) for each stanza as indicated innames(stanza_groups)
.Leading
is a Boolean vector indicating which stanza is used for estimating biomass (and therefore recruitment) for eachnames(stanza_groups)
.K
is a numeric vector with the generalized von Bertalanffy growth coefficient (in weight not length) for eachunique(stanza_groups)
d
is a numeric vector with the generalized von Bertalanffy growth coefficient (in weight not length) for eachunique(stanza_groups)
Amat
is a numeric vector with the age at 50% maturity (in years) for eachunique(stanza_groups)
. If provided, it is used to back-calculate and replace the value forWmat
Wmat
is a numeric vector with the weight at 50% maturity (relative to W_inf) for eachunique(stanza_groups)
.Wmatslope
is a numeric vector with the slope of the logistic maturity at weight (relative to W_inf) for eachunique(stanza_groups)
Usage
data(gulf_of_alaska)
Details
The data compiled come from a variety of sources:
Sablefish biomass index and age-comps are from a cooperative longline survey, and the index is treated as a relative index.
Pollock biomass index and age-comps are from a bottom trawl survey, and the index is treated as a relative index with a prior on the catchability coefficient.
Copepod biomass is from an oblique tow bongo net survey, with data provided by Dave Kimmel.
Euphausiid biomass is from the Seward Long-Term Ecological Research program, provided by Russ Hopcroft to the Ecological Status Report.
Other P_over_B, Q_over_B and Diet_proportions values are derived from Rpath models, provided by Bia Dias.
Age-structed inputs are determined based on life-history information and biological samples. See Thorson et al. (In review) for more details regarding data standardization and sources