微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

曼哈顿图未标记所有数据点

如何解决曼哈顿图未标记所有数据点

我在下面有一个名为 test 的数据。我正在尝试使用 manhattan 库中的 qqman 函数绘制此数据,但如果您在图中注意到,它并未标记 Pvalue < -log10(1e-08) 以下的所有点。我在这里错过了什么?

library(qqman)    
manhattan(test,main = "",ylim=c(0,15),col = c("blue4","orange3"),suggestiveline = -log10(1e-06),genomewideline = -log10(1e-08),annotatetop = TRUE,annotatePval = -log10(1e-08),chrlabs = as.character(sort(unique(test$CHR))))

数据:

    test <- structure(list(SNP = c("CSNK1G2","MEIS2","RP11-128A17.1","SPPL2B","ZBTB10","SDK1","1:29115501:C:G","DCBLD2","MINPP1","RN7SL411P","CSRP2BP","POFUT1","SAMHD1","ARFGEF2","PCSK1","NRCAM","FTSJ3","MTSS1","VNN2","C22orf23","TRIP11","RP11-267M23.3","RNF17","LAD1","GCAT","DPP6","MYRIP","COL14A1","THUMPD3","5:90519607:G:C","LMF2","2:130182303:G:C","RP11-378A13.1","XRCC2","nop9","PTPRR","IKZF4","SND1","PER3","REXO1","ZNF554","1:37828266:GA:G","COQ5","DAGLB","3:113869718:C:T","UBQLN1","GCN1L1","3:152340425:C:A","RP11-96C23.10","SH3BP2","ATRNL1","MAP3K9","CEP170B","SOAT1","RGSL1","RPL5","FOLH1B","AC068799.1","MSI1","TRIAP1","NCOR2","POLE","TNFAIP2","COCH","ZFYVE26","ELMSAN1","KIF7","PTX4","LA16c-OS12.2","FLJ45079","CCDC178","MUM1","AP3D1","RN7SL202P","POLQ","HLTF","AHSG","ITPR1","SHQ1","Metazoa_SRP","MYL5","ANXA6","PKHD1L1","TNKS","RP11-195F19.29","LA16c-390E6.3","PKD1","CRYBG3","ARID3C","CGN","G2E3","KIF26A","SI","VPS13A","OLFM3"),P = c(1.459e-12,3.993e-12,2.925e-11,2.932e-11,1.599e-10,5.258e-10,8.619e-10,9.423e-10,1.377e-09,1.415e-09,1.734e-09,1.759e-09,1.98e-09,2.554e-09,3.798e-09,6.174e-09,6.569e-09,1.55e-08,1.899e-08,2.218e-08,2.485e-08,2.776e-08,2.796e-08,3.455e-08,3.904e-08,5.134e-08,5.333e-08,5.357e-08,5.687e-08,6.479e-08,6.9e-08,6.901e-08,7.179e-08,7.28e-08,7.375e-08,7.408e-08,7.733e-08,7.814e-08,8.17e-08,1.019e-07,1.061e-07,1.068e-07,1.172e-07,1.189e-07,1.205e-07,1.211e-07,1.217e-07,1.218e-07,1.22e-07,1.239e-07,1.243e-07,1.245e-07,1.247e-07,1.249e-07,1.25e-07,1.252e-07,1.254e-07,1.259e-07),CHR = c(19,15,19,8,7,1,3,10,20,5,17,6,22,14,13,2,12,9,4,11,16,18,1),BP = c(1953357L,37037001L,37094404L,2334772L,2351704L,80486958L,4077071L,29115501L,98811473L,98811485L,87552289L,40210113L,18162024L,32230868L,36930902L,48989666L,96399107L,108234727L,63824712L,124704342L,80487332L,132752654L,37956695L,91993921L,94534777L,24874294L,201386508L,37816581L,154807181L,40209992L,120158188L,9384525L,90519607L,50505437L,130182303L,218253844L,152660792L,24306548L,70764560L,56002806L,128087016L,7819427L,1825726L,2834416L,37828266L,120509934L,6447649L,113869718L,83713065L,120162755L,152340425L,87008531L,2831832L,115215944L,70809013L,104884202L,179302620L,182474354L,182488987L,92832097L,89670555L,117277892L,120357950L,120444886L,124332440L,132639181L,103126688L,30885788L,67804167L,73738854L,89633246L,1486462L,190308L,77881996L,33266975L,1356531L,2115209L,3918955L,121490123L,149071684L,186620833L,4683555L,72848177L,72922141L,680597L,151122844L,109450979L,120280088L,9679900L,34655657L,1465246L,2108561L,97893057L,34625661L,151520093L,30597650L,104178576L,165063524L,77405807L,101846861L)),row.names = c(135977L,89108L,89111L,136668L,136719L,266871L,252346L,20439L,206542L,206543L,39549L,201043L,173375L,174747L,176078L,177956L,230107L,244948L,121300L,259465L,266873L,232530L,187699L,87423L,267501L,75719L,12894L,187625L,249450L,201042L,259059L,206062L,229682L,190471L,152953L,160353L,249361L,81933L,71724L,69236L,245983L,27397L,135369L,136925L,22276L,61367L,254396L,192020L,278877L,61243L,196181L,39522L,212204L,31288L,85338L,79906L,10456L,10807L,10810L,28627L,57504L,60879L,61339L,61351L,62524L,63596L,79508L,82184L,85041L,85741L,96907L,98641L,99118L,124968L,128649L,133175L,136343L,139848L,192458L,195816L,198417L,201965L,205639L,205643L,214737L,222732L,258461L,259099L,267666L,276638L,98628L,99960L,206314L,276618L,4727L,82139L,79769L,196857L,278603L,192L
),class = "data.frame")

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。