def cal_max_dup_tags ( genome_size, tags_number, p=1e-5 ):
"""Calculate the maximum duplicated tag number based on genome
size, total tag number and a p-value based on binomial
distribution. Brute force algorithm to calculate reverse CDF no
more than MAX_LAMBDA(100000).
"""
return binomial_cdf_inv(1-p,tags_number,1.0/genome_size)
def cal_max_dup_tags ( genome_size, tags_number, p=1e-5 ):
"""Calculate the maximum duplicated tag number based on genome
size, total tag number and a p-value based on binomial
distribution. Brute force algorithm to calculate reverse CDF no
more than MAX_LAMBDA(100000).
"""
return binomial_cdf_inv(1-p,tags_number,1.0/genome_size)