Should we condition to survival?

R=NULL
R=NULL
time=proc.time()
for(i in 1:10000){
  s=sim_phyl()
  R[i]=s$r
}
t2 = proc.time()-time
paste(t2[3]/60,'min')

## [1] "30.8261833333333 min"

sum(R)

## [1] 2

R=NULL
time=proc.time()
for(i in 1:10000){
  s=sim_phyl(mu0=0.2)
  R[i]=s$r
}
t2 = proc.time()-time
paste(t2[3]/60,'min')

## [1] "41.00505 min"

sum(R)

## [1] 2

R=NULL
time=proc.time()
for(i in 1:10000){
  s=sim_phyl(mu0=0.4)
  R[i]=s$r
}
t2 = proc.time()-time
paste(t2[3]/60,'min')

## [1] "54.9860666666667 min"

sum(R)

## [1] 81

Comparison. Conditioning to survival

time = proc.time()
s=sim_phyl()
t2 = proc.time()-time
paste(t2[3],'sec')

## [1] "0.203 sec"

time = proc.time()
llik(pars=c(0.8,0.1,40),n=s$n,E=s$E,t=s$wt)

## [1] 242.5088

t2 = proc.time()-time
paste(t2[3],'sec')

## [1] "0.00300000000000011 sec"

time = proc.time()
llik(pars=c(0.8,0.1,40),n=s$n,E=s$E,t=s$wt,conditionToSurvival = TRUE)

## [1] 242.7147

t2 = proc.time()-time
paste(t2[3],'sec')

## [1] "0.032 sec"

time = proc.time()
s=sim_phyl()
t2 = proc.time()-time
paste(t2[3],'sec')

## [1] "0.256 sec"

time = proc.time()
llik(pars=c(0.8,0.1,40),n=s$n,E=s$E,t=s$wt)

## [1] 330.1948

t2 = proc.time()-time
paste(t2[3],'sec')

## [1] "0.00300000000000011 sec"

time = proc.time()
llik(pars=c(0.8,0.1,40),n=s$n,E=s$E,t=s$wt,conditionToSurvival = TRUE)

## [1] 330.4473

t2 = proc.time()-time
paste(t2[3],'sec')

## [1] "0.0349999999999997 sec"