# NOT RUN {
GenoS <- SimGeno(Ped.griffin, nSnp=400, ParMis=0.4)
griffin.sex <- sapply(Ped.griffin$ID,
function(x) substr(x, start=nchar(x), stop=nchar(x)))
LH_griffin <- data.frame(ID = Ped_griffin$ID,
Sex = ifelse(griffin.sex=="F", 1, 2),
BirthYear = Ped_griffin$BY)
SeqOUT_griffin <- sequoia(GenoS, LH_griffin,
MaxSibIter = 10,
args.AP = list(Smooth = FALSE))
# }
Run the code above in your browser using DataLab