This is a test version of Biostars. For the public version, visit https://www.biostars.org.
compare two GWAS results

I have performed two GWAS with the same trait but different populations but I don't know how to compare the two results.

thank you!!!!

gwas snp
• 2,421 views • link
updated by DM • DM • written by sihambarkus • sihambarkus

2 answers

Create a mirrored Manhattan plot using GWAS data

https://github.com/anastasia-lucas/hudson

enter image description here

• 0 views • link
written by Shicheng Guo Shicheng Guo

Hi,

You can use Stephen's qqman package to create Miami plots to compare your results. 

library(qqman)
par(mfrow=c(2,1))
par(mar=c(0,5,3,3))
manhattan(gwasResults,ylim=c(0,10),cex=2.2,cex.lab=2.5,font.lab=2,font.axis=2,cex.axis=1.6,las=2,font=4)
par(mar=c(5,5,3,3))
manhattan(gwasResults,ylim=c(10,0),cex=2.2,cex.lab=2.5,font.lab=2,font.axis=2,cex.axis=1.6,las=2,font=4,xlab="",xaxt="n")
dev.off()
• 0 views • link
written by DM • DM

Log in to answer this question.

More posts like this