Post date: Sep 02, 2015 8:1:0 PM
I modified varne to write all bootstrap samples. I used this to calculate Ne with ETPIs for the L1, L2, and L3 lines. I calculated Ne as previously described, but then took quantiles over the bootstrap samples. Here is the R code:
Nmung<-scan("bootM13toM14.txt")
Nl1<-scan("bootM14toL1F100.txt")
Nl2<-scan("bootM14toL2F87.txt")
Nl3<-scan("bootM14toL3F85.txt")
t1<-231
t2<-194
t3<-189
## L1 Ne
r<-(t1-100)/t1
Ne<-(Nmung*Nl1*r-Nmung*Nl1)/(Nl1*r-Nmung)
quantile(Ne,probs=c(0.5,0.025,0.975))
# 50% 2.5% 97.5%
#253.7027 246.1102 260.9276
## L2 Ne
r<-(t2-87)/t2
Ne<-(Nmung*Nl2*r-Nmung*Nl2)/(Nl2*r-Nmung)
quantile(Ne,probs=c(0.5,0.025,0.975))
# 50% 2.5% 97.5%
#314.8783 304.5410 325.0229
## L3 Ne
r<-(t3-85)/t3
Ne<-(Nmung*Nl3*r-Nmung*Nl3)/(Nl3*r-Nmung)
quantile(Ne,probs=c(0.5,0.025,0.975))
# 50% 2.5% 97.5%
#335.4818 324.3981 347.7378