[BioC] a warining from edgeR
Gordon K Smyth
smyth at wehi.EDU.AU
Tue Jul 31 01:09:59 CEST 2012
It would appear that one gene out of 54374 has received a dispersion
estimate that has not fully converged. That is very unlikely to affect
your overall results.
In the devel version of edgeR, to become the release version in October,
this warning will no longer occur, because we now compute the maximum of
the spline interpolant exactly.
Best wishes
Gordon
> Date: Sun, 29 Jul 2012 11:13:06 -0400
> From: wang peter <wng.peter at gmail.com>
> To: bioconductor at r-project.org
> Subject: [BioC] a warining from edgeR
>
> I have 35 samples and 54374 genes
>
> such is the coding, do u think is it ok?
>
>
>
> #normalization
> dge <- calcNormFactors(dge)
>
> treatment=factor(c(rep('control',6),rep('treated',24),rep('control',5)))
> time=factor(c('0h','0h','0h','24h','24h','24h','0h','0h','0h','6h','6h','6h','6h','12h','12h','12h','12h','18h','18h','18h','18h',
> '24h','24h','24h','36h','36h','36h','48h','48h','48h','6h','12h','18h','36h','48h'))
> design <- model.matrix(~treatment*time)
>
> dge <- estimateGLMCommonDisp(dge, design)
> dge <- estimateGLMTagwiseDisp(dge, design)
>
>
>
> Warning message:
> In maximizeInterpolant(spline.pts, apl.smooth[j, ]) :
> max iterations exceeded
>
>
>
>
>
>> sessionInfo()
> R version 2.15.1 (2012-06-22)
> Platform: x86_64-pc-mingw32/x64 (64-bit)
>
> locale:
> [1] LC_COLLATE=Chinese (Simplified)_People's Republic of China.936
> [2] LC_CTYPE=Chinese (Simplified)_People's Republic of China.936
> [3] LC_MONETARY=Chinese (Simplified)_People's Republic of China.936
> [4] LC_NUMERIC=C
> [5] LC_TIME=Chinese (Simplified)_People's Republic of China.936
>
> attached base packages:
> [1] stats graphics grDevices utils datasets methods base
>
> other attached packages:
> [1] edgeR_2.4.6 limma_3.10.3
>
>
> --
> shan gao
> Room 231(Dr.Fei lab)
> Boyce Thompson Institute
> Cornell University
> Tower Road, Ithaca, NY 14853-1801
> Office phone: 1-607-254-1267(day)
> Official email:sg839 at cornell.edu
> Facebook:http://www.facebook.com/profile.php?id=100001986532253
>
______________________________________________________________________
The information in this email is confidential and intend...{{dropped:4}}
More information about the Bioconductor
mailing list