Exemplo n.º 1
0
// check for alleles not in population database
bool
AlleleSet::checkBackground(PMF<Allele> const &background) const
{
    bool ret = true;

	std::vector< PMF<Allele> >::const_iterator ip;
	for(ip = m_pmfs.begin(); ip != m_pmfs.end(); ++ip)
	{
		PMF<Allele>::const_iterator ia;
		for(ia = ip->begin(); ia != ip->end(); ++ia)
		{
			if (background.find(ia->first) == background.end())
			{
				// Allele not in database.

				warn << startl << "allele not in population database: " << ia->first.string() << std::endl;
				ret = false;
			}
		}
	}
	return ret;
}