Learn R Programming

phylodyn (version 0.9.0)

plot_BNPR: Plot a BNPR output

Description

Plot a BNPR output

Usage

plot_BNPR(BNPR_out, traj = NULL, xlim = NULL, ylim = NULL, nbreaks = 40,
  lty = 1, lwd = 2, col = "black", main = "", log = "y",
  ylab = "Effective Population Size", xlab = "Time", xmarline = 3,
  axlabs = NULL, traj_lty = 2, traj_lwd = 2, traj_col = col,
  newplot = TRUE, credible_region = TRUE, heatmaps = TRUE,
  heatmap_labels = TRUE, heatmap_labels_side = "right", heatmap_width = 7,
  yscale = 1, ...)

Arguments

BNPR_out
output of BNPR or BNPR_PS.
traj
function summarizing the true effective population size trajectory.
xlim
numeric x-axis interval.
ylim
numeric y-axis interval.
nbreaks
integer number of bins for sampling heatmap.
lty
line type for estimated trajectory.
lwd
line width for estimated trajectory.
col
color for estimated trajectory.
main
character main plot title.
log
character which axes to plot log-scale. Defaults to "y".
ylab
character y-axis label.
xlab
character x-axis label.
xmarline
numeric if not using default x-axis labels, how far to put the labels from the axis.
axlabs
character vector x-axis labels.
traj_lty, traj_lwd, traj_col
line type, line width, and line color for the true trajectory.
newplot
boolean whether to create a new plot or superimpose over a previously open plot.
credible_region
logical whether to display pointwise credible region.
heatmaps
boolean whether to display sampling and coalescent heatmaps.
heatmap_labels
boolean whether to display labels on heatmaps.
heatmap_labels_side
string which side of plot to display heatmaps.
heatmap_width
numeric how wide heatmaps should be.
yscale
numeric scaling applied to all effective population calculations.
...
additional arguments to be passed onto plot().