Description Usage Arguments Value Author(s) References Examples
Estiamtes function data for tip species curves from tip coefficients.
1 | get.aligned.function.data(tip.coefficients, ylength, ymin = 0.01, ymax = 0.99)
|
tip.coefficients |
Matrix of estimated regression coefficients of tip curves. Row names should correspond to species names. The first column should contain the logit glm intercept; the second column contains the logit glm slope. |
ylength |
How many landmarks (points on the curve) to evaluate. |
ymin |
Because 0 and 1 are undefined for inverse logit functions, the minimum and maximum values are defaulted to .01 and .99. Can be adjusted as needed. |
ymax |
Because 0 and 1 are undefined for inverse logit functions, the minimum and maximum values are defaulted to .01 and .99. Can be adjusted as needed. |
A data frame of aligned X-coordinates for function-valued traits for a given constant Y.
Eric W. Goolsby
Goolsby, E.W. 2015. "Phylogenetic comparative methods for evaluating the evolutionary history of function-valued traits." Systematic Biology. In press.
1 2 3 4 5 6 7 8 9 | # simulate evolution of a function-valued trait (glm with logit link)
sim_data <- sim.curves()
# get tip coefficients and aligned function data
tip.coefficients <- get.tip.coefficients(formula = y~x,tree = sim_data$tree,data = sim_data$data)
data <- get.aligned.function.data(tip.coefficients)
# estimate evolutionary rates
evo.model.fitted <- evo.model(sim_data$tree,data)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.