bioa.Amplicons.calc_score

Amplicons.calc_score(avg_cov, std_cov, avg_div, std_div)

Calculates the score of the amplicons set.

It is the sum of the logs of the following probabilities:

P(zscore < zscore(Min Coverage))

P(zscore < zscore(Min Diversity))

P(zscore < zscore(Min Non-Zero Diversity))

Parameters :

avg_cov : float

Average (mean) coverage.

std_cov : float

Standard deviation of coverage.

avg_div : float

Average (mean) diversity.

std_div : float

Standard deviation of diversity.

Returns :

score : float

Score is defined as: log(prob_cov_zscore) + log(prob_div_zscore) + log(prob_nzdiv_score)