Error at DEXSeqHTML
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!
• 991 views
•
link
0 answers
No answers yet.
Log in to answer this question.