Hi All,
 I am pretty new to edgeR, this is the first time I use the package, as
well as the first time I try to find differential expression using RNA-Seq
data (this is the first time I have biological replicates).
I apologize in advance for my lack of knowledge.

I have 10 samples, 5 biological replicates of each condition, and I am just
trying to define differentially expressed genes between the two conditions
I have a matrix of counts in tab delimited format, in which the first
column is the GeneID and the first row is a header detailing the samples
this is my code:

> x <- read.delim("RNAseq_AllgenesFCHC_Counts.txt",row.names="GeneID",
header=TRUE)
> group <- factor(c(1,1,1,1,1,2,2,2,2,2))
> y <- DGEList(counts=x,group=group)
> y <- estimateCommonDisp(y)

and I get the following error:
Error in if (mx < tol) { : missing value where TRUE/FALSE needed

can anyone explain to me what the error means?

Thank you very much for your help
Lucia

	[[alternative HTML version deleted]]

