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