# NOT RUN {
data(PedigWithErrors)
data(Phen)
Phen <- Phen[80:140,]
keep <- Phen$Indiv
Pedig <- prePed(PedigWithErrors, keep=keep, thisBreed="Hinterwaelder", lastNative=1970)
fA <- pedIBD(Pedig, keep.only=keep)
fD <- pedIBDatN(Pedig, thisBreed="Hinterwaelder", keep.only=keep)
x <- pedIBDorM(Pedig, thisBreed="Hinterwaelder", keep.only=keep)
Kin <- kinlist(fA=fA, fB=x$pedIBDorM, fC=x$pedIBDorMM, fD=fD)
help.opticont(Kin, Phen)
# }
Run the code above in your browser using DataLab