[BioC] buglet (?) snprma (oligo)
Benilton Carvalho
bcarvalh at jhsph.edu
Thu Sep 6 00:37:29 CEST 2007
Hi Tobias,
sorry for the late reply.
I've just committed a change which uses
preprocessCore:::normalize.quantiles.use.target instead of the code
mentioned below.
thanks,
b
On Sep 2, 2007, at 3:18 AM, Tobias Verbeke wrote:
> Dear list,
>
> I discovered an oddity in snprma (oligo,
> svn revision 26822). It appears that ranks
> are used to extract elements from vectors in
>
> for (i in 1:ncol(tmpExprs))
> tmpExprs[, i] <- reference[rank(tmpExprs[, i])]
>
>
> This might break iff there are two values that
> are exactly equal as can be seen from
>
> bar <- 2 * c(2, 4, 4, 5, 1)
> rank(bar)
> foo <- 1:5
> foo[bar]
>
> A solution would be to specify an alternative
> ties.method such as
>
> rank(bar, ties.method = "first")
>
> Any thoughts?
>
> Regards,
> Tobias
>
> --
>
> Tobias Verbeke - Consultant
> Business & Decision Benelux
> Rue de la révolution 8
> 1000 Brussels - BELGIUM
>
> +32 499 36 33 15
> tobias.verbeke at businessdecision.com
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives: http://news.gmane.org/
> gmane.science.biology.informatics.conductor
More information about the Bioconductor
mailing list