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
=
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
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")
3 matches
Mail list logo