R/set_initModelClass_kin.R

In TIMP: Fitting Separable Nonlinear Models in Spectroscopy and Microscopy

setMethod("initModelClass", signature(model="kin"),                       
          function (model) 
          {
            model@mirf <-  length(model@measured_irf) != 0
            model@irf <-  length(model@irfpar) != 0 || model@mirf
            model@dispmu <-  length(model@parmu) != 0
            model@disptau <-  length(model@partau) != 0
            model@fullk <-  length(dim(model@kmat)) != 1
            if(model@fullk) 
              model@seqmod <- FALSE
            model@wavedep <- (model@dispmu || model@disptau || model@weight || 
                                model@lclp0 || model@lclpequ || length(model@parmu) > 
                                0)
            model@getX <- !(model@dispmu || model@disptau || model@weight || 
                              length(model@parmu) > 0)
            model@clpdep <- model@wavedep
            if(model@fullk) 
              model@ncomp <- nrow(model@kmat) + length(model@kinpar2)
            else if (model@numericalintegration)
              model@ncomp <- length(model@initialvals)
            else 
              model@ncomp <- length(model@kinpar) + length(model@kinpar2)
            if(length(model@nl)==0) {
              model@ncolc <- array(model@ncomp, 1)
            } else {
              model@ncolc <- array(model@ncomp, model@nl)
            }
            
            if (length(model@cohspec$type) == 0) # TODO: find a better test
              model@cohspec$type <- ""
            if (length(model@oscspec$type) == 0) # TODO: find a better test
              model@oscspec$type <- ""
            if(length(model@speckin2$seqmod) == 0)
              model@speckin2$seqmod <- FALSE
            if(length(model@speckin2$jvec) == 0) {
              model@speckin2$fullk <- FALSE
            } else { 
              model@speckin2$fullk <- TRUE
            }
            model@usekin2 <- if( length(model@kinpar2) == 0) FALSE
            else TRUE
            if (length(model@cohspec) != 0) 
              model <- getCoh(model)
            if (length(model@oscspec) != 0) 
              model <- getOsc(model)
            model <- getAnisotropy(model)  
            
            model
          })