R/getPrel.R

In TIMP: Fitting Separable Nonlinear Models in Spectroscopy and Microscopy

"getPrel" <-
function (model) 
{
	# dependent params. are specified as list or vector
	# in prel; get their indices in the parameter group 
        # in its vectorized form 

  prelspec <- model@prelspec
  ppars <- intersect(slotNames(theta()), slotNames(model))
  pinde <- vector("list", length(ppars))
  for(p in 1:length(pinde)) pinde[[p]] <- vector()
  names(pinde) <- ppars
  for(rel in prelspec) {
    if(length(rel$ind1) == 1) 
      pinde[[rel$what1]] <- append(pinde[[rel$what1]], rel$ind1)
    else {
      indf <- rel$ind1[2]
      if(rel$ind1[1] > 1) {
        slW<-slot(model, rel$what1)
        indf <- indf + sum(unlist(lapply(slW, length))[1:(rel$ind1[1]-1)]) 
      }
      pinde[[ rel$what1 ]] <-
        append(pinde[[ rel$what1 ]], indf)
    }
  }
  pinde
}