基于亲本基因型推断杂种基因型

时间:2017-03-27 06:21:32

标签: r marker hybrid

我有杂交亲本的标记信息,并希望根据其亲本基因型推断出杂种基因型。以下是我的数据parent marker info的子集示例。在每个基因型列中:  1)0 - 一个等位基因的纯合子;  2)1 - 另一个等位基因的纯合子;  3)2 - 杂合子和  4)&​​#34; 999"表示对国家缺失或不自信。

我想以这种方式编码混合基因型hybrid marker coding method

如何编写R代码来构建杂交的基因型,这些基因型是从每个父母与其他父母交叉而不是相互交叉而衍生的,例如parent1 / parent2,parent1 / parent3和parent3 / parent5?

以下是生成子样本的代码。

Marker <- c(1001, 1002, 1003, 1004, 1005, 1006, 1007, 1008, 1009, 1010, 1011, 1012, 1013, 1014)
Parent1 <- c(0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 2, 0)
Parent2 <- c(0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0)
Parent3 <- c(0, 0, 0, 0, 0, 0, 0, 999, 0, 999, 0, 0, 0, 0)
Parent4 <- c(999, 0, 0, 0, 999, 0, 0, 0, 0, 1, 0, 0, 0, 0)
Parent5 <- c(0, 0, 0, 0, 0, 0, 2, 999, 0, 999, 0, 0, 0, 0)
Parent6 <- c(0, 0, 0, 0, 0, 2, 0, 0, 0, 1, 0, 0, 0, 0)
Parent7 <- c(999, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0)
Parent8 <- c(999, 0, 0, 0, 999, 0, 0, 0, 0, 1, 0, 0, 0, 0)
Genotypes <- cbind(Marker, Parent1, Parent2, Parent3, Parent4, Parent5, Parent6, Parent7, Parent8)

非常感谢,

0 个答案:

没有答案