Dear list members, I am having problems using the function estimateGLMTrendedDisp() in edgeR. When I run it on my dataset, it returns the error: eset_b_tmm <- estimateGLMTrendedDisp(eset_a_tmm , design) Error in dispCoxReid(y, design, offset = offset, ...) : no data rows with required number of counts If I run the estimateGLMCommonDisp() function in my dataset it works. If I run the example from the Documentation page it works. I have read the documentation but I do not find the source of my problem. My object eset_a_tmm is a "DGEList". The only difference I see between my dataset and the y from the Documentation page is that my data has samples and some genes with 0 (zero) counts. Does anyone have a clue of what I am doing wrong? Thanks very much, Fernando sessionInfo() R version 2.13.0 (2011-04-13) Platform: x86_64-pc-mingw32/x64 (64-bit) locale: [1] LC_COLLATE=English_United States.1252 LC_CTYPE=English_United States.1252 [3] LC_MONETARY=English_United States.1252 LC_NUMERIC=C [5] LC_TIME=English_United States.1252 attached base packages: [1] splines stats graphics grDevices utils datasets methods base other attached packages: [1] DESeq_1.4.0 locfit_1.5-6 lattice_0.19-23 akima_0.5-4 Biobase_2.12.0 [6] edgeR_2.2.1 loaded via a namespace (and not attached): [1] annotate_1.30.0 AnnotationDbi_1.14.0 DBI_0.2-5 [4] genefilter_1.34.0 geneplotter_1.30.0 grid_2.13.0 [7] limma_3.8.0 RColorBrewer_1.0-2 RSQLite_0.9-4 [10] survival_2.36-5 tools_2.13.0 xtable_1.5-6 [[alternative HTML version deleted]]