Properties of bias and variance of two multiallelic estimators of F(ST)

This study compares theoretical and simulated properties of two estimators of fixation indices F(ST) for multiallelic data from a single locus. The estimators are due to Weir and Cockerham (straight theta(WC)) and to Robertson and Hill (;straight theta(RH)), respectively. Both estimators are linear combinations of biallelic estimators which differ in the way frequent and rare alleles are weighted. Consequently, their sampling properties differ as far as bias and variance are concerned. In the infinite island model at migration-drift equilibrium, in the case of one multiallelic locus, we analytically approximate the bias of the two estimators and show that;straight theta(WC) is nearly unbiased, whereas;straight theta(RH) is negatively biased. Hence, we propose a correction of bias of the latter. Moreover, we reexamine the properties of variance of the initial estimators: due to their construction, their variances are minimal, each over different parameter ranges,;straight theta(RH) performing better for low differentiation and;straight theta(WC) for high differentiation. These theoretical properties are confirmed by simulations, which show that our correction of bias of;straight theta(RH) performs well and does not affect its property of minimal variance for low differentiation. Hence, we advocate the use of;straight theta(WC) for high values of differentiation, and the use of;straight theta(RH), with our correction of bias, for low or moderate differentiation.