R/addGroupOfUnusedAnimals.R

In rmsharp/nprcmanager: Genetic Tools for Colony Management

#' addGroupOfUnusedAnimals adds a group to the saved groups if needed
#'
## Copyright(c) 2017-2020 R. Mark Sharp
## This file is part of nprcgenekeepr

#' @return A list of groups, which are each lists of animal Ids that are unused
#' animals at the end of the iteration.
#'
#' @param savedGroupMembers list of groups of animals in the form of a vector
#' of animal Ids.
#' @param candidates character vector of IDs of the animals available for
#' use in the group.
#' @param ped dataframe that is the `Pedigree`. It contains pedigree
#' information including the IDs listed in \code{candidates}.
#' @param minAge integer value indicating the minimum age to consider in group
#' formation. Pairwise kinships involving an animal of this age or younger will
#'  be ignored. Default is 1 year.
#' @param harem logical variable when set to \code{TRUE}, the formed groups
#' have a single male at least \code{minAge} old.
#'
addGroupOfUnusedAnimals <- function(savedGroupMembers, candidates, ped,
                                    minAge, harem) {
  if (harem) { # Sires were added to groupMembers
    candidates <- removePotentialSires(candidates, minAge, ped)
  }

  # Adding a group for the unused animals
  n <- length(savedGroupMembers) + 1
  savedGroupMembers[[n]] <-
    ifelse(isEmpty(setdiff(candidates, unlist(savedGroupMembers))),
           c(NA), list(setdiff(candidates, unlist(savedGroupMembers))))[[1]]
  savedGroupMembers
}