R/p1.R

p1 <- function(t,l,m,rho){ rho*(l-m)^2 * exp(-(l-m)*t)/(rho*l+(l*(1-rho)-m)*exp(-(l-m)*t))^2 }
tanja819/TreePar documentation built on May 31, 2019, 2:57 a.m.