R/allGenos.R

In relMix: Relationship Inference for DNA Mixtures

#' Find all possible genotypes
#'
#' Finds all possible genotypes based on input alleles.
#' @param alleles Vector of input alleles, numeric or character
#' @return Matrix of all possible genotypes, one row per genotype
#' @author Guro Dorum
#' @export
#Reworked from paramlink's allGenotypes
allGenos <- function(alleles){

  n <- length(alleles)
  if (n < 2)
    return(matrix(nrow = 0, ncol = 2))
  v1 <- rep.int(alleles[-n],(n-1):1)
  v2 <-  NULL
  for (i in 2:n) v2 <- c(v2, alleles[i:n])
  rbind(cbind(v1=alleles, v2=alleles),cbind(v1, v2, deparse.level = 0))
}