plotPrior: Plot the prior distribution used in the MCMC analysis

In Caetanods/ratematrix: Bayesian Estimation of the Evolutionary Rate Matrix

Plot the prior distribution used in the MCMC analysis

Description

Function plots the prior distribution used in the MCMC analysis.

Usage

plotPrior(handle, n = 1000, root = FALSE, color = "black", ...)

Arguments

handle	the output object from the 'ratematrixMCMC' function.
n	number of samples from the prior to be plotted (default is 1000).
root	whether to plot the prior for the root value instead of the evolutionary rate matrix (default is FALSE).
color	color for the plot (default is "black").
...	other parameters to be passed to the function 'plotRatematrix' or 'plotRootValue'. See help page for list of possible parameters.

Details

Function will make a plot of the prior for the evolutionary rate matrix by default. One can plot the prior for the root value instead by setting 'root' to TRUE.

The prior distribution often has a different range of parameter values when compared to the posterior distribution. Depending on the prior configuration the range of the prior can be orders of magnitude larger than the posterior distribution. In this case, it is important to observe the scale of the x axis when comparing the prior and the posterior distribution. One can use the 'set.xlim' parameter to restrict the x axis for plotting the prior to be similar to the posterior distribution. However, often the region of parameter space of the posterior distribution has a low likelihood under the prior. This results in problems to take samples from that region to make the plot. This problem can be identified when the 'set.xlim' argument is changed and the plot shows only a few samples.

Value

A plot similar to 'plotRatematrix'.

Author(s)

Daniel S. Caetano and Luke J. Harmon

Caetanods/ratematrix documentation built on Jan. 4, 2023, 3:12 p.m.