library("affy")
Data<-ReadAffy();
eset<-mas5(Data,normalize=T)
if (method=='mas5')
{
if(norm=='F')
{ eset<-mas5(Data,normalize=F) } else { eset<-mas5(Data,normalize=T) }
}
if (method=='rma')
{ eset <-rma(Data) }
chip.names <- rownames( attr(attr(Data,'phenoData'),'data') );
expdata <- exprs(eset);
colnames(expdata) <- chip.names;
write.table(expdata,file="mas5.txt",sep="\t");
沒有留言:
張貼留言