ped <- list(
id = c(1, 2, 3),
father = c(0, 0, 1),
mother = c(0, 0, 2),
sex = c(1, 2, 1), # 1 = male, 2 = female
genotype = c("AA", "AB", "BB")
)
bplotped <- list(
id = c(4, 5),
father = c(0, 4),
mother = c(0, 5),
sex = c(1, 2),
genotype = c("AB", "AA")
)
Run the code above in your browser using DataLab