YinLiLin / CMplot

📊 Circular and Rectangular Manhattan Plot
520 stars 112 forks source link

Multi_tracks plotting does not work #60

Closed shanwai1234 closed 2 years ago

shanwai1234 commented 2 years ago

Hello,

I want to draw "Multi_tracks Rectangular-Manhattan plot" using your demo data and your demo code.

SNPs <- list(
    pig60K$SNP[pig60K$trait1<1e-6],
    pig60K$SNP[pig60K$trait2<1e-6],
    pig60K$SNP[pig60K$trait3<1e-6]
)
CMplot(pig60K, plot.type="m",multracks=TRUE,threshold=c(1e-6,1e-4),threshold.lty=c(1,2), 
        threshold.lwd=c(1,1), threshold.col=c("black","grey"), amplify=TRUE,bin.size=1e6,
        chr.den.col=c("darkgreen", "yellow", "red"), signal.col=c("red","green","blue"),
        signal.cex=1, file="jpg",memo="",dpi=300,file.output=TRUE,verbose=TRUE,
        highlight=SNPs, highlight.text=SNPs, highlight.text.cex=1.4)

However, it consistently gave me error message as below:

Error in CMplot(pig60K, plot.type = "m", multracks = TRUE, threshold = c(1e-06, : No shared SNPs between Pmap and highlight!

Could you give me any help on this? Thank you so much!