bias.plot.BR: Bomb Radiocarbon Plots

In ballengerj/FishyR: Functions to Facilitate Standard Analyses Conducted in Fisheries Science

Bomb Radiocarbon Plots

Description

bias.plot.BR plots output from function calc.bias

Usage

bias.plot.BR(bias.out, plottype)

Arguments

bias.out	list containing output from function calc.bias
plottype	integer, with possible values of 1 and 2 that controls what, type of plots is plotted: 1 the simulated h distribution is plotted, along with arrows showing the h values for different values of percentage bias and broken lines showing the 95% confidence interval for h 2 a cumulative distribution of percentage bias is plotted, with broken lines showing the 95% confidence interval for bias

Value

A plot.

See Also

Other Bomb Radiocarbon Analyses: calc.bias(), calc.h(), closest.pt(), get.h.from.bias(), get.nearest(), get.ref.line(), sim.hdist()

Examples

data(snapper)
snapper.red <- snapper[snapper$C14year >= 1955 & snapper$C14year <= 1972, ]
data(bluenose)
bluenose.red <- bluenose[bluenose$C14year >= 1955 &
  bluenose$C14year <= 1972, ]
bluenose.red <- bluenose.red[bluenose.red$C14 >= -506 &
  bluenose.red$C14 <= 94.0, ]
out <- calc.bias(snapper.red, bluenose.red, nsim = 5000, plot.type = 0,
  x = "C14year", y = "C14", age = "age", cy = "catch.year",
  age.err = "age.err", s.age = "samp.age", y.err = "C14.err",
  bias.range = c(-50, 50), bias.step = 1)
bias.plot.BR(out, plottype = 1)
bias.plot.BR(out, plottype = 2)

ballengerj/FishyR documentation built on June 17, 2022, 10:33 p.m.