Level 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
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.