library(sdmTMB)
pcod_mod <-
sdmTMB(
formula = density ~ 1,
data = pcod,
spatial = "off",
spatiotemporal = "off",
time = "year",
# family = delta_gamma(),
family = tweedie(),
extra_time = 2018:2020,
silent = FALSE)
# predict(pcod_mod) # broken, will fix
predict(pcod_mod, newdata = pcod) # this will work