Question: MAS5 algorithm in mas5() function of xps function when process Human Gene 1.0 ST Array or Human Exon 1.0 ST Array
0
gravatar for linweili_tj
3.0 years ago by
linweili_tj0 wrote:

Human Gene 1.0 ST array and exon arrays have only PM probesets and no MM probesets, thus generally it is not possible to perform MAS5 with these arrays. However, R package "xps" can indeed apply MAS5 to exon arrays and Human Gene 1.0 ST array. What is the detailed algorithm of MAS5 in "xps" package when it process these arrays? What are the differences between the modified MAS5 algorithm of "xps" and the regular MAS5 algorithm mentioned in sadd_whitepaper.pdf (Affymetrix website)?

ADD COMMENTlink modified 3.0 years ago • written 3.0 years ago by linweili_tj0
Answer: MAS5 algorithm in mas5() function of xps function when process Human Gene 1.0 ST
2
gravatar for cstrato
3.0 years ago by
cstrato3.9k
Austria
cstrato3.9k wrote:

Hi,

For IV arrays the MAS5 algorithm implemented in 'xps' is identical to the MAS5 
algorithm described in 'sadd_whitepaper.pdf'. Please have a look at Fig.4 of 
Chapter 2.2 of vignette 'APTvsXPS.pdf' where I compare the different implementations 
of MAS 5, i.e. Expression Console (EC), APT, xps, and affy. 

To understand why 'xps' is able to implement MAS5 for whole genome and exon arrays
you need to know how Affymetrix does calculate the background intensities for the
different arrays: 

- For IV arrays Affymetrix does not simply subtract MM from PM for 
each probe, but first calculates an 'Ideal Mismatch' (IM), which is then subtracted 
from PM.

- For whole genome and exon arrays Affymetrix does not use MM probes but uses 
special ”genomic” and ”antigenomic” probes to compute an IM.


Affymetrix explains this already in 'Expression analysis algorithm tutorial.pdf'
for the original MAS4 algorithm, see Chapter III.c. Background Calculation. As
you see the array is divided into 4x4 sectors, and the IM is calculated for each
sector. 
Please read Chapter 3.2 (sector) of my vignette 'xpsPreprocess.pdf'.

In their 'sadd_whitepaper.pdf' for MAS5 Affymetrix explains that to calculate the
Ideal Mismatch (IM) they use the one-step biweight algorithm (page 6). Once again 
the array is divided into 4x4 sectors (see the Figure on page 4). 
Please read Chapter 3.3 (weightedsector) of my vignette 'xpsPreprocess.pdf'.

For whole genome and exon arrays Affymetrix has created background probes with 
different GC-contents, see their 'exon_background_correction_whitepaper.pdf'.
However, these probes are now used to compute a background for PM probes with 
the same GC-content. 
Please read Chapter 3.4 (gccontent) of my vignette 'xpsPreprocess.pdf'.

In summary, this means that I can compute the background intensities separately
for each type of array and then use the calculated background probes for further
processing. In the case of the whole genome and exon arrays I use simply the
background intensities computed with the GC-content.

Best regards,
Christian

ADD COMMENTlink written 3.0 years ago by cstrato3.9k

Dear Professor Stratowa,

Thank you for your detailed and nice replies! I understand the modified MAS5 algorithm implemented in "xps" package. I can not find Expression analysis algorithm tutorial.pdf on Affymetrix website. Can you send me the tutorial or the web link?

Best regards,

Weili

ADD REPLYlink written 3.0 years ago by linweili_tj0
Dear Weili Please find attached the tutorial. Best regards, Christian On 08/31/16 10:28, linweili_tj [bioc] wrote: > Activity on a post you are following on support.bioconductor.org > <https: support.bioconductor.org=""> > > User linweili_tj <https: support.bioconductor.org="" u="" 11383=""/> wrote > Comment: MAS5 algorithm in mas5() function of xps function when process > Human Gene 1.0 ST Array or Human Exon 1.0 ST Array > <https: support.bioconductor.org="" p="" 86604="" #86643="">: > > Dear Professor Stratowa, > > Thank you for your detailed and nice replies! I understand the modified > MAS5 algorithm implemented in "xps" package. I can not find /*Expression > analysis algorithm tutorial.pdf */on Affymetrix website. Can you send me > the tutorial or the web link? > > Best regards, > > Weili > > ------------------------------------------------------------------------ > > Post tags: xps, MAS5 algorithm, Human Gene 1.0 ST Array, Exon Array > > You may reply via email or visit > C: MAS5 algorithm in mas5() function of xps function when process Human Gene 1.0 ST >
ADD REPLYlink written 3.0 years ago by cstrato3.9k
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: 245 users visited in the last hour