R/sim.bd.age.loop.R

In TreeSim: Simulating Phylogenetic Trees

sim.bd.age.loop <-
function(age,numbsim,lambda,mu,frac=1,mrca=FALSE,complete=TRUE,K) {
	if (mrca == FALSE) {
		phy <- sim2.bd.age(age,numbsim,lambda,mu,K)
		if (complete==FALSE) {
			phy <- reconstructed.age(phy,frac)
		}
	} else {
		if (complete == TRUE) {
			phy <- sim2.bd.mrca(age,numbsim,lambda,mu,K)
		}	else {
			phy <- sim2.bd.mrca.reconst(age,numbsim,lambda,mu,frac,K)
			}
	}
	phy
	}