Function to simulate a step in a breeding scheme
ogc.mobps(
  population,
  animallist,
  relationship.matrix.ogc,
  depth.pedigree.ogc,
  bve.pedigree.error,
  ogc.target = "min.sKin",
  ogc.uniform = NULL,
  ogc.lb = NULL,
  ogc.ub = NULL,
  ogc.ub.sKin = NULL,
  ogc.lb.BV = NULL,
  ogc.ub.BV = NULL,
  ogc.eq.BV = NULL,
  ogc.ub.sKin.increase = NULL,
  ogc.lb.BV.increase = NULL,
  bve.p_i.list = NULL,
  miraculix = FALSE,
  miraculix.mult = FALSE,
  import.position.calculation = NULL,
  decodeOriginsU = decodeOriginsR,
  nbits = NULL,
  store.sparse = FALSE,
  verbose = TRUE,
  bit.storing = FALSE
)contributions of each individual in selection
Population list
List of individuals to include in ogc
Method to calculate relationship matrix for OGC (Default: "pedigree", alt: "vanRaden", "CE", "non_stand", "CE2", "CM")
Depth of the pedigree in generations (default: 7)
Set to FALSE to ignore/correct for any pedigree errors
Target of OGC (default: "min.sKin" - minimize inbreeding; alt: "max.BV" / "min.BV" - maximize genetic gain; both under constrains selected below)
This corresponds to the uniform constrain in optiSel
This corresponds to the lb constrain in optiSel
This corresponds to the ub constrain in optiSel
This corresponds to the ub.sKin constrain in optiSel
This corresponds to the lb.BV constrain in optiSel
This corresponds to the ub.BV constrain in optiSel
This corresponds to the eq.BV constrain in optiSel
This corresponds to the upper bound (current sKin + ogc.ub.sKin.increase) as ub.sKin in optiSel
This corresponds to the lower bound (current BV + ogc.lb.BV.increase) as lb.BV in optiSel
XXX
XXX
XXX
XXX
XXX
XXX
XXX
Set to FALSE to not display any prints
Set to TRUE if the MoBPS (not-miraculix! bit-storing is used)
population <- creating.diploid(nsnp=1000, nindi=100, n.additive = 10)
population <- breeding.diploid(population, breeding.size=100, selection.size=c(25,25),
ogc = TRUE)
Run the code above in your browser using DataLab