how to do the following using justRMA in affy package?
2
0
Entering edit mode
@james-anderson-1641
Last seen 10.2 years ago
An embedded and charset-unspecified text was scrubbed... Name: not available Url: https://stat.ethz.ch/pipermail/bioconductor/attachments/20070510/ e2282ce8/attachment.pl
• 379 views
ADD COMMENT
0
Entering edit mode
Seth Falcon ★ 7.4k
@seth-falcon-992
Last seen 10.2 years ago
James Anderson <janderson_net at="" yahoo.com=""> writes: > 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? Try reading the help pages for these functions: list.files sample save paste for (i in 1:50) { Create a character vector of your 100 CEL files. Sample 70 elements. Run justRMA save to file paste("data_", i, ".txt", sep="") } + seth -- Seth Falcon | Computational Biology | Fred Hutchinson Cancer Research Center http://bioconductor.org
ADD COMMENT
0
Entering edit mode
@james-w-macdonald-5106
Last seen 12 hours ago
United States
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.
ADD COMMENT

Login before adding your answer.

Traffic: 914 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6