我有杂交亲本的标记信息,并希望根据其亲本基因型推断出杂种基因型。以下是我的数据parent marker info的子集示例。在每个基因型列中: 1)0 - 一个等位基因的纯合子; 2)1 - 另一个等位基因的纯合子; 3)2 - 杂合子和 4)" 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)
非常感谢,