plot.medusa | R Documentation |

summarizing piecewise diversification models estimated by MEDUSA

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

`x` |
an object of class |

`cex` |
text size |

`time` |
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 <phylo.jwb@gmail.com>, 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.

`medusa`

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

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

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.