Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
DEG for paired samples, biological replicates hypatia RNA Sequencing 24 05-07-2014 02:25 AM
DESeq-biological replicates differ between groups pmgr Bioinformatics 0 07-15-2013 05:52 AM
Cuffdiff with biological replicates vs individual samples not adding up. Noa Bioinformatics 0 07-02-2013 11:50 PM
Using DESeq with Biological Replicates with Different RNAseq Read Lengths jbvt Bioinformatics 3 05-21-2013 08:34 AM
DESeq without biological replicates (dataset: Marioni et al.) Andrea Apolloni Bioinformatics 3 01-13-2012 05:58 AM

Thread Tools
Old 11-05-2013, 09:51 PM   #1
Location: china

Join Date: Aug 2013
Posts: 14
Question DESeq NB test for samples without biological replicates

Hi guys,

I hv encountered problems performing negative binomial test for samples without replicates while control group has replicates with DESeq.

Experimental design:
conduct the experiment under three conditions with common control: A B C (same samples treated at 3 distinct levels -- 50%, 80%, 95%)
Condition A has two biological replicates, Condition B and C have no biological replicates, Control group has two biological replicates.
noted as A_1, A_2, B, C, Ctrl_1, Ctrl_2
Pipeline of transcriptomic RNA-seq(paired-end):
map clean reads to ref genome using Tophat -> get the raw counts using SAMtools and HTSeq -> normalized the raw counts matrix using DESeq -> negative binomial test using DESeq

# Do stats based on a negative binomial distribution from Simon's vignette
With replicates:
cds = estimateDispersions(cds)
results = nbinomTest(cds, "condition", "control")
Without replicates:
cds = estimateDispersions(cds, method = "blind")
results = nbinomTest(cds, "condition", "control")

So, my confusion goes for a reasonable strategy to conduct NB test without replicates:
I have two replicates of control group(Ctrl1,Ctrl2) while no replicates for conditon B and conditon C, which way to perform a negative binomial test fits best for this case?
My current strategy is performing NB test between Ctrl1 and B/C(ignore Ctrl2).
jason_ARGONAUTE is offline   Reply With Quote
Old 11-06-2013, 01:38 AM   #2
Devon Ryan
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480

You'll be better off using the control replicates than if you were to use the "blind" method.
dpryan is offline   Reply With Quote

deseq nb test, without replicates

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 09:07 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO