Data will eventually be in the icesdata
package
library(icesdata)
data(icesdata)
Figure 1 Time series of standardised catches from ICES Cat 1 stock assessments.
Figure 2 Standarised SSB from ICES Cat 1 stock assessments.
Figure 3 Time series of standardised catches and SSB from ICES Cat 1 stock assessments.
Figure 4 Stock recruitment relationhips.
Figure 5 Log likelihoods
Figure 6 Log likelihoods
Figure 10
Error in splinefun(initial, year/g) : zero non-NA points
Error in splinefun(initial, year/g) : zero non-NA points
Error in splinefun(initial, year/g) : zero non-NA points
Error in splinefun(initial, year/g) : zero non-NA points
Figure 11
Figure 12
Figure 13
Figure 14
Figure 15
Error in h(simpleError(msg, call)) :
error in evaluating the argument 'formula' in selecting a method for function 'model.frame': f.lower = f(lower) is NA
Error in h(simpleError(msg, call)) :
error in evaluating the argument 'formula' in selecting a method for function 'model.frame': f.lower = f(lower) is NA
Error in h(simpleError(msg, call)) :
error in evaluating the argument 'formula' in selecting a method for function 'model.frame': f.lower = f(lower) is NA
Error in h(simpleError(msg, call)) :
error in evaluating the argument 'formula' in selecting a method for function 'model.frame': f.lower = f(lower) is NA
Error in h(simpleError(msg, call)) :
error in evaluating the argument 'formula' in selecting a method for function 'model.frame': f.lower = f(lower) is NA
Error in h(simpleError(msg, call)) :
error in evaluating the argument 'formula' in selecting a method for function 'model.frame': f.lower = f(lower) is NA
Error in h(simpleError(msg, call)) :
error in evaluating the argument 'formula' in selecting a method for function 'model.frame': f.lower = f(lower) is NA
Error in h(simpleError(msg, call)) :
error in evaluating the argument 'formula' in selecting a method for function 'model.frame': f.lower = f(lower) is NA
Error in h(simpleError(msg, call)) :
error in evaluating the argument 'formula' in selecting a method for function 'model.frame': f.lower = f(lower) is NA
Figure 16
Figure 17
.id SRR msy emsy bmsy k.eb k.ssb
1 anf.27.3a46 Beverton & Holt 22434.239 152742.16 51115.10 326975.7 208133.9
2 anf.27.3a46 Ricker 27910.156 228489.04 85899.50 345069.2 216817.0
3 ank.27.78abd Beverton & Holt 13101.669 147559.09 48453.97 339069.0 164540.2
4 ank.27.78abd Ricker 18335.838 204051.14 66672.03 344143.0 164543.0
5 aru.27.5a14 Ricker 10468.198 52303.51 57380.23 134525.6 143642.4
6 aru.27.5a14 Beverton & Holt 8212.342 37167.19 41039.42 126867.0 148535.8
hmsy fmsy r r.eb m.spwn harvest.spwn wts
1 0.15885100 0.5778499 0.2888444 0.13908683 0 0 0.10731522
2 0.13028064 0.3929187 0.2742502 0.16452605 0 0 0.10731522
3 0.09298113 0.3152508 0.1934339 0.07981245 0 0 0.64602319
4 0.09415578 0.3216049 0.2355227 0.10574192 0 0 0.64602319
5 0.22332270 0.2014256 0.1548028 0.16693530 0 0 0.05857032
6 0.24968876 0.2232793 0.1388462 0.15766754 0 0 0.05857032
ks.D matsel
1 0.4516129 0.1077140
2 0.4516129 0.1077140
3 0.8125000 -0.3727192
4 0.8125000 -0.3727192
5 0.7727273 -0.4367806
6 0.7727273 -0.4367806