R/one_cluster_or_more_cata.R

In ClustBlock: Clustering of Datasets

.one_cluster_or_more_cata <- function(Data, nblo, nperm = 100, alpha = 0.05) {
  n <- nrow(Data)
  p <- ncol(Data)
  nvar <- p / nblo
  Blocks <- rep(nvar, nblo)
  J <- rep(1:nblo, times = Blocks)
  vp <- .vp_s(Data, nblo)[2]
  lambda_per <- NULL
  for (i in 1:nperm)
  {
    Data_per <- Data
    for (j in 1:ncol(Data_per))
    {
      Data_per[, j] <- sample(Data[, j])
    }
    lambda_per[i] <- .vp_s(Data_per, nblo)[2]
  }

  pval <- 0
  for (i in 1:nperm)
  {
    if (vp <= lambda_per[i]) {
      pval <- pval + 1
    }
  }
  pval <- pval / nperm
  decision <- pval > alpha
  return(list(decision = decision, pvalue = pval))
}