[BioC] RE: Specifying Sample Names in ReadAffy (Correction)
Paul Boutros
Paul.Boutros at utoronto.ca
Fri Jul 4 23:15:04 MEST 2003
Rafael pointed out that there was a typo in my example. The correct version
is below:
> samples.files <- c(
"RAE230A_043003_IM01T_LH.CEL",
"RAE230A_043003_IM02T_LH.CEL",
"RAE230A_043003_IM03T_LH.CEL",
"RAE230A_043003_IM04T_LH.CEL",
);
> samples.names <- c(
"LET",
"LEC",
"WWT",
"WWC",
);
> Data <- ReadAffy(
filenames = samples.files,
sampleNames = samples.names
);
> -----Original Message-----
> From: Paul Boutros [mailto:pcboutro at engmail.uwaterloo.ca]
> Sent: Friday, July 04, 2003 5:37 PM
> To: bioconductor at stat.math.ethz.ch
> Subject: Specifying Sample Names in ReadAffy;
>
>
> Hello all,
>
> I would like to be able to specify the sampleNames parameter to
> ReadAffy to
> ease later parsing and analysis of the data somewhat. The documentation
> (Textual description of affy vignete and ?ReadAffy) indicates
> that I need to
> set up a character array with both the file-names and the desired sample-
> names. This is what I have used (limited data-set for testing):
>
> > samples.files <- c(
> "RAE230A_043003_IM01T_LH.CEL",
> "RAE230A_043003_IM02T_LH.CEL",
> "RAE230A_043003_IM03T_LH.CEL",
> "RAE230A_043003_IM04T_LH.CEL",
> );
>
> > samples.names <- c(
> "LET",
> "LEC",
> "WWT",
> "WWC",
> );
>
> > Data <- ReadAffy(
> filenames = samples.files);
> sampleNames = samples.animals,
> );
>
> But this fails with the error:
>
> Error in sub(pattern, replacement, x, ignore.case, extended) :
> invalid argument
> In addition: Warning message:
> sampleNames not same length as filenames. Using filenames as
> sampleNames instead
> in: ReadAffy(filenames = samples.files, sampleNames = samples.animals,
>
> So I tested:
> > length(samples.files);
> [1] 4
> > length(samples.names);
> [1] 4
>
> Have I misunderstood the purpose of the sampleNames argument? Or
> ss there
> another way of changing the "column" names in an eSet and in
> write.exprs()
> output?
>
> Any suggestions much appreciated!
> Paul
>
> ----------------------------------------
> This mail sent through www.mywaterloo.ca
More information about the Bioconductor
mailing list