hello,
I have started working with DESeq for RNASeq data from plants.
Among the first plots I have got, there is the dispersion (D) vs counts mean (M). This plot is fitted with the curve D=a+b/M (which in a log-log scale is a decreasing curve that reaches a plateau for high value of M, see attachment).
However, my data show a linear correlation (still in the log-log plot).
Is my analysis wrong? can I change the fitting expression? why the fitting procedure does not set a=0 which would give a linear relationship between Log(D) and Log(M)?
thanks!
I have started working with DESeq for RNASeq data from plants.
Among the first plots I have got, there is the dispersion (D) vs counts mean (M). This plot is fitted with the curve D=a+b/M (which in a log-log scale is a decreasing curve that reaches a plateau for high value of M, see attachment).
However, my data show a linear correlation (still in the log-log plot).
Is my analysis wrong? can I change the fitting expression? why the fitting procedure does not set a=0 which would give a linear relationship between Log(D) and Log(M)?
thanks!
Comment