data(bird.orders)
plot(compute.brlen(bird.orders, 1))
plot(compute.brlen(bird.orders, runif, min = 0, max = 5))
split.screen(c(2,2))
screen(1)
plot(compute.brlen(bird.orders, power=1), main=expression(rho==1))
screen(2)
plot(compute.brlen(bird.orders, power=3), main=expression(rho==3))
screen(3)
plot(compute.brlen(bird.orders, power=0.5), main=expression(rho==0.5))
screen(4)
plot(compute.brlen(bird.orders, power=0.1), main=expression(rho==0.1))
close.screen(all=TRUE)
Run the code above in your browser using DataLab