主要参数:
M:为0-1-2编码的基因型数据,正常来说,都为0或者2,不应该有杂合的位点
ped:需要杂交的系谱,ID,亲本1,亲本2
注意事项:
亲本不允许有缺失值,不允许有杂合位点,这样生成的F1...),each=2), P2 = rep(c("ID4","ID5"),3))
ped
library(tidyverse)
ped = ped %>% mutate(ID = str_c(P1,"_"...,P2)) %>% select(ID,P1,P2)
ped
生成F1的代码:
synthetic.cross(M012,ped = ped,indiv = "ID",mother = "P1",father...这里为“exact”,即无法准确推测的设置为缺失,还有其它的选项可以选择。...结果文件: