SEQanswers

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

Reply
 
Thread Tools
Old 11-05-2013, 09:51 PM   #1
jason_ARGONAUTE
Member
 
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
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

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

Tags
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