plot_seasonality: Plot Seasonality
In mdkarcher/phylodyn: Statistical Tools for Phylodynamics

Plot Seasonality

plot_seasonality(BNPR_out, zero_date, start = 0, years = NULL,
  period = 1, ylim = NULL, nbreaks = 40, lty = 1, lwd = 2,
  col_years = "gray", col_mean = "black", main = "", axlabs = NULL,
  xlab = "Time", xmarline = 3, ylab = "Effective Population Size",
  log_y = TRUE, heatmaps = TRUE, heatmap_labels = TRUE,
  heatmap_labels_side = "right", heatmap_width = 7, legend = NULL,
  yscale = 1, ...)

`BNPR_out`	output of BNPR or BNPR_PS.
`zero_date`	numeric which time to present as 0 on graph.
`start`	numeric when to start displaying data.
`years`	numeric how many years to count back in time.
`period`	numeric how long of a period is a "year".
`ylim`	numeric y-axis interval.
`nbreaks`	integer number of bins for sampling heatmap.
`lty`	numeric line type for the estimated trajectory paths.
`lwd`	numeric line width for the estimated trajectory paths.
`col_years`	color of year paths.
`col_mean`	color of mean path.
`main`	character main plot title.
`axlabs`	character vector x-axis labels.
`xlab`	character x-axis label.
`xmarline`	numeric if not using default x-axis labels, how far to put the labels from the axis.
`ylab`	character y-axis label.
`log_y`	logical whether to log-scale the y-axis.
`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.
`legend`	character legend text. The defauly, NULL, disables.
`yscale`	numeric scaling applied to all effective population calculations.
`...`	additional arguments to be passed onto plot().