times

computationaltimes

s = sim.tree(seed=3)
j=1
m = seq(10,1000,by = 10)
times = vector('numeric',length = length(m))
for(i in m){
  time=proc.time()
  sim.srt(tree = s$tree.extant,pars = c(0.8,0.1,40),parallel = TRUE,n_trees = i)
  times[j] = get.time(time)
  j=j+1
}

qplot(m,times)

m1=m
times1=times
sum(times)

## [1] 4993.968

m = seq(2000,10000,by = 2000)
times = vector('numeric',length = length(m))
for(i in m){
  time=proc.time()
  sim.srt(tree = s$tree.extant,pars = c(0.8,0.1,40),parallel = TRUE,n_trees = i)
  times[j] = get.time(time)
  j=j+1
}

qplot(c(m1,m),c(times1,times[101:105]))

so,to…