tseries: Extract Time Series Statistics from FLStock Objects

In laurieKell/FLCandy: Candidate methods for FLR

Extract Time Series Statistics from FLStock Objects

Description

Creates a data frame of time series statistics from FLStock objects, including catch, ebiomass, SSB, fishing mortality, harvest rate, and mean natural mortality.

Calculates the surplus production and expected yield etc for the estinates of SSB and biomass

Usage

tseries(object, ...)

## S4 method for signature 'FLBRP'
tseries(object)

## S4 method for signature 'FLStock'
tseries(
  object,
  flqs = list(catch = function(object) catch(object), eb = function(object)
    ebiomass(object), ssb = function(object) ssb(object), f = function(object)
    fbar(object), h = function(object) catch(object)/ebiomass(object), m =
    function(object)
    FLQuant(aaply(m(object)[ac(range(object)["minfbar"]:range(object)["maxfbar"])], 2,
    mean), dimnames = dimnames(fbar(object))))
)

## S4 method for signature 'FLStocks'
tseries(
  object,
  flqs = list(catch = function(object) catch(object), eb = function(object)
    ebiomass(object), ssb = function(object) ssb(object), f = function(object)
    fbar(object), h = function(object) catch(object)/ebiomass(object), m =
    function(object)
    FLQuant(aaply(m(object)[ac(range(object)["minfbar"]:range(object)["maxfbar"])], 2,
    mean), dimnames = dimnames(fbar(object))))
)

Arguments

object	an FLBRP object
...	Additional arguments (not currently used)
seasons	a numeric with seasons

Value

A data frame containing time series of:

• catch - Catch values

• eb - Exploitable biomass

• ssb - Spawning stock biomass

• f - Fishing mortality (Fbar)

• h - Harvest rate (catch/ebiomass)

• m - Mean natural mortality

FLQuants object

See Also

expand

Examples

## Not run: 
# For single stock
data(ple4)
ts1 <- tseries(ple4)

# For multiple stocks
stocks <- FLStocks(stock1=ple4, stock2=ple4)
ts2 <- tseries(stocks)

## End(Not run)

## Not run: 

laurieKell/FLCandy documentation built on April 17, 2025, 5:23 p.m.