R/sim.rateshift.taxa.R

sim.rateshift.taxa <-function(n,numbsim,lambda,mu,frac,times,complete=TRUE,K=0,norm=TRUE){
	out<-lapply(1:numbsim,sim.rateshift.taxa.help,n=n,lambda=lambda,mu=mu,frac=frac,times=times,complete=complete,K=K,norm=norm)
	out1<-lapply(out,function(x){ x[[1]][[1]]})
	out2<-sapply(out,function(x){ x[[2]]})
	for (i in 1:length(out1)){
		out1[[i]]$root.edge<-out2[i]-max(getx(out1[[i]],sersampling=1)[,1])
	}
	#out3<-list(out1,out2)
	#out3
	out1
	}

Try the TreeSim package in your browser

Any scripts or data that you put into this service are public.

TreeSim documentation built on May 2, 2019, 3:23 a.m.