This is a test version of Biostars. For the public version, visit https://www.biostars.org.
R library/package for pairwise differential expression for microarray data?

Does anyone know of a reliable package for doing pairwise DE analysis on microarray data in the form of an ExpressionSet object? I've used edgeR before for doing pairwise comparison with rna-seq, but I'm a little new to this and can't really find anything to do it on affy data specifically...

microarray affy differential expression
• 3,865 views • link
updated by Ram Ram • written by -5¢ • -5¢

1 answer

I use samr:

R. Tibshirani, G. Chu, Balasubramanian Narasimhan and Jun Li (2011). samr: SAM: Significance Analysis of Microarrays. R package version 2.0. http://CRAN.R-project.org/package=samr

Expression data are a matrix E, with identifiers PROBE_IDENTIFIERS and symbols GENE_SYMBOLS, and labels indicating pairs. Labels are integers starting at 1, with positive or negative indicating paired samples. Works for any array (or indeed any matrix of numbers).

data = list(x=data.matrix(E),
            y=labels,
            geneid=PROBE_IDENTIFIERS
            genenames = GENE_SYMBOLS,
            logged2=TRUE)
samr.obj = samr(data, resp.type="Two class paired", nperms=100)
• 1 views • link
updated by Ram Ram • written by David Quigley David Quigley

thanks very much!

• 1 views • link
written by -5¢ • -5¢

You're very welcome. If that answers your question, please mark the question as answered.

• 1 views • link
updated by Ram Ram • written by David Quigley David Quigley

Log in to answer this question.

More posts like this