Questions tagged [deseq2]

Deseq2 is an R package for analyzing RNA sequencing data

Filter by
Sorted by
Tagged with
2 votes
0 answers
43 views

How to incorporate negative controls in DESeq2

We are doing a comparison between two outcomes (positive and negative). We could not have any positive controls as we do not have any "control" data to set as baseline, either from ...
Karthik Nair's user avatar
0 votes
0 answers
38 views

how is log2foldchange calculated in DESeq2

I am trying to calculate log2FoldChange in DESeq2 by hand because I am trying to reproduce my results and understand how DESeq2 works. What I want to know is: Are the counts first fitted to a ...
Guri Giaever's user avatar
1 vote
0 answers
36 views

enrichment of DE genes within wide regions not pathways, hypergeometric test or GSEA?

I haven't found any relevant posts, so here is my question (it's more a question on statistics I believe): I have the deseq2 output of DE genes( around 25000 genes) , and among them 1000 are ...
dew's user avatar
  • 11
2 votes
2 answers
67 views

Understanding DESeq2 output for resultNames

I'm following along with some pre written code to learn how to do differential expression. I am using the same data set and have gone back through the code and copy and pasted what they used and the ...
Sabrina's user avatar
  • 73
3 votes
1 answer
65 views

Issue with Differential Gene expression analysis with Deseq2

I am using the Deseq2 package to perform a DE analysis on multiple factors. So, basically to perform the analysis I consider patients characterized by the Brugada Syndrome type 1 and other patients ...
Alessia Avon's user avatar
2 votes
1 answer
183 views

Correlation heatmap of RNA-seq clusters all samples together leading to very low no. of DEGs

