R/sample_a_d_a.R

In fcampelo/TaxonSampling: Hierarchical Taxonomic Sampling for Maximized Diversity at All Levels

sample_a_d_a <- function(m, m_i, childrenCount, childrenCountSpp, ...){
  # First round
  # If there's enough to sample at least one entry per lineage,
  # do it to increase diversity
  if (m > 0 & sum(childrenCount > m_i) <= m){
    child       <- names(childrenCount)[childrenCount > m_i]
    m_i[child]  <- m_i[child] + 1
    m           <- m - length(child)
  }

  # Subsequent rounds
  m_i <- sample_a_d_y(m, m_i, childrenCount, childrenCountSpp, ...)

  return(m_i)
}