Recombination and new ms.to.bases in compare_f_estimators.r
code to plot expected D values (fig 1)