Affymetrix Power Tools To Measure Whether The Gene Is Expressed Or Not?
1
0
Entering edit mode
12.1 years ago
GR ▴ 400

Dear All,

I have few genes that I want to know are expressed or not in few tissues and conditions (microarray data). I was wondering if anyone has used affymetrix power tools to detect the signal intensities to say whether a gene is expressed in the given condition/tissue etc. I got the normalized signal intensities but I don't know what threshold to use to say a particular gene is expressed or not.

Any other way/steps to do this will also be helpful.

Thanks, RT

gene expression microarray • 2.6k views
ADD COMMENT
0
Entering edit mode
12.1 years ago
Ido Tamir 5.2k
  • mas5:

In mas5 there is the stratification into a (absent), m (marginally), p(present). I think the mas5 normalization is not very good (underestimates the low range (or was it overestimates ?, I forgot)). The a,m,p classification however is rather accurate.

  • GCRMA:

You could normalize with GCRMA, then plot the density of expression. You normally see a huge peak that represents the not expressed genes, then expressed genes come up. 2 distributions: low end = background intensities of not expressed genes which fluctuates higher = intensities from expressed genes. Choose a cutoff somewhere where you think most of the data comes from the expressed genes. In the data shown I would say that between intensities 4-5 probes (not genes!!) might be detected as expressed (~ marginally in mas5) and above 5 they are reliably expressed and detected.

http://bioinformatics.knowledgeblog.org/files/2011/06/Screen-shot-2011-06-15-at-15.44.11.png

ADD COMMENT
0
Entering edit mode

Hi Ido,

Do you know if there is any way to extract p-values for the present/absence calls. I could just get the P/A/M labels using the following R code. I want to add p-values to these results. Any suggestions?

Thanks, Reetu

affy.data <- ReadAffy()

data.mas5calls = mas5calls(affy.data)

data.mas5calls.calls = exprs(data.mas5calls)

ADD REPLY

Login before adding your answer.

Traffic: 1452 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6