Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Add boundaries argument to run_sdmtmb() and add new area to boundaries_data #32

Open
wants to merge 2 commits into
base: main
Choose a base branch
from

Conversation

iantaylor-NOAA
Copy link
Collaborator

What was changed/addressed

I don't know what the current workflow is, but for exploring the fraction of Petrale Sole north of 46 degrees 53.3 minutes, I found it helpful to be able to add a boundary argument to run_sdmtmb(). Perhaps this is unnecessary if you typically run the model once and then separately call calc_index_areas() after the initial estimation. If so, go ahead and delete this PR.

I also added new "North of 46deg 53.3min" and "South of 46deg 53.3min" areas to boundaries_data. If you'd rather not have that file include stuff like this, but think the other commit is useful, please cherry pick as you wish.

These changes allowed me to get everything I wanted in a single call to run_sdmtmb():

run_sdmtmb(
    data = data$data_filtered[[1]],
    formula = "catch_weight ~ 0 + fyear + pass_scaled",
    family = sdmTMB::delta_lognormal(),
    anisotropy = TRUE,
    n_knots = 500,
    boundaries = boundaries_data[c("Coastwide", "North of 46deg 53.3min", "South of 46deg 53.3min")],
    spatiotemporal = list("iid", "iid")
)

Related to issue(s)

No issue. Need for petrale calculation discussed by email.

- move Roxygen description to template and
- run devtools::document() which modified sanity_data.Rd
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

1 participant