volcano plot display selected labels
1
0
Entering edit mode
7 weeks ago
edus_bioinfo ▴ 20

Hi all, With the code below I try to draw volcano plot but there is a small problem. The labels that I selected is not shown in the graph.

volc <- EnhancedVolcano(dgeList_cpm_filtered_gtf,
                        lab = dgeList_cpm_filtered_gtf$gene_name,
                        x = 'logFC',
                        y = 'PValue',
                        xlab = bquote(~Log[2]~ 'FC'),
                        selectLab = c('VCAM1','NCAM2','FGFR1',
                                      'NLGN1','FGF5','UGT2B7','AKR1C2','SPP1','NR1H4','ALDH1A1'),
                        pCutoff = 0.05,
                        FCcutoff = 2,
                        pointSize = 1,
                        labSize = 4.0,
                        xlim = c(-20, 20),

                        axisLabSize = 18,
                        colAlpha = 1,
                        legendLabels = c("NS", expression(log[2] ~ FC), "p-value", expression(p - value ~ and~ log[2] ~ FC)),
                        legendPosition = 'right',
                        legendLabSize = 12,
                        legendIconSize = 4.0,
                        drawConnectors = FALSE,
                        widthConnectors = 0.75,
                        col=c('gray59', 'springgreen3', 'deepskyblue1', 'tomato1'),
                        title = "",
                        subtitle = "")

enter image description here

How can I show all 10 selected points in the graph? Thanks in advance

volcanoplot enhancedvolcano rna-seq • 163 views
ADD COMMENT
1
Entering edit mode
7 weeks ago

Hi, please can you do 2 things to resolve this:

  • verify that all of the desired genes are in dgeList_cpm_filtered_gtf$gene_name

  • use drawConnectors = TRUE

Kevin



Olá, por favor pode fazer 2 coisas para resolver isto:

  • verificar que todos os genes ficam parte de dgeList_cpm_filtered_gtf$gene_name
  • use drawConnectors = TRUE

Kevin

ADD COMMENT

Login before adding your answer.

Traffic: 1812 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6