diff options
| -rw-r--r-- | NEWS.md | 2 | ||||
| -rw-r--r-- | R/mhmkin.R | 16 |
2 files changed, 18 insertions, 0 deletions
@@ -2,6 +2,8 @@ - 'R/mhmkin.R': New method for performing multiple hierarchical mkin fits in one function call, optionally in parallel. +- 'R/mhmkin.R': 'anova.mhmkin' for conveniently comparing the resulting fits. + - 'R/illparms.R': New generic to show ill-defined parameters with methods for 'mkinfit', 'mmkin', 'saem.mmkin' and 'mhmkin' objects. - 'R/multistart.R': New method for testing multiple start parameters for hierarchical model fits, with function 'llhist' and new generic 'parplot' for diagnostics, and new generics 'which.best' and 'best' for extracting the fit with the highest likelihood @@ -35,6 +35,12 @@ mhmkin <- function(objects, backend = "saemix", algorithm = "saem", ...) { #' @export #' @rdname mhmkin +mhmkin.mmkin <- function(object, ...) { + mhmkin(list(object), ...) +} + +#' @export +#' @rdname mhmkin mhmkin.list <- function(objects, backend = "saemix", ..., cores = if (Sys.info()["sysname"] == "Windows") 1 else parallel::detectCores(), cluster = NULL) @@ -190,3 +196,13 @@ update.mhmkin <- function(object, ..., evaluate = TRUE) { if(evaluate) eval(call, parent.frame()) else call } + +#' export +anova.mhmkin <- function(object, ..., + method = c("is", "lin", "gq"), test = FALSE, model.names = "auto") { + if (identical(model.names, "auto")) { + model.names <- paste(rownames(object), "-", colnames(object)) + } + rlang::inject(anova(!!!(object), method = method, test = test, model.names = model.names)) +} + |