R/getAniA.R

In TIMP: Fitting Separable Nonlinear Models in Spectroscopy and Microscopy

"getAniA" <- function(A, dataset, ani, anipar) {
  if (ani$angle[dataset] != "MA" && !ani$useparperp) {
    dA <- nrow(A)
    if (ani$angle[dataset] == "PAR") {
      gamma <- 2
    }
    if (ani$angle[dataset] == "PERP") {
      gamma <- -1
    }
    for (i in 1:dA) {
      for (j in 1:dA) {
        if (ani$anifunc[i] == "const" || ani$anifunc[i] == "exp") {
          A[j, i] <- A[j, i] * (1 + (gamma *
            anipar[ani$parperm[[i]][1]]))
        }
        if (ani$anifunc[i] == "expvar") {
          A[j, i] <- A[j, i] * (1 + (gamma * (anipar[ani$parperm[[i]][1]]
          - anipar[ani$parperm[[i]][2]])))
        }
      }
    }
  }
  A
}