Create a labelled genealogy from an epidemiological record

getLabGenealogy(epirecord, epsilon = 0.001234)

Arguments

epirecord

an epidemiological record, as output from the function makeEpiRecord. It must be a matrix of at least two rows and with five columns named "Infectee", "Infector", "InfnTime", "RecTime", "DoneFlag".

epsilon

an optional small number to be used for branch lengths which would otherwise be zero.

Value

An object of class phylo representing the transmission tree from infectors to infectees.

See also

Author

Caroline Colijn c.colijn@imperial.ac.uk

Michelle Kendall michelle.louise.kendall@gmail.com

Examples

## Generate an epidemiological record:
myepirecord <- makeEpiRecord(c(1,2,3,4))
## make the corresponding genealogy from this record:
mygenealogy <- getLabGenealogy(myepirecord)
plot(mygenealogy)