Question: tutorial/documentation to use oligo/affy for mta1.0 Affymetrix arrays
0
gravatar for iankoudinova
3.3 years ago by
iankoudinova10
iankoudinova10 wrote:

I want to run Bioconductor packages on mta1.0 affymetrix arrays. I see that there are pd.mta.1.0, mta10stprobeset.db and mta10sttranscriptcluster.db available. There is no workflow for this type of arrays yet availbale. I tried to use affy, and received the message that

The affy package is not designed for this array type.
Please use either the oligo or xps package.

Can anyone give me some instructions what should I use, what kinds of files should be present in the current directory? Do I need to upload something from Affymetrix? I do already have mta10 cdf environment.

I built cdf environment for this type of arrays and run the local QC script which requires simpleaffy, affyPLM, affyQCReport, arrayQualityMetrics. It does work fine with mta10cdf environment I provided.  It helps me to choose outliers  Now  I want to run differential expression analysis on the same CEL files of mta1.0 arrays. 

 

Thank you

affy mta10 analysis • 778 views
ADD COMMENTlink modified 3.3 years ago • written 3.3 years ago by iankoudinova10
Answer: tutorial/documentation to use oligo/affy for mta1.0 Affymetrix arrays
0
gravatar for James W. MacDonald
3.3 years ago by
United States
James W. MacDonald51k wrote:

It depends on what you intend to do with the MTA array. The main selling point for these arrays is that you can test for differential transcript abundance. If you want to do that sort of thing, then the only recourse I know of is to use Affy's software (Expression Console followed by Transcriptome Analysis Console).

If instead you want to use the array like a Gene ST array, and just summarize all probes for each transcript, then the workflow isn't different than any other Affy microarray.

library(oligo)
dat <- read.celfiles(list.celfiles()) 
## assuming you are running R in the dir with the celfiles
## otherwise you need to do
dat <- read.celfiles(filenames = <vector of file names with path>)
eset <- rma(dat)

Then make comparisons using something like limma.

ADD COMMENTlink written 3.3 years ago by James W. MacDonald51k
Answer: tutorial/documentation to use oligo/affy for mta1.0 Affymetrix arrays
0
gravatar for iankoudinova
3.3 years ago by
iankoudinova10
iankoudinova10 wrote:

Ok. Thank you. It is very helpful. 

Now is the question how to convert mta1.0 probe set ids to gene symbols and/or entrez?

ADD COMMENTlink written 3.3 years ago by iankoudinova10

If you have a comment, or additional question, click on the ADD COMMENT link and type in the box that appears. The Add  your answer box below is intended for people to use if they are actually answering a question.

To annotate any microarray, you use the corresponding ChipDb package. In this case it's the mta10sttranscriptcluster.db package, and you would use the select function.

symbols <- select(mta10sttranscriptcluster.db, <vector of IDs>, "SYMBOL")

or perhaps mapIds

symbols <- mapIds(mta10sttranscriptcluster.db, <vector of IDs>, "SYMBOL", "PROBEID")

There are vignettes for the AnnotationDbi package that explain things, and I gave a workshop at BioC2016 that may be useful as well.

 

ADD REPLYlink written 3.3 years ago by James W. MacDonald51k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 237 users visited in the last hour