from PhyNetPy.BirthDeath import CBDP
# a birth rate of .05, an extinction rate of .01,
# 6 extant taxa, and the default 100% sampling rate
model = CBDP(.05, .01, 6)
model.generateTree() #to see the network, call .printGraph() on the result
model.generateTrees(100) #to generate 100 trees
#Use of the optional sampling rate parameter
model2 = CBDP(.05, .01, 6, .5)
#to generate trees under the Yule model, set extinction rate to 0
yuleModel = CBDP(.05, 0, 6)