R/diffChange.R

In TIMP: Fitting Separable Nonlinear Models in Spectroscopy and Microscopy

"diffChange" <-
  function(modellist, diffsadd) {
    ## diffschange has structure
    ## list(list(what, dataset, spec), ...)
    ## e.g.,  list(list(what="kinpar", datatset=2,
    ##             spec=c(1,2,3), type="multifree"), ...)
    ## dataset can be a vector of indices into dataset list
    ## type is optional; if multifree then parameters in
    ## what are free per-dataset

    for (diffs in diffsadd) {
      if (length(diffs$type) == 0) {
        diffs$type <- "perd"
      }
      for (i in 1:length(diffs$dataset)) {
        slot(
          modellist[[diffs$dataset[i]]],
          diffs$what
        ) <- diffs$spec
      }
    }
    newl <- list()
    for (i in 1:length(diffsadd)) {
      if (length(diffsadd[[i]]$type) == 0) {
        diffsadd[[i]]$type <- "perdataset"
      }
      if (diffsadd[[i]]$type == "multifree") {
        for (j in 1:length(diffsadd[[i]]$dataset)) {
          newl[[length(newl) + 1]] <- diffsadd[[i]]
          newl[[length(newl)]]$dataset <- diffsadd[[i]]$dataset[j]
        }
      } else {
        newl[[length(newl) + 1]] <- diffsadd[[i]]
      }
    }


    list(modellist = modellist, diffsadd = newl)
  }