compare two GWAS results

Question

[sihambarkus]

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!!!!

Answer 1

[Shicheng Guo]

Create a mirrored Manhattan plot using GWAS data

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

Answer 2

[DM]

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()