An Approximate Bayesian Estimator Suggests Strong, Recurrent Selective Sweeps in Drosophila
Figure 4
Approximate Bayesian estimation of the strength and rate of selection as well as the neutral θ, when estimation is based upon the means and SDs of π, S, θH and ZnS.
The model is one in which s and 2Nλ are fixed. For the strong selection case s = 1.0E−02, and 2Nλ = 2.0E−05, for weak selection s = 1.0E−04, and 2Nλ = 2.0E−03. ρ = 0.1/site and θ = 0.01/site. Shown are the distributions of 1000 MAP estimates. The dotted lines indicate the true values. The distributions for 10 50 kb region datasets are given in black, and for 1000 500 bp datasets in gray. As shown, the use of these multiple summary statistics improves estimation relative to π alone (Figure S1), reducing the RMSEs (Table S1).