# NOT RUN {
data(PedigWithErrors)
data(Phen)
keep <- Phen$Indiv
Pedig <- prePed(PedigWithErrors, keep=keep, thisBreed="Hinterwaelder", lastNative=1970)
cand <- candes(phen=Phen, fA=pedIBD(Pedig, keep.only=keep))
con <- list(ub.fA=0.03, uniform="female")
Offspring <- opticont("max.BV", cand, con, trace = FALSE)
N <- 150
Candidate <- Offspring$parent
Candidate$nOff <- noffspring(Candidate, N)$nOff
sum(Candidate$nOff[Candidate$Sex=="male"])
#[1] 150
sum(Candidate$nOff[Candidate$Sex=="female"])
#[1] 150
round(2*N*Candidate$oc-Candidate$nOff, 2)
# }
Run the code above in your browser using DataLab