Skip to content

Instantly share code, notes, and snippets.

Avatar

Giulio Valentino Dalla Riva gvdr

  • School of Mathematics and Statistics, University of Canterbury
  • Christchurch, New Zealand
View GitHub Profile
@gvdr
gvdr / do_the_mfob
Last active Dec 19, 2015
A Markovian trait based fields of bullets method. Based on dendropy and random, mfob is a small module that permits a fast implementation of a fields of bullets model where the probability of extinction depends on a binary trait evolved from the root as a markov process.
View do_the_mfob
"""
Created on Fri Jun 21 10:52:11 2013
@author: Giulio Valentino Dalla Riva
@email: gvd16@uclive.ac.nz
"""
#generate a tree list following a birth and death process
#the process is conditioned to N leaves (with success)
#the list is long rang
def generate_trees_bdN(birth=0.57721,death=0.130357,N=49,rang=314):