X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=man%2Fcompar.ou.Rd;h=84be45483a7ad914fbaf47cac673ca7e1d4f4d0b;hb=d1546ec66ff1a8ea123adefebe14f6316c23705f;hp=b6f4b5ceea48114ff38b7c7b1bf6a57abcce9268;hpb=c827059eeafc8cbe41c812b26979543ab287803e;p=ape.git diff --git a/man/compar.ou.Rd b/man/compar.ou.Rd index b6f4b5c..84be454 100644 --- a/man/compar.ou.Rd +++ b/man/compar.ou.Rd @@ -9,8 +9,9 @@ compar.ou(x, phy, node = NULL, alpha = NULL) character.} \item{phy}{an object of class \code{"phylo"}.} \item{node}{a vector giving the number(s) of the node(s) where the - parameter `theta' (the character optimum) is assumed to change. By - default there is no change (same optimum thoughout lineages).} + parameter `theta' (the trait optimum) is assumed to change. The + node(s) can be specified with their labels if \code{phy} has node + labels. By default there is no change (same optimum thoughout lineages).} \item{alpha}{the value of \eqn{\alpha}{alpha} to be used when fitting the model. By default, this parameter is estimated (see details).} } @@ -76,28 +77,26 @@ compar.ou(x, phy, node = NULL, alpha = NULL) Hansen, T. F. (1997) Stabilizing selection and the comparative analysis of adaptation. \emph{Evolution}, \bold{51}, 1341--1351. } -\author{Emmanuel Paradis \email{Emmanuel.Paradis@mpl.ird.fr}} +\author{Emmanuel Paradis} \seealso{ \code{\link{ace}}, \code{\link{compar.lynch}}, \code{\link{corBrownian}}, \code{\link{corMartins}}, \code{\link{pic}} } \examples{ -\dontrun{ data(bird.orders) ### This is likely to give you estimates close to 0, 1, and 0 ### for alpha, sigma^2, and theta, respectively: -compar.ou(rnorm(23), bird.orders) +compar.ou(x <- rnorm(23), bird.orders) ### Much better with a fixed alpha: -compar.ou(rnorm(23), bird.orders, alpha = 0.1) +compar.ou(x, bird.orders, alpha = 0.1) ### Let us 'mimick' the effect of different optima ### for the two clades of birds... x <- c(rnorm(5, 0), rnorm(18, 5)) ### ... the model with two optima: -compar.ou(x, bird.orders, node = -2, alpha = .1) +compar.ou(x, bird.orders, node = 25, alpha = .1) ### ... and the model with a single optimum: compar.ou(x, bird.orders, node = NULL, alpha = .1) ### => Compare both models with the difference in deviances -## with follows a chi^2 with df = 1. -} +## which follows a chi^2 with df = 1. } \keyword{models}