Description Usage Arguments Details Value Author(s) References See Also Examples

View source: R/plot.timeSlice.ML.R

Summarises phenotypic rate variation on phylogenies through

1 2 3 4 5 6 7 8 9 10 11 12 13 |

`x` |
Output of a timeSlice analysis in |

`...` |
Other functions to pass to |

`cutoff` |
Value for differences in AIC scores when comparing models. More complex models with an AIC score more than this number of units lower than simpler models are retained (as per runMedusa in the geiger package). |

`AICc` |
If |

`lowerBound` |
Minimum value for parameter estimates. |

`upperBound` |
Maximum value for parameter estimates. |

`phylo.plot` |
Logical. If |

`colour.ramp` |
The colours signifying different rates from low (first colour) to high (second colour) |

`cex.plot` |
Character expansion for the plot of rates through time |

`model.average` |
Logical only applicable to "timeSlice" models. Will return the model averaged timeSlice for models in which multiple shifts were considered (i.e, when splitTime is NULL). If TRUE, the function returns a plot showing the relative weights of each shift time and the model-averaged rates through time that are weighted by their relative weights. If TRUE, plot.phylo is ignored. |

This functions summarises the output of a "timeSlice" model in `transformPhylo.ML`

(see below). The best overall model is chosen based on AIC (or AICc if AICc=TRUE). The cut-off point for improvement in AIC score between successively more complex models can be defined using cutoff. The default cutoff is 4 but this is somewhat arbitrary and a "good" cut-off may well vary between data sets so it may well be worth exploring different cutoffs.

ModelFit Summary of the best optimal rate shift model or the model average of each split time (if model averaging was used).

Rates Summary of the rate parameters from the best rate shift model or the model averaged rates through time.

optimalTree A phylo object with branch lengths scaled relative to rate and a plot of estimated rates through time with their associated CIs.

Mark Puttick

To Add

1 2 3 4 5 6 7 8 9 10 11 12 13 14 | ```
data(anolis.tree)
data(anolis.data)
attach(anolis.data)
male.length <- matrix(Male_SVL, dimnames=list(rownames(anolis.data)))
sortedData <- sortTraitData(anolis.tree, male.length)
phy <- sortedData$phy
male.length <- sortedData$trait
phy.clade <- extract.clade(phy, 182)
male.length.clade <- as.matrix(male.length[match(phy.clade$tip.label,
rownames(male.length)),])
timeSlice.10.ml <- transformPhylo.ML(y=male.length.clade, phy=phy.clade, model="timeSlice",
splitTime=c(10))
outputSummary <- plot(timeSlice.10.ml, cutoff=0.001, cex=0.5,
colour.ramp=c("blue", "red"))
``` |

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.