data(minnbreast)
bpeds <- with(minnbreast,
pedigree(id, fatherid, motherid, sex, affected=proband, famid=famid))
## pedigree with id=8
bped.id8 <- bpeds['8']
print(bped.id8)
## to plot:
#plot(bped.id8)
## show this pedigree with mixed zygosity quadruplets
rel8 <- data.frame(id1=c(137,138,139), id2=c(138,139,140), code=c(1,2,2))
bped.id8 <- with(minnbreast[minnbreast$famid==8,],
pedigree(id, fatherid, motherid, sex, affected=proband,
relation=rel8))
## to plot:
#plot(bped.id8)
## the 8th pedigree
bped8 <- bpeds[8]
print(bped8)
## to plot:
#plot(bped8, cex=.5)
Run the code above in your browser using DataLab