R/sim2.bd.reverse.R

In tanja819/TreeSim: Simulating Phylogenetic Trees

sim2.bd.reverse <- function(n,numbsim,lambda,mu,fast){	phy <- list()
	time<-vector()
	frac<-1
	for (j in 1:numbsim){
		if (fast)
		  temp <- fast.tree(n,lambda,mu,frac)
		else
	    temp <- sim2.bd.reverse.single(n,lambda,mu,frac)
		phy <- c(phy, list(temp[[1]]))
		time<-c(time,temp[[2]])
		}
	phy2<-list(phy,time)
	phy2
}