BNPdensity (version 2019.9.11)

comp1: Ties function: univariate

Description

This function computes the distinct observations and their frequencies in a numeric vector.

Usage

comp1(y)

Details

For internal use.

Examples

# NOT RUN {
## The function is currently defined as
function(y) {
  n <- length(y)
  mat <- outer(y, y, "==")
  jstar <- led <- rep(FALSE, n)
  for (j in seq(n)) {
    if (!led[j]) {
      jstar[j] <- TRUE
      if (j == n) {
        break
      }
      ji <- seq(j + 1, n)
      tt <- mat[ji, j] %in% TRUE
      led[ji] <- led[ji] | tt
    }
    if (all(led[-seq(j)])) {
      break
    }
  }
  ystar <- y[jstar]
  nstar <- apply(mat[, jstar], 2, sum)
  r <- length(nstar)
  idx <- match(y, ystar)
  return(list(ystar = ystar, nstar = nstar, r = r, idx = idx))
}
# }