[BioC] normalize.quantiles.robust: how to implement weights for different arrays?

James W. MacDonald jmacdon at med.umich.edu
Mon Sep 18 20:39:05 CEST 2006


k. brand wrote:
> Dear BioCers,
> 
> Regarding-
> 
> normalize.quantiles.robust(x,weights=NULL,
> 
> Can some one illustrate how to implement weights for my different arrays?
> 
> Im trying to 'up weight' a good hyb against a not so good hyb, in an 
> otherwise standard RMA approach. Trying to follow the package 
> description i failed  to effect this (see failed attempt below). Any 
> suggestions greatly appreciated.
> 
> thanks in advance,
> 
> Karl
> 
> 
>  > library(affyPLM)
> Loading required package: affydata
> Loading required package: gcrma
> Loading required package: matchprobes
>  > dat <- ReadAffy()
>  > datbg <- bg.correct.rma(dat)
>  >normalize.quantiles.robust(datbg,weights="Tco1A.CEL"=1,"Tmi1A.CEL"=1,"Tsh1A.CEL"=1,"Tco2A.CEL"=1,

 From the help for normalize.quantiles.robust()

weights: A vector of weights, one for each chip


What you have passed is not a vector of weights. As far as R is 
concerned, you have passed a character string "Tco1A.CEL", followed by a 
bunch of other arguments that don't match up to the available arguments 
in the function call. You need something like

weights = rep(1:2, each = 6)

Which will be a vector of weights.

You might consider perusing "An Introduction to R".

Best,

Jim


> Error: syntax error in 
> normalize.quantiles.robust(datbg,weights="Tco1A.CEL"="
>  > "Tmi2A.CEL"=1,"Tsh2A.CEL"=1,"Tco3B.CEL"=2,"Tmi3B.CEL"=2,
> Error: syntax error in ""Tmi2A.CEL"=1,"
>  > "Tsh3B.CEL"=2,"Tco4B.CEL"=2,"Tmi4B.CEL"=2,"Tsh4B.CEL"=2)
> Error: syntax error in ""Tsh3B.CEL"=2,"
> 
> 
>  > sessionInfo()
> Version 2.3.0 (2006-04-24)
> i386-pc-mingw32
> 
> attached base packages:
> [1] "tools"     "methods"   "stats"     "graphics"  "grDevices" "utils"
> [7] "datasets"  "base"
> 
> other attached packages:
> mouse4302cdf      affyPLM        gcrma  matchprobes     affydata 
>   affy
>      "1.12.0"      "1.8.0"      "2.4.1"      "1.4.0"      "1.8.0" 
> "1.10.0"
>        affyio      Biobase
>       "1.0.0"     "1.10.0"
> 
> 


-- 
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.



More information about the Bioconductor mailing list