Re: [R] text annotation on Manhattn plot in qqman

2020-05-20 Thread Michael Dewey
Dear Ana That looks like something is hard coded in manhattan(). The simplest thing might be to contact the maintainer of the package and ask. You can make a copy of manhattan or textxy and modify them but I think the maintainer is the simplest course of action. Michael On 20/05/2020 16:10,

Re: [R] text annotation on Manhattn plot in qqman

2020-05-20 Thread Ana Marija
HI Michael, Thank you so much! That worked!!! Now I am just trying to increase the size of text of SNP and GENE on plot I tried this: a$newname <- paste(a$SNP,"\n", a$GENE) manhattan(a, chr="CHR", bp="BP", snp="newname", p="P",cex = 0.5,annotatePval = 0.0001) but I am getting this error: Error

Re: [R] text annotation on Manhattn plot in qqman

2020-05-20 Thread Michael Dewey
a$newname <- paste(a$SNP, a$GENE) manhattan(a, chr="CHR", bp="BP", snp="newname", p="P",annotatePval = 0.0001) However note that I do not use manhattan() so you may need to alter the parameters as I am assuming a is where it finds the remaining parameters. You may also need to play with the se

Re: [R] text annotation on Manhattn plot in qqman

2020-05-19 Thread Michael Dewey
Dear Ana Perhaps paste together SNP and GENE using paste() and then supply that as the snp parameter. Michael On 19/05/2020 17:12, Ana Marija wrote: Hello, I am making manhattan plot with: library(qqman) manhattan(a, chr="CHR", bp="BP", snp="SNP", p="P",annotatePval = 0.0001) and I would l

Re: [R] text annotation on Manhattn plot in qqman

2020-05-19 Thread Ana Marija
Hi Michael, can you please send me code how that would be done? Thanks Ana On Tue, May 19, 2020 at 11:18 AM Michael Dewey wrote: > > Dear Ana > > Perhaps paste together SNP and GENE using paste() and then supply that > as the snp parameter. > > Michael > > On 19/05/2020 17:12, Ana Marija wrote:

[R] text annotation on Manhattn plot in qqman

2020-05-19 Thread Ana Marija
Hello, I am making manhattan plot with: library(qqman) manhattan(a, chr="CHR", bp="BP", snp="SNP", p="P",annotatePval = 0.0001) and I would like to annotate these two SNPs which are above the threshold so that they have GENE name beside them: > a[a$SNP=="rs4081570",] SNPP CHR