Error at DEXSeqHTML

Question

[nitandressa]

Hi guys!

So, I'm trying to use DEXseq to get the differential "bin" expression and all analysis goes well until the point of plotting the HTML report. That's when I get the following error:

> DEXSeqHTML(ad2r, FDR=0.05, fitExpToVar="area")
Error in relevel.default(mf$exon, "others") : 
  'relevel' only for (unordered) factors
In addition: There were 12 warnings (use warnings() to see them)

I checked my design table file and it's all ok there. I also checked the counts file list and they all correspond to the design table. The code before:

> sTad2 #Design table
> 
cFad2 #List of counts file

> formulaFullModel = ~ samples + exon + batch:exon + area:exon
> formulaReducedModel = ~ samples + exon + batch:exon
> ad2 = DEXSeqDataSetFromHTSeq(cFad2, sampleData=sTad2, flattenedfile=flattenedFile, design = formulaFullModel)
> ad2=estimateSizeFactors(ad2)
> ad2=estimateDispersions(ad2, formula = formulaFullModel)
> ad2=testForDEU(ad2, reducedModel = formulaReducedModel, fullModel = formulaFullModel)
> ad2=estimateExonFoldChanges(ad2, fitExpToVar="area")
> ad2r = DEXSeqResults(ad2)
> DEXSeqHTML(ad2r, FDR=0.05, fitExpToVar="area")

Thanks in advance!