I am using DEXSeq so far without any problem. But I am interested in normalized counts by length in dexseq plot. Therefore, I assume i will get considerably less fluctuated expressions across exons. To do that, I have changed the countData in before plotting, (i know the analyze is based on count data, therefore i did not do that at the beginning of the analyze). Like that;
dxr2$countData <- count.norm.data plotDEXSeq( dxr2, "ENSG00000107771", legend=TRUE, cex.axis=1.2, cex=1.3, lwd=2 )
Well, it worked for some cases, but for some cases i have following error;
Fit for gene/exon ENSG00000182718 failed, coefficients for this gene won't show up. NULL
Should I change something from source code ?
So, any comments and help will be much appreciated.