> dDA

function (x, cll, pool = TRUE) 

{

  x <- data.matrix(x)

  n <- nrow(x)

  p <- ncol(x)

  cl0 <- as.integer(min(cll, na.rm = TRUE) - 1)

  cll <- as.integer(cll) - cl0

  inaC <- is.na(cll)

  clL <- cll[!inaC]

  K <- max(clL)

  if (K != length(unique(clL))) 

    stop(sQuote("cll"), " did not contain *consecutive* integers")

  nk <- integer(K)

  m <- v <- matrix(0, p, K)

  colVars <- function(x, means = colMeans(x, na.rm = na.rm), 

                      na.rm = FALSE) {

    x <- sweep(x, 2, means)

    colSums(x * x, na.rm = na.rm)/(nrow(x) - 1)

  }

  sum.na <- function(x) sum(x, na.rm = TRUE)

  for (k in 1:K) {

    which <- (cll == k)

    nk[k] <- sum.na(which)

    lsk <- x[which, , drop = FALSE]

    m[, k] <- colMeans(lsk, na.rm = TRUE)

    if (nk[k] > 1) 

      v[, k] <- colVars(lsk, na.rm = TRUE, means = m[, 

                                                     k])

  }

  structure(list(call = match.call(), cl0 = cl0, n = n, p = p, 

                 K = K, means = m, vars = v, nk = nk, pool = pool), class = "dDA")

}