Forums

Back

CPUE plot

YL
Yevhen Leonchyk, modified 3 Years ago.

CPUE plot

Youngling Posts: 36 Join Date: 8/1/17 Recent Posts

Dear SS3 team,

Please could you tell me how SS3 normalizes different abundance indices for building that stand_cpue_all plot? Is it explained in the tech manual?

Thank you in advance.

YL
Yevhen Leonchyk, modified 3 Years ago.

RE: CPUE plot

Youngling Posts: 36 Join Date: 8/1/17 Recent Posts

It seems to me that ss3 uses mean value, doesn't it?

IS
Ian Stewart, modified 3 Years ago.

RE: CPUE plot

Youngling Posts: 6 Join Date: 3/2/16 Recent Posts

Hello Yevhen,

Yes, this plot just uses the simple mean of the series to standardize each index and is created in the r4ss code rather than in SS directly.

YL
Yevhen Leonchyk, modified 3 Years ago.

RE: CPUE plot

Youngling Posts: 36 Join Date: 8/1/17 Recent Posts

Thank you for the reply.

How about such method

(Index-mean(Index))/sd(Index),

sd – standard deviation.

Might it be better for a wide variation case? But it seems to be "bad" when SD is close to zero.

KD
Kathryn Doering, modified 3 Years ago.

RE: CPUE plot

Youngling Posts: 57 Join Date: 2/19/19 Recent Posts

Hi Yevhen,

 

To follow up on Ian's response, the place in the r4ss development branch code this standardization occurs: https://github.com/r4ss/r4ss/blob/development/R/SSplotIndices.R#L564

And most of the code for making this plot is also in the same function: https://github.com/r4ss/r4ss/blob/development/R/SSplotIndices.R#L779-L811

 

Thanks,

Kathryn

YL
Yevhen Leonchyk, modified 3 Years ago.

RE: CPUE plot

Youngling Posts: 36 Join Date: 8/1/17 Recent Posts
Dear Kathryn,

thank you for the links!

U
Anonymous, modified 3 Years ago.

RE: CPUE plot

Hi Yevhen,

The initial intention of for this figure was to be able to identify consistent trends across indices on a single figure, hence the scaling around 1.0 .  This allows indices that have very different scales in relative abundance to be easily comparable.  The alternative calculation you suggested could be a useful additional figure for r4ss to include that could be complimentary to the current figure.  

KD
Kathryn Doering, modified 3 Years ago.

RE: CPUE plot

Youngling Posts: 57 Join Date: 2/19/19 Recent Posts

Hi Yevhen,

 

Sorry, I missed your response about z-scoring indices instead of dividing by the mean to standardize. I can see why that would be useful. Is there any particular reason that plot would be more useful than the current one (or in addition to the current one)? r4ss welcomes contributions and it is always possible to add/change plots depending on the needs of the r4ss users.

 

Thanks,

Kathryn

YL
Yevhen Leonchyk, modified 3 Years ago.

RE: CPUE plot

Youngling Posts: 36 Join Date: 8/1/17 Recent Posts

I use z-scores for scaling in Cluster analysis. But the mean-centered method is "one-size-fits-all".
​​​​​​​Thank you for the explanation.

IT
Ian Taylor, modified 3 Years ago.

RE: CPUE plot

Youngling Posts: 117 Join Date: 12/8/14 Recent Posts
Thanks Yevhen Leonchyk for pointing out the lack of information on the standardization approach for that figure. And thanks Ian S. and Kathryn for the explanations.

I just expanded the caption slightly so it now reads "Standardized indices overlaid. Each index is rescaled to have mean observation = 1.0."
The revised caption is in the development branch of r4ss available at 
remotes::install_github('r4ss/r4ss', ref = 'development')
and will be merged into the main branch eventually.

I should also note that Allan Hicks made good progress a few years ago on a revised version of this figure which adds uncertainty intervals and (optionally) model expectations as shown at https://github.com/r4ss/r4ss/issues/67. This conversation can serve as inspiration to get back to that effort as time allows. Feel free to make further suggestions on that github issue or here on VLab.
-Ian

On Fri, May 29, 2020 at 10:18 AM Yevhen Leonchyk <VLab.Notifications@noaa.gov> wrote:

I use z-scores for scaling in Cluster analysis. But the mean-centered method is "one-size-fits-all".
​​​​​​​Thank you for the explanation.


--
Yevhen Leonchyk Stock Synthesis Virtual Lab Forum https://vlab.noaa.gov/web/stock-synthesis/public-forums/-/message_boards/view_message/10298233 VLab.Notifications@noaa.gov