11 #ifndef ST_QUALITY_HXX
12 #define ST_QUALITY_HXX
18 #define st_assert(cond) arb_assert(cond)
26 void add(
double val) {
52 void add(
double val) {
54 square_sum += val*val;
59 square_sum += other.square_sum;
80 void add(
size_t range,
double probability) {
82 likelihood[range].
add(probability);
85 st_assert(seq_rel_pos>=0.0 && seq_rel_pos<1.0);
86 add(seq_rel_pos*ranges, probability);
120 #error st_quality.hxx included twice
121 #endif // ST_QUALITY_HXX
LikelihoodRanges stat_cons
size_t overall_range_count()
LikelihoodRanges stat_five
void add_relative(double seq_rel_pos, double probability)
ColumnQualityInfo(int seq_len, int bucket_size)
static char * alignment_name
double get_median() const
size_t get_range_count() const
GB_ERROR st_ml_check_sequence_quality(GBDATA *gb_main, const char *tree_name, const char *alignment_name, ColumnStat *colstat, const WeightedFilter *weighted_filter, int bucket_size, int marked_only, st_report_enum report, const char *dest_field)
void add(const Sampler &other)
void add(const VarianceSampler &other)
double get_variance(double median) const
char * generate_string(Sampler &t_values)
LikelihoodRanges stat_half
double get_variance() const
LikelihoodRanges stat_user
void add(size_t range, double probability)
VarianceSampler summarize_all_ranges()
LikelihoodRanges(size_t ranges)