GWAS做关联区域分析时,根据p值获取关联区域中,为什么要用如下语句以合并P值(此语句是2023.4.10日添加的) cat pvalue.bed region.bed |sort -k1,1 -k2,2n >all.region.bed 合并P值在获取关联区域的作用是什么,可否像之前录像讲的那样不用这个语句吗?

请先 登录 后评论

2 个回答

这一步不是在合并p值,是在合并snp位点和上下游10kb区域的位置信息,需要做

请先 登录 后评论
omicsgene - 生物信息
擅长:重测序,遗传进化,转录组,GWAS

bedtools merge 加上P值是因为,上一步flank 取的是p值两边的region,如果不加pvalue.bed 会出现空档: https://bedtools.readthedocs.io/en/latest/content/tools/merge.html

attachments-2025-03-c3BG0j0P67d379c95314c.png

请先 登录 后评论