MLEmid.MLEbin.hist: One figure with eight histograms for each of MLEmid and...

View source: R/plotting.R

MLEmid.MLEbin.histR Documentation

One figure with eight histograms for each of MLEmid and MLEbin methods and four binning types

Description

The default plot here reproduces Figure 4 of MEPS paper, showing the estimated exponent $b$ for 10,000 simulated data sets, binned using four methods and fitted using the MLEmid and MLEbin methods.

Usage

MLEmid.MLEbin.hist(
  results.list,
  vertCol = "red",
  vertThick = 1,
  xrange = c(-2.4, -1.1),
  xbigticks.by = 0.4,
  xsmallticks.by = 0.1,
  ylimA = c(0, 7000),
  yBigTickLab.by = 3000,
  yBigTickNoLab.by = 1000,
  ySmallTick.by = 500,
  binwidth = NA,
  omi = c(0.12, 0.05, 0.2, 0),
  mfrow = c(4, 2),
  mai = c(0.5, 0.5, 0, 0.3),
  xaxs = "i",
  yaxs = "i",
  mgp = c(2, 0.5, 0),
  cex = 0.8,
  cexAxis = 0.9,
  inset = c(0, -0.04),
  legLabMid = c("(a)", "(c)", "(e)", "(g)"),
  legLabBin = c("(b)", "(d)", "(f)", "(h)")
)

Arguments

results.list

output list from MLEbin.simulate(); see ?MLEbin.simulate() for details

vertCol

colour for vertical line at true value of b

vertThick

thickness of vertical line at true value of b

xrange

range of x values (should change to xlimA for consistency)

xbigticks.by

increment between big tick marks on x-axis (all labelled)

xsmallticks.by

increment between small tick marks on x-axis

ylimA

range of y values

yBigTickLab.by

increment between labelled big tick marks on y-axis

yBigTickNoLab.by

increment between all big tick marks on y-axis

ySmallTick.by

increment between small tick marks on y-axis

binwidth

bin width for histograms; if NA then gets calculated such that b.known is a midpoint of a bin

omi, mfrow, mai, xaxs, yaxs, mgp, cex, inset

standard options for par(), defaults are for Figure 4

cexAxis

font size for axis labels

legLabMid

character vector of length 4 for labelling each panel in the MLEmid column

legLabBin

character vector of length 4 for labelling each panel in the MLEBin column

Value

figure with eight histograms, one for each combination of binning type and fitting method.

Author(s)

Andrew Edwards


andrew-edwards/sizeSpectra documentation built on June 28, 2023, 7:09 p.m.