def multilocus_Gst_est(Ht_est, Hs_est):
    """Averages Ht_est and Hs_est across loci before calculating Gst."""
    Ht_est = sum(Ht_est) / float(len(Ht_est))
    Hs_est = sum(Hs_est) / float(len(Hs_est))
    if Ht_est == 0.0:
        return 0.0
    else:
        Gst_est = (Ht_est - Hs_est) / Ht_est
        return Gst_est
