c.phybreak: Combine mcmc-chains of phybreak-objects.

In donkeyshot/phybreak: Analysis of Outbreaks with Sequence Data

Description

The function combines the sample slots of its arguments. The variable slot (current state) is taken from the first argument

Usage

## S3 method for class 'phybreak'
c(...)

Arguments

...	Objects of class phybreak to be concatenated.

Value

The first phybreak-object with mcmc samples of all objects.

Author(s)

Don Klinkenberg [e-mail]

References

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

Examples

#First create a phybreak-object
simulation <- sim_phybreak(obsize = 5)
MCMCstate1 <- phybreak(dataset = simulation)
MCMCstate1 <- burnin_phybreak(MCMCstate1, ncycles = 20)
MCMCstate1 <- sample_phybreak(MCMCstate1, nsample = 50, thin = 2)

MCMCstate2 <- phybreak(dataset = simulation)
MCMCstate2 <- burnin_phybreak(MCMCstate2, ncycles = 20)
MCMCstate2 <- sample_phybreak(MCMCstate2, nsample = 50, thin = 2)

MCMCstate <- c(MCMCstate1, MCMCstate2)

donkeyshot/phybreak documentation built on Sept. 17, 2021, 9:32 p.m.