This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Choose between pvalue & pvalue.adj

Hi everyone,

I am currently carrying out a DE analysis on single cell, I have just arrived at the volcano plot and I do not know what to use for the y axis: the pval or the pval.adj.

From what I understood, the pval.adj are calculated because of the multiple test problem, but I would like to know whether to use the pval.adj every time or only if the histogram of the pval is not good (or other reason).

Thank you in advance for your help.

volcanoplot
• 752 views • link
updated by Istvan Albert Istvan Albert • written by Guitout • Guitout

1 answer

You should use the correct p-value depending on what you are trying to test for.

If the effect size you are testing for can be affected by multiple comparisons then you should use .adj. If the effect size that you are testing for is not affected by multiple comparisons then you should use the unadjusted p-value.

• 0 views • link
written by Istvan Albert Istvan Albert

Log in to answer this question.

More posts like this
  • not centered volcano plot, Seurat
    written by sooni • sooni
    Hello. I performed DE analysis on single cell data using Seurat, following the pseudobulk tutorial on [this page][1]. After DE analysis, I drew a volcano …
  • Histogram of normalized counts for each cell type
    written by odi • odi
    So, I am visualizing wanting to make a histogram in R using my normalized counts (from SCTransform). Specifically, i want the y-axis to be normalized …
  • Adjust Y axis scale in ggplot2
    written by santiagoboccardo31 • santiagoboccardo31
    Hi, I am trying to create a volcano plot from RNAseq data. The problem is that the adjusted p-values (plotted in Y-axis) have a very …
  • Bad volcano plot
    written by firestar firestar
    I have run differential gene expression on a 10X single-cell dataset in Seurat. See code below. m <- FindMarkers(obj_dge,ident.1=test,ident.2=ref,group.by="cell_type_condition_cluster", test.use="MAST",only.pos=FALSE,min.pct=0.25, logfc.threshold=0,max.cells.per.ident=200, min.cells.group=10,random.seed=100, assay="RNA",slot="data", latent.vars=c("nCount_RNA","percent_mito","S.Score","G2M.Score")) I …
  • why my volcano plot looks strange?
    written by hellocita • hellocita
    Hi all, I am dealing with an LC-MS Metabolomics data. I perform log transform, quantile normalization and use linear model to regress out batch. And …
  • volcano plot help code
    written by priyadhawan11234 • priyadhawan11234
    I want to make a volcano plot with 4 different cell types and these cell types should have different shape and color. But i want …
  • Plotting .tab/.bam file
    written by zack.henning • zack.henning
    Hey all, I am trying to find a way to create a scatter plot and histogram using matplotlib for an alignment I generated. I aligned …
  • Histogram plot of padj values
    written by cado • cado
    Hi, I'm carrying out DE analysis on RNASeq data in DESeq2. I have a question about the histogram plots I used to plot the p-adj …
  • How to interpret and trim this plots?
    written by eyonesi • eyonesi
    hello everyone I am running denovo rnaseq experiment at quality control step. I can’t interpret and trim two plot duplication level and GC content at …
  • Associating Gc Content With Genomic Features
    written by Gvj Gvj
    <p>Hey all,</p> <p>When I have identified certain events in genome and would like to check whether it is associated with GC content, histone mark etc, …