大家好,我是邓飞,今天介绍一下,知道群体的单倍型数据,知道群体的表型数据,如何计算单倍型和表型数据的显著性分析,以及做一个小提琴图?1,准备单倍型群体的vcf文件将单倍型,按照区间,提取群体的vcf文件,如下:2,整理表型数据的txt文件如下:3,使用R包geneHapR如下:4,读取vcf文件如下:5,读取表型数据如下:6,计算群体的单倍型(具体到个体)如下:7,将群体的单倍型数据保存到本地如下:8,单倍型和表型的显著性分析可视化如下:图表:有结果可以得出,H001的单倍型的y1表现比H002的单倍型要好,H001是优势单倍型。群体中共有119个H001,9个H002。示例数据和代码:数据和代码放到星球上面了: