R/Functions_L0_auxilary.R
isCsdHeterozygous.Rd
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.
isCsdHeterozygous(pop, simParamBee = NULL)
Pop-class
SimParamBee
, global simulation parameters
logical
We could expand isCsdHeterozygous
to work also with
Colony-class
and MultiColony-class
if needed
founderGenomes <- quickHaplo(nInd = 8, nChr = 1, segSites = 100)
SP <- SimParamBee$new(founderGenomes)
basePop <- createVirginQueens(founderGenomes)
#> Error in get(x = "SP", envir = .GlobalEnv): object 'SP' not found
drones <- createDrones(x = basePop[1], nInd = 1000)
#> Error in get(x = "SP", envir = .GlobalEnv): object 'SP' not found
droneGroups <- pullDroneGroupsFromDCA(drones, n = 10, nDrones = nFathersPoisson)
#> Error in get(x = "SP", envir = .GlobalEnv): object 'SP' not found
# Create a Colony class
colony <- createColony(x = basePop[2])
#> Error in get(x = "SP", envir = .GlobalEnv): object 'SP' not found
colony <- cross(colony, drones = droneGroups[[1]])
#> Error in get(x = "SP", envir = .GlobalEnv): object 'SP' not found
colony <- buildUp(x = colony, nWorkers = 6, nDrones = 3)
#> Error in get(x = "SP", envir = .GlobalEnv): object 'SP' not found
colony <- addVirginQueens(x = colony, nInd = 4)
#> Error in get(x = "SP", envir = .GlobalEnv): object 'SP' not found
# Use isCsdHeterozygous on a Population
isCsdHeterozygous(getQueen(colony))
#> Error in get(x = "SP", envir = .GlobalEnv): object 'SP' not found
isCsdHeterozygous(getWorkers(colony))
#> Error in get(x = "SP", envir = .GlobalEnv): object 'SP' not found