X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=man%2Fcompar.ou.Rd;h=84be45483a7ad914fbaf47cac673ca7e1d4f4d0b;hb=80d1c453d63d6aec18f0b731d59918b99e189d86;hp=4c1c02bca4edd037b41bc86e56ac2350bde83ffd;hpb=90f18c75d642f56b020bc6e0cdd0c5949c1d9a1d;p=ape.git diff --git a/man/compar.ou.Rd b/man/compar.ou.Rd index 4c1c02b..84be454 100644 --- a/man/compar.ou.Rd +++ b/man/compar.ou.Rd @@ -86,9 +86,9 @@ compar.ou(x, phy, node = NULL, alpha = NULL) 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)) @@ -97,6 +97,6 @@ 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 -## wicth follows a chi^2 with df = 1. +## which follows a chi^2 with df = 1. } \keyword{models}