test_ibdPlot <- function() {
  k0 <- c(0, 0, 0.25, 0.5, 0.75, 1)
  k1 <- c(0, 1, 0.5, 0.5, 0.25, 0)
  relation <- c("Dup", "PO", "FS", "HS", "FC", "U")
  color <- c("magenta", "cyan", "red", "blue", "lightgreen", "black")
  ibdPlot(k0, k1, relation=relation, color=color)
  ibdPlot(k0, k1, color=color)
  ibdPlot(k0, k1, relation=relation)
  ibdPlot(k0, k1)
  ibdAreasDraw()
  checkIdentical(relation, ibdAssignRelatedness(k0,k1))
}
