[BioC] MAS on Bioconductor
    Rafael A. Irizarry 
    ririzarr at jhsph.edu
       
    Mon Mar 24 14:53:29 MET 2003
    
    
  
for the answer to your mas 5 question read this:
http://stat-www.berkeley.edu/users/bolstad/MAS5diff/Mas5difference.html
for writing pms to delimeted use the "pm" function to extract the pms 
out of your affybatch, for example:
Data <- ReadAffy()
pms <- pm(Data)
and read the help file on write.table.
On 
Mon, 24 Mar 2003, Kelvin 
Yen wrote:
> Hi,
> 
> I'm a new user to Bioconductor and I've been using the Affy package to 
> analyze the .Cel files I have.  I tried replicating the analysis I did 
> using Affy's MAS program with bioconductor's MAS algorithm but the 
> results I get are completely different than the ones I get using the 
> original MAS software.  I'm analyzing six .Cel files and using MAS and 
> a subsequent t-test on the values, I get many(~30) values that are 
> e-05.  But using Bioconductor's MAS algorithm, I only get one.  Here 
> are the commands that I'm using to generate the data.
> 
> Data <- ReadAffy()
> eset <- expresso(Data, normalize=FALSE, bgcorrect.method="mas", 
> pmcorrect.method="mas, summary.method="mas")
> eset <- affy.scalevalue.exprSet(eset)
> write.exprs(eset, file="mas.xls")
> 
> I'm using the MGU74aV2 chip.  Any help would be greatly appreciated.
> 
> Kelvin Yen
> PS I eventually wanted to run the RMA algorithm on my data but first 
> wanted to make sure I could replicate data that I already knew the 
> value to.
> 
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
> https://www.stat.math.ethz.ch/mailman/listinfo/bioconductor
>
    
    
More information about the Bioconductor
mailing list