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

compare changes in likelihoods by fleet? #148

Open
iantaylor-NOAA opened this issue Jul 14, 2021 · 3 comments
Open

compare changes in likelihoods by fleet? #148

iantaylor-NOAA opened this issue Jul 14, 2021 · 3 comments

Comments

@iantaylor-NOAA
Copy link
Collaborator

The r4ss::SStableComparisons() function is useful for producing the tables showing the changes in likelihood by data type that we've shown so far, but it would often be useful to look at the finer scale changes in likelihood by fleet (the table equivalent of Piner Plot vs the main profile plot showing the sum across all fleets).

The list returned by r4ss::SSsummarize() function (as used in https://github.com/iantaylor-NOAA/Lingcod_2021/blob/main/models/model_STAR_requests.R#L31-L37 with model files on Google Drive) includes a data.frame likelihoods_by_fleet but the format is terrible (it's just rbinding the likelihoods_by_fleet table from SS_output):

     model           Label          ALL X1_Comm_Trawl  X2_Comm_Fix    X3_Rec_WA
194      1    Catch_lambda           NA   1.00000e+00  1.00000e+00  1.00000e+00
195      1      Catch_like  6.99642e-10   5.80240e-10  4.13638e-11  2.12453e-11
196      1 Init_equ_lambda           NA   1.00000e+00  1.00000e+00  1.00000e+00
197      1   Init_equ_like  0.00000e+00   0.00000e+00  0.00000e+00  0.00000e+00
198      1     Surv_lambda           NA   1.00000e+00  1.00000e+00  1.00000e+00
199      1       Surv_like -1.62875e+02  -1.61206e+01 -1.48465e+01 -5.23057e+01
200      1      Surv_N_use           NA   1.70000e+01  1.70000e+01  4.00000e+01
201      1     Surv_N_skip           NA   0.00000e+00  0.00000e+00  0.00000e+00

@kellijohnson-NOAA, if you have time at some point, could you work some tidyverse magic to provide a table of changes in likelihood for each fleet for the Survey (really index), Length, and Age data. For these models I think we can ignore all the rows with labels "lambda", "N_use", "N_skip").

This would help us understand how something like removing fixed-gear ages improves the fit to the fixed-gear length as opposed to all length data.

@iantaylor-NOAA iantaylor-NOAA added status: help wanted Extra attention is needed topic: package function topic: presentation Related to presentation labels Jul 14, 2021
@kellijohnson-NOAA
Copy link
Collaborator

I can look at this now

@kellijohnson-NOAA kellijohnson-NOAA added this to the 2021assessment milestone Nov 16, 2021
@kellijohnson-NOAA
Copy link
Collaborator

@iantaylor-NOAA do we need this for the document or can we label it as an r4ss post mortem issue?

@kellijohnson-NOAA kellijohnson-NOAA self-assigned this Nov 16, 2021
@iantaylor-NOAA iantaylor-NOAA added postmortem: r4ss and removed status: help wanted Extra attention is needed topic: presentation Related to presentation labels Nov 16, 2021
@iantaylor-NOAA
Copy link
Collaborator Author

Labeled "postmortem: r4ss" and added to project

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

No branches or pull requests

2 participants