Entering edit mode
4.5 years ago
Kumar
▴
170
I need to add p < 0.01 cut off value in my manhattan plot. Here is the code. Please suggest, how I need to calculate and add for drawing manhattan plot.
CMplot(manhat,plot.type="c",r=0.6,
col=c("grey30","grey60"),chr.labels=c("1","2","3","4","5","6","7","8","9","10","11","12","13","14"),
LOG10=TRUE, cir.chr.h=0.5,amplify=TRUE,threshold.lty=c(1,2),threshold.col=c("red",
"blue"),signal.line=1,signal.col=c("red","green"),chr.den.col=c("darkgreen","yellow","red"),
outward=FALSE,file="jpg",dpi=600,file.output=TRUE,verbose=TRUE,width=10,height=10)
INPUT File:
SNP CHR BP P
bio1 1077947 1.00E-05
bio2 849080 5.08E-05
bio3 387297 9.48E-04
bio4 643714 9.48E-04