testmanhattan

library(qqman)

## 
## For example usage please run: vignette('qqman')
## 
## Citation appreciated but not required:
## Turner, S.D. qqman: an R package for visualizing GWAS results using Q-Q and manhattan plots. biorXiv DOI: 10.1101/005165 (2014).

d <- read.csv("~/Downloads/try1.csv", row.names=1)
head(d)

##   SNP CHR BP      P
## 1 rs1   1  1  0.842
## 2 rs2   1  2  1.262
## 3 rs3   1  3 10.000
## 4 rs4   1  4  0.001
## 5 rs5   1  5 10.000
## 6 rs6   1  6  0.001

manhattan(d)

manhattan(subset(d, CHR==3))