R/Functions_L1_Pop.R
    combineBeeGametesHaploDiploid.RdLevel 1 function that produces diploid offspring from a mated queen.
  Drones are haploid, while the queen is diploid, so we first generate gametes
  (with recombination) from her and merge them with drone genomes (=gametes),
  where we randomly re-sample drones to get the desired number of progeny.
  This is an utility function, and you most likely want to use the
  cross function.
combineBeeGametesHaploDiploid(queen, drones, nProgeny = 1, simParamBee = NULL)Pop-class, with a single diploid individual
Pop-class, with one or more haploid individual(s)
integer, number of progeny to create per cross
SimParamBee, global simulation parameters
Pop-class with diploid individuals
This would be the right approach to handle haplo-diploid inheritance in bees, but it causes a raft of downstream issues, since AlphaSimR assumes that individuals have the same ploidy. Hence, we don't use this function.