commfish / seak_sablefish

NSEI sablefish stock assessment
8 stars 5 forks source link

Bootstrap confidence intervals for CPUE #23

Closed jysullivan closed 4 years ago

jysullivan commented 5 years ago

@ben-williams recommended using bootstrap CIs instead of standard deviations to characterize variability in survey and fishery cpue.

e.g. stat_summary(fun.data = mean_cl_boot, geom = ‘smooth’)

jysullivan commented 5 years ago

@ben-williams I used Ben Bolker's suggestion: https://stackoverflow.com/questions/38554383/bootstrapped-confidence-intervals-with-dplyr

The code still has the alternative figs using +/- 1 sd b/c I'm not sold on using a confidence interval. I'm really interested in characterizing the variability in the data, not so much in displaying how well the mean is estimated.

A comparison of the two methods for Fishery CPUE:

Top: +/- 1 sd fshcpue_1997_2018

Bottom: 95% bootstrap CI fshcpue_bootCI_1997_2018

ben-williams commented 5 years ago

If you are trying to explore the variability in the data, that is fine, but be sure to write this up in the document as there is typically more focus on understanding the degree of certainty of an estimate. Further, I would encourage you to think about what will be the better implementation for an asa (likely inverse var?)