Closed richelbilderbeek closed 6 years ago
Let's first put it in raket :+1:
Added it on the develop
branch, emailed @J-Damhuis and @Tomdkkr that they may complete the function. Would be a good exercise for them :angel:
In practice, this already works by using a small difference between speciation-initiation and extinction rate, instead of using pbd_expected_n_extant
.
AFAIK, there is no equation of the expected number of extant species for the PBD model [*], i.e.:
I suggest to write this function for the PBD package (hence the
pbd_
prefix) by simulating PBD trees and averaging. Clumsy, but AFAIK, the best we can do.[*] To be precise, [1] claims it is too hard, and I could not find it in later papers.
References
See also the cer_raket issue