Run back-calculation in Bayesian framework

bcalc(data, ...)

Arguments

data

data.frame containing results from otolith reading plus information on the length of hatching (column l0p) of a certain species. Data should contain:

  • id: unique fish id per individual.

  • radi: measurements of otolith growth rings (in mm).

  • agei: age estimation of fish.

  • lencap: length at capture (in mm).

  • radcap: radius of otolith at capture (in mm).

  • l0p: length of fish at hatching (in mm).

...

arguments to add to sampling.

Details

Returns a list of model fit and a dataset called lengths with the estimated lengths at age:

  • l_m: mean length in mm.

  • l_sd: sd of length.

  • l_lb: lower 95% CI quantile.

  • l_ub: upper 95% CI quantile.

Input data should include radi at age 0 measurements as well, but can handle missing data (NA) for the cases where it is not possible to measure the radius at hatching.

Author

Nina M. D. Schiettekatte