[BioC] buglet (?) snprma (oligo)

[Tobias Verbeke]

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