R/seqM.R

In CenterForStatistics-UGent/combi: Compositional omics model based visual integration

#' A small auxiliary function for the indices of the lagrange multipliers
#'
#' @param y an integer, the current dimension
#' @param nLambda1s the number of centering restrictions
#' @param normal a logical, is there a normalization restriction?
#'
#' @return a vector containing the ranks of the current lagrangian multipliers
seqM = function(y, normal = TRUE, nLambda1s = 1) {
    (y - 1) * (normal + nLambda1s + (y - 2)/2) +
        seq_len(y + nLambda1s - 1 + normal)
}