R/multiLin.R

In TIMP: Fitting Separable Nonlinear Models in Spectroscopy and Microscopy

"multiLin" <-
  function(thetaClass, diffs, parvec) {
    cnt <- 1
    newpar <- 0
    if (is.list(diffs$ind2)) {
      for (i in 1:length(parvec)) {
        newpar <- newpar + (parvec[[i]] * slot(thetaClass, diffs$what2[[i]])[[diffs$ind2[[i]]]])
      }
      # }
      # else {
      #  ## this is not working right when what2 indexes a list
      #  if(is.list(slot(thetaClass, diffs$what2))) {
      #    for(i in 1:length(diffs$ind2)) {
      #    newpar <- newpar + ( parvec[i] * slot(thetaClass, diffs$what2[i])[diffs$ind2[i]])
      #    }
    } else {
      for (i in 1:length(diffs$ind2)) {
        newpar <- newpar + (parvec[cnt] * slot(
          thetaClass,
          diffs$what2[i]
        )[diffs$ind2[i]])
        cnt <- cnt + 1
      }
    }
    # }
    newpar
  }