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: 5ef1cd270ed31fd31f8b864b511024aec5a9a8f6 [formerly d16c9a1 [formerly d52e5bb6c3137e90cac96753ec81791ce74246f0] [formerly ccedaefc3b38bd99f8b497892ee51141267af7f0 [formerly 094daba021335dc5efcd158ae9f99ebd931c2f04]] [formerly a7ef685610e2281e5f21c3ae75cdf17b02ab56a1 [formerly 2ca6f6ee216df4c52b835e8eccc6fd1d7bda381c] [formerly 5ef1cd270ed31fd31f8b864b511024aec5a9a8f6 [formerly 1c9d08b]]]]
Former-commit-id: 6cb86b6
Former-commit-id: 3fbf388d9ebfb6e1d94f8f0fd2009ea127a4d07d [formerly d5d11a2607160bfd48d2db6847d55a4be4ca3470] [formerly f5ce05828667e7f5a89b0d3e6cdc367e23911e14 [formerly 7c6f24cf028e870fb3d813eceff14dfa7139a8cf]]
Former-commit-id: 71334b7ce4137d624a5b15f0441e238817e4bf9d [formerly 4f69e3e4ec28788b831f633ccd07b028abe7e6df]
Former-commit-id: cd83817379f76689427bdb39f43793ef594136c4
Former-commit-id: 5132b4d060524d26364414d772985d99b48b33fc [formerly 7a0244f4db1393ad559033cb5b281ea1d0e7ff32] [formerly 0bb903e985f042f874c31ce572261c9b2b5175ec [formerly c5e5c3e]]
Former-commit-id: d9ce8bc9ac6bbab27ba398c1eec59486fad9f4d7 [formerly 50a6056cd0a6c77d4a06034876f317dded85bd6e]
Former-commit-id: 2aebfc425f459f392f7824a247a072ddc422cec7
  • Loading branch information
kellijohnson-NOAA committed May 20, 2022
1 parent d392d14 commit a51a873
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 a51a873

Please sign in to comment.