[BioC] anyway to speed up anova on every row of a 280k-row matrix
Robert Gentleman
rgentlem at fhcrc.org
Tue Apr 7 18:57:31 CEST 2009
if the design matrix is constant then there are all sorts of way to speed things
up. Using apply is probably the worst thing you could do...
have a look at GSEAlm for some ways to compute fits, residuals and other
diagnostics quickly.
best wishes
Robert
Kasper Daniel Hansen wrote:
> Use limma, as you are supposedly fitting the same model to every row.
>
> Kaspr
>
> On Apr 7, 2009, at 8:12 , Zheng, Xin (NIH) [C] wrote:
>
>> Hi there,
>>
>> I was doing anova on an array file. I know Partek can do it in
>> minutes, yet it cost half an hour in R with apply(matrix, 1, function
>> of anova). Any hint? Thanks.
>>
>> Xin Zheng
>>
>>
>> [[alternative HTML version deleted]]
>>
>> _______________________________________________
>> Bioconductor mailing list
>> Bioconductor at stat.math.ethz.ch
>> https://stat.ethz.ch/mailman/listinfo/bioconductor
>> Search the archives:
>> http://news.gmane.org/gmane.science.biology.informatics.conductor
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives:
> http://news.gmane.org/gmane.science.biology.informatics.conductor
>
--
Robert Gentleman, PhD
Program in Computational Biology
Division of Public Health Sciences
Fred Hutchinson Cancer Research Center
1100 Fairview Ave. N, M1-B514
PO Box 19024
Seattle, Washington 98109-1024
206-667-7700
rgentlem at fhcrc.org
More information about the Bioconductor
mailing list