# NOT RUN {
data(PedigWithErrors)
data(Phen)
Phen[1:5,]
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)
Kin2 <- kinlist(fA=fA, fD=fD)
data(Kin)
identical(Kin, Kin2)
#[1] TRUE
# }
Run the code above in your browser using DataLab