### a simple example with some random data
data(bird.orders)
x <- rnorm(45) # the tree has 23 tips and 22 nodes
### the standard-error for x should be as large as
### the estimated parameter
yule.cov(bird.orders, ~ x)
### compare with the simple Yule model, eventually
### with a likelihood ratio test
yule(bird.orders)
### another example with a tree that has a multichotomy
### but we cannot run yule() because of this!
data(bird.families)
y <- rnorm(272) # 137 tips + 135 nodes
yule.cov(bird.families, ~ y)
Run the code above in your browser using DataCamp Workspace