R/CD_mA.R

In CGE: Computing General Equilibrium

#' @export
CD_mA <- function(alpha, Beta, p) {
  # computing Cobb-Douglas demand structure matrix in a monetary economy
  nonnegative_Beta <- Beta
  nonnegative_Beta[Beta < 0] <- 0
  A <- CD_A(alpha, nonnegative_Beta, p)
  tmpA <- A

  Indx <- which(Beta < 0, arr.ind = T)
  for (k in 1:nrow(Indx)) {
    A[Indx[k, 1], Indx[k, 2]] <- t(p) %*% tmpA[, Indx[k, 2]] / (-Beta[Indx[k, 1], Indx[k, 2]])
  }
  A
}