plot: plot

In laurieKell/biodyn: Biomass Dynamic Model

Description

Creates a ggplot2 object that plots time series of biomass, harvest rate and catch. The basic object can then be modified by adding ggplot2 layers.

Usage

## S4 method for signature 'biodyn,missing'
plot(x, y, probs = c(0.95, 0.75, 0.5, 0.25, 0.05),
  na.rm = FALSE, type = 7, worm = NULL, fn = list(Stock = function(x)
  stock(x), Harvest = function(x) harvest(x), Yield = function(x) catch(x)),
  facet = facet_wrap(~qname, scales = "free", ncol = 1), ...)

Arguments

x	an object of class biodyn
y	second argument
probs	numeric vector of probabilities with values in [0,1]. (Values up to 2e-14 outside that range are accepted and moved to the nearby endpoint.)
na.rm	FALSE
type	an integer between 1 and 9 selecting one of the nine quantile algorithms detailed below to be used.
worm	iters
fn	functions
facet	facet for panels

Value

an ggplot2 object

See Also

plotPrd, plotEql

Examples

## Not run: 

x =sim()
x =window(x,end=49)
bd1=fwd(x,harvest=rlnorm(200,log(harvest(x)[,-1]),.2))
bd2=fwd(x,harvest=rlnorm(200,log(harvest(x)[,-1])*1.5,.2))
plot(biodyns("1"=bd1,"2"=bd2))

## End(Not run)

laurieKell/biodyn documentation built on May 20, 2019, 7:58 p.m.