Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
How to deal with biological replicates of nucleosome MNase-seq data? wisense Epigenetics 6 02-11-2014 12:09 PM
Help Using DESeq without replicates aprilw Bioinformatics 6 12-19-2013 07:07 AM
rRNA depletion samples partially degradated ValentinadM RNA Sequencing 1 02-01-2013 03:26 AM
DESeq: question about with replicates and without any replicates. nb509 RNA Sequencing 2 10-25-2011 06:04 AM
Can DESeq and edgeR deal with in-balanced RNA-seq data? asiangg Bioinformatics 3 05-26-2010 04:45 AM

Thread Tools
Old 09-12-2013, 10:51 AM   #1
Junior Member
Location: us

Join Date: May 2013
Posts: 8
Default How to deal with samples partially without replicates in DESeq?

I have four groups of samples, 2 groups with replicates and 2 groups without. The data frame is like this:
+ condition=c("GF","SPF","9343","9343","SPT","SPT"))
Error: unexpected symbol in:
> expt_design<-data.frame(row.names=colnames(colon),
+ condition=c("GF","SPF","9343","9343","SPT","SPT"))
> expt_design
X9343.1 9343
X9343.2 9343
> conditions=expt_design$condition
> conditions
[1] GF SPF 9343 9343 SPT SPT
Levels: 9343 GF SPF SPT

When I tried to run it gave me this error message:
Error in .local(object, ...) :
None of your conditions is replicated. Use method='blind' to estimate across conditions, or 'pooled-CR', if you have crossed factors.

Plz help me to identify the problem and tell me how to deal with partially replicated samples
zatoichi888 is offline   Reply With Quote
Old 09-12-2013, 11:08 AM   #2
Devon Ryan
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480

Did you read the error message, the fix is contained within it. Also, you should read the DESeq tutorial.
dpryan is offline   Reply With Quote
Old 09-12-2013, 11:12 PM   #3
Simon Anders
Senior Member
Location: Heidelberg, Germany

Join Date: Feb 2010
Posts: 994

You might have confused DESeq by passing your "condition" data as a data.frame rather than as a factor. Try again using

cds <- newCountDataSet( countTable, exptDesign$condition)

And consider switching to DESeq2. It not only has improved power; we have also improved the interface to avoid problems like this one.
Simon Anders is offline   Reply With Quote

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 01:28 AM.

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