plot.medusa {geiger}

MEDUSA: modeling evolutionary diversification using stepwise AIC

Description

summarizing piecewise diversification models estimated by MEDUSA

Usage

	## S3 method for class 'medusa'
plot(x, cex = 0.5, time = TRUE, ...)

Arguments

x

an object of class medusa for plotting

cex

text size

time

logical. should a time axis be plotted. default = TRUE.

...

additional arguments to be passed to internal functions

Details

The medusa model returns a raw list object. This function is used to generate a modified edge matrix (see read.tree for details on the edge matrix), giving all relevant information about the estimated diversification process. The returned z-matrix includes: the ancestor (anc) and descendant (dec) relationships between nodes of the tree (using ape indices); the beginning (t.0) and ending (t.1) times and length (t.len) of each branch; the diversities at the start (n.0) and end (n.t) of each branch; the piecewise model assigned partition to the branch; whether the branch is associated with a shift; the timing of the shift (t.shift); the net-diversification rate (r) and relative-extinction rate (epsilon) associated with the branch as well as for the direct ancestor of the branch (ancestral.r and ancestral.epsilon). The z-matrix also includes a summary attribute that shows which model is chosen and associated information on model fit (see Examples).

The raw output of medusa contains an optimized model. The summary output may then be sent to a plotting function which will display the location on the tree where breakpoints have been placed. Note that the first piecewise model corresponds to the root and all descendants (until another breakpoint is encountered).

Author(s)

JW Brown <[e-mail]>, RG FitzJohn, ME Alfaro, LJ Harmon, and JM Eastman

References

Alfaro, ME, F Santini, C Brock, H Alamillo, A Dornburg, DL Rabosky, G Carnevale, and LJ Harmon. 2009. Nine exceptional radiations plus high turnover explain species diversity in jawed vertebrates. Proceedings of the National Academy of Sciences 106: 13410-13414.

See Also

medusa

Examples

    
    dat=get(data(whales))
    phy=dat$phy
    richness=dat$richness
        
    res <- medusa(phy, richness)
    
    # select best model based on AICc (showing the third model as best)
    plot(res, cex=0.5, label.offset=1) # using plot.medusa()
    title("AICc-chosen model")