aboutsummaryrefslogtreecommitdiff
path: root/R/parms.mkinfit.R
diff options
context:
space:
mode:
authorJohannes Ranke <jranke@uni-bremen.de>2022-09-28 16:34:57 +0200
committerJohannes Ranke <jranke@uni-bremen.de>2022-09-28 16:34:57 +0200
commit3529f5ff498d7d054c7b1911ddfc4b242902b85d (patch)
tree4c642bfddcc68e353fe75e8037d39ad8f269d56e /R/parms.mkinfit.R
parent75f361bed527b91bec205c5452add13247760d61 (diff)
Fix handling of multistart fits with failures
Diffstat (limited to 'R/parms.mkinfit.R')
-rw-r--r--R/parms.mkinfit.R69
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)
-}

Contact - Imprint