I am writing to you to take an input or may be you can provide a different perspective. I am at wits end :( So I have nearly 200 samples. These are separated into two groups (group I treated with ...
Angelo's user avatar
  • 227
2 votes
2 answers
210 views

Integrate bulk RNA and ATAC-seq genes

I have up and down-regulated genes from bulk RNA result from DeSeq2. Genes that are differential in chromatin accessibility between two conditions from DiffBind. I would like to find the genes overlap ...
Chris's user avatar
  • 101
1 vote
1 answer
126 views

DESeq2 throwing error while normalizing raw microarray expression data due to presence of negative values

This question was also asked on Biostars I am trying to download and analyze a miRNA expression dataset from NCBI GEO (GSE25631). I specifically want non-normalized ...
accibio's user avatar
  • 145
2 votes
2 answers
51 views

RNAseq alignment: best practices for aligning to multiple isoforms?

I have Illumina RNAseq data and would like to maximize my power to find candidate genes that are differentially expressed genes between experimental conditions. Many of my (de novo assembled and ...
DavidR's user avatar
  • 23
1 vote
1 answer
342 views

Differential gene expression: how do I correctly compare three groups to each other?

I'm trying to figure out the right way to do differential gene expression analysis with a discrete variable with three groups. The context is, I want to assess differential expression as a function of ...
mrz123456's user avatar
0 votes
1 answer
82 views

Batch correction in differential expression analysis

I have sent two sets (two batches in matter of sending for sequencing) of different samples (plasma) to small RNA-seq to Qiagen company This is how my meta data look ...
Angel's user avatar
  • 1,981
1 vote
0 answers
787 views

Manually calculating log2 fold change values from DESeq2 normalized counts [closed]

This question was also asked on Biostars and Bioconductor Support My aim is to perform gene clustering based on RNA seq datasets (raw expression values) from NCBI database. After clustering I realized ...
jayesh kumar's user avatar
3 votes
2 answers
350 views

DESeq2 for large number of samples takes too much RAM

I am trying to run a very large number of transposase-accessible chromatin (ATAC)-seq samples through DESeq2 to find peaks of differential chromatin accessible across my study genome. I have about 100 ...
Shashank Nagaraja's user avatar
1 vote
1 answer
81 views

DESeq2 and EdgeR

I am new to using both DESeq2 and EdgeR in Bioconductor used for transforming my RNA expression data. However, I am struggling to understand their specific purpose, differences between them and ...
noor fatimah's user avatar
1 vote
1 answer
77 views

Filtering reads greater than 5 from HT-seq count files

I have some raw counts from HTSeq after aligning with hg38 human reference genome. I want to do filtering in a way that the filtered count files should have the same number of lines. The reason behind ...
Aranyak Goswami's user avatar
1 vote
0 answers
41 views

Creating multiple phenotype datasets using bootstrap method "Bootstrap-samples-by-column-of-a-data-frame-in-r" for DEG analysis

I am working on a datasets and after some discussion with my group, we doubt that maybe one or more of our controls are different than the other controls. The motivation is to see if one or more ...
salman1490's user avatar
2 votes
4 answers
219 views

nested samples in design DEseq2

I'm running a DEG analysis with Deseq2 by specifing the following design that includes 12 samples from 6 pts in 4 different conditions: ...
GuidoL's user avatar
  • 21
1 vote
2 answers
105 views

Basic RNA Differential Expression in R

I have two matrices, one for individuals before treatment and one for the same individuals after treatment. Both matrices are raw read counts of RNA expression. ...
aurelius_37770's user avatar
1 vote
1 answer
71 views

deseq2 single factor design output interpretation

This is regarding the single factor design For example if I have Age or other continuous numerical variable how to provide that into the design formula. For this post do i need to 'You could ...
kcm's user avatar
  • 1,774
1 vote
1 answer
98 views

deseq2 model design : Different gene output

Here I have two designs: Design 1 ...
kcm's user avatar
  • 1,774
3 votes
1 answer
191 views

deseq2 makeContrasts function

This article talked about various deseq2 designs etc. One of the designs I would like to use is explained as this: Control versus treatment average ...
kcm's user avatar
  • 1,774
1 vote
1 answer
528 views

deseq2 full and reduced model interpretation

I would like to know how to interpret the output of the formula although i been thorough loads of literature but I'm not confident yet. So this is my full model I'm running ...
kcm's user avatar
  • 1,774
0 votes
0 answers
68 views

Different results of DESeq2 in different ways of comparison

I have different treatments (infected and control ones). When using DESeq2, I have different results in two ways of comparison. first way: factor: treatment, factor level1: Moribund Infected, factor ...
Zohreh's user avatar
  • 11
2 votes
1 answer
417 views

What is the difference between Normalized Expression in EdgeR vs DESeq2?

I am trying to access the normalized expression in both edgeR and DESeq2, yet the results are different. Does anyone know why? How to get normalized expression using edgeR: ...
Nova's user avatar
  • 165
2 votes
1 answer
274 views

baseMean threshold

I have an RNA-seq dataset and I am using DEseq2 to find differentially expressed genes between the two groups. I used pre-filtering to remove any genes that have ...
suzie321's user avatar
1 vote
1 answer
125 views

Error in DESeq2 analysis with dimnames

I'm getting the following error in a DESeq2 analysis and I can't seem to figure out the issue: ...
Nova's user avatar
  • 165
3 votes
0 answers
65 views

LRT test in deseq2 on multiple condition

I have some doubts when i'm running LRT This is my condition data-frame small subset ...
kcm's user avatar
  • 1,774
1 vote
2 answers
366 views

How to input data and metadata from NCBI for RNA-Seq analysis in R

So, I have downloaded some .txt and .csv files containing the raw count data (I am supposed to find out differentially expressed genes and see which genes and up- and down-regulated). I have also ...
Addy's user avatar
  • 11
0 votes
1 answer
555 views

How does DESeq2 "collapseReplicates()" function work on read counts data?

Comparing read counts from an RNA-seq experiment for two select genes before and after using DESeq2's collapseReplicates() and ...
Gawain's user avatar
  • 115
1 vote
0 answers
42 views

Batch Effects that are Close to Collinear with Main Effect

Two years ago, I started working on a project uses both RNAseq and ATACseq. It's supposed to be a simple Healthy Control (HC) vs disease study. The sample collection was done in 2 phases, and it was ...
Jeff's user avatar
  • 177
0 votes
1 answer
234 views

How to remedy a DESeq2 collapsing technical replicates error?

Goal: To ensure "the sum of the counts for [my samples] is the same as the counts in the [samples] columns in ddsColl" after collapsing replicates using ...
Gawain's user avatar
  • 115
0 votes
1 answer
1k views

Use of "baseMean" in DeSEQ2

The basemean is described as the "mean of normalized counts of all samples". My question is, how would we interpret differences in basemean between genes, pertaining to reliability and such? ...
CodeMonke's user avatar
-1 votes
1 answer
59 views

Visualizing top expressed genes

I have a spreadsheet of CPM values I want to visualize the top 20 genes expressed across my samples This is how my data look like ...
Angel's user avatar
  • 1,981
1 vote
1 answer
129 views

Best R program to use for RNA seq analysis

We have an RNA seq taken from a the same tissue in two different developmental stages. We have one gene, with the number of reads from stage A and stage B, and we want to see which genes are ...
Merav Gold's user avatar
1 vote
0 answers
43 views

The correct design for time series experiment

I have a treatment and control in two time points like this ...
Angel's user avatar
  • 1,981
1 vote
0 answers
37 views

How to simulate replicates for DGE analysis?

I am prototyping with data visualization of DGE results, and I would like to work on the analysis pipeline before the real data is available. Currently, I only have 3 samples for wild type and 1 ...
pietro_molina's user avatar
0 votes
2 answers
297 views

Help with performing PCA analysis on data from 3 different datasets of normalized counts data for RNA-seq experiment

So, I have limited knowledge of R but I need to do a PCA analysis of 3 different datasets of gene expression as a result of combined growth or mono-culture growth. The 3 different datasets I performed ...
Justin1609's user avatar
0 votes
1 answer
96 views

Understanding the different designs in DESeq2

I am using DEseq2 and trying to understand the results obtained using different models. I have a data design with 2 genotypes and 2 time points. ...
yathrakaaran's user avatar
1 vote
1 answer
218 views

Different contrasts in DESeq2

I have a treatment and control in two time points like this ...
Angel's user avatar
  • 1,981
1 vote
2 answers
291 views

DEseq2 contrast in two time points

I need help with some clarifications please I have control versus treatment in two time points like ...
Angel's user avatar
  • 1,981
0 votes
0 answers
186 views

Differential miRNA expression using RPM

I have microRNA (miRNA) expression data in RPM. I would like to do differential gene expression between two groups. How can I do this? I have considered edgeR and DESeq2 in R, but it looks like they ...
Sylvia Rodriguez's user avatar
-1 votes
1 answer
2k views

Adjusted P-value from DESeq2

What does the padj really mean? I believe it's a BH method, where it attempts to quantify the likelihood of a "false positive" but I am not sure. Also, what does it mean if the resulting ...
CodeMonke's user avatar
1 vote
1 answer
84 views

DESEQ2 Question about results()

I am currently learning to perform Differential Analysis via DESEQ2 R Package, and I believe I've made progress, able to format the data correctly [maybe] for DDS(). When I run the results function to ...
CodeMonke's user avatar
0 votes
2 answers
591 views

Beginner with DESeq2 having issue with analysis

I am trying to run DESeq2 analysis on my featureCounts output counts table. I ran the following code: ...
Justin1609's user avatar
1 vote
1 answer
365 views

DESeq2. Which is better: using the altHypothesis argument in the results function or manually filtering for your desired P-value and Log2FoldChange?

In DESeq2 you can identify your differentially expressed (DE) genes using the results function. I noticed people in my lab using the results() function with the minimum number of arguments supplied (...
Sam's user avatar
  • 13
2 votes
1 answer
166 views

How does FDRtool work?

I have a question about using FDRtool. In the below code (on RNA seq data whose p values were acquired using Deseq2), the FDRtool was first used and thereafter p.adjust using the benjamini hochberg ...
Emma's user avatar
  • 31
3 votes
2 answers
294 views

Calculating most abundant transcript from RNA-Seq data

vcf2maf uses VEP to annotate variants, and I believe selects the default Ensembl transcript to use for annotation. Sometimes the ...
Tomas Bencomo's user avatar
2 votes
2 answers
862 views

DESeq2 multiple treatments, multiple time points, multiple cell lines

Note: this question has also been asked on Bioconductor Support I know there are a lot of questions asking similar things here on this forum and I checked the vingette and did a lot of other research ...
tnocs's user avatar
  • 23
1 vote
2 answers
571 views

Un-normalize DESeq2 counts

I obtained a matrix of RNA-seq count data that has been normalized by DESeq2's median of ratio method. I know that DESeq2 wants to take in un-normalized counts, but I do not have access to those data. ...
geom_na's user avatar
  • 237
1 vote
2 answers
81 views

Bulk RNA-Seq Read Length Normalization across different samples

I have 20 samples out of which 14 are 100 bp in length and 6 are 150 bp. Is there a way to normalize the read length for cross-sample differential expression comparison? What would be the best way to ...
so_close_yet's user avatar