SIMplyBee (version 0.4.1)

isCsdHeterozygous: Test if individuals are heterozygous at the csd locus

Description

Level 0 function that returns if individuals of a population are heterozygous at the csd locus. See SimParamBee for more information about the csd locus.

Usage

isCsdHeterozygous(pop, simParamBee = NULL)

Value

logical

Arguments

pop

Pop-class

simParamBee

SimParamBee, global simulation parameters

Details

We could expand isCsdHeterozygous to work also with Colony-class and MultiColony-class if needed

Examples

founderGenomes <- quickHaplo(nInd = 8, nChr = 1, segSites = 100)
SP <- SimParamBee$new(founderGenomes)
SP$nThreads = 1L
basePop <- createVirginQueens(founderGenomes)

drones <- createDrones(x = basePop[1], nInd = 1000)
droneGroups <- pullDroneGroupsFromDCA(drones, n = 10, nDrones = nFathersPoisson)

# Create a Colony class
colony <- createColony(x = basePop[2])
colony <- cross(colony, drones = droneGroups[[1]])
colony <- buildUp(x = colony, nWorkers = 6, nDrones = 3)
colony <- addVirginQueens(x = colony, nInd = 4)

# Use isCsdHeterozygous on a Population
isCsdHeterozygous(getQueen(colony))
isCsdHeterozygous(getWorkers(colony))