diff options
Diffstat (limited to 'R/parms.mkinfit.R')
-rw-r--r-- | R/parms.mkinfit.R | 62 |
1 files changed, 0 insertions, 62 deletions
diff --git a/R/parms.mkinfit.R b/R/parms.mkinfit.R deleted file mode 100644 index a1f2d209..00000000 --- a/R/parms.mkinfit.R +++ /dev/null @@ -1,62 +0,0 @@ -#' Extract model parameters from mkinfit models -#' -#' This function always returns degradation model parameters as well as error -#' model parameters, 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. Methods are implemented for -#' [mkinfit()] objects and for [mmkin()] objects. -#' @param \dots Not used -#' @return For mkinfit objects, a numeric vector of fitted model parameters. -#' For mmkin row objects, a matrix with the parameters with a -#' row for each dataset. If the mmkin object has more than one row, a list of -#' such matrices is returned. -#' @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? -#' @rdname parms -#' @export -parms.mkinfit <- function(object, transformed = FALSE, ...) -{ - if (transformed) object$par - else c(object$bparms.optim, object$errparms) -} - -#' @rdname parms -#' @export -parms.mmkin <- function(object, transformed = FALSE, ...) -{ - if (nrow(object) == 1) { - res <- sapply(object, parms, transformed = transformed, ...) - colnames(res) <- colnames(object) - } else { - res <- list() - for (i in 1:nrow(object)) { - res[[i]] <- parms(object[i, ], transformed = transformed, ...) - } - names(res) <- rownames(object) - } - return(res) -} |