plot_BEAST: Plot trajectory estimate and credible intervals
In mdkarcher/phylodyn: Statistical Tools for Phylodynamics

Plot trajectory estimate and credible intervals

plot_BEAST(midpts, med, up, low, cutoff = NULL, ylim = NULL,
  xlab = "Time", ylab = "Effective Population Size", axlabs = NULL,
  samp_times = NULL, n_sampled = NULL, nbreaks = 40,
  med_col = grDevices::rgb(0.33, 0.484, 0.828),
  cred_col = grDevices::rgb(0.33, 0.484, 0.828, 0.4),
  heatmap_width = 7, heatmap_legend_cex = 1, ...)

`midpts`	numeric vector of x-values at which to plot estimate and credible interval bounds.
`med`	numeric vector of estimated trajectory values.
`up`	numeric vector of trajectory pointwise credible interval upper bounds.
`low`	numeric vector of trajectory pointwise credible interval lower bounds.
`cutoff`	numeric maximum time (in past) plotted.
`ylim`	numeric y-axis interval.
`xlab`	character x-axis label.
`ylab`	character y-axis label.
`axlabs`	character vector x-axis labels.
`samp_times`	numeric vector of sampling times.
`n_sampled`	integer vector of number of samples taken at each element of samp_times.
`nbreaks`	integer number of bins in the sampling time histogram heatmap.
`med_col`	color of trajectory estimate line.
`cred_col`	color of trajectory pointwise credible interval region.
`heatmap_width`	numeric vertical width of sampling time histogram heatmap.
`heatmap_legend_cex`	numeric plotting expansion for heatmap legend.
`...`	additional parameters to be passed to plot().