gwas脚本中的tassel计算变异位点p值的命令中,结果会自动过滤掉indel,如何保留indel,并计算其p值?

咱们gwas分析的脚本中的tassel的run_pipeline.pl工具为什么在计算变异位点p值的时候,会自动过滤掉indel。我想同时计算indel的p值,需要怎么做?代码如下:

run_pipeline.pl -Xmx30G -fork1 -h hapmap.ordered.hmp.txt \

    -fork2 -t $tassel_trait -fork3 -q $structure -excludeLastTrait \

   -combine5 -input1 -input2 -input3 -intersect -FixedEffectLMPlugin \

   -endPlugin -export hd_GLM_structure

是否可以通过在这条命令中添加什么参数,让它不过滤indel,同时计算indel的p值呢?
请先 登录 后评论

1 个回答

omicsgene - 生物信息
擅长:重测序,遗传进化,转录组,GWAS

目前课程的流程不支持indel数据做gwas关联分析。由于高通量测序得到的SNP数量本身是过量的;增加INDEL数据并不一定会增加关联效果;

indel数据需要特殊处理一下转换成类似SNP的数据才可以合并分析;

请先 登录 后评论
  • 1 关注
  • 0 收藏,556 浏览
  • 提出于 2024-09-03 17:57

相似问题