Tabulates summary statistics for mean length at age from MCMC simulations using Reca.
MCMC simulations are typically obtained with RunRecaModels
.
Summary statistics are obtained from the posterior distribution, and
the interval is reported as equal-tailed credible intervals are reported.
Mean length for plus-groups are a weighted by the relative catch-at-age in each composite age group. For iterations where all of the plus-group ages have a zero catch at age, this weight is not defined, and summary statistics are obtained from the remaining iterations.
If 'RecaCatchAtAge' contains estimate for a set of aggregation variables, such as area, gear, stock, etc., summary statistics will be presented similarly.
Rounding of numbers according to the argument 'Decimals' is done with round
,
so that negative numbers specify rounding to powers of ten, and rounding of the digit 5 is towards the even digit.
The units considered valid for mean lengths are those listed for quantity 'length' in StoxUnits
Arguments
- RecaCatchAtAge
Results from MCMC simulations (
RecaCatchAtAge
).- PlusGroup
If given, ages 'PlusGroup' or older are included in a plus group.
- IntervalWidth
The width of the reported credible interval. A value of 0.9 gives 90 per cent credible intervals. Defaults to `r stoxFunctionAttributes$ReportRecaLengthAtAge$functionParameterDefaults$IntervalWidth`.
- Decimals
integer specifying the number of decimals to report for 'MeanIndividualLength', 'SD', 'Low' and 'High'. Defaults to `r stoxFunctionAttributes$ReportRecaLengthAtAge$functionParameterDefaults$Decimals`.
- Threshold
threshold for reporting mean weight. Rows with an estimated Catch At Age (number of individuals) lower than this will have NA reported for their mean length Defaults to `r stoxFunctionAttributes$ReportRecaLengthAtAge$functionParameterDefaults$Threshold`.
- Unit
unit for 'MeanIndividualLength', 'SD', 'Low' and 'High'
See also
RunRecaModels
for running Reca-analysis
Examples
lengthAtAge <- ReportRecaLengthAtAge(RstoxFDA::RecaCatchAtAgeExample,
PlusGroup = 13, Unit="cm", Decimals = 0)
lengthAtAge$MeanLengthByAge
#> AgeGroup Age MeanIndividualLength SD Low High
#> 1: Age 1 1 28 2 25 32
#> 2: Age 2 2 38 4 35 41
#> 3: Age 3 3 47 1 45 48
#> 4: Age 4 4 51 1 50 52
#> 5: Age 5 5 56 1 55 57
#> 6: Age 6 6 61 1 60 62
#> 7: Age 7 7 67 1 65 67
#> 8: Age 8 8 71 1 70 72
#> 9: Age 9 9 75 1 73 76
#> 10: Age 10 10 79 1 77 80
#> 11: Age 11 11 82 1 81 83
#> 12: Age 12 12 85 1 84 87
#> 13: Age 13+ 13 95 2 92 98