Volcano plot
0
0
Entering edit mode
2.8 years ago
fahim ▴ 20

I want to draw a volcano plot of my DGE.So according to my data analysis in R studio, I found 15521 DGE.So from this I sort top sign DGE by giving adj.p.value <0.01.Where I found 185 DGE. But now I am confused about the drawing of volcano plot. I want to know the upregulated and downregulated genes among them. So from which count I should see the plot(15521 or 185). Moreover I want to use below code

EnhancedVolcano(res,
    lab = rownames(res),
    x = 'logFC',
    y = 'adj.p.value',
    pCutoff = ,
    FCcutoff = ,
    cutoffLineType = 'twodash',
    cutoffLineWidth = 0.8,
    pointSize = 4.0,
    labSize = 6.0,
    colAlpha = 1,
    legend Labels=c('Not sig.','Log (base 2 ) FC','p-value',
      'p-value & Log (base 2) FC'),
    legendPosition = 'right',
    legendLabSize = 16,
    legendIconSize = 5.0)

But what is the range of pcutoff and FCcutoff I should give there.And I want to clarify one more thing,that is if my x and y label was logFC & adj.p.value, so what should be the lebel legend there - Log (base 2) FC or Log (base) FC & p-value or adj.p.value?

EnhancedVolcano • 712 views
ADD COMMENT
0
Entering edit mode

I found 15521 DGE

...you may want to reconsider this analysis. This is absurdly much. How did analyse that, and which species/setup is the experiment?

ADD REPLY
0
Entering edit mode

15521 this is the overall gene expression after normalization. From them I have to draw plot by giving adj.p.value and logFC cut off. In this 15521, there is sign and non sign gene.

ADD REPLY

Login before adding your answer.

Traffic: 1240 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