[BioC] decideTests

Jason Skelton jps at sanger.ac.uk
Thu Nov 4 11:45:30 CET 2004


Hi Gordon

>Well, the error message is pretty self-explanatory!  You need to give it an object of class
>"MArrayLM", i.e., you need to give it output from lmFit() and eBayes(), as in all the examples in
>the User's Guide.
>
>You didn't say how you created FIT and EB.  I am guessing now that you used lm.series() or
>gls.series().  I am been phasing these functions out of the user-interface for a long time -- they
>are lower-level functions that produce unclassed list objects.  There is no need for you to use
>them.  So far I have been writing functions that are tolerant of whether the input arguments are
>formally classed or not, but in future I may not do that.
>
>So, please, use lmFit(), not lm.series().
>  
>
Sorry, but i used lmFit()

 > COR <- duplicateCorrelation(nwaMA, ndups=2, spacing=840, DESIGN)
 > COR$cor
[1] 0.3662469
 > FIT <- lmFit(nwaMA, ndups=2, spacing=840, DESIGN, 
correlation=0.3662469, method="ls")
 > class(FIT)
[1] "MArrayLM"
attr(,"package")
[1] "limma"

I've been using this for a long time now and stopped using lm.series() a 
long time ago....

I've now managed to get this to work however
I was using the wrong version of ebayes.......?
I was using ebayes and NOT eBayes (mainly from habit)
apologies for the confusion

"you need to give it output from lmFit() and eBayes()"


Well just the output from eBayes ?
As the former gives the following message:

Error in array(x, c(length(x), 1), if (!is.null(names(x))) list(names(x),  :
        attempt to set an attribute on NULL


Thanks for your help and advice
Cheers

Jason

||||



More information about the Bioconductor mailing list