R/perturbation.R

In Compositional: Compositional Data Analysis

Defines functions perturbation

Documented in perturbation

perturbation <- function(x, y, oper = "+") {

  if ( is.vector(y) ) {
    a <- Rfast::eachrow(x, y, oper = oper)
  } else {
    if (oper == "+")  a <- x * y
    if (oper == "-")  a <- x / y
  }
  
  a / Rfast::rowsums(a)
}

Try the Compositional package in your browser