最近有一个老师在飞哥的学习圈子中问了一个问题,恰好是我擅长的领域,我还在我编写的包(learnasreml)中定义了几个相关的函数,就解答了一下,顺便写篇博客。
邓老师,想问一下如何根据猪场子代父母代数据找出它们之间的系谱关系图呢,最好是方便统计的那种
这个问题,可以从以下三个方面解答:
数据介绍:
library(data.table)
library(tidyverse)
library(nadiv)
library(learnasreml)
ped = fread("d:/pedigree_test.csv")
head(ped)
dim(ped)
下面,介绍一下这三种方法的具体实施方法。
首先,这是三列系谱数据,我们需要将其转为15列系谱数据,包括每个个体的三个世代的数据。用的是我写的包learnasreml
中的pedigree_3_to_15
函数:
如果我们想查看1092个体的三代系谱记录,可以找到ID列是1092的行,分别给出个体的三代系谱数据:
如果我们想要筛选个体记录比较完整的个体,可以以DDD为条件,筛选一下:
然后,将其变为3列系谱。这里用到的是pedigree_15_to_3
上一步中,得到的系谱还有937条,我们对着937条系谱进行可视化。
首先,计算亲缘关系A矩阵:
Amat = makeA(pped1) %>% as.matrix()
Amat[1:10,1:10]
对其进行可视化:
或者单独提取出聚类分析图:
也可以做成圈图的聚类可视化:
当然,上面品种较多,可视化不太直观,如果挑选出感兴趣的少数个体(比如20~50个),用上面的方法是非常方便的。
系谱可视化,可以使用visPedigree
包中的visPedigree
library(visPedigree)
pped2 = visPedigree::tidyped(ped1)
visped(pped2)
可以非常直观看出个体间的关系。
以上就是根据系谱查看个体间关系的三种方法,希望对大家有所帮助。