R/diffFree.R

In TIMP: Fitting Separable Nonlinear Models in Spectroscopy and Microscopy

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

    for (diffs in diffsfree) {
      for (i in 1:length(diffs$dataset)) {
        if (length(diffs$start) != 0) {
          if (length(diffs$ind) == 2) {
            slot(
              modellist[[diffs$dataset[i]]],
              diffs$what
            )[[diffs$ind[1]]][diffs$ind[2]] <-
              diffs$start
          } else {
            slot(
              modellist[[diffs$dataset[i]]],
              diffs$what
            )[diffs$ind] <- diffs$start
          }
        }
      }
    }
    modellist
  }