# NOT RUN {
data(PedigWithErrors)
data(Phen)
Phen <- Phen
keep <- Phen$Indiv
Pedig <- prePed(PedigWithErrors, keep=keep, thisBreed="Hinterwaelder", lastNative=1970)
pKin <- pedIBD(Pedig, keep.only=keep)
pKinatN <- pedIBDatN(Pedig, thisBreed="Hinterwaelder", keep.only=keep)
cand <- candes(Phen, pKin=pKin, pKinatN=pKinatN)
cand$current[,c("Type", "Var", "Breed", "Val", "Name")]
# Type Var Breed Val Name
#BV trait BV missing -4.903304e-17 BV
#NC trait NC missing 4.235493e-01 NC
#pKin kinship pKin missing 2.197712e-02 pKin
#pKinatN nat. kin. pKinatN missing 7.531245e-02 pKinatN
### Shorthand:
cand$mean
# BV NC pKin pKinatN
#1 -4.903304e-17 0.4235493 0.02197712 0.07531245
cand$mean$pKin
#[1] 0.02197712
#Object cand can now be used as an argument of function opticont.
# }
Run the code above in your browser using DataLab