R/specparF.R

In TIMP: Fitting Separable Nonlinear Models in Spectroscopy and Microscopy

"specparF" <- 
  function (specpar, xi, i, specref, specdispindex, specdisppar, 
            parmufunc = "") 
  {
    disppar <- specpar
    for (j in 1:length(specdispindex)) {
      if (parmufunc == "exp") {
        if (length(specdisppar[[j]]) == 1) 
          specdisppar[[j]][2] <- specdisppar[[1]][2]
        disppar[[specdispindex[[j]][1]]][specdispindex[[j]][2]] <- simpExp(specpar[[specdispindex[[j]][1]]][specdispindex[[j]][2]], 
                                                                           specdisppar[[j]], xi, specref)
      }
      if (parmufunc == "multiexp") {
        disppar[[specdispindex[[j]][1]]][specdispindex[[j]][2]] <- 
          simpExp(specpar[[specdispindex[[j]][1]]][specdispindex[[j]][2]], 
                  specdisppar[[j]], xi, specref)
        
        
      }
      if (parmufunc == "" || parmufunc == "poly") {
        disppar[[specdispindex[[j]][1]]][specdispindex[[j]][2]] <- simpPol(specpar[[specdispindex[[j]][1]]][specdispindex[[j]][2]], 
                                                                           specdisppar[[j]], xi, specref)
      }
    }
    disppar
  }