CEL.list
2
0
Entering edit mode
@denrique-escobar-espinoza-1303
Last seen 10.3 years ago
i have decided to read cel files from a list containing full path to each cel files, one per line: &$>head /atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/CEL.list /atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R0.CEL /atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R1.CEL /atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03A20M1-R0.CEL /atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03C20M1-R0.CEL /atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03D20M1-R0.CEL /atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03E20M1-R0.CEL /atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03F20M1-R0.CEL /atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03G20M1-R0.CEL in R console i ve done: > (cels =celslist ); [[1]] [1] "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R0.CEL" [2] "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R1.CEL" [3] "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03A20M1-R0.CEL" [4] "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03C20M1-R0.CEL" [5] "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03D20M1-R0.CEL" [6] "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03E20M1-R0.CEL" [7] "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03F20M1-R0.CEL" [8] "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03G20M1-R0.CEL" > (cels1=celslist[[1]]); [1] "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R0.CEL" [2] "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R1.CEL" [3] "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03A20M1-R0.CEL" [4] "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03C20M1-R0.CEL" [5] "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03D20M1-R0.CEL" [6] "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03E20M1-R0.CEL" [7] "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03F20M1-R0.CEL" [8] "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03G20M1-R0.CEL" in the fisrt caSE <(cels =celslist)> i have this error: Error in just.rma(filenames = l$filenames, phenoData = l$phenoData, description = l$description, : Could not open file /atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/c("/atlas/affy/proto_GQ/GQ00 5/M/20/MG_U74Av2/GQ005M02D20M1-R0.CEL", "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R1.CEL", "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03A20M1-R0.CEL", "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03C20M1-R0.CEL", "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03D20M1-R0.CEL", "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03E20M1-R0.CEL", "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03F20M1-R0.CEL", "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03G20M1-R0.CEL" in the last case <(cels1=celslist[[1]])> i ve this error: Error in just.rma(filenames = l$filenames, phenoData = l$phenoData, description = l$description, : Could not open file /atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2//atlas/affy/proto_GQ/GQ005/M /20/MG_U74Av2/GQ005M02D20M1-R0.CEL can someone tell me how to edit my CEL.list file in order to put 1 absolute path to each *.CEL and force R to normalize in the current directory *.CEL that are in other directories?
• 1.6k views
ADD COMMENT
0
Entering edit mode
@james-w-macdonald-5106
Last seen 2 days ago
United States
Hi Enrique, There is probably a way to use a file like that, but wouldn't it be easier to just pass the path to your call to justRMA()? eset <- justRMA(celfile.path = "/atlas/affy/proto_GQ/GQ005/M/20/") Best, Jim D.Enrique ESCOBAR ESPINOZA wrote: > i have decided to read cel files from a list containing full path to > each cel files, > one per line: > &$>head /atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/CEL.list > /atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R0.CEL > /atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R1.CEL > /atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03A20M1-R0.CEL > /atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03C20M1-R0.CEL > /atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03D20M1-R0.CEL > /atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03E20M1-R0.CEL > /atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03F20M1-R0.CEL > /atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03G20M1-R0.CEL > in R console i ve done: > >>(cels =celslist ); > > [[1]] > [1] "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R0.CEL" > [2] "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R1.CEL" > [3] "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03A20M1-R0.CEL" > [4] "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03C20M1-R0.CEL" > [5] "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03D20M1-R0.CEL" > [6] "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03E20M1-R0.CEL" > [7] "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03F20M1-R0.CEL" > [8] "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03G20M1-R0.CEL" > > >>(cels1=celslist[[1]]); > > [1] "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R0.CEL" > [2] "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R1.CEL" > [3] "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03A20M1-R0.CEL" > [4] "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03C20M1-R0.CEL" > [5] "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03D20M1-R0.CEL" > [6] "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03E20M1-R0.CEL" > [7] "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03F20M1-R0.CEL" > [8] "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03G20M1-R0.CEL" > in the fisrt caSE <(cels =celslist)> i have this error: > Error in just.rma(filenames = l$filenames, phenoData = l$phenoData, > description = l$description, : > Could not open file > /atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/c("/atlas/affy/proto_GQ/GQ 005/M/20/MG_U74Av2/GQ005M02D20M1-R0.CEL", > "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R1.CEL", > "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03A20M1-R0.CEL", > "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03C20M1-R0.CEL", > "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03D20M1-R0.CEL", > "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03E20M1-R0.CEL", > "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03F20M1-R0.CEL", > "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03G20M1-R0.CEL" > in the last case <(cels1=celslist[[1]])> i ve this error: > Error in just.rma(filenames = l$filenames, phenoData = l$phenoData, > description = l$description, : > Could not open file > /atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2//atlas/affy/proto_GQ/GQ005 /M/20/MG_U74Av2/GQ005M02D20M1-R0.CEL > > can someone tell me how to edit my CEL.list file in order to put 1 > absolute path to each *.CEL and force R to normalize in the current > directory *.CEL that are in other directories? > > _______________________________________________ > 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
0
Entering edit mode
Hi list, I wonder if anybody actually did get around to provide a cdf package for the Zea mays test chip (cornchip0)? Even though this package is superseded now, there are some publicly available experiments, using this platform. Last year there was a quick email exchange but I could never find the actual package. In the case that it was not released: Through some channels I was able to get the actual "cdf" file from Affymetrix and make a cdf package, but does someone have a step by step instructions how to make a windows binary (the links I found are mostly outdated: Brian Ripley's page points to another page...) or does anybody have everything necessary set up ;-) Kind regards, Bj?rn -+-+-+-+-+-+-+-+-+-+-+- Bj?rn Usadel, PhD Max Planck Institute of Molecular Plant Physiology System Regulation Group Am M?hlenberg 1 D-14476 Golm Germany Tel (+49 331) 567-8114 Email usadel at mpimp-golm.mpg.de WWW mapman.mpimp-golm.mpg.de
ADD REPLY
0
Entering edit mode
@james-w-macdonald-5106
Last seen 2 days ago
United States
Please don't reply just to me. The idea of a list is to keep the questions and answers available for others as well. I don't understand your response. How does this duplicate data? Also, what is your 'norm'? My understanding of your original question was that you wanted to operate on celfiles that are not in your working directory, which is what the celfile.path argument is for. Am I missing something? Best, Jim D.Enrique ESCOBAR ESPINOZA wrote: > i don t want to duplicate data, > my norm should be in another place than my cel files > > --- "James W. MacDonald" <jmacdon at="" med.umich.edu=""> wrote: > > >>Hi Enrique, >> >>There is probably a way to use a file like that, but wouldn't it be >> >>easier to just pass the path to your call to justRMA()? >> >>eset <- justRMA(celfile.path = "/atlas/affy/proto_GQ/GQ005/M/20/") >> >>Best, >> >>Jim >> >> >>D.Enrique ESCOBAR ESPINOZA wrote: >> >>>i have decided to read cel files from a list containing full path >> >>to >> >>>each cel files, >>>one per line: >>>&$>head /atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/CEL.list >>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R0.CEL >>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R1.CEL >>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03A20M1-R0.CEL >>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03C20M1-R0.CEL >>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03D20M1-R0.CEL >>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03E20M1-R0.CEL >>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03F20M1-R0.CEL >>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03G20M1-R0.CEL >>>in R console i ve done: >>> >>> >>>>(cels =celslist ); >>> >>>[[1]] >>>[1] >> >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R0.CEL" >> >>>[2] >> >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R1.CEL" >> >>>[3] >> >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03A20M1-R0.CEL" >> >>>[4] >> >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03C20M1-R0.CEL" >> >>>[5] >> >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03D20M1-R0.CEL" >> >>>[6] >> >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03E20M1-R0.CEL" >> >>>[7] >> >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03F20M1-R0.CEL" >> >>>[8] >> >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03G20M1-R0.CEL" >> >>> >>>>(cels1=celslist[[1]]); >>> >>>[1] >> >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R0.CEL" >> >>>[2] >> >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R1.CEL" >> >>>[3] >> >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03A20M1-R0.CEL" >> >>>[4] >> >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03C20M1-R0.CEL" >> >>>[5] >> >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03D20M1-R0.CEL" >> >>>[6] >> >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03E20M1-R0.CEL" >> >>>[7] >> >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03F20M1-R0.CEL" >> >>>[8] >> >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03G20M1-R0.CEL" >> >>>in the fisrt caSE <(cels =celslist)> i have this error: >>>Error in just.rma(filenames = l$filenames, phenoData = >> >>l$phenoData, >> >>>description = l$description, : >>> Could not open file >>> >> > /atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/c("/atlas/affy/proto_GQ/GQ 005/M/20/MG_U74Av2/GQ005M02D20M1-R0.CEL", > >>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R1.CEL", >>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03A20M1-R0.CEL", >>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03C20M1-R0.CEL", >>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03D20M1-R0.CEL", >>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03E20M1-R0.CEL", >>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03F20M1-R0.CEL", >>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03G20M1-R0.CEL" >>>in the last case <(cels1=celslist[[1]])> i ve this error: >>>Error in just.rma(filenames = l$filenames, phenoData = >> >>l$phenoData, >> >>>description = l$description, : >>> Could not open file >>> >> > /atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2//atlas/affy/proto_GQ/GQ005 /M/20/MG_U74Av2/GQ005M02D20M1-R0.CEL > >>>can someone tell me how to edit my CEL.list file in order to put >> >>1 >> >>>absolute path to each *.CEL and force R to normalize in the >> >>current >> >>>directory *.CEL that are in other directories? >>> >>>_______________________________________________ >>>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. >> > > > > -------------------------------------------------- > D.Enrique ESCOBAR ESPINOZA > (DESS Bioinfomatique > B.Sc. Biologie Moleculaire > Certificat Informatique) > http://adn.bioinfo.uqam.ca/~escd07097301/ > http://spaces.msn.com/members/escobarebio/ > ICQ#: 201778618 > ------------------------------------------------- > Tel: (514) 690-8398 > Montreal QC Canada > > __________________________________________________ > Do You Yahoo!? > Tired of spam? Yahoo! Mail has the best spam protection around > http://mail.yahoo.com -- 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
0
Entering edit mode
may i haven t explained the things correctly may be because my mother tongue is spanish & my second one is french... you had a good idea, but the cel files i want can be in diferent directories, maybe i want to normalize (nrom) 10 hgu133a celfiles that are in 10 different paths or directories... Thanks --- "James W. MacDonald" <jmacdon at="" med.umich.edu=""> wrote: > Please don't reply just to me. The idea of a list is to keep the > questions and answers available for others as well. > > I don't understand your response. How does this duplicate data? > Also, > what is your 'norm'? > > My understanding of your original question was that you wanted to > operate on celfiles that are not in your working directory, which > is > what the celfile.path argument is for. Am I missing something? > > Best, > > Jim > > D.Enrique ESCOBAR ESPINOZA wrote: > > i don t want to duplicate data, > > my norm should be in another place than my cel files > > > > --- "James W. MacDonald" <jmacdon at="" med.umich.edu=""> wrote: > > > > > >>Hi Enrique, > >> > >>There is probably a way to use a file like that, but wouldn't it > be > >> > >>easier to just pass the path to your call to justRMA()? > >> > >>eset <- justRMA(celfile.path = > "/atlas/affy/proto_GQ/GQ005/M/20/") > >> > >>Best, > >> > >>Jim > >> > >> > >>D.Enrique ESCOBAR ESPINOZA wrote: > >> > >>>i have decided to read cel files from a list containing full > path > >> > >>to > >> > >>>each cel files, > >>>one per line: > >>>&$>head /atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/CEL.list > >>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R0.CEL > >>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R1.CEL > >>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03A20M1-R0.CEL > >>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03C20M1-R0.CEL > >>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03D20M1-R0.CEL > >>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03E20M1-R0.CEL > >>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03F20M1-R0.CEL > >>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03G20M1-R0.CEL > >>>in R console i ve done: > >>> > >>> > >>>>(cels =celslist ); > >>> > >>>[[1]] > >>>[1] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R0.CEL" > >> > >>>[2] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R1.CEL" > >> > >>>[3] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03A20M1-R0.CEL" > >> > >>>[4] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03C20M1-R0.CEL" > >> > >>>[5] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03D20M1-R0.CEL" > >> > >>>[6] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03E20M1-R0.CEL" > >> > >>>[7] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03F20M1-R0.CEL" > >> > >>>[8] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03G20M1-R0.CEL" > >> > >>> > >>>>(cels1=celslist[[1]]); > >>> > >>>[1] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R0.CEL" > >> > >>>[2] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R1.CEL" > >> > >>>[3] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03A20M1-R0.CEL" > >> > >>>[4] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03C20M1-R0.CEL" > >> > >>>[5] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03D20M1-R0.CEL" > >> > >>>[6] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03E20M1-R0.CEL" > >> > >>>[7] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03F20M1-R0.CEL" > >> > >>>[8] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03G20M1-R0.CEL" > >> > >>>in the fisrt caSE <(cels =celslist)> i have this error: > >>>Error in just.rma(filenames = l$filenames, phenoData = > >> > >>l$phenoData, > >> > >>>description = l$description, : > >>> Could not open file > >>> > >> > > > /atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/c("/atlas/affy/proto_GQ/GQ00 5/M/20/MG_U74Av2/GQ005M02D20M1-R0.CEL", > > > >>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R1.CEL", > >>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03A20M1-R0.CEL", > >>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03C20M1-R0.CEL", > >>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03D20M1-R0.CEL", > >>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03E20M1-R0.CEL", > >>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03F20M1-R0.CEL", > >>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03G20M1-R0.CEL" > >>>in the last case <(cels1=celslist[[1]])> i ve this error: > >>>Error in just.rma(filenames = l$filenames, phenoData = > >> > >>l$phenoData, > >> > >>>description = l$description, : > >>> Could not open file > >>> > >> > > > /atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2//atlas/affy/proto_GQ/GQ005/M /20/MG_U74Av2/GQ005M02D20M1-R0.CEL > > > >>>can someone tell me how to edit my CEL.list file in order to put > >> > >>1 > >> > >>>absolute path to each *.CEL and force R to normalize in the > >> > >>current > >> > >>>directory *.CEL that are in other directories? > >>> > >>>_______________________________________________ > >>>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. > >> > > > > > > > > -------------------------------------------------- > > D.Enrique ESCOBAR ESPINOZA > > (DESS Bioinfomatique > > B.Sc. Biologie Moleculaire > === message truncated ===
ADD REPLY
0
Entering edit mode
If you give the full path as part of the filenames to justRMA then set celfile.path="" Best, Ben On Tue, 2006-08-22 at 07:22 -0700, D.Enrique ESCOBAR ESPINOZA wrote: > may i haven t explained the things correctly > may be because my mother tongue is spanish & my second one is > french... > you had a good idea, but the cel files i want can be in diferent > directories, > maybe i want to normalize (nrom) 10 hgu133a celfiles that are in 10 > different paths or directories... > > Thanks > --- "James W. MacDonald" <jmacdon at="" med.umich.edu=""> wrote: > > > Please don't reply just to me. The idea of a list is to keep the > > questions and answers available for others as well. > > > > I don't understand your response. How does this duplicate data? > > Also, > > what is your 'norm'? > > > > My understanding of your original question was that you wanted to > > operate on celfiles that are not in your working directory, which > > is > > what the celfile.path argument is for. Am I missing something? > > > > Best, > > > > Jim > > > > D.Enrique ESCOBAR ESPINOZA wrote: > > > i don t want to duplicate data, > > > my norm should be in another place than my cel files > > > > > > --- "James W. MacDonald" <jmacdon at="" med.umich.edu=""> wrote: > > > > > > > > >>Hi Enrique, > > >> > > >>There is probably a way to use a file like that, but wouldn't it > > be > > >> > > >>easier to just pass the path to your call to justRMA()? > > >> > > >>eset <- justRMA(celfile.path = > > "/atlas/affy/proto_GQ/GQ005/M/20/") > > >> > > >>Best, > > >> > > >>Jim > > >> > > >> > > >>D.Enrique ESCOBAR ESPINOZA wrote: > > >> > > >>>i have decided to read cel files from a list containing full > > path > > >> > > >>to > > >> > > >>>each cel files, > > >>>one per line: > > >>>&$>head /atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/CEL.list > > >>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R0.CEL > > >>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R1.CEL > > >>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03A20M1-R0.CEL > > >>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03C20M1-R0.CEL > > >>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03D20M1-R0.CEL > > >>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03E20M1-R0.CEL > > >>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03F20M1-R0.CEL > > >>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03G20M1-R0.CEL > > >>>in R console i ve done: > > >>> > > >>> > > >>>>(cels =celslist ); > > >>> > > >>>[[1]] > > >>>[1] > > >> > > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R0.CEL" > > >> > > >>>[2] > > >> > > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R1.CEL" > > >> > > >>>[3] > > >> > > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03A20M1-R0.CEL" > > >> > > >>>[4] > > >> > > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03C20M1-R0.CEL" > > >> > > >>>[5] > > >> > > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03D20M1-R0.CEL" > > >> > > >>>[6] > > >> > > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03E20M1-R0.CEL" > > >> > > >>>[7] > > >> > > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03F20M1-R0.CEL" > > >> > > >>>[8] > > >> > > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03G20M1-R0.CEL" > > >> > > >>> > > >>>>(cels1=celslist[[1]]); > > >>> > > >>>[1] > > >> > > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R0.CEL" > > >> > > >>>[2] > > >> > > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R1.CEL" > > >> > > >>>[3] > > >> > > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03A20M1-R0.CEL" > > >> > > >>>[4] > > >> > > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03C20M1-R0.CEL" > > >> > > >>>[5] > > >> > > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03D20M1-R0.CEL" > > >> > > >>>[6] > > >> > > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03E20M1-R0.CEL" > > >> > > >>>[7] > > >> > > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03F20M1-R0.CEL" > > >> > > >>>[8] > > >> > > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03G20M1-R0.CEL" > > >> > > >>>in the fisrt caSE <(cels =celslist)> i have this error: > > >>>Error in just.rma(filenames = l$filenames, phenoData = > > >> > > >>l$phenoData, > > >> > > >>>description = l$description, : > > >>> Could not open file > > >>> > > >> > > > > > > /atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/c("/atlas/affy/proto_GQ/GQ 005/M/20/MG_U74Av2/GQ005M02D20M1-R0.CEL", > > > > > > >>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R1.CEL", > > > >>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03A20M1-R0.CEL", > > > >>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03C20M1-R0.CEL", > > > >>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03D20M1-R0.CEL", > > > >>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03E20M1-R0.CEL", > > > >>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03F20M1-R0.CEL", > > >>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03G20M1-R0.CEL" > > >>>in the last case <(cels1=celslist[[1]])> i ve this error: > > >>>Error in just.rma(filenames = l$filenames, phenoData = > > >> > > >>l$phenoData, > > >> > > >>>description = l$description, : > > >>> Could not open file > > >>> > > >> > > > > > > /atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2//atlas/affy/proto_GQ/GQ005 /M/20/MG_U74Av2/GQ005M02D20M1-R0.CEL > > > > > >>>can someone tell me how to edit my CEL.list file in order to put > > >> > > >>1 > > >> > > >>>absolute path to each *.CEL and force R to normalize in the > > >> > > >>current > > >> > > >>>directory *.CEL that are in other directories? > > >>> > > >>>_______________________________________________ > > >>>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. > > >> > > > > > > > > > > > > -------------------------------------------------- > > > D.Enrique ESCOBAR ESPINOZA > > > (DESS Bioinfomatique > > > B.Sc. Biologie Moleculaire > > > === message truncated === > > _______________________________________________ > 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
ADD REPLY
0
Entering edit mode
You can specify several paths as such: cel.locations<-c("/directory1","/directory2") rma.Eset<-justRMA(celfile.path=cel.locations) Which saves you the trouble of itemizing the locations of particular cel files. Regards, Duane Hassane _________________________________________________ Duane C. Hassane, Ph.D. Center for Pediatric Biomedical Research University of Rochester Medical Center NOTICE OF CONFIDENTIALITY This message contains information which may be confidential and privileged. Unless you are the addressee (or authorized to receive for the addressee), you may not use, copy or disclose to anyone the message or any information contained in the message (including attachments). If you have received this message in error, please advise the sender by reply e-mail and delete the message. -----Original Message----- From: bioconductor-bounces@stat.math.ethz.ch [mailto:bioconductor-bounces at stat.math.ethz.ch]On Behalf Of D.Enrique ESCOBAR ESPINOZA Sent: Tuesday, August 22, 2006 9:22 AM To: James W. MacDonald; bioconductor at stat.math.ethz.ch Subject: Re: [BioC] CEL.list may i haven t explained the things correctly may be because my mother tongue is spanish & my second one is french... you had a good idea, but the cel files i want can be in diferent directories, maybe i want to normalize (nrom) 10 hgu133a celfiles that are in 10 different paths or directories... Thanks --- "James W. MacDonald" <jmacdon at="" med.umich.edu=""> wrote: > Please don't reply just to me. The idea of a list is to keep the > questions and answers available for others as well. > > I don't understand your response. How does this duplicate data? > Also, > what is your 'norm'? > > My understanding of your original question was that you wanted to > operate on celfiles that are not in your working directory, which > is > what the celfile.path argument is for. Am I missing something? > > Best, > > Jim > > D.Enrique ESCOBAR ESPINOZA wrote: > > i don t want to duplicate data, > > my norm should be in another place than my cel files > > > > --- "James W. MacDonald" <jmacdon at="" med.umich.edu=""> wrote: > > > > > >>Hi Enrique, > >> > >>There is probably a way to use a file like that, but wouldn't it > be > >> > >>easier to just pass the path to your call to justRMA()? > >> > >>eset <- justRMA(celfile.path = > "/atlas/affy/proto_GQ/GQ005/M/20/") > >> > >>Best, > >> > >>Jim > >> > >> > >>D.Enrique ESCOBAR ESPINOZA wrote: > >> > >>>i have decided to read cel files from a list containing full > path > >> > >>to > >> > >>>each cel files, > >>>one per line: > >>>&$>head /atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/CEL.list > >>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R0.CEL > >>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R1.CEL > >>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03A20M1-R0.CEL > >>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03C20M1-R0.CEL > >>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03D20M1-R0.CEL > >>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03E20M1-R0.CEL > >>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03F20M1-R0.CEL > >>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03G20M1-R0.CEL > >>>in R console i ve done: > >>> > >>> > >>>>(cels =celslist ); > >>> > >>>[[1]] > >>>[1] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R0.CEL" > >> > >>>[2] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R1.CEL" > >> > >>>[3] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03A20M1-R0.CEL" > >> > >>>[4] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03C20M1-R0.CEL" > >> > >>>[5] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03D20M1-R0.CEL" > >> > >>>[6] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03E20M1-R0.CEL" > >> > >>>[7] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03F20M1-R0.CEL" > >> > >>>[8] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03G20M1-R0.CEL" > >> > >>> > >>>>(cels1=celslist[[1]]); > >>> > >>>[1] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R0.CEL" > >> > >>>[2] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R1.CEL" > >> > >>>[3] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03A20M1-R0.CEL" > >> > >>>[4] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03C20M1-R0.CEL" > >> > >>>[5] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03D20M1-R0.CEL" > >> > >>>[6] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03E20M1-R0.CEL" > >> > >>>[7] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03F20M1-R0.CEL" > >> > >>>[8] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03G20M1-R0.CEL" > >> > >>>in the fisrt caSE <(cels =celslist)> i have this error: > >>>Error in just.rma(filenames = l$filenames, phenoData = > >> > >>l$phenoData, > >> > >>>description = l$description, : > >>> Could not open file > >>> > >> > > > /atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/c("/atlas/affy/proto_GQ/GQ00 5/M/20/MG_U74Av2/GQ005M02D20M1-R0.CEL", > > > >>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R1.CEL", > >>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03A20M1-R0.CEL", > >>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03C20M1-R0.CEL", > >>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03D20M1-R0.CEL", > >>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03E20M1-R0.CEL", > >>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03F20M1-R0.CEL", > >>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03G20M1-R0.CEL" > >>>in the last case <(cels1=celslist[[1]])> i ve this error: > >>>Error in just.rma(filenames = l$filenames, phenoData = > >> > >>l$phenoData, > >> > >>>description = l$description, : > >>> Could not open file > >>> > >> > > > /atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2//atlas/affy/proto_GQ/GQ005/M /20/MG_U74Av2/GQ005M02D20M1-R0.CEL > > > >>>can someone tell me how to edit my CEL.list file in order to put > >> > >>1 > >> > >>>absolute path to each *.CEL and force R to normalize in the > >> > >>current > >> > >>>directory *.CEL that are in other directories? > >>> > >>>_______________________________________________ > >>>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. > >> > > > > > > > > -------------------------------------------------- > > D.Enrique ESCOBAR ESPINOZA > > (DESS Bioinfomatique > > B.Sc. Biologie Moleculaire > === message truncated === _______________________________________________ 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 -----Original Message----- From: bioconductor-bounces@stat.math.ethz.ch [mailto:bioconductor-bounces at stat.math.ethz.ch]On Behalf Of D.Enrique ESCOBAR ESPINOZA Sent: Tuesday, August 22, 2006 9:22 AM To: James W. MacDonald; bioconductor at stat.math.ethz.ch Subject: Re: [BioC] CEL.list may i haven t explained the things correctly may be because my mother tongue is spanish & my second one is french... you had a good idea, but the cel files i want can be in diferent directories, maybe i want to normalize (nrom) 10 hgu133a celfiles that are in 10 different paths or directories... Thanks --- "James W. MacDonald" <jmacdon at="" med.umich.edu=""> wrote: > Please don't reply just to me. The idea of a list is to keep the > questions and answers available for others as well. > > I don't understand your response. How does this duplicate data? > Also, > what is your 'norm'? > > My understanding of your original question was that you wanted to > operate on celfiles that are not in your working directory, which > is > what the celfile.path argument is for. Am I missing something? > > Best, > > Jim > > D.Enrique ESCOBAR ESPINOZA wrote: > > i don t want to duplicate data, > > my norm should be in another place than my cel files > > > > --- "James W. MacDonald" <jmacdon at="" med.umich.edu=""> wrote: > > > > > >>Hi Enrique, > >> > >>There is probably a way to use a file like that, but wouldn't it > be > >> > >>easier to just pass the path to your call to justRMA()? > >> > >>eset <- justRMA(celfile.path = > "/atlas/affy/proto_GQ/GQ005/M/20/") > >> > >>Best, > >> > >>Jim > >> > >> > >>D.Enrique ESCOBAR ESPINOZA wrote: > >> > >>>i have decided to read cel files from a list containing full > path > >> > >>to > >> > >>>each cel files, > >>>one per line: > >>>&$>head /atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/CEL.list > >>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R0.CEL > >>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R1.CEL > >>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03A20M1-R0.CEL > >>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03C20M1-R0.CEL > >>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03D20M1-R0.CEL > >>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03E20M1-R0.CEL > >>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03F20M1-R0.CEL > >>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03G20M1-R0.CEL > >>>in R console i ve done: > >>> > >>> > >>>>(cels =celslist ); > >>> > >>>[[1]] > >>>[1] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R0.CEL" > >> > >>>[2] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R1.CEL" > >> > >>>[3] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03A20M1-R0.CEL" > >> > >>>[4] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03C20M1-R0.CEL" > >> > >>>[5] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03D20M1-R0.CEL" > >> > >>>[6] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03E20M1-R0.CEL" > >> > >>>[7] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03F20M1-R0.CEL" > >> > >>>[8] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03G20M1-R0.CEL" > >> > >>> > >>>>(cels1=celslist[[1]]); > >>> > >>>[1] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R0.CEL" > >> > >>>[2] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R1.CEL" > >> > >>>[3] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03A20M1-R0.CEL" > >> > >>>[4] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03C20M1-R0.CEL" > >> > >>>[5] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03D20M1-R0.CEL" > >> > >>>[6] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03E20M1-R0.CEL" > >> > >>>[7] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03F20M1-R0.CEL" > >> > >>>[8] > >> > >>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03G20M1-R0.CEL" > >> > >>>in the fisrt caSE <(cels =celslist)> i have this error: > >>>Error in just.rma(filenames = l$filenames, phenoData = > >> > >>l$phenoData, > >> > >>>description = l$description, : > >>> Could not open file > >>> > >> > > > /atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/c("/atlas/affy/proto_GQ/GQ00 5/M/20/MG_U74Av2/GQ005M02D20M1-R0.CEL", > > > >>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R1.CEL", > >>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03A20M1-R0.CEL", > >>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03C20M1-R0.CEL", > >>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03D20M1-R0.CEL", > >>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03E20M1-R0.CEL", > >>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03F20M1-R0.CEL", > >>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03G20M1-R0.CEL" > >>>in the last case <(cels1=celslist[[1]])> i ve this error: > >>>Error in just.rma(filenames = l$filenames, phenoData = > >> > >>l$phenoData, > >> > >>>description = l$description, : > >>> Could not open file > >>> > >> > > > /atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2//atlas/affy/proto_GQ/GQ005/M /20/MG_U74Av2/GQ005M02D20M1-R0.CEL > > > >>>can someone tell me how to edit my CEL.list file in order to put > >> > >>1 > >> > >>>absolute path to each *.CEL and force R to normalize in the > >> > >>current > >> > >>>directory *.CEL that are in other directories? > >>> > >>>_______________________________________________ > >>>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. > >> > > > > > > > > -------------------------------------------------- > > D.Enrique ESCOBAR ESPINOZA > > (DESS Bioinfomatique > > B.Sc. Biologie Moleculaire > === message truncated === _______________________________________________ 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
ADD REPLY
0
Entering edit mode
Ben may have a trick that will work, but I don't know of any way to get read.affybatch() to use an arbitrary set of directories for the celfiles. You could probably hack the code to do something like that pretty easily if you really wanted to. Best, Jim D.Enrique ESCOBAR ESPINOZA wrote: > may i haven t explained the things correctly > may be because my mother tongue is spanish & my second one is > french... > you had a good idea, but the cel files i want can be in diferent > directories, > maybe i want to normalize (nrom) 10 hgu133a celfiles that are in 10 > different paths or directories... > > Thanks > --- "James W. MacDonald" <jmacdon at="" med.umich.edu=""> wrote: > > >>Please don't reply just to me. The idea of a list is to keep the >>questions and answers available for others as well. >> >>I don't understand your response. How does this duplicate data? >>Also, >>what is your 'norm'? >> >>My understanding of your original question was that you wanted to >>operate on celfiles that are not in your working directory, which >>is >>what the celfile.path argument is for. Am I missing something? >> >>Best, >> >>Jim >> >>D.Enrique ESCOBAR ESPINOZA wrote: >> >>>i don t want to duplicate data, >>>my norm should be in another place than my cel files >>> >>>--- "James W. MacDonald" <jmacdon at="" med.umich.edu=""> wrote: >>> >>> >>> >>>>Hi Enrique, >>>> >>>>There is probably a way to use a file like that, but wouldn't it >> >>be >> >>>>easier to just pass the path to your call to justRMA()? >>>> >>>>eset <- justRMA(celfile.path = >> >>"/atlas/affy/proto_GQ/GQ005/M/20/") >> >>>>Best, >>>> >>>>Jim >>>> >>>> >>>>D.Enrique ESCOBAR ESPINOZA wrote: >>>> >>>> >>>>>i have decided to read cel files from a list containing full >> >>path >> >>>>to >>>> >>>> >>>>>each cel files, >>>>>one per line: >>>>>&$>head /atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/CEL.list >>>>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R0.CEL >>>>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R1.CEL >>>>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03A20M1-R0.CEL >>>>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03C20M1-R0.CEL >>>>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03D20M1-R0.CEL >>>>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03E20M1-R0.CEL >>>>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03F20M1-R0.CEL >>>>>/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03G20M1-R0.CEL >>>>>in R console i ve done: >>>>> >>>>> >>>>> >>>>>>(cels =celslist ); >>>>> >>>>>[[1]] >>>>>[1] >>>> >>>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R0.CEL" >>>> >>>> >>>>>[2] >>>> >>>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R1.CEL" >>>> >>>> >>>>>[3] >>>> >>>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03A20M1-R0.CEL" >>>> >>>> >>>>>[4] >>>> >>>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03C20M1-R0.CEL" >>>> >>>> >>>>>[5] >>>> >>>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03D20M1-R0.CEL" >>>> >>>> >>>>>[6] >>>> >>>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03E20M1-R0.CEL" >>>> >>>> >>>>>[7] >>>> >>>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03F20M1-R0.CEL" >>>> >>>> >>>>>[8] >>>> >>>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03G20M1-R0.CEL" >>>> >>>> >>>>>>(cels1=celslist[[1]]); >>>>> >>>>>[1] >>>> >>>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R0.CEL" >>>> >>>> >>>>>[2] >>>> >>>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R1.CEL" >>>> >>>> >>>>>[3] >>>> >>>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03A20M1-R0.CEL" >>>> >>>> >>>>>[4] >>>> >>>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03C20M1-R0.CEL" >>>> >>>> >>>>>[5] >>>> >>>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03D20M1-R0.CEL" >>>> >>>> >>>>>[6] >>>> >>>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03E20M1-R0.CEL" >>>> >>>> >>>>>[7] >>>> >>>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03F20M1-R0.CEL" >>>> >>>> >>>>>[8] >>>> >>>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03G20M1-R0.CEL" >>>> >>>> >>>>>in the fisrt caSE <(cels =celslist)> i have this error: >>>>>Error in just.rma(filenames = l$filenames, phenoData = >>>> >>>>l$phenoData, >>>> >>>> >>>>>description = l$description, : >>>>> Could not open file >>>>> >>>> > /atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/c("/atlas/affy/proto_GQ/GQ 005/M/20/MG_U74Av2/GQ005M02D20M1-R0.CEL", > >>>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M02D20M1-R1.CEL", >> >>>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03A20M1-R0.CEL", >> >>>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03C20M1-R0.CEL", >> >>>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03D20M1-R0.CEL", >> >>>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03E20M1-R0.CEL", >> >>>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03F20M1-R0.CEL", >>>> >>>>>"/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03G20M1-R0.CEL" >>>>>in the last case <(cels1=celslist[[1]])> i ve this error: >>>>>Error in just.rma(filenames = l$filenames, phenoData = >>>> >>>>l$phenoData, >>>> >>>> >>>>>description = l$description, : >>>>> Could not open file >>>>> >>>> > /atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2//atlas/affy/proto_GQ/GQ005 /M/20/MG_U74Av2/GQ005M02D20M1-R0.CEL > >>>>>can someone tell me how to edit my CEL.list file in order to put >>>> >>>>1 >>>> >>>> >>>>>absolute path to each *.CEL and force R to normalize in the >>>> >>>>current >>>> >>>> >>>>>directory *.CEL that are in other directories? >>>>> >>>>>_______________________________________________ >>>>>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. >>>> >>> >>> >>> >>>-------------------------------------------------- >>>D.Enrique ESCOBAR ESPINOZA >>>(DESS Bioinfomatique >>>B.Sc. Biologie Moleculaire >> > === message truncated === > > > __________________________________________________ > Do You Yahoo!? > Tired of spam? Yahoo! Mail has the best spam protection around > http://mail.yahoo.com -- 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 REPLY
0
Entering edit mode
In general you should not be calling read.affybatch directly. Instead users should be calling ReadAffy(). If you use ReadAffy() the previously mentioned trick eg cel.locations<-c("/directory1","/directory2") my.abatch <-ReadAffy(celfile.path=cel.locations) will work fine. On Tue, 2006-08-22 at 11:26 -0400, James W. MacDonald wrote: > Ben may have a trick that will work, but I don't know of any way to get > read.affybatch() to use an arbitrary set of directories for the > celfiles. You could probably hack the code to do something like that > pretty easily if you really wanted to. > > Best, > > Jim > -- Ben Bolstad <bmb at="" bmbolstad.com=""> http://bmbolstad.com
ADD REPLY
0
Entering edit mode
I guess there is always the possibility to wrap a function with an other function with a different signature (or a processing of the parameters) that looks more convenient, but using 'list.files' seems to be doing the trick, and would introduce a nifty function that can used in other contexts E.g: read.affybatch(filenames=list.files(c("directory1", "directory2"))) L. > In general you should not be calling read.affybatch directly. Instead > users should be calling ReadAffy(). > > If you use ReadAffy() the previously mentioned trick eg > > cel.locations<-c("/directory1","/directory2") > my.abatch <-ReadAffy(celfile.path=cel.locations) > > will work fine. > > > > On Tue, 2006-08-22 at 11:26 -0400, James W. MacDonald wrote: >> Ben may have a trick that will work, but I don't know of any way to get >> read.affybatch() to use an arbitrary set of directories for the >> celfiles. You could probably hack the code to do something like that >> pretty easily if you really wanted to. >> >> Best, >> >> Jim >> > > > -- > Ben Bolstad <bmb at="" bmbolstad.com=""> > http://bmbolstad.com > > _______________________________________________ > 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 >
ADD REPLY

Login before adding your answer.

Traffic: 606 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