# comp2: Ties function: bivariate In BNPdensity: Ferguson-Klass Type Algorithm for Posterior Normalized Random Measures

## Description

This function computes the distinct observations (couples) and their frequencies in a bivariate numeric vetor.

## Usage

 `1` ```comp2(y, z) ```

## Details

For internal use.

## Examples

 ``` 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30``` ```## The function is currently defined as function (y, z) { if (length(y) != length(z)) stop("Vectors y and z should have equal length!") n <- length(y) matY <- outer(y, y, "==") matZ <- outer(z, z, "==") mat <- matY & matZ 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] zstar <- z[jstar] nstar <- apply(mat[, jstar], 2, sum) rstar <- length(nstar) idx <- match(y, ystar) return(list(ystar = ystar, zstar = zstar, nstar = nstar, rstar = rstar, idx = idx)) } ```

BNPdensity documentation built on May 29, 2017, 9:33 p.m.