[BioC] limma Normalization question

Gordon Smyth smyth at wehi.edu.au
Tue Nov 22 02:46:08 CET 2005


>[BioC] limma Normalization question
>Cecilia McGregor cmcgre1 at lsu.edu
>Tue Nov 22 00:54:53 CET 2005
>
>Dr Smyth
>
>I have some control spots that are not differentially expressed on the array.
>I tried to use the modifyWeights as on page 21 of the Users Guide (6 Oct 
>2005).
>However I realized that it doesn't work, cause I don't use weights in my
>read.maimages. I used a image analysis program, not recognized by limma, 
>so used
>the following to read the data.
>
>RG <- read.maimages(targets$FileName, 
>columns=list(R="Rf",G="Gf",Rb="Rb",Gb="Gb"))
>
>I tried to add "phony" column for weights (all values=1) in the intensity 
>files,
>so I can modify it, but have not been able to get it to work (yet).
>
>Cecilia

Dear Cecilia,

There's no need to change the data files (a last resort!), you could just 
use for example

RG$weights <- array(1, dim(RG))

Best wishes
Gordon



More information about the Bioconductor mailing list