R/get_delta_mat.R

In tystan/deltacomp: Functions to Analyse Compositional Data and Produce Confidence Intervals for Relative Increases and Decreases in the Compositional Components

get_delta_mat <- function(deltas, comparisons, comps, mean_comps) {

  n_delta <- length(deltas)
  poss_comps0 <- create_comparison_matrix(comparisons, comps, mean_comps)
  
  delta_mat <- NULL
  for(d in 1:n_delta) {
    delta_mat <- rbind(delta_mat, deltas[d] * poss_comps0)
  }

  return(delta_mat)
  
}