Start to implement likelihood in speciation model