[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?

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)



> 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.
Affymetrix and cDNA Microarray Core
University of Michigan Cancer Center
1500 E. Medical Center Drive
7410 CCGC
Ann Arbor MI 48109

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