h_process_recon: Prepare HiSSE marginal reconstruction for plotting

Description Usage Arguments Value Examples

View source: R/utilhisse.R

Description

A function to process the output from hisse::MarginRecon. First we get model averaged rates for tips and nodes of the phylogeny. Then convert the tree from phylo format into a treedata class and assign the data.frame of mode-averaged rates to the @data slot. We output the node and tip rates tables and the reformated tree object for downstream plotting with ggtree and ggplot.

Usage

1
h_process_recon(hisse_recon)

Arguments

hisse_recon

An object produced by hisse::MarginRecon, or a list of such objects over which we can model-average using hisse::GetModelAveRates

Value

A list with three components:

Examples

1
2
3
4
5
data("diatoms")
processed_hisse <- h_process_recon(hisse_recon=diatoms$cid4_recon)
processed_hisse$tip_rates
processed_hisse$node_rates
processed_hisse$tree_data

teofiln/utilhisse documentation built on Oct. 4, 2019, 2:44 p.m.