calc_historical_probs: Calculate probabilities from an MCMC model to compare with...

In pacific-hake/hake-assessment: Stock Assessment for Pacific Hake

Calculate probabilities from an MCMC model to compare with estimates from past assessments

Description

Calculate probabilities from an MCMC model to compare with estimates from past assessments

Usage

calc_historical_probs(model, start = 2012, end)

Arguments

model	A model as loaded by create_rds_file()
start	First assessment year to do comparisons
end	Final assessment year to do comparisons

Details

Calculate a model's estimate of P(B_t+1 < B_t) and P(B_t+1 < B_40%) to then compare with the estimates of those quantities in year t's assessment, in combine_historical_probs().

Value

Data frame with columns

• year - year (lower case to check when binding columnsin combine_historical_probs()

• P_decline_curr - current estimate (from model) of the probability that the spawning biomass declined from year to year + 1

• P_below_B40_curr - current estimate (from model) of the probability that the spawning biomass was below B_40 in year + 1

pacific-hake/hake-assessment documentation built on Jan. 14, 2025, 9:12 p.m.