[BioC] Combine function in matchprobes
w.huber at dkfz-heidelberg.de
w.huber at dkfz-heidelberg.de
Sat Jan 17 11:23:07 MET 2004
Hi Min-Han Tan,
You specified CDF packages instead of probe packages. The combine
function will produce a more intelligible error message rather than the
one you reported in version 1.0.3 which will be on the web shortly.
Best wishes
Wolfgang
-------------------------------------
Wolfgang Huber
Division of Molecular Genome Analysis
German Cancer Research Center
Heidelberg, Germany
Phone: +49 6221 424709
Fax: +49 6221 42524709
Http: www.dkfz.de/abt0840/whuber
-------------------------------------
On Fri, 16 Jan 2004, Robert Gentleman wrote:
> Well, the man page for combine says:
> Usage:
>
> combine(batch, probepkg, newcdf, verbose=TRUE)
>
> Arguments:
>
> batch: A 'list' of AffyBatches.
>
> probepkg: A character vector of the same length as 'batch', containing
> the names of the probe sequences libraries (for example,
> 'hgu133aprobe') that go with the AffyBatches.
>
> newcdf: Character. The name of the name CDF environment that is to
> be
> created.
>
> verbose: Logical. If TRUE, messages are printed to the console.
>
>
> And I do not see where you have indicated the probe packages?
> You do need them (and they are available on the Bioconductor site).
>
> Robert
>
>
> On Fri, Jan 16, 2004 at 03:03:27PM -0500, Tan, MinHan wrote:
> > Hi,
> >
> > I tried combining 2 Affybatches (HGU133A and HGU133Plus 2.0) using
> > combine in matchprobes, but ran into the following problems.
> >
> > I used both R 1.9.0 and R 1.8.1 with the latest downloads from
> > Bioconductor.
> >
> > > merged.data <- combine(list(DataHGU133A,DataHGU133P), c("hgu133acdf",
> > "hgu133plus2cdf"), "mycdf", verbose=TRUE)
> > package:hgu133acdf hgu133acdf
> > package:hgu133plus2cdf hgu133plus2cdf
> > 0 probes in common between hgu133acdf and hgu133acdf
> > Error in sequs[[i]] : subscript out of bounds
> >
> > Looking at the combine function, I seem to be running into problems at
> > this point?
> >
> > > tmp <- get(upp[1], uppname[1])
> > > str(tmp)
> > length 0 <environment>
> > > upp[1]
> > [1] "hgu133acdf"
> > > uppname[1]
> > [1] "package:hgu133acdf"
> >
> >
> > Your advice is very much appreciated!
> >
> >
> > Min-Han Tan
> >
> > This email message, including any attachments, is for the so...{{dropped}}
> >
> > _______________________________________________
> > Bioconductor mailing list
> > Bioconductor at stat.math.ethz.ch
> > https://www.stat.math.ethz.ch/mailman/listinfo/bioconductor
>
> --
> +---------------------------------------------------------------------------+
> | Robert Gentleman phone : (617) 632-5250 |
> | Associate Professor fax: (617) 632-2444 |
> | Department of Biostatistics office: M1B20 |
> | Harvard School of Public Health email: rgentlem at jimmy.harvard.edu |
> +---------------------------------------------------------------------------+
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
> https://www.stat.math.ethz.ch/mailman/listinfo/bioconductor
>
More information about the Bioconductor
mailing list