R/getDiffThetaClChange.R

In TIMP: Fitting Separable Nonlinear Models in Spectroscopy and Microscopy

"getDiffThetaClChange" <-
  function(th, parorder, modellist, thetaClasslist, diffchange) {
    pcnt <- 1
    for (diffs in diffchange) {
      d <- diffs$dataset[1]
      parchange <- modellist[[d]]@parnames

      if (diffs$what %in% parchange) {
        parslot <- getPar(
          modellist[[d]],
          parorder[[pcnt]], th, thetaClasslist[[d]]
        )
        pcnt <- pcnt + 1
        if (is.list(slot(modellist[[d]], diffs$what))) {
          tlist <- list()
          cnt <- 1

          for (i in 1:length(slot(modellist[[d]], diffs$what))) {
            tlist[[i]] <- parslot[cnt:(cnt +
              length(slot(
                modellist[[d]],
                diffs$what
              )[[i]]) - 1)]
            cnt <- cnt + length(slot(
              modellist[[d]],
              diffs$what
            )[[i]])
          }
          parslot <- tlist
        }
        for (datai in diffs$dataset) {
          slot(thetaClasslist[[datai]], diffs$what) <- parslot
        }
      }
      if (diffs$what == "prelspec") {
        thL <- addPrelCl(
          thetaClasslist[[d]], modellist[[d]],
          th, parorder[[pcnt]]
        )
        for (datai in diffs$dataset) {
          thetaClasslist[[datai]] <- thL
        }
        pcnt <- pcnt + 1
      }
    }

    thetaClasslist
  }