Warning in lfproc when using rlog or varianceStabilisingTransformation from DESeq2
1
0
Entering edit mode
katja • 0
@katja-7627
Last seen 8.9 years ago
Sweden

Hi,

during the analysis of sRNA sequencing data I encountered a warning that I do not understand.

From counts of sRNA sequences I built a DESeqDataSet (7 groups of sRNA sequences x 28 samples) and performed normalisation with rlog function.

rld <- rlog(dds, blind = TRUE)

After that I got a warning:

NOTE: fitType='parametric', but the dispersion trend was not well captured by the
  function: y = a/x + b, and a local regression fit was automatically substituted.
  specify fitType='local' or 'mean' to avoid this message next time.

Warning message:
In lfproc(x, y, weights = weights, cens = cens, base = base, geth = geth,  :
  Estimated rdf < 1.0; not estimating variance

When I used fitType = "mean", there was no warning, but it appeared when I chose fitType = "local":

Warning message:
In lfproc(x, y, weights = weights, cens = cens, base = base, geth = geth,  :
  Estimated rdf < 1.0; not estimating variance

The same thing also happend when I used a function varianceStabilisingTransformation. I also used blind = TRUE and tried different settings of fitType.

I read the vignette from DESeq2 package, documentation for rlog and varianceStabilisingTransformation, but could not find the answer, what the warning means. I also could not find the explanation of lfproc function or rdf. Did anyone encounter the same problem or does anyone know what it means? Could it be that there is not enough datapoints for the methods to work?

Thanks for your help,

Katja

DESeq2 rlog varianceStabilizingTransformation lfproc • 2.8k views
ADD COMMENT
0
Entering edit mode
@mikelove
Last seen 1 minute ago
United States
I'll take a look at writing a better warning. Yes, with only 7 rows, there are too few points to fit a smooth trend line of dispersion ~ mean. So the mean fit type is all I would use.
ADD COMMENT

Login before adding your answer.

Traffic: 891 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6