R/baseCorlambda.R

In TIMP: Fitting Separable Nonlinear Models in Spectroscopy and Microscopy

"baseCorlambda" <-
  function(data, basev) {
    psisim <- data@psi.df
    tmin <- max(basev[1], 1)
    tmax <- min(basev[2], data@nt)
    lmin <- max(basev[3], 1)
    lmax <- min(basev[4], data@nl)

    for (j in lmin:lmax) {
      baseline <- sum(psisim[tmin:tmax, j]) / (length(tmin:tmax))
      psisim[, j] <- psisim[, j] - baseline
    }
    data@psi.df <- psisim
    data
  }