如何使用Kinship2软件包将多个配偶添加到我的血统书中?

时间:2018-11-24 15:15:15

标签: r inheritance bioinformatics genetics family-tree

我正在尝试使用Kinship2 R软件包绘制非常复杂的血统书。大多数家庭表现出很高的血缘关系和配偶关系,我该如何在我的血统书上证明这一点。

这是我当前的情节

##data
df <- data.frame(id=c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29),
             sex=c(1,2,1,2,2,2,1,2,2,1,1,1,2,2,2,2,2,2,2,2,2,1,2,2,1,2,1,1,1),
             dadid=c(0,0,1,1,1,1,1,0,0,0,0,0,0,3,3,10,10,10,11,11,11,11,12,12,7,7,7,0,28),
             momid=c(0,0,2,2,2,2,2,0,0,0,0,0,0,8,8,4,4,4,5,5,5,5,6,6,13,13,13,0,16),
             status=c(0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0))


ped <- pedigree(id = df$id, dadid = df$dadid, momid = df$momid, sex = df$sex, status=df$status)
plot(ped)

我想在28和16之间添加近亲双喜,将第二个配偶添加到3,以及为14添加一个配偶。

我该怎么做?

Current plot

0 个答案:

没有答案