## Not run: 
#    require(kinship2)
#    pedmus <- pedigree(musculus$id, musculus$dadid, musculus$momid, musculus$sex)
#    plot(pedmus, affected=musculus$diet)
# 
#    fam  <- makefamid(musculus$id, musculus$dadid, musculus$momid)
#    kmatrix <- makekinship(fam, musculus$id, musculus$dadid, musculus$momid)
#    kmatrix[1:5,1:15]
# ## End(Not run)
Run the code above in your browser using DataLab