R/aspic-control.R

In laurieKell/mpb: Management Procedures

setMethod('setControl<-', signature(object='aspic',value="FLPar"), 
          function(object,value,min=0.1,max=10.0) {
  
  fit=control(object)[,"fit"]
  nms=dimnames(object@params)$params
  object@control=FLPar(array(rep(c(1,NA,NA,NA,1),each=length(nms)), 
                             dim     =c(length(nms),5,dims(object@params)$iter), 
                             dimnames=list(params=nms,option=c("fit","min","val","max","lambda"),
                                           iter  =dimnames(object@params)$iter)))
  object@control[nms,"val"]=value[nms,]
  object@control[nms,"min"]=value[nms,]*min
  object@control[nms,"max"]=value[nms,]*max
  object@control[,"fit"]=fit
  
  return(object)})