R/fstats.R

In heavywatal/rwtl: Personal R package

#' Fst estimators
#'
#' @description
#' `fst_HBK()` calculates Kst by Hudson, Boos, and Kaplan (1992).
#' @param within,between mean branch length or diversity
#' @param n number of subpopulations
#' @rdname fstats
#' @export
fst_HBK = function(within, between, n = 2) {
  (between - within) / (between + within / (n - 1))
}

#' @description
#' `fst_HSM()` calculates Fst by Hudson, Slatkin, and Maddison (1992).
#' @rdname fstats
#' @export
fst_HSM = function(within, between) {
  1.0 - within / between
}