[BioC] Error DEXSeq

[Roberta Carriero [guest]]

I have a warning message when I estimate the dipersion parameter. I can't carry on because of some other error messages in the following steps, as I attached below.
Thanks in advance

 -- output of sessionInfo(): 

sizeFactors (ecs)
 7A31.counts  7A32.counts  7A33.counts 46BR1.counts 46BR2.counts 46BR3.counts 
   1.0177858    1.0753635    1.1738656    0.9085247    0.9827212    0.9198726 
> ecs<- estimateDispersions ( ecs )
Dispersion estimation. (Progress report: one dot per 100 genes)
...................................................................................................................................
Warning messages:
1: In .local(object, ...) :
  Exons with less than 11 counts will be discarded. For more details read the documentation, parameter minCount
2: In .local(object, ...) :
  Genes with more than 70 testable exons will be kicked out of the analysis. For more details read the documentation, parameter maxExon
> ecs <- fitDispersionFunction(ecs)
Error in if (sum(log(coefs/oldcoefs)^2) < 0.005) break : 
  missing value where TRUE/FALSE needed
In addition: Warning messages:
1: In glmgam.fit(mm, disps[good], start = coefs) :
  Too much damping - convergence tolerance not achievable
2: In log(coefs/oldcoefs) : NaNs produced
> head( fData (ecs)$dispBeforeSharing)
[1] 0.000000e+00 6.003931e-03 4.852941e-03 0.000000e+00 2.771527e-10
[6] 0.000000e+00
> ecs at dispFitCoefs
[1] NA NA
> head(fData(ecs)$dispFitted)
[1] NA NA NA NA NA NA
> plotDispEsts(ecs)
Error: could not find function "plotDispEsts"
> ecs<-testForDEU(ecs)
Error in testForDEU(ecs) : 
  No dispersion values found, call function fitDispersionFunction first.
> 


--
Sent via the guest posting facility at bioconductor.org.