Tag: log2FoldChange

Strangely too low P-value and Adjusted P-value(FDR) DESeq2 and edgeR

My data is experimental data that has been overexpressed for a specific gene. Data samples are divided into 3 groups according to the over-expression time and each group has 3 samples. (total 9 samples) I conducted DGE analysis on the control group and one case group with DESeq2. cts <-…

Continue Reading Strangely too low P-value and Adjusted P-value(FDR) DESeq2 and edgeR

Change log2FoldChange range – plotMA

You can use base R graphics to make these plots. The data is sitting there in columns of the res object, so you can filter it directly, and use boolean vectors to pick out the things you need: # make sure there are no NA values sum(is.na(res$log2FoldChange)) # choose some…

Continue Reading Change log2FoldChange range – plotMA

sigh of log2FC values in DESeq2

sigh of log2FC values in DESeq2 0 @194b0276 Last seen 10 hours ago United States Can someone explain to me how is the sign of log2FoldChange is set in the results of DEseq2? I was pretty sure that it is calculated as log2(Counts_treatment/Count_reference), where reference is determined alphabetically (unless specified…

Continue Reading sigh of log2FC values in DESeq2

Understanding the output of Negative Binomial in DESeq2

Are zero is a default option of log2fold change to be considered as up and down? Yes, by default, the null hypothesis is that the log2FoldChange is zero. Does this mean that there are seven genes in total that will have significant adjusted p-value? Yes, you are interpreting your results…

Continue Reading Understanding the output of Negative Binomial in DESeq2

DESeq2 log2FC and p-value filtering

DESeq2 log2FC and p-value filtering 1 Dear all, I want to filter DE genes by log2FC > 2/-2 and p-value 0.0001. Why this two versions of lof2FC and p-value filtering differ in up- and down-regulated genes numbers (from 800 to 1500) and which one is the right one ? 1….

Continue Reading DESeq2 log2FC and p-value filtering

Correct way to make multiple comparisons on DESeq2?

I have a project where I have done RNA-seq (paired-end sequencing on Illumina HiSeq) of a worm at different days of development i.e. Ages 0-12. For each age, I have sequenced 3 replicate specimens. I’m new to DESeq2 and I was wondering if what I did below is correct. library(DESeq2)…

Continue Reading Correct way to make multiple comparisons on DESeq2?

Outliers on DESEq2 Results

I have an RNAseq dataset, where one of the genes I intend to analyze has hundreds of counts ranging from 10 to 12, with a few counts > 9000. I process this data in Deseq2 and get that the gene is differentially expressed across several samples of interest. What can…

Continue Reading Outliers on DESEq2 Results

Need help to remove NA values from data frame

Need help to remove NA values from data frame 2 I have this data frame : and I want to remove those rows which contain NA values from the log2fold change column How can I do this through R? DeSEQ2 R • 256 views Hi Anas, If your data frame…

Continue Reading Need help to remove NA values from data frame

Volcano plot not centered at origin

I am making a volcano plot for my phosphoproteomics data and seemingly the data are not centered around the origin. The reason that I know there is something wrong is because before my code was producing the correct volcano plot, and now it is shifted. What are some reasons why…

Continue Reading Volcano plot not centered at origin

Handle inflated log2FC while using interaction term in DESeq2

Hi guys, I’m working with a 8 samples experiment (lower than 3x replicates, I know..) with a design like > colData(dds) DataFrame with 8 rows and 2 columns condition traitment <factor> <factor> 4200-JS-1 norm ctrl 4200-JS-2 norm ctrl 4200-JS-3 norm trt 4200-JS-4 norm trt 4200-JS-5 hyper ctrl 4200-JS-6 hyper ctrl…

Continue Reading Handle inflated log2FC while using interaction term in DESeq2

Is it a good idea to plot a volcano plot starting at position x =1

Is it a good idea to plot a volcano plot starting at position x =1 1 Hello everyone! I was wandering if it makes sense to start a volcano plot (x=log2FoldChange) starting at position one, so basically discarding all values below 1 and above -1 leaving only values outside this…

Continue Reading Is it a good idea to plot a volcano plot starting at position x =1

Circos plot with logfold change RNA seq data

I am new to circos plot analysis and have been trying to use the cyclize package. I want to display mRNA differential gene expression data based on data analyses of 8 libraries and links between their respective target genes. The dataset I am working with looks like this geneid baseMean…

Continue Reading Circos plot with logfold change RNA seq data

lfcShrink probelm in many 0 count genes RNA-seq data

Hi, Dr love. I post a question about weird MAplot or volcano plot of DESeq2 diff result and also in biostar. ATpoint give a useful answer about too many 0 count genes and prefiltering. It seems that too many 0 count genes makes lfc shrink have a probelm. And I…

Continue Reading lfcShrink probelm in many 0 count genes RNA-seq data

weird MAplot or volcano plot of DESeq2 diff result

Hi, every one. I find a werid MAplot or volcano plot of DESeq reuslt. I am wondering whether you can give me some advice. This diff result is from two cell type bulk RNA-seq. I use two specific marker to get these two cell type using Flow cytometer. I alreadly…

Continue Reading weird MAplot or volcano plot of DESeq2 diff result

same padj for all the genes after DEseq analysis

Hi everyone, I have done a pair comparison with DEseq2 to find differentially expressed genes between two samples with 6 replicates, for the DEseq2 result, i got exactly same padj value for all the genes and it is not significant, is this normal ? I don’t think the padj should…

Continue Reading same padj for all the genes after DEseq analysis