variant_filter_numerical.hpp

Go to the documentation of this file.

#ifndef GENESIS_POPULATION_FILTER_VARIANT_FILTER_NUMERICAL_H_
#define GENESIS_POPULATION_FILTER_VARIANT_FILTER_NUMERICAL_H_
 
/*
    Genesis - A toolkit for working with phylogenetic data.
    Copyright (C) 2014-2024 Lucas Czech
 
    This program is free software: you can redistribute it and/or modify
    it under the terms of the GNU General Public License as published by
    the Free Software Foundation, either version 3 of the License, or
    (at your option) any later version.
 
    This program is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    GNU General Public License for more details.
 
    You should have received a copy of the GNU General Public License
    along with this program.  If not, see <http://www.gnu.org/licenses/>.
 
    Contact:
    Lucas Czech <lucas.czech@sund.ku.dk>
    University of Copenhagen, Globe Institute, Section for GeoGenetics
    Oster Voldgade 5-7, 1350 Copenhagen K, Denmark
*/
 
#include "genesis/population/sample_counts.hpp"
#include "genesis/population/filter/sample_counts_filter_numerical.hpp"
#include "genesis/population/filter/sample_counts_filter.hpp"
#include "genesis/population/filter/variant_filter.hpp"
#include "genesis/population/variant.hpp"
 
#include <functional>
#include <iosfwd>
#include <memory>
#include <stdexcept>
#include <string>
#include <utility>
#include <vector>
 
namespace genesis {
namespace population {
 
// =================================================================================================
//     Variant Filter Numerical Params
// =================================================================================================
 
struct VariantFilterNumericalParams
{
    // -------------------------------------------
    //     Numeric
    // -------------------------------------------
 
    size_t min_read_depth = 0;
 
    size_t max_read_depth = 0;
 
    size_t deletions_count_limit = 0;
 
    // -------------------------------------------
    //     SNP vs Invariant
    // -------------------------------------------
 
    bool only_snps = false;
 
    bool only_biallelic_snps = false;
 
    size_t snp_min_count = 0;
 
    size_t snp_max_count = 0;
 
    double snp_min_allele_frequency = 0.0;
};
 
// =================================================================================================
//     Variant Filter Numerical Functions
// =================================================================================================
 
// --------------------------------------------------------------------------------------
//     apply_variant_filter_numerical
// --------------------------------------------------------------------------------------
 
bool apply_variant_filter_numerical(
    Variant& variant,
    VariantFilterNumericalParams const& params,
    VariantFilterStats& stats
);
 
bool apply_variant_filter_numerical(
    Variant& variant,
    VariantFilterNumericalParams const& params
);
 
// --------------------------------------------------------------------------------------
//     make_variant_filter_numerical_excluding
// --------------------------------------------------------------------------------------
 
inline std::function<bool(Variant&)> make_variant_filter_numerical_excluding(
    VariantFilterNumericalParams const& params
) {
    return [params]( Variant& variant ){
        return apply_variant_filter_numerical( variant, params );
    };
}
 
inline std::function<bool(Variant&)> make_variant_filter_numerical_excluding(
    VariantFilterNumericalParams const& params,
    VariantFilterStats& stats
) {
    return [params, &stats]( Variant& variant ){
        return apply_variant_filter_numerical( variant, params, stats );
    };
}
 
// --------------------------------------------------------------------------------------
//     make_variant_filter_numerical_tagging
// --------------------------------------------------------------------------------------
 
inline std::function<void(Variant&)> make_variant_filter_numerical_tagging(
    VariantFilterNumericalParams const& params
) {
    return [params]( Variant& variant ){
        apply_variant_filter_numerical( variant, params );
    };
}
 
inline std::function<void(Variant&)> make_variant_filter_numerical_tagging(
    VariantFilterNumericalParams const& params,
    VariantFilterStats& stats
) {
    return [params, &stats]( Variant& variant ){
        apply_variant_filter_numerical( variant, params, stats );
    };
}
 
inline std::function<void(Variant&)> make_variant_filter_numerical_tagging(
    SampleCountsFilterNumericalParams const& sample_count_params,
    VariantFilterNumericalParams const& variant_params,
    bool all_need_pass = false
) {
    return [sample_count_params, variant_params, all_need_pass]( Variant& variant ){
        apply_sample_counts_filter_numerical(
            variant, sample_count_params, all_need_pass
        );
        apply_variant_filter_numerical( variant, variant_params );
    };
}
 
inline std::function<void(Variant&)> make_variant_filter_numerical_tagging(
    SampleCountsFilterNumericalParams const& sample_count_params,
    VariantFilterNumericalParams const& variant_params,
    VariantFilterStats& variant_stats,
    SampleCountsFilterStats& sample_count_stats,
    bool all_need_pass = false
) {
    return [variant_params, sample_count_params, &variant_stats, &sample_count_stats, all_need_pass](
        Variant& variant
    ){
        apply_sample_counts_filter_numerical(
            variant, sample_count_params, variant_stats, sample_count_stats, all_need_pass
        );
        apply_variant_filter_numerical( variant, variant_params, variant_stats );
    };
}
 
} // namespace population
} // namespace genesis
 
#endif // include guard