While running DESeq2 on a Human dataset (commands given below), I have noticed the following error. How I can handle my data with this case? Is that the same case of Error in estimateDispersions of DEXSeq? Could you tell me how to solve the problem, please?
> sampleFiles =grep("count",list.files("/folder/to/path/"),value=TRUE)
> sampleTable<-data.frame(sampleName=sampleFiles, fileName=sampleFiles, condition=sampleCondition)
> ddsHTSeq<-DESeqDataSetFromHTSeqCount(sampleTable=sampleTable, directory="/folder/to/path/", design=~condition)
# estimating size factors
# estimating dispersions
# gene-wise dispersion estimates
# mean-dispersion relationship
# Error in estimateDispersionsFit(object, fitType = fitType, quiet = quiet) :
# all gene-wise dispersion estimates are within 2 orders of magnitude
# from the minimum value, and so the standard curve fitting techniques will not work.
# One can instead use the gene-wise estimates as final estimates:
# dds <- estimateDispersionsGeneEst(dds)
# dispersions(dds) <- mcols(dds)$dispGeneEst