R/getConToPlot.R

In TIMP: Fitting Separable Nonlinear Models in Spectroscopy and Microscopy

"getConToPlot" <- function (C, cohspec, cohcol) 
{
    cohmax <- vector()
	  if(!identical(cohcol, 0)) {	
		  for (i in 1:length(cohcol)) { 

		        cohmax <- append(max(abs( C[,ncol(C)-i+1])), cohmax)
		  	cvec <- C[, ncol(C)-i+1 ] / 
					 max(abs( C[, ncol(C)-i+1 ] ))
			cvec[is.nan(cvec)]<-0 
			C[, ncol(C)-i+1 ] <- cvec
			
	          }       
	  }
	  attr(C, "max") <- cohmax

	  C
}