R/sim.genetree.R

sim.genetree <- function(n,numbsim){
statistics<-vector()
for (i in 1:numbsim){
	statisticstemp<-c(1:4)*0
	lineages<-c(1:n)*0+1
	while (length(lineages)>1){
			merge<-sort(sample(1:length(lineages),2))
			statisticstemp<-stats(lineages,merge[1],merge[2],statisticstemp)
			lineages[merge[1]]<-lineages[merge[1]]+lineages[merge[2]]
			lineages<-lineages[-merge[2]]	
	}
	statistics<-rbind(statistics, statisticstemp)
}
rownames(statistics)<-c(1:length(statistics[,1]))
colnames(statistics)<-c("Colless","s","Sackin","cherries")
statistics
}

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.