Fit an FLBRP object

This method carries out the necessary calculation to obtain the reference points for an object of class FLBRP. Results of the calculation are stored in @refpts slot of the object, and the object is then returned.

To directly obtain the recalculated refpts slot of an FLBRP object, the computeRefpts method can be used. This is equivalent to fitting the object using brp and then extracting the @refpts slot.

The plot method for the FLBRP class will show the relationship between equilibrium SSB, Yield and Profit against levels of fishing mortality and that of recruitment, Yield and profit against SSB.

# S4 method for FLBRP
brp(object)

# S4 method for FLBRP
computeRefpts(object)

# S4 method for FLBRP,missing
plot(
  x,
  refpts = dimnames(x@refpts)$refpt,
  obs = FALSE,
  labels = TRUE,
  shapes = "missing",
  colours = "missing",
  panels = NULL,
  ncol = 2,
  ...
)

Arguments

object: An object of class FLBRP
x: An object of class FLBRP
refpts: Reference points to include, defaults to all in standard refpts slot bu virgin. A vector of class character
obs: Should observations be plotted? Defaults to FALSE.

Value

An object of class FLBRP with reference points estimated

An object of class ggplot2.

Examples

data(ple4brp)

ple4brp <- brp(ple4brp)
refpts(ple4brp)
#> An object of class "FLPar"
#>         quant
#> refpt    harvest   yield     rec       ssb       biomass   revenue   cost     
#>   virgin  0.00e+00  0.00e+00  5.09e+05  1.74e+06  1.79e+06  0.00e+00  3.50e+07
#>   msy     2.75e-01  8.20e+04  1.13e+06  7.33e+05  8.32e+05  1.14e+08  4.33e+07
#>   crash   5.25e-01  7.38e-07  1.53e-05  2.95e-06  4.07e-06  9.80e-04  5.08e+07
#>   f0.1    1.55e-01  6.47e+04  8.92e+05  1.15e+06  1.23e+06  9.14e+07  3.97e+07
#>   fmax    2.08e-01  7.64e+04  1.02e+06  9.61e+05  1.06e+06  1.07e+08  4.12e+07
#>   spr.30  1.94e-01  7.38e+04  9.86e+05  1.01e+06  1.10e+06  1.04e+08  4.08e+07
#>   mey     2.61e-01  8.18e+04  1.12e+06  7.79e+05  8.78e+05  1.13e+08  4.28e+07
#>         quant
#> refpt    profit   
#>   virgin -3.50e+07
#>   msy     7.03e+07
#>   crash  -5.08e+07
#>   f0.1    5.18e+07
#>   fmax    6.58e+07
#>   spr.30  6.29e+07
#>   mey     7.07e+07
#> units:  NA 
m(ple4brp)[1:3,] <- 0.2
computeRefpts(ple4brp)
#> An object of class "FLPar"
#>         quant
#> refpt    harvest   yield     rec       ssb       biomass   revenue   cost     
#>   virgin  0.00e+00  0.00e+00  6.14e+05  1.56e+06  1.62e+06  0.00e+00  3.50e+07
#>   msy     2.42e-01  6.38e+04  1.15e+06  6.69e+05  7.62e+05  8.87e+07  4.23e+07
#>   crash   4.63e-01  1.06e-06  2.59e-05  4.99e-06  6.80e-06  1.42e-03  4.89e+07
#>   f0.1    1.56e-01  5.48e+04  1.01e+06  9.70e+05  1.06e+06  7.74e+07  3.97e+07
#>   fmax    2.09e-01  6.25e+04  1.11e+06  7.82e+05  8.74e+05  8.73e+07  4.13e+07
#>   spr.30  1.95e-01  6.10e+04  1.09e+06  8.32e+05  9.23e+05  8.56e+07  4.08e+07
#>   mey     2.29e-01  6.36e+04  1.14e+06  7.14e+05  8.07e+05  8.86e+07  4.19e+07
#>         quant
#> refpt    profit   
#>   virgin -3.50e+07
#>   msy     4.64e+07
#>   crash  -4.89e+07
#>   f0.1    3.77e+07
#>   fmax    4.61e+07
#>   spr.30  4.47e+07
#>   mey     4.67e+07
#> units:  NA 
data(ple4brp)
plot(ple4brp)

# ADD observations
plot(ple4brp, obs=TRUE)

# SELECT which refpts to plot
plot(ple4brp, refpts=c("msy", "crash", "f0.1"))

plot(ple4brp, refpts=c("msy", "crash", "f0.1"), colours=c("green", "red", "blue"))

plot(ple4brp, refpts=c("msy", "crash", "f0.1"),
  colours=c("green", "red", "blue"), shapes=c(21, 22, 23))

# method will only plot existing referenced points
refpts(ple4brp) <- refpts(ple4brp)[c("msy", "fmax", "spr.30"),]
plot(ple4brp)

# Select which panels to plot
plot(ple4brp, panels=1:4)

# and in which format
plot(ple4brp, panels=1:3, ncol=1)

Arguments

Value

See also

Examples