[BioC] CEL.list vol 2

D.Enrique ESCOBAR ESPINOZA escobarebio at yahoo.com
Tue Aug 22 18:02:47 CEST 2006


lets say that i want to work with *.CEls files that are par of an
experiment:
>getwd()
>/atlas/affy/experiments/telomerase/

lets say that i m gathering .CELs in differents paths like this:
/atlas/affy/proto_GQ/GQ005/M/14/MG_U74Av2/GQ005M03C14M1-R0.CEL
/atlas/affy/proto_GQ/GQ005/M/16/MG_U74Av2/GQ005M03D16M1-R0.CEL
/atlas/affy/proto_GQ/GQ005/M/18/MG_U74Av2/GQ005M03E18M1-R0.CEL
/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03F20M1-R0.CEL
/atlas/affy/proto_GQ/GQ006/F/14/MG_U74Av2/GQ006M03C14F1-R0.CEL
/atlas/affy/proto_GQ/GQ006/F/16/MG_U74Av2/GQ006M03D16F1-R0.CEL
/atlas/affy/proto_GQ/GQ006/F/18/MG_U74Av2/GQ006M03E18F1-R0.CEL
/atlas/affy/proto_GQ/GQ006/F/20/MG_U74Av2/GQ006M03F20F1-R0.CEL

and i decide to put all the paths in a CEL.list
in the curretn directory:
/atlas/affy/experiments/telomerase/
and its absolute path is 
/atlas/affy/experiments/telomerase/CEL.list

I want to save my Rdatas and all analysis in 
/atlas/affy/experiments/telomerase/

> # load my cels
> celslist<-
scan(file="CEL.list",what=list("c"),sep="\n",quiet=TRUE);
> (cels =celslist     );
[[1]]
[1] "/atlas/affy/proto_GQ/GQ005/M/14/MG_U74Av2/GQ005M03C14M1-R0.CEL"
[2] "/atlas/affy/proto_GQ/GQ005/M/16/MG_U74Av2/GQ005M03D16M1-R0.CEL"
[3] "/atlas/affy/proto_GQ/GQ005/M/18/MG_U74Av2/GQ005M03E18M1-R0.CEL"
[4] "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03F20M1-R0.CEL"
[5] "/atlas/affy/proto_GQ/GQ006/F/14/MG_U74Av2/GQ006M03C14F1-R0.CEL"
[6] "/atlas/affy/proto_GQ/GQ006/F/16/MG_U74Av2/GQ006M03D16F1-R0.CEL"
[7] "/atlas/affy/proto_GQ/GQ006/F/18/MG_U74Av2/GQ006M03E18F1-R0.CEL"
[8] "/atlas/affy/proto_GQ/GQ006/F/20/MG_U74Av2/GQ006M03F20F1-R0.CEL"
> (cels1=celslist[[1]]);
[1] "/atlas/affy/proto_GQ/GQ005/M/14/MG_U74Av2/GQ005M03C14M1-R0.CEL"
[2] "/atlas/affy/proto_GQ/GQ005/M/16/MG_U74Av2/GQ005M03D16M1-R0.CEL"
[3] "/atlas/affy/proto_GQ/GQ005/M/18/MG_U74Av2/GQ005M03E18M1-R0.CEL"
[4] "/atlas/affy/proto_GQ/GQ006/F/14/MG_U74Av2/GQ006M03C14F1-R0.CEL"
[5] "/atlas/affy/proto_GQ/GQ005/M/20/MG_U74Av2/GQ005M03D20M1-R0.CEL"
[6] "/atlas/affy/proto_GQ/GQ006/F/16/MG_U74Av2/GQ006M03D16F1-R0.CEL"
[7] "/atlas/affy/proto_GQ/GQ006/F/18/MG_U74Av2/GQ006M03E18F1-R0.CEL"
[8] "/atlas/affy/proto_GQ/GQ006/F/20/MG_U74Av2/GQ006M03F20F1-R0.CEL"

when i do
eset  <-just.rma(filenames = cels , phenoData =NULL)
eset1 <-just.rma(filenames = cels1, phenoData =NULL)

i always receive
>Error in just.rma(filenames = l$filenames, phenoData =
>l$phenoData,
>description = l$description,  :
>        Could not open file


what can i do for putting in filenames
the correct different path foreach different .CEL?



More information about the Bioconductor mailing list