hi Hui Zhao,

This code looks correct for the exploratory analysis without replicates.
Can you explain and paste the content of the errors you are experiencing?

It looks like you are rounding the values of the readcount object. Note
that our software is designed to work on raw, unnormalized counts of reads
in features.

Mike


On Sun, Dec 8, 2013 at 10:41 PM, Hui Zhao [guest] <guest@bioconductor.org>wrote:

>
> I want to know how to deal with data without replicates using DESeq2,can
> you tell me how to deal with it.You can see my codes below,but it have
> errors when deal with sample without replicates.
>
>  -- output of sessionInfo():
>
>      library('DESeq2')
>      readcount<-read.delim("readcount.xls",row.names=1)
>      readcount<-round(readcount)
>      readcount<-readcount[,c(1,3)]
>      colData<-data.frame(condition=c("A","B"))
>      row.names(colData)<-c("A","B")
>      colData
>
>  dds<-DESeqDataSetFromMatrix(countData=readcount,colData=colData,design=~condition)
>
>  colData(dds)$condition<-factor(colData(dds)$condition,levels=c("A","B"))
>      dds <- DESeq(dds)
>      res<-results(dds)
>      res<-res[order(res$padj),]
>      res<-as.data.frame(res)
>
>
> --
> Sent via the guest posting facility at bioconductor.org.
>

	[[alternative HTML version deleted]]

