sumstat_jsfs(name = "jsfs", populations = c(1, 2), per_locus = FALSE,
transformation = identity)
TRUE
, the JSFS is return for each locus instead
of globally. In this case, the result is a list, where each entry is the
JSFS for the corresponding locus.