search

MomentsLD / moments

MIT License
10 stars 3 forks source link

Add F-statistic functions to Spectrum_mod #183

Closed apragsdale closed 6 months ago

apragsdale commented 7 months ago

An LDstats object has class functions to compute f2, f3 and f4 from the heterozygosity statistics. These are not currently available as SFS class functions. This PR adds these and tests them against the output of those same statistics as computed from the LD functions.

apragsdale commented 7 months ago

With the checks for valid input and corner cases, these functions are a bit bulky. They could be put in their own functions in Util and called from the class function.

We also still need a bunch of tests of valid input, handling the index order of populations if they are "out of order", etc.