[BioC] enrichment packages that accept t-stat (or related stat) as input
Robert Castelo
robert.castelo at upf.edu
Fri Aug 30 12:14:57 CEST 2013
hi Juliet,
On 08/30/2013 02:43 AM, Juliet Hannah wrote:
> Hi Robert,
>
> Thanks for your response. I will look into it.
>
> Also is it correct GSVA always requires an expression matrix. It
> seems that it integrates with limma, so if I have done an analysis in
> limma does this mean that I should be
> able to use GSVA for an enrichment analysis.
the output of GSVA can be used as input to limma but the converse is not
true. GSVA calculates sample-wise gene set summaries of gene expression,
and therefore, it starts from a gene-by-sample expression data matrix
and produces a geneset-by-sample expression data matrix.
limma moderated t-statistics give you a measure of differential
expression for every gene over the entire set of samples, i.e., not
sample-wise, so it can be seen as one-dimensional vector whose values
represent summaries of differential expression for every gene.
to get *geneset* summaries of limma genewise t-statistics you may use
the approach i mentioned below (which you can find in the vignette of
the Category package). to gather more understand on what this approach
is doing, as well as seeing an alternative for detecting not only shifts
in mean expression but also in scale, you may consult the following
paper, others in the list may give you more appropriate bibliographic
references:
Irizarry, Rafael A.; Wang, Chi; Zhou, Yun; and Speed, Terence P., "GENE
SET ENRICHMENT ANALYSIS MADE SIMPLE" (April 2009). Johns Hopkins
University, Dept. of Biostatistics Working Papers. Working Paper 185.
http://biostats.bepress.com/jhubiostat/paper185
i do not know of what R package implements these calculations in a
user-friendly way. if you don't find any, i can paste for you here the
few lines of code that make the calculations described in this paper.
cheers,
robert.
> Thanks,
>
> Juliet
>
>
> On Thu, Aug 29, 2013 at 2:43 AM, Robert Castelo <robert.castelo at upf.edu
> <mailto:robert.castelo at upf.edu>> wrote:
>
> Juliet,
>
> i think the first 5 pages in the vignette entitled "Using Categories
> to Analyze Microarray Data" from the Category package:
>
> http://www.bioconductor.org/__packages/release/bioc/html/__Category.html
> <http://www.bioconductor.org/packages/release/bioc/html/Category.html>
>
> may be doing what you are looking for.
>
> cheers,
> robert.
>
>
> On 08/28/2013 08:04 PM, Juliet Hannah wrote:
>
> All,
>
> I am looking for an Bioconductor enrichment package that does
> something
> similar to GSEA for pre-computed test statistics. This method
> would not
> rely on a cutoff. That is, rather than passing an expression
> matrix, one
> can compute summarizes outside of the package (such as a limma
> t), and then
> pass these. Any suggestions?
>
> Thanks,
>
> Juliet
>
> [[alternative HTML version deleted]]
>
> _________________________________________________
> Bioconductor mailing list
> Bioconductor at r-project.org <mailto:Bioconductor at r-project.org>
> https://stat.ethz.ch/mailman/__listinfo/bioconductor
> <https://stat.ethz.ch/mailman/listinfo/bioconductor>
> Search the archives:
> http://news.gmane.org/gmane.__science.biology.informatics.__conductor
> <http://news.gmane.org/gmane.science.biology.informatics.conductor>
> .
>
>
> --
> Robert Castelo, PhD
> Associate Professor
> Dept. of Experimental and Health Sciences
> Universitat Pompeu Fabra (UPF)
> Barcelona Biomedical Research Park (PRBB)
> Dr Aiguader 88
> E-08003 Barcelona, Spain
> telf: +34.933.160.514 <tel:%2B34.933.160.514>
> fax: +34.933.160.550 <tel:%2B34.933.160.550>
>
>
--
Robert Castelo, PhD
Associate Professor
Dept. of Experimental and Health Sciences
Universitat Pompeu Fabra (UPF)
Barcelona Biomedical Research Park (PRBB)
Dr Aiguader 88
E-08003 Barcelona, Spain
telf: +34.933.160.514
fax: +34.933.160.550
More information about the Bioconductor
mailing list