plot.medusa: MEDUSA: modeling evolutionary diversification using stepwise...

plot.medusaR Documentation

MEDUSA: modeling evolutionary diversification using stepwise AIC


summarizing piecewise diversification models estimated by MEDUSA


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



an object of class medusa for plotting


text size


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


additional arguments to be passed to internal functions


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).


JW Brown <>, RG FitzJohn, ME Alfaro, LJ Harmon, and JM Eastman


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



    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")

geiger documentation built on April 4, 2023, 1:07 a.m.