Forums

Back

SS 3.30.13.00 released Mar 13, 2019

Richard Methot, modified 4 Years ago.

SS 3.30.13.00 released Mar 13, 2019

Youngling Posts: 219 Join Date: 11/24/14 Recent Posts

Latest SS News:  March 13, 2019

Greetings Stock Synthesis (SS) Users,

In this release of Stock Synthesis 3.30.13.00 we have included some fixes, augmentations and new features.  Thanks to some users who were pushing the envelope on time-varying parameters for finding some needed fixes to the code.  You can find all the changes listed in the Change Log and below.   The fixes were minor and mostly in the output, but there is one big fix affecting the calculation of the log likelihood for the parameter deviations.  We will put out a fuller report on that later and encourage you to contact us if you need info on this change sooner.

Please remember to check the change log frequently so you are aware of anything relevant for your SS project.  It also is a good idea to look over the *.ss_new files created by the newest version of SS so you can see listed all available input options.  Of course, the user manual, r4ss, and graphical Stock Synthesis Interface (SSI) have been updated.   The SSI now provides previews for the shape of most selectivity functions.

We have a new SS team member, Kathryn Doering.  She will be involved with VLAB maintenance, creation of an example library, and developing more R-based routines to enhance SS.  Welcome Kathryn.

If you have any SS usage questions or issues please post in the FORUM.

Problems using the SS VLAB site can be addressed to nmfs.stock.synthesis@noaa.gov

Please contact VLAB administrators at vlab.admin@noaa.gov to change your username or password. 

 

ADMB 12.0

