4 #include <Sequence/summstats/ld.hpp> 5 #include <Sequence/VariantMatrix.hpp> 6 #include <Sequence/VariantMatrixViews.hpp> 10 TwoLocusCounts::TwoLocusCounts(std::int8_t i_, std::int8_t j_,
int n_)
11 : i{ i_ }, j{ j_ }, n{ n_ }
15 std::vector<TwoLocusCounts>
16 two_locus_haplotype_counts(
const VariantMatrix& m, std::size_t sitei,
17 const std::size_t sitej,
18 const bool skip_missing)
22 std::vector<TwoLocusCounts> rv;
23 for (
auto i = ri.begin(), j = rj.begin(); i < ri.end(); ++i, ++j)
25 if (!skip_missing || ((*i >= 0 || *j >= 0) && skip_missing))
28 = std::find_if(rv.begin(), rv.end(),
29 [i, j](
const TwoLocusCounts& t) {
30 return t.i == *i && t.j == *j;
32 if (exists == rv.end())
34 rv.emplace_back(*i, *j, 1);
The namespace in which this library resides.
ConstRowView get_ConstRowView(const VariantMatrix &m, const std::size_t row)
Return a ConstRowView from VariantMatrix m at row row. std::out_of_range is thrown if row is out of r...