[BioC] "normalize.quantiles.robust" produces the SAME intensities from different weights?!
k. brand
k.brand at erasmusmc.nl
Thu Sep 21 19:35:46 CEST 2006
Ben, thank you for the fast help- it worked!
Karl
on 9/21/2006 5:16 AM Ben Bolstad said the following:
> Thanks for this report. It turned out that there was a minor bug in
> normalize.quantiles.robust() and incorrect parameters were being passed
> down to the C code. This is due to some changes I made many months ago
> now.
>
> A corrected version of the function will appear in the next version of
> affy. But briefly the final line in that function should have been:
>
>
> .Call("R_qnorm_robust_c",x,copy,weights,as.integer(use.median),as.integer(use.log2),as.integer(use.huber),PACKAGE="affy")
>
> Best,
>
> Ben
>
> On Wed, 2006-09-20 at 17:29 +0200, k. brand wrote:
>> Dear BioCers,
>>
>> My interpretation of the function "normalize.quantiles.robust" is the
>> ability to selectively increase or reduce the contribution of specific
>> chips during normalization.
>>
>> Although i get different results from equally weighted chips (which is
>> the same as "justRMA" as a verification), all variations of unequal
>> weightings produce do the SAME intensities?! See 2 examples below.
>>
>> Can some one find the error in my script, or better yet, give an example
>> of a better way to achieve my objective-
>>
>> to selectively 'pull' low intensities of duplicate chips form a poor hyb
>> (on expired arrays), up to the levels of duplicate chips of a good
>> hyb(on unexpired arrays).
>>
>> TIA,
>>
>> Karl
>>
>>
>
>
--
Karl Brand <k.brand at erasmusmc.nl>
Department of Genetics
Erasmus MC
Dr Molewaterplein 50
3015 GE Rotterdam
lab +31 (0)10 408 7409 fax +31 (0)10 408 9468
More information about the Bioconductor
mailing list