The dataset consists of sexual size-dimorphism data for 38 species of anoles from Cuba, Hispaniola, Jamaica, and Puerto Rico (Butler, Schoener, and Losos 2000). Each of these species belongs to one of six microhabitat types, or "ecomorphs" (sensu Williams, 1972): trunk-ground, grass-bush, trunk, trunk-crown, twig, and crown-giant. The data were used to demonstrate an evolutionary association between habitat type and degree of sexual size dimorphism.
A data frame with 38 observations on the following 6 variables.
node: Labels for the nodes.
species: Names of extant species.
log.SSD: Log sexual size dimorphism of extant species.
ancestor: Name of ancestor node.
time: Time of node.
OU.1: a factor with one level,
OU.7: a factor with levels corresponding to ecomorph (
Size dimorphism was calcuated as the log-ratio of male snout-to-vent length to female snout-to-vent length (males are larger).
In this example, we tested three models of evolution: Brownian motion, Ornstein-Uhlenbeck with one global optimum, and Ornstein-Uhlenbeck with seven optima (one for each ecomorph type plus an additional one for an "unknown" type).
For the seven-optima model, we assigned each terminal branch to an optimum according to the ecomorph type of the extant species. Because we had no information to help guide hypotheses about internal branches, we assigned internal branches to the "unknown" selective regime. The phylogeny of these species is consistent with and adaptive radiation, with a burst of speciation events early in the evolutionary history of this clade (see phylogeny in Butler & King (2004) or example below).
Marguerite A. Butler, Aaron A. King
## Analysis of sexual size dimorphism data ## Save time for CRAN tree <- with(anolis.ssd,ouchtree(node,ancestor,time/max(time),species)) plot(tree,node.names=TRUE) h1 <- brown(anolis.ssd['log.SSD'],tree) h1 plot(h1) h2 <- hansen(anolis.ssd['log.SSD'],tree,anolis.ssd['OU.1'],sqrt.alpha=1,sigma=1) h2 plot(h2) h3 <- hansen(anolis.ssd['log.SSD'],tree,anolis.ssd['OU.7'],sqrt.alpha=1,sigma=1) h3 plot(h3)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.