get_expanded: get_expanded is a helper function to return a matrix of...

In ericward-noaa/bycatch: Fit Bayesian Models to Bycatch Data

get_expanded is a helper function to return a matrix of posterior predictive values for unobserved bycatch

Description

get_expanded is a helper function to return a matrix of posterior predictive values for unobserved bycatch

Usage

get_expanded(fitted_model)

Arguments

fitted_model

Data and fitted model returned from estimation

Value

matrix (MCMC draws x time steps) of posterior predictive values for unobserved bycatch

Examples

d <- data.frame(
  "Year" = 2002:2014,
  "Takes" = c(0, 0, 0, 0, 0, 0, 0, 0, 1, 3, 0, 0, 0),
  "expansionRate" = c(24, 22, 14, 32, 28, 25, 30, 7, 26, 21, 22, 23, 27),
  "Sets" = c(391, 340, 330, 660, 470, 500, 330, 287, 756, 673, 532, 351, 486)
)
fit <- fit_bycatch(Takes ~ 1,
  data = d, time = "Year",
  effort = "Sets",
  family = "poisson",
  expansion_rate = "expansionRate",
  time_varying = FALSE
)
expanded <- get_expanded(fit)

ericward-noaa/bycatch documentation built on July 19, 2023, 6:10 p.m.