I just realized I'm only seeing this in my smallRNA libraries (adapter removed, t/rRNA removed, size-selected for 20-25nt in sRNA Workbench).
I map the reads using bowtie (-v 0).
I generate read counts with htseq-count, then build my count table(s).
I run DESeq following along with the vignette section 3.1.
So far, every 1v1 count table I've looked (7 of the possible 15) at has called more significantly changed loci (padj < 0.05) than I get when looking at the exact same comparison using a count table that includes all 30 of my bio-reps.
The biggest 'jump' was from 76 to 372 loci for one comparison.
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
Originally posted by rndouglas View PostThis is what I had been doing all along (running everything in one count table), but then this morning I thought I'd try A vs. B in a separate count table. I was surprised to find almost double the loci with padj < 0.05 compared to when I ran everything in one count table (and hence my new-found concern).
Leave a comment:
-
Have a look at the size factors. If one of them from the full dataset is very different than the others, that can cause this sort of result.
Leave a comment:
-
This is what I had been doing all along (running everything in one count table), but then this morning I thought I'd try A vs. B in a separate count table.
I was surprised to find almost double the loci with padj < 0.05 compared to when I ran everything in one count table (and hence my new-found concern).
Leave a comment:
-
It is better to have a single count table. This shall lead to increased statistical power. You are likely to get somewhat higher number of differentially expressed genes as DESeq would be able to tease out more signal from the noise.
Leave a comment:
-
What to include in my count table(s) for DESeq
I have five experimental groups (A, B, C, D, and E; 3-6 bioreps each) that I want to compare to each other and a control (Y; 3 bioreps) using DESeq.
When I run DESeq is it better to use a count table that includes A, B, C, D, E, and Y, then run each comparison, or is it better to make a count table for each comparison (ie: a table for A and Y, another table for B and Y, another table for A and B, etc.) and end up with ~15 different count tables?
I ask because I end up with different lists of significantly changed loci depending on how I run the analysis and I'm not sure which is more 'correct.'
Latest Articles
Collapse
-
by seqadmin
The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...-
Channel: Articles
04-22-2024, 07:01 AM -
-
by seqadmin
Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...-
Channel: Articles
04-04-2024, 04:25 PM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, Yesterday, 11:49 AM
|
0 responses
13 views
0 likes
|
Last Post
by seqadmin
Yesterday, 11:49 AM
|
||
Started by seqadmin, 04-24-2024, 08:47 AM
|
0 responses
16 views
0 likes
|
Last Post
by seqadmin
04-24-2024, 08:47 AM
|
||
Started by seqadmin, 04-11-2024, 12:08 PM
|
0 responses
61 views
0 likes
|
Last Post
by seqadmin
04-11-2024, 12:08 PM
|
||
Started by seqadmin, 04-10-2024, 10:19 PM
|
0 responses
60 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 10:19 PM
|
Leave a comment: