R/graycode_to_partition.R

In senhu/FCBMA: Factor collapsing with Bayesian model averaging for regression models

#' Convert graycode to listing of partition
#'
#' Convert a graycode to listing of partition (equivalence class), to illustrate which elements are grouped together.
#'
#' @param x a graycode (in any format, not necessarily in canonical format)
#'
#' @return a character string listing equivalence class
#'
#' @seealso \code{\link[partitions]{listParts}}
#'
#' @examples
#' graycode_to_partition(c(1,2,3,1,2))
#' graycode_to_partition(c(2,1,3,2,1))
#'
#' @export graycode_to_partition

graycode_to_partition <- function(x) {
  y = split(seq_along(x), x)
  z = lapply(y, paste, collapse=",")
  res <- paste(c("(",paste(z, collapse = ")(" ), ")"), collapse = "")
  return(res)
}

# graycode.to.partition <- function(x) {
#   x <- convert.cano.graycode(x)
#   res <- split(seq_along(x), x)
#   class(res) <- c(class(res), "equivalence")
#   return(res)
# }