phybreak (version 0.2.0)

plot.phybreakdata: Plotting a phybreakdata object.

Description

Plots a phybreakdata-object twice: (1) as transmission tree and (2) as phylogenetic tree, using the default graphical parameters of plotTrans and plotPhylo. The default is to plot the current state, but any posterior sample can be chosen, as well as various consensus trees. Consensus tree "edmonds" plots only a transmission tree, consensus tree "mcc" only a phylogenetic tree.

Usage

# S3 method for phybreakdata
plot(x, ...)

Arguments

x

An object of class phybreakdata.

...

Some methods for this generic require additional arguments. None are used in this method.

References

Klinkenberg et al. (2017) Simultaneous inference of phylogenetic and transmission trees in infectious disease outbreaks. PLoS Comput Biol, 13(5): e1005495.

Examples

# NOT RUN {
#First build a phybreak-object containing samples.
simulation <- sim.phybreak(obsize = 5)
MCMCstate <- phybreak(data = simulation$sequences, times = simulation$sample.times)
MCMCstate <- burnin.phybreak(MCMCstate, ncycles = 20)
MCMCstate <- sample.phybreak(MCMCstate, nsample = 50, thin = 2)

plot(MCMCstate, plot.which = "mpc")
# }