Open swvanderlaan opened 6 years ago
library(hudson)
two files with 4 columns : SNP, CHR, POS, pvalue
gmirror(top=female_hudson, bottom=male_hudson, tline=0.05/nrow(female_hudson), bline=0.05/nrow(male_hudson), toptitle = "GWAS Comparison: Female samples", bottomtitle = "GWAS Comparison: Male samples", highlight_p = c(0.05/nrow(female_hudson),0.05/nrow(male_hudson)), highlighter="dodgerblue")
Add in option to choose for a stratified Manhattan (bottom vs. upper for instance male vs. female), i.e. a MIAMI-plot. There is a specific package that has this function --> borrow it from there... 🔶