diff options
author | Johannes Ranke <jranke@uni-bremen.de> | 2022-09-28 16:34:57 +0200 |
---|---|---|
committer | Johannes Ranke <jranke@uni-bremen.de> | 2022-09-28 16:34:57 +0200 |
commit | 3529f5ff498d7d054c7b1911ddfc4b242902b85d (patch) | |
tree | 4c642bfddcc68e353fe75e8037d39ad8f269d56e /R/parms.mkinfit.R | |
parent | 75f361bed527b91bec205c5452add13247760d61 (diff) |
Fix handling of multistart fits with failures
Diffstat (limited to 'R/parms.mkinfit.R')
-rw-r--r-- | R/parms.mkinfit.R | 69 |
1 files changed, 0 insertions, 69 deletions
diff --git a/R/parms.mkinfit.R b/R/parms.mkinfit.R deleted file mode 100644 index 83766355..00000000 --- a/R/parms.mkinfit.R +++ /dev/null @@ -1,69 +0,0 @@ -#' Extract model parameters -#' -#' This function returns degradation model parameters as well as error -#' model parameters per default, in order to avoid working with a fitted model -#' without considering the error structure that was assumed for the fit. -#' -#' @param object A fitted model object. -#' @param \dots Not used -#' @return Depending on the object, a numeric vector of fitted model parameters, -#' a matrix (e.g. for mmkin row objects), or a list of matrices (e.g. for -#' mmkin objects with more than one row). -#' @seealso [saem], [multistart] -#' @examples -#' # mkinfit objects -#' fit <- mkinfit("SFO", FOCUS_2006_C, quiet = TRUE) -#' parms(fit) -#' parms(fit, transformed = TRUE) -#' -#' # mmkin objects -#' ds <- lapply(experimental_data_for_UBA_2019[6:10], -#' function(x) subset(x$data[c("name", "time", "value")])) -#' names(ds) <- paste("Dataset", 6:10) -#' \dontrun{ -#' fits <- mmkin(c("SFO", "FOMC", "DFOP"), ds, quiet = TRUE, cores = 1) -#' parms(fits["SFO", ]) -#' parms(fits[, 2]) -#' parms(fits) -#' parms(fits, transformed = TRUE) -#' } -#' @export -parms <- function(object, ...) -{ - UseMethod("parms", object) -} - -#' @param transformed Should the parameters be returned as used internally -#' during the optimisation? -#' @param errparms Should the error model parameters be returned -#' in addition to the degradation parameters? -#' @rdname parms -#' @export -parms.mkinfit <- function(object, transformed = FALSE, errparms = TRUE, ...) -{ - res <- if (transformed) object$par - else c(object$bparms.optim, object$errparms) - if (!errparms) { - res[setdiff(names(res), names(object$errparms))] - } - else return(res) -} - -#' @rdname parms -#' @export -parms.mmkin <- function(object, transformed = FALSE, errparms = TRUE, ...) -{ - if (nrow(object) == 1) { - res <- sapply(object, parms, transformed = transformed, - errparms = errparms, ...) - colnames(res) <- colnames(object) - } else { - res <- list() - for (i in 1:nrow(object)) { - res[[i]] <- parms(object[i, ], transformed = transformed, - errparms = errparms, ...) - } - names(res) <- rownames(object) - } - return(res) -} |