That doesn't help me in this context. I don't know how to select specific genes from an Affymetrix-type dataset. Once selected, sure, a basic t-test or regional ANOVA will be fine, but I can do neither of those things until I can specify the genes I'm interested in.
Comparing Two Genes In An Affymetrix Microarray Dataset In R
I'm working with an Affymetrix rat transcriptome in R and would like to compare, say, expression levels of MAOA and MAOA versus a housekeeping gene like ACTB. I've gone through the basic initial stages:
rawdata = ReadAffy()
eset = rma(rawdata)
and am aware of the affy and limma packages but not entirely sure how to go about specifying only two or three genes for analysis/comparison. Help?
1 answer
Oh ok.As far as I know the AffyBatch objects do not contain the gene symbol annotations.You will need to use some annotation file to map TranscriptClusterIDs to gene symbols. You could use Bioconductor packages to do that or download annotation data directly from Affy.
In your case exprs(eset) should give you an expression matrix where the rows are the probesets and the row names are the TranscriptClusterIDs.
Log in to answer this question.
More posts like this
-
Question on performing RMA normalization
written by Natasha •
-
Microarray analysis using limma multiple groups
written by Paul
-
lesser fold change value using RMA background correction
written by Paul
-
differential expression analysis from microarray data
written by au.rinki.bio •
-
Analysis of microarray data: "eset" problem
written by orzech_mag
-
Microarray Data Analysis Using Bioconductor
written by H@rry •
-
What Is The Best Way To Calculate Gene Expression From Affymetrix Data
written by Sirus
-
Extract Pre-Normalisation Expression Data From Probe Set Values
written by munch •
-
Affymetrix Data Using Mas5 Giving Large Fc Values
written by Bade •
-
Header Row In Write.Exprs Is File Names Instead Of Designated Sample Names
written by Gotgenes
I dont know if I understand correctly. I think you first must to annotate. Maybe http://www.bioconductor.org/help/workflows/annotation-data/ could help. After that, I think you can use http://bioconductor.org/packages/release/bioc/html/affycomp.html.