Skip to content

Commit

Permalink
Changes upper bounds of logkappa2
Browse files Browse the repository at this point in the history
using VAST_3.9.0 the upper bound of logkappa2 is reached
using the default settings b/c of the small number of mesh points.
This commit changes the upper bound to 10 using
run_model = FALSE to get the correct map.
This trick can be found in
James-Thorson-NOAA/VAST#300
which was started by @kellijohnson-NOAA b/c of similar issues
when running the model for West Coast species.

Former-commit-id: 1c9d08b
  • Loading branch information
kellijohnson-NOAA committed May 20, 2022
1 parent afdc4ed commit 6cb86b6
Showing 1 changed file with 15 additions and 1 deletion.
16 changes: 15 additions & 1 deletion sandbox/pcod_delta_gamma.R
Original file line number Diff line number Diff line change
Expand Up @@ -67,7 +67,21 @@ fit_vast <- fit_model(
b_i = pcod[, "density"],
a_i = pcod[, "effort"],
input_grid = input_grid,
working_dir = paste0(here("doc", "appendix-VAST"), "/")
working_dir = paste0(here("doc", "appendix-VAST"), "/"),
run_model = FALSE
)
upperbounds <- fit_vast$Upper
upperbounds["logkappa2"] <- 10
fit_vast <- fit_model(
settings = settings,
Lat_i = pcod[, "lat"],
Lon_i = pcod[, "lon"],
t_i = pcod[, "year"],
b_i = pcod[, "density"],
a_i = pcod[, "effort"],
input_grid = input_grid,
working_dir = paste0(here("doc", "appendix-VAST"), "/"),
upper = upperbounds
)
# saveRDS(fit_vast, file = f)
# } else {
Expand Down

0 comments on commit 6cb86b6

Please sign in to comment.