[BioC] how to do the following using justRMA in affy package?

James W. MacDonald jmacdon at med.umich.edu
Fri May 11 15:53:48 CEST 2007


James Anderson wrote:
> Hi,
> 
> I have 100 cel files (U133A) in a folder, I want to randomly select
> 70 cel files from them and use justRMA to get the normalized
> intensity. I need to do this 50 times (every time, the 70 cel files
> are not necessarily the same). More imporantly, I want to save the
> expression value obtained by justRMA to file name Data_1.txt
> Data_2.txt Data_3.txt, ..., Data_50.txt in an automatic fashion. I
> tried myself first and could not figure out how to do it, could
> anybody offer me any help or suggestions?

library(affy)
for(i in 1:50){
	ind <- sample(1:100, 70)
	eset <- justRMA(filenames=list.celfiles()[ind])
	write.exprs(eset, paste("Data_", i, ".txt", sep=""), 					col.names=NA)
}

Best,

Jim


> 
> Many thanks,
> 
> James
> 
> 
> ---------------------------------
> 
> 
> [[alternative HTML version deleted]]
> 
> _______________________________________________ Bioconductor mailing
> list Bioconductor at stat.math.ethz.ch 
> https://stat.ethz.ch/mailman/listinfo/bioconductor Search the
> archives:
> http://news.gmane.org/gmane.science.biology.informatics.conductor


-- 
James W. MacDonald, M.S.
Biostatistician
Affymetrix and cDNA Microarray Core
University of Michigan Cancer Center
1500 E. Medical Center Drive
7410 CCGC
Ann Arbor MI 48109
734-647-5623


**********************************************************
Electronic Mail is not secure, may not be read every day, and should not be used for urgent or sensitive issues.



More information about the Bioconductor mailing list