Re: [R] how to change manhattan plot code to get a different color per chromosome

2020-06-16 Thread Ana Marija
Hi Jim, thank you so much for this elaborate code I will definitely use some hacks from it in the future. For now, for the purpose of compactness of the code I just set manually colors which I want to be there: manhattan(results_log,chr="CHR",bp="POS",p="META_pval",snp="MARKER",suggestiveline =

Re: [R] how to change manhattan plot code to get a different color per chromosome

2020-06-16 Thread Jim Lemon
Hi Ana, Your attached image seems to have bailed out before landing in the list. Here's how to do it using a simple manhattan plot with the data from: https://reneshbedre.github.io//assets/posts/mhat/gwas_res_sim.csv gwas<-read.csv("gwas_res_sim.csv",stringsAsFactors=FALSE) # get the data into

[R] how to change manhattan plot code to get a different color per chromosome

2020-06-15 Thread Ana Marija
Hello, Is there is a way to set colors in this plot to look like this one in attach (different color for each CHR-there is 22 of them)? library(qqman) results_log <- read.table("meta_p_pos_chr.F", head=TRUE,stringsAsFactors=FALSE) png("META.png")