[Bioc-devel] Bug in DEXSeq, me thinks

[Steve Lianoglou]

And then add names back to muhats, of course:

On Mon, Mar 12, 2012 at 8:32 PM, Steve Lianoglou
<mailinglist.honeypot at gmail.com> wrote:
> Hi,
>
> It just so happened that DEXSeq, for reasons I couldn't explain, kept
> throwing errors during a call to `estimateDispersions`. I couldn't
> explain it because I really thought I didn't change much, but I've got
> some magic combination of samples I'm currently trying to
> estimateDispersion on that none of them throw any errors while
> estimate `muhats` (line 245 of core_functions.R):
>
> muhats <- sapply(testablegenes, function(gn) {
> ...
> muhat
> })
[snip]

names(muhats) <- testablegenes

-steve
-- 
Steve Lianoglou
Graduate Student: Computational Systems Biology
 | Memorial Sloan-Kettering Cancer Center
 | Weill Medical College of Cornell University
Contact Info: http://cbio.mskcc.org/~lianos/contact