### An example from Paradis (2003) using the avian orders:
data(bird.orders)
### Number of species in each order from Sibley and Monroe (1990):
S <- c(10, 47, 69, 214, 161, 17, 355, 51, 56, 10, 39, 152, 6, 143, 358, 103, 319, 23, 291, 313, 196, 1027, 5712)
xcut<-95*0.207407
out<-bd.groups.optim(bird.orders, S, xcut)
### Not accounting for the cutoff when defining the bird taxa estimates 0 extinction:
out1<-bd.groups.optim(bird.orders, S, 0)
### xcut=0 produces almost the same results as the following ape function
bd.ext(bird.orders, S)
Run the code above in your browser using DataLab