R/getStartTri.R

In TIMP: Fitting Separable Nonlinear Models in Spectroscopy and Microscopy

"getStartTri" <- function(xx, nexp, ncomp){
  y <- as.matrix(rowSums(xx))
  dim(y) <- c(ncomp, nexp)
  ampList <- vector("list", length=nexp)
  for(i in 1:nexp) ampList[[i]] <- rep(0,ncomp)
  fixedamps <- vector("list", length=ncomp) 
  for(i in 1:ncomp) {
    ycol <- y[i,]
    mind <- which(ycol==max(ycol))
    ycol <- ycol /  max(ycol)
    fixedamps[[i]] <- c(mind,i)
    for(j in 1:nexp) 
      ampList[[j]][i] <- ycol[j]
  }
  list(fixedamps=fixedamps, ampList=ampList)
}