R/diffRemove.R

In TIMP: Fitting Separable Nonlinear Models in Spectroscopy and Microscopy

"diffRemove" <-
  function(modellist, diffsremove) {
    ## diffsremove has structure
    ## list(list(what, ind, dataset), ...)
    ## e.g.,  list(list(what="kinpar", ind=c(1,2), dataset=2), ...)
    ## dataset can be a vector of indices into dataset list

    for (diffs in diffsremove) {
      for (i in 1:length(diffs$dataset)) {
        slot(
          modellist[[diffs$dataset[i]]],
          diffs$what
        ) <-
          slot(
            modellist[[diffs$dataset[i]]],
            diffs$what
          )[[diffs$ind[1]]][-diffs$ind[2]]
      }
    }
    modellist
  }