老师您好,针对这个问题https://www.omicsclass.com/question/7725。你们提到:你这个数据里面是不是P值有0的情况,取了对数就无穷大了,导致出问题。但是,我发现里面没有P值为0的情况,那是不是脚本gwas_manhattan_plot.r有问题呢?需要修改gwas_manhattan_plot.r脚本的什么地方才可以解决顶部只能画半圆的问题呢?还有,我发现P-Value的原始文件有NaN的情况,是什么原因呢?是否对GWAS有影响呢?

此图是问题https://www.omicsclass.com/question/7725的GWAS图P-value的原始输入文件,红圈是最大值(但是对这个值进行曼哈顿图可视化后只有半个圆形):

可视化的曼哈顿图:attachments-2025-03-AfBMOmXn67cfadb0c8a1f.png

请先 登录 后评论

1 个回答

“发现P-Value的原始文件有NaN的情况”,应该是你有些SNP位点缺失率过高之类的情况,大概是你前期的数据过滤不是很严格,你画图的时候直接删掉这样的数据就可以

关于你提出的曼哈顿图上面的点显示不全的问题,你可以对脚本进行如下的修改:

scale_y_continuous(expand = c(0,0.05)

即把脚本中的y轴坐标轴上限增加5%

请先 登录 后评论
  • 1 关注
  • 0 收藏,147 浏览
  • John 提出于 6天前

相似问题