binary distances R ( tanimoto)

binaryDistances.R

# http://www.sequentix.de/gelquest/help/distance_measures.htm

bdist <- function( x, y , method ){
  a <- sum( x == 1 & y == 1)
  b <- sum( x == 1 & y == 0)
  c <- sum( x == 0 & y == 1)
  d <- sum( x == 0 & y == 0)
  
  switch(method,
         tanimoto = (a+d)/(a+2*(b+c)+d),
         stop("Method not implemented")
  )
}