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.

See also

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)