[BioC] help w/ bg.correct.gcrma
Zhijin Wu
zwu at jhsph.edu
Sat Nov 1 05:54:57 MET 2003
Thanks, James. I will update the "return" line in gcrma to be
compatible with R-1.8.
Jean
----- Original Message -----
From: James MacDonald <jmacdon at med.umich.edu>
Date: Friday, October 31, 2003 9:17 pm
Subject: Re: [BioC] help w/ bg.correct.gcrma
> The first try should read like this:
>
> bgaffy
> <-
> bg.correct.gcrma(rawaffy,gcgroup=getGroupInfo
(rawaffy),estimate="eb",triple.goal=TRUE)
>
> Note the error in your message. You can also set
> triple.goal=FALSE, but
> you have to do one or the other.
>
> As for the multi-argument returns message, that is simply a change
> as of
> R-1.8.0 that really only concerns the package maintainer. For now the
> functions will work fine (except for the annoying warnings).
>
> HTH,
>
> Jim
>
>
>
> James W. MacDonald
> Affymetrix and cDNA Microarray Core
> University of Michigan Cancer Center
> 1500 E. Medical Center Drive
> 7410 CCGC
> Ann Arbor MI 48109
> 734-647-5623
> >>> Jenny Drnevich <drnevich at uiuc.edu> 10/31/03 16:29 PM >>>
> Hi,
>
> I want to explore the GC-based background correction in gcrma. I'd
> liketo
> get the bg corrected PM values, not the summary values. I think I
> can do
>
> this with bg.correct.gcrma, but I'm having problems getting it to
> run
> correctly. I'm looking at some of the Affy spike-in data, and
> running R
> 1.8.0 and Bioconductor 1.3. I followed the format I found from
> ?bg.correct.gcrma:
>
> > bgaffy
> <-
> bg.correct.gcrma(rawaffy,gcgroup=getGroupInfo(rawaffy),estimate="eb")
Loading required package: hgu95aprobe
> Loading required package: matchprobes
> Error in FUN(array(newX[, i], d.call, dn.call), ...) :
> Argument "triple.goal" is missing, with no default
> In addition: Warning messages:
> 1: multi-argument returns are deprecated in: return(y = yhat, wt)
> 2: multi-argument returns are deprecated in: return(y = yhat, wt)
> 3: multi-argument returns are deprecated in: return(y = yhat, wt)
> 4: multi-argument returns are deprecated in: return(y = yhat, wt)
> 5: multi-argument returns are deprecated in: return(y = yhat, wt)
>
> Additionally, I tried:
>
> bgaffy
> <-
> bg.correct.gcrma(rawaffy,gcgroup=getGroupInfo(rawaffy),estimate="mle")
> and got over 50 of the "multi-argument returns" warning messages. Any
> ideas
> about what I need to do?
>
> Thanks,
> Jenny
>
>
> Jenny Drnevich, Ph.D.
> Department of Animal Biology
> 515 Morrill Hall
> 505 S Goodwin Ave
> Urbana, IL 61801
> USA
>
> ph: 217-244-6826
> fax: 217-244-4565
> e-mail: drnevich at uiuc.edu
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
> https://www.stat.math.ethz.ch/mailman/listinfo/bioconductor
>
> _______________________________________________
> 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