R/F_heq_nb_jac.R

In CenterForStatistics-UGent/RCM: Fit row-column association models with the negative binomial distribution for the microbiome

#' The jacobian of the linear equality constraints
#'
#' @param Alpha the current estimate of the environmental gradient
#' @param alphaK a matrix with the environmental gradients
#' of the lower dimensions
#' @param d an integer, the number of environmental variables,
#'  including dummies
#' @param k an integer, the current dimension
#' @param centMat a centering matrix
#' @param ... further arguments for other methods, not needed in this one
#'
#' @return The jacobian matrix
heq_nb_jac = function(Alpha, alphaK, d, k, 
    centMat, ...) {
    if (k == 1) {
        return(rbind(centMat, 2 * Alpha))
    } else {
        rbind(centMat, 2 * Alpha, t(alphaK))
    }
}