makeplot.asdsf: Plot the Standard Deviation of Split Frequencies over the...

In rwty: R We There Yet? Visualizing MCMC Convergence in Phylogenetics

Description

This function takes two or more rwty.chain ojects and returns a plot of ASDSF as the run progresses. The solid line with points shows the Average Standard Deviation of Split Frequences at the current generation The grey ribbon shows the upper and lower 95

Usage

makeplot.asdsf(chains, burnin = 0, window.size = 20, min.freq = 0)

Arguments

chains	A list of rwty.chain objects.
burnin	The number of trees to eliminate as burnin. Defaults to zero.
window.size	The number of trees between each point at which the ASDSFs is calculated (note, specified as a number of sampled trees, not a number of generations)
min.freq	The minimum frequency for a node to be used for calculating ASDSF.

Value

output A cumulative plot of ASDSF across all chains

Examples

## Not run: 
data(fungus)
p <- makeplot.asdsf(fungus, burnin = 20)
p

## End(Not run)

Example output

Loading required package: ape
Loading required package: ggplot2
[1] "Creating ASDSF plot"
$asdsf.plot

Warning messages:
1: Transformation introduced infinite values in continuous y-axis 
2: Transformation introduced infinite values in continuous y-axis 
3: Transformation introduced infinite values in continuous y-axis 
4: Transformation introduced infinite values in continuous y-axis 
5: Transformation introduced infinite values in continuous y-axis 
6: Transformation introduced infinite values in continuous y-axis 

rwty documentation built on May 2, 2019, 4 p.m.