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

[Feature]: Log scale outputs #547

Open
e-perl-NOAA opened this issue Jan 2, 2024 · 2 comments
Open

[Feature]: Log scale outputs #547

e-perl-NOAA opened this issue Jan 2, 2024 · 2 comments
Assignees
Labels
misc. output statistics related to logL wishlist request new feature; bigger than revision; OK to remove after adding to Milestone

Comments

@e-perl-NOAA
Copy link
Collaborator

e-perl-NOAA commented Jan 2, 2024

Describe the solution you would like.

Producing log scale outputs for approaches like MVLN is a feature that some users would find useful. See
APPROXIMATION OF KOBE POSTERIORS FROM STOCK SYNTHESIS FOR NORTH ATLANTIC BLUE SHARK and ss3diags use of MVLN (download to view plots or go back to Vignettes folder and go to the .rmd version of ss3diags_handbook to view markdown version).

The requested option already exists. See below from a starter file:

102 # Depletion basis:  denom is: 0=skip; 1=X*SPBvirgin; 2=X*SPBmsy; 3=X*SPB_styr; 4=X*SPB_endyr; 5=X*dyn_Bzero;  values>=11 invoke N multiyr (up to 9!) with 10's digit; >100 invokes log(ratio)
1 # Fraction (X) for Depletion denominator (e.g. 0.4)
1 # SPR_report_basis:  0=skip; 1=(1-SPR)/(1-SPR_tgt); 2=(1-SPR)/(1-SPR_MSY); 3=(1-SPR)/(1-SPR_Btarget); 4=rawSPR
3 # F_std_reporting_units: 0=skip; 1=exploitation(Bio); 2=exploitation(Num); 3=sum(Apical_F's); 4=mean F for range of ages (numbers weighted); 5=unweighted mean F for range of ages
#COND 10 15 #_min and max age over which mean F will be calculated with F_reporting=4 or 5
102 # F_std_scaling: 0=no scaling; 1=F/Fspr; 2=F/Fmsy; 3=F/Fbtgt; where F means annual F_std, Fmsy means F_std@msy; values >=11 invoke N multiyr (up to 9!) using 10's digit; >100 invokes log(ratio)

This request also was in #79 . It was implemented in ~March 2022 probably with 3.30.19. However, it was not added to the change log or issue tracker.

@e-perl-NOAA Please check the User Manual to see if this in documented in the starter.ss section.

Describe alternatives you have considered

Having users approximate/calculate log scale outputs from normal outputs.

Statistical validity, if applicable

A comparison of MVLN to MCMC seems like an avenue people are looking into.

Describe if this is needed for a management application

No.

Additional context

No response

@e-perl-NOAA e-perl-NOAA added wishlist request new feature; bigger than revision; OK to remove after adding to Milestone misc. output labels Jan 2, 2024
@Rick-Methot-NOAA Rick-Methot-NOAA added the statistics related to logL label Jan 3, 2024
@Henning-Winker
Copy link

I note that the existing features:

102 # Depletion basis: denom is: >100 invokes log(ratio)
102 # F_std_scaling: >100 invokes log(ratio)

Can be readily implemented for the production of Kobe plots and Kobe-2-strategy matrices in tuna RFMOs, where log(SSB/SSBmsy) and log(F/Fmsy) tend to be quantities of interest (but not always), assuming these are the log(ratio).

To gain full flexibility for producing decision tables, all possible settings for Depletion basis, F_std_scaling as well as log(SSB), log(SSB0), log(R0), log(F_Btgt) should be generated on log-scale.

This is noting that

0 # F_std_scaling: 0=no scaling

Is the preferred option as it allows to work out the covariance for any F target (Fspr, Fsb, Fmsy, F0.1), see:

https://github.com/jabbamodel/ss3diags/blob/master/Vignette/ss3diags_handbook.pdf

@e-perl-NOAA
Copy link
Collaborator Author

@Rick-Methot-NOAA In the section on the starter file, on pages 15 and 17 it describes those options
SS330_User_Manual.pdf

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
misc. output statistics related to logL wishlist request new feature; bigger than revision; OK to remove after adding to Milestone
Projects
None yet
Development

No branches or pull requests

3 participants