Download BET 2023 assessment report:
- Stock assessment of bigeye tuna in the western and central Pacific Ocean: 2023
WCPFC-SC19-2023/SA-WP-05
Download BET 2023 diagnostic model:
- Clone the bet-2023-diagnostic repository or download as main.zip file
Download BET 2023 grid results:
- The bet-2023-grid repository includes a bet-2023-grid-results.zip file
The BET 2023 assessment used a structural uncertainty grid with 54 models:
Axis | Levels | Option |
---|---|---|
Tag mixing | 2 | 1, 2* quarters |
Size data weighting | 3 | Sample sizes divided by 10, 20*, 40 |
Age data weighting | 3 | 0.5, 0.75*, 1 |
Steepness | 3 | 0.65, 0.80*, 0.95 |
The bet-2023-grid-results.zip file contains all files necessary to run or browse the BET 2023 grid models.
The grid models are run from a par file, as described in the corresponding doitall.sh
script. This starting par file is the best of 20 jittered par files from the pre-grid analysis.
The final par and rep files are consistently named final.par
and plot-final.par.rep
to facilitate harvesting results from across the 54 grid member models.
Preview of zip file contents:
bet-2023-grid-results.zip
├── bin
│ └── mfclo64
└── grid
├── m1_s10_a050_h65
│ ├── 14.par
│ ├── 15.par
│ ├── bet.age_length
│ ├── bet.frq
│ ├── bet_hess_inv_diag
│ ├── bet_pos_hess_cor
│ ├── bet.tag
│ ├── bet.var
│ ├── doitall.sh
│ ├── final.par
│ ├── mfcl.cfg
│ ├── neigenvalues
│ ├── plot-final.par.rep
│ ├── test_plot_output
│ └── xinit.rpt
├── m1_s10_a050_h80
│ ├── ...
The estimation_uncertainty.R script uses Monte Carlo simulations to add estimation uncertainty to the structural uncertainty grid estimates of reference points. The resulting means and quantiles are found in estimation_uncertainty.csv.
See also Section 6.2.3 and Table 5 in the BET 2023 stock assessment report.
The script requires the FLCore and FLR4MFCL packages:
install_github("flr/FLCore")
install_github("PacificCommunity/FLR4MFCL")