Performs data checks and data conversions,
and stores some data-related parameters in preparation for running
ParameterizeRecaModels
.
Usage
PrepareRecaEstimate(
StoxBioticData,
StoxLandingData,
FixedEffects = character(),
RandomEffects = character(),
UseCarEffect = FALSE,
CarEffect = character(),
CarNeighbours,
UseAgingError = FALSE,
AgeErrorMatrix,
UseStockSplitting = FALSE,
UseStockSplittingError = FALSE,
StockSplittingParameters,
CellEffect = c("Off", "All"),
MinAge = integer(),
MaxAge = integer(),
MaxLength = numeric(),
LengthResolution = numeric(),
HatchDay = integer()
)
Arguments
- StoxBioticData
StoxBioticData
data with samples from fisheries and approriate columns appended for specifying cells. See details.- StoxLandingData
StoxLandingData
data with landings from fisheries and approriate columns appended for specifying cells. See details.- FixedEffects
optional, vector identifying column names that should be treated as fixed effects. See details
- RandomEffects
optional, vector identifying column names that should be treated as random effects. See details
- UseCarEffect
if TRUE, CAR-effect (conditional autoregressive effect) will be used.
- CarEffect
optional, character identifying the column name that should be treated as CAR-effect. Mandatory if 'UseCarEffect' is TRUE.
- CarNeighbours
CarNeighbours
, mandatory if 'carEffect' is given. Identifies which values of the carEffect are to be considered as neighbours.- UseAgingError
If TRUE, aging, error parameters will be incorporated in the models.
- AgeErrorMatrix
AgeErrorMatrix
, optional, specifies the probabilities of misreading ages. mandatory if UseAgingError is TRUE.- UseStockSplitting
If TRUE, models will be condigured to provide estimates for each of two stocks based on otholitt-typing. See
StockSplittingParameters
.- UseStockSplittingError
If TRUE, the model for error in stock classification (otolithtypes) will be applied.
- StockSplittingParameters
Parameters for stock splitting. Mandatory if 'UseStockSplitting' is TRUE. May be obtained with
DefineStockSplittingParameters
.- CellEffect
Configures the cell effect. Defaults to Off. If 'All', an interaction term will be added with all covariates that in the cell (whether they are fixed or random effects). Any CAR-effect is always included in the cell effect.
- MinAge
optional, must match dimensions of any 'AgeErrorMatrix'. If not provided it will be derived from data.
- MaxAge
optional, must match dimensions of any 'AgeErrorMatrix'. If not provided it will be derived from data.
- MaxLength
optional, maximal fish length in data in cm. If not provided it will be derived from data.
- LengthResolution
optional, resolution for length measurements in cm. If not provided modal value from data is used.
- HatchDay
defaults to Defaults to 1. encoding the day of the year when fish is consider to transition from one age to the next. 1 represents Jan 1st.
Value
RecaData
Data prepared for running Reca.
Details
Parameters are obtained for 'cells' using ParameterizeRecaModels
and and applied to landings-data when predicting total catch-at-age using RunRecaModels
.
Cells
Reca conceptually decomposes a fishery into cells. Parameters are estimated for each cell, and later total catch at age is predicted for each cell. Cells must therefore be defined for both the census of landings, and for the samples.
'Cells' are defined by all effects / covariates that are specified for both 'StoxLandingData' and 'StoxBioticData' (all effects are specified by one of the arguments: 'FixedEffects', 'RandomEffects', and 'CarEffect'). Any variable that are to be included in the cell definition must always be observed (not containing NAs) and be coded coherently (the same coding system) in both 'StoxBioticData' and 'StoxLandingData'.
Covariates
Only effects that are part of the cell definition may be specified as 'Fixed effect' or 'CAR effect'. For fixed effect all values / levels of the effect must be sampled. For CAR effect, the cell or neighbouring cells (as defined by 'CarNeighbours') must be sampled
Any variable in samples 'StoxBioticData' may be included as a random effect as long as it is always observed (not containing NAs).
A random effect identifying haul / landing is always included in Reca. Do not add haul-identifiers as covariates.
Stock splitting
The option UseStockSplitting requires that a Variable 'otolithtype' is added to the table 'Individual' in StoxBioticData. This is developed for splitting coastal cod from atlantic cod, and that is reflected in the Reca documentation, but this is abstracted in StoX, so that other stocks may be encoded to fit the analysis. For some cases it may even be extended to other domains defined for individual fish.
See also
ReportFdaSampling
for inspecting the data availability for potential
cell definitions.
ParameterizeRecaModels
,
and RunRecaModels
for parameterising and running the Reca-analysis.
DefineStratumPolygon
,
AddStratumStoxBiotic
,
and AddStratumStoxLanding
for making 'Stratum' part of the cell definition,
and #' DefineCarNeighbours
, for maing 'Stratum' available as a CAR-effect
DefinePeriod
,
AddPeriodStoxBiotic
,
and AddPeriodStoxLanding
for making 'Period' part of the cell definition.
#' DefineTranslation
,
AddGearGroupStoxBiotic
,
and AddGearGroupStoxLanding
for making 'GearGroup' part of the cell definition.
DefineAgeErrorMatrix
for making an age-error matrix available for the Reca-analysis.
DefineStockSplittingParameters
for configuring stock-splitting parameters,
and AddToStoxBiotic
for adding otolith-type to samples for stock-splitting analysis.