[BioC] justMAS reliability

Eitan Halper-Stromberg eitan.halper-stromberg at nist.gov
Fri Jun 24 21:50:54 CEST 2005


Hi,
   I am trying to speed up my mas5 processing by using justMAS on my  
affybatch. My AffyBatch is made of 8 chips and the cdf for these  
chips is hgu133a. I have tested justMAS by running it several times  
on this affybatch and I am obtaining different expression measure  
results even though I am doing the exact same command. My affybatch  
is called 'F'  and an array of geneNames is called 'gn', and I type:

 > j<-justMAS(F)

 > j at exprs[gn=="1007_s_at"]
1007_s_at      <NA>      <NA>      <NA>      <NA>      <NA>       
<NA>      <NA>
6.603113  6.171309  6.399581  6.171309  6.171309  6.171309  6.171309   
6.036457

 > i<-justMAS(F)

 > i at exprs[gn=="1007_s_at"]
1007_s_at      <NA>      <NA>      <NA>      <NA>      <NA>       
<NA>      <NA>
6.603113       NaN  6.171309  5.482231  6.171309  5.438713  5.387300   
7.311159

Does anybody know why I am getting this discrepancy?

-Eitan



More information about the Bioconductor mailing list