2021-06-20

scaffold水平GWAS分析

自己记录一下学习生信过程中的坑。

使用plink软件进行分析,得到结果之后,由于SCAFFOLD组装条数多,在画曼哈顿图是无法使用所有SCAFFOLD进行,故可选根据阈值提取显著位点


提取gwas分析的显著位点

在根据显著位点提取其对应SCAFFOLD


第一个命令获得显著位点所在SCAFFOLD号,第二个命令生成提取脚本,最后提取

合并SCAFFOLD为一个文件


由于作图使用的R包要求要求染色体\scaffold号为数字,故提取编号进行处理


生成作图所需文件

作图


©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容