From b1740ade9a1746ccdb325b95915ef88872489f03 Mon Sep 17 00:00:00 2001 From: Johannes Ranke Date: Fri, 28 Oct 2022 11:59:54 +0200 Subject: Export status method for mhmkin, move to status.R --- R/mhmkin.R | 39 --------------------------------------- 1 file changed, 39 deletions(-) (limited to 'R/mhmkin.R') diff --git a/R/mhmkin.R b/R/mhmkin.R index 47b84b9d..7f3ff9fa 100644 --- a/R/mhmkin.R +++ b/R/mhmkin.R @@ -136,45 +136,6 @@ print.mhmkin <- function(x, ...) { print(status(x)) } -status.mhmkin <- function(object, ...) { - if (inherits(object[[1]], "saem.mmkin")) { - test_func <- function(fit) { - if (inherits(fit$so, "try-error")) { - return("E") - } else { - if (!is.null(fit$FIM_failed)) { - return_values <- c("fixed effects" = "Fth", - "random effects and error model parameters" = "FO") - return(paste(return_values[fit$FIM_failed], collapse = ", ")) - } else { - return("OK") - } - } - } - } else { - stop("Only mhmkin objects containing saem.mmkin objects currently supported") - } - result <- lapply(object, test_func) - result <- unlist(result) - dim(result) <- dim(object) - dimnames(result) <- dimnames(object) - - class(result) <- "status.mhmkin" - return(result) -} - -#' @export -print.status.mhmkin <- function(x, ...) { - class(x) <- NULL - print(x, quote = FALSE) - cat("\n") - if (any(x == "OK")) cat("OK: Fit terminated successfully\n") - if (any(x == "Fth")) cat("Fth: Could not invert FIM for fixed effects\n") - if (any(x == "FO")) cat("FO: Could not invert FIM for random effects and error model parameters\n") - if (any(x == "Fth, FO")) cat("Fth, FO: Could not invert FIM for fixed effects, nor for random effects and error model parameters\n") - if (any(x == "E")) cat("E: Error\n") -} - #' @export AIC.mhmkin <- function(object, ..., k = 2) { if (inherits(object[[1]], "saem.mmkin")) { -- cgit v1.2.1