Uses simulate_data()
, and fits a negative binomial gam model for response fish_count
using
predictor variables mark_stage
("Jack", "AD Adult", "UM Adult"), yearfac
(factor version of year), and a smooth across doy
(numeric for day of year). Output is suitable model for
plot_seasonal_gam_panels()
.
Examples
if (FALSE) { # \dontrun{
model = simulate_gam()
model_panels = plot_seasonal_gam_panels(model)
wrap_panels(model_panels)
} # }