SS is now being compiled exclusively with ADMB 12.0 (except the 32 bit versions).  An important new feature in ADMB 12.0 improves the MCMC algorithm.  See below from the admb-project.org website:

  • Dr. Cole Monnahan has developed a new MCMC functionality into ADMB. The current Metropolis algorithm was updated to have more meaningful console output and not require an sdreport object. The old method remains, and the new one is accessible with argument ‘-rwm’. A new algorithm, the no-U-turn sampler (a self-tuning extension of the existing ‘hybrid’ method was also added, and available with ‘-nuts’. The functionality matches the software platform Stan closely. The user can specify whether to use the MLE covariance, an arbitrary matrix, or adapt one during the warmup phase (recommended). Both these additions work much with the R package adnuts, which streamlines workflow, including parallel chains, diagnostic checking, and inference. See adnuts package vignette for more information.

 

Future big features on the wish list:

  • Option to replace input of growth parameters with input of empirical length-at-age
  • Change input of tagging parameters to use a much more concise design matrix approach
  • Adding otolith weight as a proxy for fish age as a new data type
  • Area-specific spawner-recruitment
  • Area-specific growth in addition to morph specific growth

 

 

 

SS VLAB Release 3.30.13.00 changes:

Type

Category

Description

 Input_Change?

fix

tag

fix formatting problem for tag output in compreport.sso

 

fix

survey

fix issue with float Q that may have been source of some poor convergences

 

fix

controls

fix implementation of seasons as pseudo-years

 

fix

survey

fix ss_trans to handle mirror q better

 

fix

output

fix reporting when using seasonal MGparms

 

fix

output

fix output of Dirichlet parameters

 

fix

input

fix reading of sizecomp obs before styr

 

fix

output

fix output of equilibrium F when using F_reporting=5

 

fix

output

fix issues with F_reporting=4 and 5 for multiseason models

 

fix

statistics

Important: corrected the -logL calculation for parameter deviations.  This affects only a constant (N*log(s) is corrected to be N*log(1)) in the logL of the devs, it does not affect how the devs are used.

 

new

survey

create new Q option so that mirror Q uses an offset to ln(Q)

Yes, if used

new

growth

add growth cessation model per Maunder et al (2018)

Yes, if used

new

time_vary

add dev vector options so last dev can persist to end of time series (codes 21, 22, 23, 24)

Yes, if used

new

recruitment

add new recdev approach based on dev being relative to R0, not SRR; add column to spawn_recr in report.sso

Yes, if used

new

forecast

add new forecast feature:  year-specific ABC buffer

Yes, if used

new

output

add table for F in report.sso

 

revise

input

revise warning for recrdist options

 

revise

tag

move tag recapture code into a separate file for ease in future updating:  SS_tagrep.tpl

 

revise

output

re-order and organize the list of tables appearing at top of report.sso

 

revise

output

small fix to posteriors.sso for Q parms

 

 

 

Stephan Smith, modified 5 Years ago.

Fwd: [Stock Synthesis - Release updates] SS 3.30.13.00 released Mar 13, 201

Youngling Posts: 2 Join Date: 5/15/12 Recent Posts
Sharing this with VLST:  Great example of the collaboration.

---------- Forwarded message ---------
From: Richard Methot <VLab.Notifications@noaa.gov>
Date: Wed, Mar 13, 2019 at 6:51 PM
Subject: [Stock Synthesis - Release updates] SS 3.30.13.00 released Mar 13, 2019
To: VLab.Notifications@noaa.gov <VLab.Notifications@noaa.gov>


Latest SS News:  March 13, 2019

Greetings Stock Synthesis (SS) Users,

In this release of Stock Synthesis 3.30.13.00 we have included some fixes, augmentations and new features.  Thanks to some users who were pushing the envelope on time-varying parameters for finding some needed fixes to the code.  You can find all the changes listed in the Change Log and below.   The fixes were minor and mostly in the output, but there is one big fix affecting the calculation of the log likelihood for the parameter deviations.  We will put out a fuller report on that later and encourage you to contact us if you need info on this change sooner.

Please remember to check the change log frequently so you are aware of anything relevant for your SS project.  It also is a good idea to look over the *.ss_new files created by the newest version of SS so you can see listed all available input options.  Of course, the user manual, r4ss, and graphical Stock Synthesis Interface (SSI) have been updated.   The SSI now provides previews for the shape of most selectivity functions.

We have a new SS team member, Kathryn Doering.  She will be involved with VLAB maintenance, creation of an example library, and developing more R-based routines to enhance SS.  Welcome Kathryn.

If you have any SS usage questions or issues please post in the FORUM.

Problems using the SS VLAB site can be addressed to nmfs.stock.synthesis@noaa.gov

Please contact VLAB administrators at vlab.admin@noaa.gov to change your username or password. 

 

ADMB 12.0

SS is now being compiled exclusively with ADMB 12.0 (except the 32 bit versions).  An important new feature in ADMB 12.0 improves the MCMC algorithm.  See below from the admb-project.org website:

  • Dr. Cole Monnahan has developed a new MCMC functionality into ADMB. The current Metropolis algorithm was updated to have more meaningful console output and not require an sdreport object. The old method remains, and the new one is accessible with argument ‘-rwm’. A new algorithm, the no-U-turn sampler (a self-tuning extension of the existing ‘hybrid’ method was also added, and available with ‘-nuts’. The functionality matches the software platform Stan closely. The user can specify whether to use the MLE covariance, an arbitrary matrix, or adapt one during the warmup phase (recommended). Both these additions work much with the R package adnuts, which streamlines workflow, including parallel chains, diagnostic checking, and inference. See adnuts package vignette for more information.

 

Future big features on the wish list:

  • Option to replace input of growth parameters with input of empirical length-at-age
  • Change input of tagging parameters to use a much more concise design matrix approach
  • Adding otolith weight as a proxy for fish age as a new data type
  • Area-specific spawner-recruitment
  • Area-specific growth in addition to morph specific growth

 

 

 

SS VLAB Release 3.30.13.00 changes:

Type

Category

Description

 Input_Change?

fix

tag

fix formatting problem for tag output in compreport.sso

 

fix

survey

fix issue with float Q that may have been source of some poor convergences

 

fix

controls

fix implementation of seasons as pseudo-years

 

fix

survey

fix ss_trans to handle mirror q better

 

fix

output

fix reporting when using seasonal MGparms

 

fix

output

fix output of Dirichlet parameters

 

fix

input

fix reading of sizecomp obs before styr

 

fix

output

fix output of equilibrium F when using F_reporting=5

 

fix

output

fix issues with F_reporting=4 and 5 for multiseason models

 

fix

statistics

Important: corrected the -logL calculation for parameter deviations.  This affects only a constant (N*log(s) is corrected to be N*log(1)) in the logL of the devs, it does not affect how the devs are used.

 

new

survey

create new Q option so that mirror Q uses an offset to ln(Q)

Yes, if used

new

growth

add growth cessation model per Maunder et al (2018)

Yes, if used

new

time_vary

add dev vector options so last dev can persist to end of time series (codes 21, 22, 23, 24)

Yes, if used

new

recruitment

add new recdev approach based on dev being relative to R0, not SRR; add column to spawn_recr in report.sso

Yes, if used

new

forecast

add new forecast feature:  year-specific ABC buffer

Yes, if used

new

output

add table for F in report.sso

 

revise

input

revise warning for recrdist options

 

revise

tag

move tag recapture code into a separate file for ease in future updating:  SS_tagrep.tpl

 

revise

output

re-order and organize the list of tables appearing at top of report.sso

 

revise

output

small fix to posteriors.sso for Q parms

 

 

 


--
Richard Methot Stock Synthesis Virtual Lab Forum https://vlab.noaa.gov/group/stock-synthesis/discussions-forums-/-/message_boards/view_message/6161546VLab.Notifications@noaa.gov
XH
Xi He, modified 5 Years ago.

RE: SS 3.30.13.00 released Mar 13, 2019

Youngling Posts: 1 Join Date: 5/11/16 Recent Posts
Hi Rick, I am helping E.J. to set up cowcod SS model.  The model has no recr dev estimates.  So I set "do_recdev" to 0 (as shown below), and command out the rest settings for recr devs.  It seems to me that ss still needs some parameters to run.  I tried a few settings but still could not get it to work.  The SS manual does not have clear instructions what to do.  What did I do wrong in the following settings? Thank you again for help. Xi # 3 #_Spawner-Recruitment; Options: 2=Ricker; 3=std_B-H; 4=SCAA; 5=Hockey; 6=B-H_flattop; 7=survival_3Parm; 8=Shepherd_3Parm; 9=RickerPower_3parm 0  # 0/1 to use steepness in initial equ recruitment calculation 0  #  future feature:  0/1 to make realized sigmaR a function of SR curvature #_          LO            HI          INIT         PRIOR PR_SD       PR_type      PHASE    env-var    use_dev   dev_mnyr dev_mxyr     dev_PH      Block    Blk_Fxn #  parm_name              6            35           8.6 8.6            99             0         1          0 0          0          0          0          0          0 # SR_LN(R0)           0.21          0.99           0.6 0.6            99             0         -2          0 0          0          0          0          0          0 # SR_BH_steep              0             2             1 1            99             0         -4          0 0          0          0          0          0          0 # SR_sigmaR             -5             5             0 0            99             0         -4          0 0          0          0          0          0          0 # SR_regime              0           0.5             0 0            99             0         -3          0 0          0          0          0          0          0 # SR_autocorr 0 #do_recdev:  0=none; 1=devvector; 2=simple deviations # 1900# first year of main recr_devs; early devs can preceed this era # 2012 # last year of main recr_devs; forecast devs start in following year # -2 #_recdev phase # 1 # (0/1) to read 13 advanced options  # 0 #_recdev_early_start (0=none; neg value makes relative to recdev_start)  # -2 #_recdev_early_phase  # 1 #_forecast_recruitment phase (incl. late recr) (0 value resets to maxphase+1)  # 1 #_lambda for Fcast_recr_like occurring before endyr+1  # 1519 #_last_yr_nobias_adj_in_MPD; begin of ramp  # 1519 #_first_yr_fullbias_adj_in_MPD; begin of plateau  # 2014.1 #_last_yr_fullbias_adj_in_MPD  # 2018 #_end_yr_for_ramp_in_MPD (can be in forecast to shape ramp, but SS sets bias_adj to 0.0 for fcast yrs)  # 0.9999 #_max_bias_adj_in_MPD (-1 to override ramp and set biasadj=1.0 for all estimated recdevs)  # 0 #_period of cycles in recruitment (N parms read below)  # -5 #min rec_dev  # 5 #max rec_dev  # 0 #_read_recdevs #_end of advanced SR options