R/ranpar.hist.R

In xpose4: Diagnostics for Nonlinear Mixed-Effect Models

#' @describeIn  par_cov_hist random parameter distributions
#' @export
ranpar.hist <-
  function(object,
           onlyfirst=TRUE,
           main="Default",
           
           ...) {
    

    if(any(is.null(xvardef("ranpar",object)))) {
      return(cat("No ETAs are defined in the current database!\n"))
    }
     
    ## create enpty list for plots
    number.of.plots <- 0
    for (i in xvardef("ranpar", object)) {
      number.of.plots <- number.of.plots + 1
    }    
    plotList <- vector("list",number.of.plots)
    plot.num <- 0 # initialize plot number
    
    ## loop 
    for (i in xvardef("ranpar", object)) {      
      
      xplot <- xpose.plot.histogram(i,
                                    object,
                                    main=NULL,
                                    onlyfirst = onlyfirst,
                                    pass.plot.list=TRUE,
                                    ...)

      
      plot.num <- plot.num+1
      plotList[[plot.num]] <- xplot
    }

    
    default.plot.title <- "Distribution of ETAs"
    plotTitle <- xpose.multiple.plot.title(object=object,
                                           plot.text = default.plot.title,
                                           main=main,
                                           ...)
    obj <- xpose.multiple.plot(plotList,plotTitle,...)
    return(obj)
    
  }